Go Back   SEQanswers > Bioinformatics > Bioinformatics
Similar Threads
Thread Thread Starter Forum Replies Last Post
No reference/alternate (DP4) read depth in VCF? Tally Bioinformatics 0 02-20-2012 03:33 PM
Maximum read Depth in Samtools Anjali Bioinformatics 2 01-16-2012 02:15 AM
How to determine maximum read depth in Samtools VarFilter Muhua wang Bioinformatics 0 10-06-2011 03:09 PM
How to decide maximum read depth in Samtools VarFilter Muhua wang Illumina/Solexa 0 10-06-2011 03:02 PM
What is read depth -D100 means in samtools? ketan_bnf Bioinformatics 1 07-28-2011 06:54 PM

Thread Tools
Old 04-16-2012, 05:22 AM   #1
Junior Member
Location: Berlin

Join Date: Jul 2011
Posts: 8
Default Read Depth in vcf (samtools / bcftools)

I generated one vcf file from 4 different bam-files (4 samples) and had a look at some variants in the vcf and in the IGViewer (IGV).

I don't understand why the value for DP in the formatstring (i.e. for one sample) often differs (is less) from what I can see in the IGV. Meaning, vcf tells me a certain variant in a certain sample has read depth (DP) of 3 but I see more reads covering that position in the IGV. I thought that maybe bases of bad quality where left out, but reads and base qualities are good. Is there any other measure that I don't know yet which filters out certain reads / bases from being reported in the DP field in the format string (and consequently is not used for genotype assignment)?

Any hint is appreciated.
Marie_Noir is offline   Reply With Quote
Old 04-17-2012, 06:48 AM   #2
Location: USA

Join Date: Apr 2012
Posts: 12

I also had same question when I started my analysis. Then I figured out that DP (Depth of coverage) is not exactly the "counts of reads" at that position. Please see following links and I am sure you will get the answer:

Hope it helps
baika is offline   Reply With Quote

read depth, samtools bcftools, vcf

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 06:54 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2022, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO