Notice:
If you happen to see a question you know the answer to, please do chime in and help your fellow community members. 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.

SelectVariants

tony743tony743 canadaMember

I am trying to select heterozygotes from the VCF file. I use this command many times and it works all the time. But something error pops up this time, do not know what is the problem. I check the documentary and seems like syntax is right because I use it all the time

1) partial data line information is like this

CHROM POS ID REF ALT QUAL FILTER INFO FORMAT 4751_42R

chr10 3101362 rs33880920 G C 189.60 . ABHom=1.00;AC=0;
AF=0.00;AN=0;BaseQRankSum=0.742;DB;DP=7;Dels=0.00;FS=0.000;HaplotypeScore=0.0000
;MQ=57.19;MQ0=0;MQRankSum=-0.742;QD=31.60;ReadPosRankSum=0.742 GT ./.
chr10 3101529 rs29314734 C T 116.56 . ABHet=0.500;ABHo
m=0.933;AC=0;AF=0.00;AN=0;BaseQRankSum=-0.212;DB;DP=14;Dels=0.00;FS=5.332;Haplot
ypeScore=0.0000;MQ=60.00;MQ0=0;MQRankSum=0.000;OND=0.056;QD=16.65;ReadPosRankSum
=-0.271 GT ./.

2) when I run
$ java -Xmx20g -jar /hpf/tools/centos6/gatk/3.6.0/GenomeAnalysisTK.jar -T SelectVariants -R genome.fa --variant 4751_42map.vcf -select 'vc.getGenotype("4751_42R”).isHet()' -o 4751_42maph.vcf

3) Error is

ERROR A GATK RUNTIME ERROR has occurred (version 3.6-0-g89b7209):
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 https://www.broadinstitute.org/gatk
ERROR
ERROR MESSAGE: Argument select-0has a bad value. Invalid expression used (vc.getGenotype("4751_42R”).isHet()). Please see the JEXL docs for correct syntax.
ERROR ------------------------------------------------------------------------------------------
Tagged:

Answers

Sign In or Register to comment.