Notice:
If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We encourage our fourm members to be more involved, jump in and help out your fellow researchers with their questions. GATK forum is a community forum and helping each other with using GATK tools and research is the cornerstone of our success as a genomics research community.We appreciate your help!

Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

CombineVariants error: Allele in genotype not in variant context ...

nchuangnchuang Member
edited July 2015 in Ask the GATK team

Hi guys,

I am trying to find private variants that are not shared with 1000G or dbSNP vcfs. I followed the video tutorial method of using CombineVariants to combine my 3 vcfs before SelectVariants for the unique set. However, running CombineVariants throws this error below. I believe this might be related to using VQSR in INDEL mode because when I tried this without calibrating for Indels I did not get this error.

INFO 17:35:28,751 HelpFormatter - --------------------------------------------------------------------------------
INFO 17:35:28,754 HelpFormatter - The Genome Analysis Toolkit (GATK) v3.4-0-g7e26428, Compiled 2015/05/15 03:25:
41
INFO 17:35:28,754 HelpFormatter - Copyright (c) 2010 The Broad Institute
INFO 17:35:28,754 HelpFormatter - For support and documentation go to http://www.broadinstitute.org/gatk
INFO 17:35:28,759 HelpFormatter - Program Args: -T CombineVariants -R /local/projects-t3/1000G/hu
man_g1k_v37_decoy.fasta -V:GBR /local/projects-t3/1000G/gbr.tranche2.master.pass.vcf -V
:ATHENS /local/projects-t3/1000G/athens.tranche2.master.pass.vcf -V:DBSNP /local/projec
ts-t3/1000G/annovar/humandb/dbsnp142.vcf -V:1kG /local/projects-t3/1000G/gatk.resource/1k
G_phase3.snps.b37.vcf -o gbr.vs.athens.tranche2.dbsnp.1kG.combined.pass.vcf -nt 8
INFO 17:35:28,764 HelpFormatter - Executing on Linux 2.6.18-371.4.1.el5 amd64;
Java HotSpot(TM) 64-Bit Server VM 1.8.0_40-b26.
INFO 17:35:28,765 HelpFormatter - Date/Time: 2015/07/27 17:35:28
INFO 17:35:28,765 HelpFormatter - --------------------------------------------------------------------------------
INFO 17:35:28,765 HelpFormatter - --------------------------------------------------------------------------------
INFO 17:35:29,759 GenomeAnalysisEngine - Strictness is SILENT
INFO 17:35:29,992 GenomeAnalysisEngine - Downsampling Settings: Method: BY_SAMPLE, Target Coverage: 1000
INFO 17:35:31,990 MicroScheduler - Running the GATK in parallel mode with 8 total threads, 1 CPU thread(s) for each of 8 data thread(s), of 64 processors available on this machine
INFO 17:35:32,075 GenomeAnalysisEngine - Preparing for traversal
INFO 17:35:32,080 GenomeAnalysisEngine - Done preparing for traversal
INFO 17:35:32,081 ProgressMeter - [INITIALIZATION COMPLETE; STARTING PROCESSING]
INFO 17:35:32,081 ProgressMeter - | processed | time | per 1M | | total | remaining
INFO 17:35:32,082 ProgressMeter - Location | sites | elapsed | sites | completed | runtime | runtime
INFO 17:36:04,522 ProgressMeter - 1:18016781 2792752.0 32.0 s 11.0 s 0.6% 92.9 m 92.3 m
INFO 17:36:10,985 GATKRunReport - Uploaded run statistics report to AWS S3

ERROR ------------------------------------------------------------------------------------------
ERROR stack trace

java.lang.IllegalStateException: Allele in genotype <:DEL> not in the variant context [AGAAG, AAGGAAG, A, AAGGA
, AAGGAAGGAAGGA]
at htsjdk.variant.variantcontext.VariantContext.validateGenotypes(VariantContext.java:1320)
at htsjdk.variant.variantcontext.VariantContext.validate(VariantContext.java:1258)
at htsjdk.variant.variantcontext.VariantContext.(VariantContext.java:395)
at htsjdk.variant.variantcontext.VariantContextBuilder.make(VariantContextBuilder.java:496)
at htsjdk.variant.variantcontext.VariantContextBuilder.make(VariantContextBuilder.java:490)
at org.broadinstitute.gatk.utils.variant.GATKVariantContextUtils.simpleMerge(GATKVariantContextUtils.java
:1139)
at org.broadinstitute.gatk.tools.walkers.variantutils.CombineVariants.map(CombineVariants.java:339)
at org.broadinstitute.gatk.tools.walkers.variantutils.CombineVariants.map(CombineVariants.java:136)
at org.broadinstitute.gatk.engine.traversals.TraverseLociNano$TraverseLociMap.apply(TraverseLociNano.java
:267)
at org.broadinstitute.gatk.engine.traversals.TraverseLociNano$TraverseLociMap.apply(TraverseLociNano.java
:255)
at org.broadinstitute.gatk.utils.nanoScheduler.NanoScheduler.executeSingleThreaded(NanoScheduler.java:274
)
at org.broadinstitute.gatk.utils.nanoScheduler.NanoScheduler.execute(NanoScheduler.java:245)
at org.broadinstitute.gatk.engine.traversals.TraverseLociNano.traverse(TraverseLociNano.java:144)
at org.broadinstitute.gatk.engine.traversals.TraverseLociNano.traverse(TraverseLociNano.java:92)
at org.broadinstitute.gatk.engine.traversals.TraverseLociNano.traverse(TraverseLociNano.java:48)
at org.broadinstitute.gatk.engine.executive.ShardTraverser.call(ShardTraverser.java:98)
at java.util.concurrent.FutureTask.run(FutureTask.java:266)
at java.util.concurrent.ThreadPoolExecutor.runWorker(ThreadPoolExecutor.java:1142)
at java.util.concurrent.ThreadPoolExecutor$Worker.run(ThreadPoolExecutor.java:617)
at java.lang.Thread.run(Thread.java:745)

ERROR ------------------------------------------------------------------------------------------
ERROR A GATK RUNTIME ERROR has occurred (version 3.4-0-g7e26428):
ERROR
ERROR This might be a bug. Please check the documentation guide to see if this is a known problem.
ERROR If not, please post the error message, with stack trace, to the GATK forum.
ERROR Visit our website and forum for extensive documentation and answers to
ERROR commonly asked questions http://www.broadinstitute.org/gatk
ERROR
ERROR MESSAGE: Allele in genotype <:DEL> not in the variant context [AGAAG, AAGGAAG, A, AAGGA, AAGGAAGGAAGGA]
ERROR ------------------------------------------------------------------------------------------

Should be < *:DEL>

Answers

Sign In or Register to comment.