Bug Bulletin: we have identified a bug that affects indexing when producing gzipped VCFs. This will be fixed in the upcoming 3.2 release; in the meantime you need to reindex gzipped VCFs using Tabix.

Error in BeagleOutputToVCF

hirokomatsuihirokomatsui Posts: 8Member
edited January 2013 in Ask the team

I'm having a problem running BeagleOutputToVCF, showing an error message:

##### ERROR stack trace
java.lang.IllegalStateException: Allele in genotype AT not in the variant context [A*, G]
        at org.broadinstitute.sting.utils.variantcontext.VariantContext.validateGenotypes(VariantContext.java:1197)
        at org.broadinstitute.sting.utils.variantcontext.VariantContext.validate(VariantContext.java:1137)

Is this because my Beagle output file include INDELs?

I'm using GenomeAnalysisTK-2.2-10-gbbafb72, and the command line is like:

> java -Xmx4g -jar /raid/software/src/GenomeAnalysisTK-2.3-5-g49ed93c/GenomeAnalysisTK.jar \
-R hg19_lite.fa \
-T BeagleOutputToVCF \
-V chr1.01.vcf \
-beagleR2:BEAGLE beagle/chr1.01.r2 \
-beaglePhased:BEAGLE chr1.01.phased \
-beagleProbs:BEAGLE chr1.01.gprobs \
-o beagle_vcf/chr1.01.vcf \
-L 1:1-10000000

Thanks, Hiroko Matsui

Post edited by Geraldine_VdAuwera on


Sign In or Register to comment.