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.

GATK not calling varients found using samtools

kohrnbkohrnb Portland State UniversityMember

I'm trying to compare my custom SNP-calling pipeline to GATK, and have found that GATK doesn't seem to be calling some SNPs that my pipeline found. For my pipeline, I'm looking at proportion of samples with each identity at a position (requirement of a depth of at least 1000 in all samples at the position before attempting calling), as determined by a file from samtools mpileup (with depth adjusted to avoid downsampling), and calling a SNP if they aren't all the same. For GATK, I'm using HaplotypeCaller with no downsampling, --forceActive, and a ploidy of 1 (with the same 18 samples in each case). My question is: why might GATK not find some SNPs that my pipeline finds?

Best Answer

Answers

Sign In or Register to comment.