7

Is the MUMmer suite capable of calculating reference sequence coverage statistics for all query sequences collectively? It would be possible to achieve by parsing the output of nucmer / show-coords / show-tiling but it seems like there should be a better way. I currently do this using a sensitive read mapper, samtools depth and some scripting.

To clarify, I'd like to know the reference coverage achieved using all of the query sequences (i.e. the whole de novo assembly, in my case).

bedeabc
  • 248
  • 1
  • 6

1 Answers1

4

I believe all you need to do is to run dnadiff from MUMmer. That will run a comparison and output a number of useful metrics.

Devon Ryan
  • 19,602
  • 2
  • 29
  • 60
5heikki
  • 156
  • 2
  • Many thanks Devon. The QRY field of the AlignedBases record in dnadiff's out.report seems to be what I was after. Although clearly a powerful tool, I was surprised how slow dnadiff was to run. – bedeabc Aug 11 '17 at 18:31