This is a test version of Biostars. For the public version, visit https://www.biostars.org.
bamtools stats output

I have merged 200 bam files using bamtools. When I use bamtools stats options on the merged bam file, I get 117% each for forward strands and reverser strands. why would that happen? it should be 50% for both strands.

**********************************************
Stats for BAM file(s):
**********************************************

Total reads:       3158957612
Mapped reads:      3158957612   (100%)
Forward strand:    3726962454   (117.981%)
Reverse strand:    3726962454   (117.981%)
Failed QC:         0    (0%)
Duplicates:        0    (0%)
Paired-end reads:  3158957612   (100%)
'Proper-pairs':    3158957612   (100%)
Both pairs mapped: 3158957612   (100%)
Read 1:            3726962454
Read 2:            3726962454
Singletons:        0    (0%)
alignment sequencing

Looked weird. bamtools looked a bit old. Could you please check the stats with samtools flagstat ?

fin swimmer

0 answers

No answers yet.

Log in to answer this question.