Question: Tool for de novo assembly evaluation ??
0
gravatar for Picasa
17 months ago by
Picasa310
Picasa310 wrote:

Hello,

I've ran Soapdenovo, which produce a statistic file (scafStatistics).

However the N50 is different from the one produced by QUAST ...

Do you have an idea what is wrong ?

By the way, do you know a fast tool for computing assembler statistics such as N50, contig number etc.

quast soap • 702 views
ADD COMMENTlink modified 17 months ago • written 17 months ago by Picasa310
0
gravatar for Chris Fields
17 months ago by
Chris Fields1.8k
University of Illinois Urbana-Champaign
Chris Fields1.8k wrote:

We've used the assemblathon script for this for quite a while, but there are likely better ones out there now.

http://korflab.ucdavis.edu/datasets/Assemblathon/Assemblathon2/Basic_metrics/assemblathon_stats.pl

Re: variation in N50, it may be dependent on whether the tool is reporting the stats on raw vs. a filtered (removing seqs less than a minimum cutoff). We generally perform both, and if at all possible use an estimated genome size so the metrics (NG50 in this case) can be compared more readily across assemblies.

ADD COMMENTlink written 17 months ago by Chris Fields1.8k
0
gravatar for Picasa
17 months ago by
Picasa310
Picasa310 wrote:

Is assemblathon script permors filtering before computing N50 ?

Cause I have different result compare to Soapdenovo and QUAST.

ADD COMMENTlink written 17 months ago by Picasa310

For future reference, you should add this as a reply to my answer, not add this as another answer.

The script processes the sequences as is and performs no filtering; the other tools may give different results depending on whether they pre-filter or not. TBH I'm not sure.

EDIT: Also, make sure the N50 you're comparing is contig or scaffold-based (and whether there are any minimal N counts needed to split scaffolds into contigs).

ADD REPLYlink modified 17 months ago • written 17 months ago by Chris Fields1.8k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1209 users visited in the last hour