To celebrate the release of GATK 4.0, we are giving away free credits for running the GATK4 Best Practices pipelines in FireCloud, our secure online analysis portal. It’s first come first serve, so sign up now to claim your free credits worth $250. Sponsored by Google Cloud. Learn more at https://software.broadinstitute.org/firecloud/documentation/freecredits

SNP calling filtering

Hi,

I am using RNA-seq data for SNP calling. I have used the Unified genotyper tool from GATK with this command:

java -Xmx22g -jar /bubo/sw/apps/bioinfo/GATK/2.3.6/GenomeAnalysisTK.jar -T UnifiedGenotyper -R path/to/myreference.fa -I path/to/mybamfile.bam -stand_call_conf 30 -stand_emit_conf 30 -o snps.raw.vcf

I have obtained SNPs mostly with very high phred scale but the filter field is "." for all my SNPs, there is no PASS or LowQual in the field.
I read in the GATK forum that:
" If the FILTER field is a ".", then no filtering has been applied to the records, meaning that all of the records will be used for analysis but without explicitly saying that any PASS. You should avoid such a situation by always filtering raw variant calls before analysis. "

I was wondering about the filtering process, is it a parameter that I should add to my command or I must do it in another way?

Thanks a lot for your help.

Best Answer

Answers

Sign In or Register to comment.