This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Extracting information from gatk VariantFiltration output

Hi

I am filtering by GATK VariantFiltration option and I need variant corresponding to PASS column in my output file , how I can extract only variant passed the filter?

This is my command 

-bash-4.1$ java -jar /home/local/software/GATK/3.7/source/GenomeAnalysisTK.jar -T VariantFiltration -R ./hs37d5.fa -o somatic.snvs.vcf -V /somatic.snvs.vcf --filterExpression "AB < 0.2 || MQ0 > 50" --filterName ″HARD_TO_VALIDATE″^C
-bash-4.1$

Any help?

gatk wgs
• 1,910 views • link
written by zizigolu zizigolu

1 answer 

java -jar ./GenomeAnalysisTK.jar -T SelectVariants -R ./genome.fa -V ./snvs.vcf -ef -o ./selected.vcf
• 0 views • link
written by zizigolu zizigolu

Log in to answer this question.

More posts like this