Learn how easy it is to sync an existing GitHub or Google Code repo to a SourceForge project! See Demo

Close

Confusing Output of Mate/paired Mapping

2011-12-15
2013-06-05
  • I liked the output stats regarding the number of mapped, multi-mapped and unmapped reads when mapping single-end/shotgun reads. However, It took me a while figure out the numbers for mates I was mapping as the output is somewhat more complicated and convoluted. Could I suggest that the output is changed to something like:

    xxx reads; of these:
      xxx (x%) were discarded on input
      xxx (x%) were paired; of these pairs:
        xxx (x%) aligned concordantly exactly 1 time
        xxx (x%) aligned concordantly >1 times
        xxx (x%) aligned concordantly exactly 1 time
        xxx (x%) aligned discordantly >1 times
        xxx (x%) had one or both mates unaligned; of these:
         xxx (x%) aligned 0 times
         xxx (x%) aligned exactly 1 time
         xxx (x%) aligned >1 times

    I think this makes it clearer how many were mappable as pairs, of those which were concordant and discordant, those which were not mapped as pairs and those which were not mapped at all.

    My $0.02 worth :)
    Nathan

     
  • spot the typo!

     
  • HI,
    i was trying to align paired end reads to human genome, i was  not able to get the correct sam alignment as i dint even receive the sam file summary at the end.