Samtools 0.1.19 and Samtools 1.2 Read depth values for forward/reverse strand do not match
0
0
Entering edit mode
10.1 years ago

Hi All,

I have used Samtools 0.1.19 and 1.2 versions mpileup to view positions with mismatches to the reference genome.

However with exactly the same settings the first 4 values in the I16 tag (reference bases forward, reference bases reverse strand, alt forward, alternative reverse) do not match between the two versions:

samtools mpileup -g -B -Q 30 -f genome.fa -l my.bed my.bam1 my.bam2 my.bam3 > my.bcf

Compare

bcftools/1.0

bcftools view my.vcf | grep 64753439

file generated with Samtools 1.2

chr7    64753439    .    A    G,<X>    0    .    DP=792;I16=122,200,36,285,15279,812733,11867,477401,6440,128800,6420,128400,2991,47307,93,105;QS=2.57949,3.42051,0;VDB=0;SGB=-104.385;RPB=0;MQB=1;MQSB=1;BQB=0;MQ0F=0    PL    0,37,123,255,247,212    0,55,155,255,233,209    114,0,21,150,159,190    128,0,77,213,255,240    113,0,39,173,214,203    112,0,16,197,255,215

file generated with Samtools 0.1.19

chr7    64753439    .    A    G,X    0    .    DP=792;I16=119,299,38,321,15322,563512,11935,397431,8360,167200,7180,143600,3909,60917,133,147;QS=0.537967,0.462033,0.000000,0.000000;VDB=0;RPB=23.5676    PL    0,113,162,255,255,239    0,106,176,255,255,226    113,0,40,161,197,206    131,0,102,242,255,255    125,0,79,212,255,245    122,0,28,237,255,231

Can someone help to identify the reason for this discrepancy?

Thank you very much!

Dessi

DP4 read-depth samtools • 3.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 2361 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6