Question: gatk v4.14.0 does not recognize matematical operators
0
gravatar for dago
9 months ago by
dago2.6k
Germany
dago2.6k wrote:

Dear all I have problem with filtering the SNPs called with gatk v4.1.4.0. I have a vcf file that looks like this:

    contig00001     8227    .       A       G       40.68   .       AC=1;AF=0.100;AN=10;BaseQRankSum=-0.861;ClippingRankSum=0.000;DP=88;FS=13.035;MLEAC=1;MLEAF=0.100;MQ=60.00;MQRankSum=0.000;QD=0.46;ReadPosRankSum=-0.300;SOR=2.792      GT:AD:DP:GQ:PL  0/0/0/0/0/0/0/0/0/1:85,3:88:34:66,0,34,78,132,196,276,380,528,782,6199
    contig00001     8244    .       T       C       43.68   .       AC=1;AF=0.100;AN=10;BaseQRankSum=-0.838;ClippingRankSum=0.000;DP=80;FS=15.529;MLEAC=1;MLEAF=0.100;MQ=60.00;MQRankSum=0.000;QD=0.55;ReadPosRankSum=-1.716;SOR=2.783      GT:AD:DP:GQ:PL  0/0/0/0/0/0/0/0/0/1:77,3:80:31:69,0,31,71,119,178,251,347,482,716,5872
    contig00001     8846    .       C       T       72.68   .       AC=1;AF=0.100;AN=10;BaseQRankSum=0.659;ClippingRankSum=0.000;DP=88;FS=2.385;MLEAC=1;MLEAF=0.100;MQ=59.93;MQRankSum=0.273;QD=0.83;ReadPosRankSum=-4.696;SOR=1.102        GT:AD:DP:GQ:PL  0/0/0/0/0/0/0/0/0/1:84,4:88:31:98,0,31,72,124,186,264,366,510,759,3392
    contig00001     9854    .       A       G       42.69   .       AC=1;AF=0.100;AN=10;BaseQRankSum=0.463;ClippingRankSum=0.000;DP=79;FS=11.687;MLEAC=1;MLEAF=0.100;MQ=60.00;MQRankSum=0.000;QD=0.55;ReadPosRankSum=1.267;SOR=2.799        GT:AD:DP:GQ:PL  0/0/0/0/0/0/0/0/0/1:74,3:77:29:68,0,29,66,111,166,235,324,450,665,2970
    contig00001     19796   .       A       G       34.70   .       AC=1;AF=0.100;AN=10;BaseQRankSum=-2.543;ClippingRankSum=0.000;DP=66;FS=10.825;MLEAC=1;MLEAF=0.100;MQ=60.00;MQRankSum=0.000;QD=0.53;ReadPosRankSum=-0.600;SOR=0.829      GT:AD:DP:GQ:PL  0/0/0/0/0/0/0/0/0/1:63,3:66:23:60,0,23,54,93,140,199,275,384,572,2745
    contig00001     20699   .       T       C       47.70   .       AC=1;AF=0.100;AN=10;BaseQRankSum=0.326;ClippingRankSum=0.000;DP=66;FS=2.442;MLEAC=1;MLEAF=0.1

I would like to filter the SNPs now, e.g. gatk VariantFiltration -V GenomeA.rawSNPs.vcf -filter "QD < 2" --filter-name "snp_def" -O test.vcf

or

gatk VariantFiltration -R ../GenomeA_contigs.fa -V GenomeA.rawSNPs.vcf -O test.vcf --filter-name "snp" --filter-expression "MQ < 60.0"

However, I keep getting the following error:

A USER ERROR has occurred: Invalid argument '<'.

I believe there is a problem in recognizing the mathematical operators or with "/' I tried various combionation with or without "/' but I am still getting the same error. I cannto quite figure out where the problem is. Can you guys spot the problem?

Thank you

snp java • 215 views
ADD COMMENTlink modified 9 months ago • written 9 months ago by dago2.6k

I also posted this issue on the gatk github, sorry for the double posting

ADD REPLYlink written 9 months ago by dago2.6k
0
gravatar for Mark
9 months ago by
Mark800
Mark800 wrote:

I think the syntax is --filter-name then -filter. You specify filter first then the filter name.

See the documentation here: https://software.broadinstitute.org/gatk/documentation/tooldocs/current/org_broadinstitute_hellbender_tools_walkers_filters_VariantFiltration.php

ADD COMMENTlink written 9 months ago by Mark800

Thanks for your comment. I jsut tried and it gives me the same error.

ADD REPLYlink written 9 months ago by dago2.6k

Oh! you're missing the reference sequence so include -R seq.fa at the start. Follow the manual precisely. I think the order of flag input is very important (and spaces too).

ADD REPLYlink written 9 months ago by Mark800

Sorry, I did that as well, but I still get the same error. I edit my post now

ADD REPLYlink written 9 months ago by dago2.6k

Okay this is the last thing I can think of. Try matching exactly the type of values. In the first example you gave QD is a float so make sure its "QD < 2.0" and in the second MQ is two decimals (shoudn't make a diff but try) MQ < 60.00. I didn't think MQ was a float I thought it was an int so its weird to see it as a float in your file. Not sure whats happening

See the caveats mentioned here: https://software.broadinstitute.org/gatk/documentation/article.php?id=1255

ADD REPLYlink modified 9 months ago • written 9 months ago by Mark800
1

Thank you @Amar, I appreciate your suggestions. I tried what you said as well 20.0, 20, 20.00. Always the same. I personally believe that there is an issue with Java.

ADD REPLYlink written 9 months ago by dago2.6k

Hello, dago! If you finally solve this problem could you share how you fix it?

ADD REPLYlink written 6 months ago by fireprophet.io10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 688 users visited in the last hour