Question: why allelic count on forward and reverse strand do not add up to depth in the results of samtools mpileup?
0
gravatar for mt1022
3.8 years ago by
mt1022220
China
mt1022220 wrote:
dmel_mitochondrion_genome       2878    .       A       <*>     0       .       DP=15;I16=0,9,0,0,347,13419,0,0,180,3600,0,0,140,2570,0,0;QS=1,0;MQ0F=0 PL:ADF:ADR      0,27,86:0,0:9,0

See the above line for example. Here the depth is DP=15, however allelic count on forward strand and reverse stand is 0, 0 and 9, 0, and sum up to 9 < 15. Why this happens?

Thanks in advance.

mpileup samtools vcf • 1.2k views
ADD COMMENTlink modified 3.8 years ago by Tom_L340 • written 3.8 years ago by mt1022220
0
gravatar for Tom_L
3.8 years ago by
Tom_L340
Tom_L340 wrote:

The first four values of the I16 tags are depths at Q13 threshold, corresponding to 10^(-13/10) base calling error, which is 5% (documentation: http://samtools.sourceforge.net/mpileup.shtml).

Would it be possible that 15 is the total depth without taking into account any base calling threshold whereas 9 passed the given threshold? Can you check that into the corresponding BAM file?

Cheers.

ADD COMMENTlink written 3.8 years ago by Tom_L340

the quality of base calling for those reads is ~30 at that position. I guess this is not the main problem.

ADD REPLYlink modified 3.7 years ago • written 3.7 years ago by mt1022220
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1146 users visited in the last hour