bamtools stats output
0
0
Entering edit mode
6.2 years ago
GK1610 ▴ 120

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 • 3.6k views
ADD COMMENT
2
Entering edit mode

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

fin swimmer

ADD REPLY

Login before adding your answer.

Traffic: 1701 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