GATK reported AD and DP

spreeth84 · 03-27-2012, 01:51 PM

I am using Unified Genotyper for SNP calling. I find that the reported Allele depths and read depths seem much lower than expected (with and without the base quality recalibration). I did the alignment using CLC bio and have a SNP file from CLC bio which reports the coverage at the SNP position of reads that passed the quality filter. I understand the reported coverage values may not be exact but for equivalent quality scores, I find the DP from GATK to be 10% or 20% of the CLC bio coverage. Has anyone else seen this? Does GATK scale down the coverage?

spreeth84 · 03-28-2012, 08:43 AM

it turned out I had to set Downsampling to None. My AD and DP were under reported due to the downsampling.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Assemblathon:Unable to reproduce results reported for SOAPdenovo	Mahtab	Bioinformatics	0	12-21-2011 05:01 PM
GATK / DepthofCoverage	nguyendofx	Bioinformatics	0	11-07-2011 11:21 AM
FS tag in GATK	ulz_peter	Bioinformatics	0	10-10-2011 01:42 AM
GATK help!	adaptivegenome	Bioinformatics	0	01-17-2011 08:01 PM
TopHat: ambiguous hits are only reported for one of the two paired end reads	GoneSouth	RNA Sequencing	0	01-19-2010 10:34 AM

03-27-2012, 01:51 PM	#1
spreeth84 Junior Member Location: Boston Join Date: Jan 2011 Posts: 9	GATK reported AD and DP I am using Unified Genotyper for SNP calling. I find that the reported Allele depths and read depths seem much lower than expected (with and without the base quality recalibration). I did the alignment using CLC bio and have a SNP file from CLC bio which reports the coverage at the SNP position of reads that passed the quality filter. I understand the reported coverage values may not be exact but for equivalent quality scores, I find the DP from GATK to be 10% or 20% of the CLC bio coverage. Has anyone else seen this? Does GATK scale down the coverage?

03-28-2012, 08:43 AM	#2
spreeth84 Junior Member Location: Boston Join Date: Jan 2011 Posts: 9	it turned out I had to set Downsampling to None. My AD and DP were under reported due to the downsampling.