BenLangmead / bowtie2

A fast and sensitive gapped read aligner
GNU General Public License v3.0
664 stars 158 forks source link

cannot validate overall alignment rate #496

Closed shl129 closed 3 days ago

shl129 commented 3 days ago

Hi,

I got the output as below. However, based on the number of reads, the overall alignment rate is (18245046+5157319+143922)/24911033 = 93.94%. Not 95.63% in the output. Did I miss something?

24911033 reads; of these: 24911033 (100.00%) were paired; of these: 1508668 (6.06%) aligned concordantly 0 times 18245046 (73.24%) aligned concordantly exactly 1 time 5157319 (20.70%) aligned concordantly >1 times

1508668 pairs aligned concordantly 0 times; of these:
  143922 (9.54%) aligned discordantly 1 time
----
1364746 pairs aligned 0 times concordantly or discordantly; of these:
  2729492 mates make up the pairs; of these:
    2175150 (79.69%) aligned 0 times
    335437 (12.29%) aligned exactly 1 time
    218905 (8.02%) aligned >1 times

95.63% overall alignment rate