I developed own DNA assembly pipeline. Input is set of reads and output is set of contigs. Many papers measure own algorithms and compare it against each other. There are basic metrics like:
- N50, N70..,
- number of contigs,
- genom coverage,
- error rate
I download tool QUAST and there I can see metric N50, L50, where is easy to develop own computation. But how I can measure error rate of my output contigs? Is there any tool which can do it? I searched in many papers but I cannot find way how they do it. I would like to compare with some assemblers with some general way. I found only hints that it should be done by aligning contigs to reference genomes but I do not know how and if it is right way.