The current GATK version is 3.7-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.
GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.
Register now for the upcoming GATK Best Practices workshop, Feb 20-22 in Leuven, Belgium. Open to all comers! More info and signup at

SNP calling filtering

HomaHoma Member Posts: 25


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


Sign In or Register to comment.