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