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).
How to calculate overall reference coverage with MUMmer?
assembly;coverage
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.