Entering edit mode
7.3 years ago
mt1022
▴
310
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.
the quality of base calling for those reads is ~30 at that position. I guess this is not the main problem.