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).

  • 1
    $\begingroup$ Would QUAST work on your data? It runs MUMmer and reports coverage. $\endgroup$
    – user172818
    Commented Aug 4, 2017 at 15:21
  • $\begingroup$ Good suggestion @user172818. I'll investigate. $\endgroup$
    – bedeabc
    Commented Aug 4, 2017 at 19:17

1 Answer 1


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.

  • $\begingroup$ 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. $\endgroup$
    – bedeabc
    Commented Aug 11, 2017 at 18:31

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.