The current GATK version is 3.6-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Powered by Vanilla. Made with Bootstrap.
Last chance to register for the GATK workshop next week in Basel, Switzerland! http://www.sib.swiss/training/upcoming-training-events/training/gatk-workshop-lecture

SelectVariants - Error

tinutinu Posts: 48Member
edited November 2013 in Ask the GATK team

Hi,

I am trying to subset a few samples from a VCF using the following command

java -Xmx8g -jar GenomeAnalysisTK-2.6-4-g3e5ff60/GenomeAnalysisTK.jar -T SelectVariants -R Homo_sapiens_assembly19.fasta -V Input.vcf -sf samples.txt-o Out.vcf

Getting the error

ERROR MESSAGE: Key Indel_FS,Indel_QD found in VariantContext field FILTER at 1:985458 but this key isn't defined in the VCFHeader.  We require all VCFs to have complete VCF headers by default.

I checked the Input.vcf and can find that the following lines exists in the VCF

##FILTER=<ID=Indel_FS,Description="FS>200.0">
##FILTER=<ID=Indel_QD,Description="QD<2.0">

Not sure why the error ?

Any help much appreciated..
Many Thanks,
Tinu

Post edited by Geraldine_VdAuwera on
Tagged:

Best Answer

Answers

  • Geraldine_VdAuweraGeraldine_VdAuwera Posts: 10,469Administrator, Dev admin

    Hi Tinu,

    For the program, Indel_FS is not equivalent to FS. What tool did you annotate you vcf with to produce the Indel_FS annotation?

    Geraldine Van der Auwera, PhD

  • tinutinu Posts: 48Member

    Hi Geraldine,

    GATK was used to annotate the VCF and a particular variant the FILTER field has values 'Indel_FS,Indel_QD'.
    So would it work if I change ',' to ';' in the VCF ?

  • tinutinu Posts: 48Member

    Yes, I tried changed the FILTER field value from 'Indel_FS,Indel_QD' to 'Indel_FS;Indel_QD' and it worked

    Many Thanks

Sign In or Register to comment.