This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to assess the quality of metagenomics assembly?

Hello everyone! I'm a beginner of metagenomics research. Presently we use soapdenovo package to perform metagenomics assembly then choose the optimal Scaffold according to N50 length (the longest one). I'm wondering that how should we assess the quality of metagenomics assembly except according to the N50 length? Thanks very much!

assembly metagenomics

OK, thanks all of you. I'll try QUAST later.

2 answers

Quast is, indeed, very useful. However, it is most useful when you have a reference, where it can calculate the number of misassemblies and error rates. For metagenomes, another crucial metric - arguably, the most important metric - is the percentage of reads that map back to the assembly. If only 50% of your reads map to the assembly... it is not very complete. But if 95% of your reads map to the assembly, then even if it is somewhat fragmented, that's probably very good.

It might also help to look at the percentage properly paired reads to detect any chimeras, something that seems especially relevant in a metagenome assembly.

Try quast

It gives you a report with the quality of your assembly.

Log in to answer this question.