This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Read depth per base DNA seq

Hello everyone,

I am doing DNA seq germ line mutation analysis. I have used fastqc, cutadapt , bwa , samtools and gatk for most of the analysis.

The problem I am facing is when I calculate per base read depth by "samtools depth" it is different than the read depth values in the vcf file ( DP).

For example for the same sample the values are like that (C = read depth from samtools depth)

 A                     B                    C               D
Chr no.1        2939189             716             DP=71
Chr no.1        3140625                11            DP=1
Chr no.1        5983221             128             DP=50

I have used following commands

samtools depth -b sorted.bam > coveragefile

for vcf

gatk HaplotypeCaller -R ref.fa -I file.bam -O output.vcf

I would really appreciate if someone can tell me why I am getting this difference and which read depth values I should take.

Thanks

read-depth gatk samtools
• 1,304 views • link
updated by Ram Ram • written by Sabeen • Sabeen

This is probably related with how both samtools and GATK deal with duplicate reads, multi-mapping reads... etc. Most likely GATK is a way more restrictive than samtools. In addition, HaplotypeCaller performs a re-assembly of the reads before calling the variants, therefore the "mapping" of the reads that samtools is using and GATK might not be 100% equal.

• 0 views • link
written by iraun iraun

1 answer

A few possibilities:

  • it is unclear if you are using exactly the same bam file in both analysis since they have different names (file.bam and sorted.bam). This would be a trivial explanation to the discrepancy. For instance If you marked duplicates in file.bam before using gatk (as advised in most pipelines), then the duplicate reads are not counted.
  • gatk filters out low quality reads while samtools does not. If you have a lot of low quality reads, it could also explain the discrepancy.
  • other differences between gatk and samtools are discussed here: Difference Between Samtools And Gatk Algorithms
• 0 views • link
written by Carlo Yague Carlo Yague

Thanks everyone i understand now.

• 0 views • link
written by Sabeen • Sabeen

Log in to answer this question.

More posts like this