| Line 65: |
| |
: [http://code.google.com/p/ngopt/wiki/How_To_Score_Genome_Assemblies_with_Mauve How to score genome assemblies using the Mauve system] |
|
: [http://code.google.com/p/ngopt/wiki/How_To_Score_Genome_Assemblies_with_Mauve How to score genome assemblies using the Mauve system] |
| |
|
|
|
| |
|
+ |
* '''Evaluated by Mauve Assembly Metrics to calculate the values for the left columns of "N50, Blast_IntactCDS"'''
|
| |
|
|
|
| |
|
+ |
: [http://code.google.com/p/ngopt/wiki/How_To_Score_Genome_Assemblies_with_Mauve How to score genome assemblies using the Mauve system] (mauve_linux_snapshot_2011-08-31)
|
| |
|
|
|
| |
* '''Evaluated by Blast with Features''' |
|
* '''Evaluated by Blast with Features''' |
| |
|
|
|
| Line 71: |
| |
: [http://gage.cbcb.umd.edu/index.html Gage] |
|
: [http://gage.cbcb.umd.edu/index.html Gage] |
| |
|
|
|
| - |
*'''Score with Mauve metrics:'''
|
+ |
*'''Score with Mauve Assembly Metrics, N50, Blast and GAGE:'''
|
| |
|
+ |
'''''Set1'''''
|
| |
{| {{table}} border="1" |
|
{| {{table}} border="1" |
| |
| align="center" style="background:#f0f0f0;"|'''Name''' |
|
| align="center" style="background:#f0f0f0;"|'''Name''' |
| Line 87: |
| |
| align="center" style="background:#f0f0f0;"|'''BrokenCDS''' |
|
| align="center" style="background:#f0f0f0;"|'''BrokenCDS''' |
| |
| align="center" style="background:#f0f0f0;"|'''IntactCDS''' |
|
| align="center" style="background:#f0f0f0;"|'''IntactCDS''' |
| - |
| align="center" style="background:#f0f0f0;"|'''ContigN50'''
|
+ |
| align="center" style="background:#f0f0f0;"|'''ContigN50^'''
|
| |
|
+ |
| align="center" style="background:#f0f0f0;"|'''ContigN90'''
|
| |
| align="center" style="background:#f0f0f0;"|'''MaxContigLength''' |
|
| align="center" style="background:#f0f0f0;"|'''MaxContigLength''' |
| - |
| align="center" style="background:#f0f0f0;"|'''!50'''
|
+ |
| align="center" style="background:#eee685;"|'''N50^'''
|
| - |
| align="center" style="background:#eee685;;"|'''Blast_IntactCDS'''
|
+ |
| align="center" style="background:#eee685;"|'''Blast_IntactCDS'''
|
| |
| align="center" style="background:#f0f0f0;"|'''Units(>200)''' |
|
| align="center" style="background:#f0f0f0;"|'''Units(>200)''' |
| - |
| align="center" style="background:#f0f0f0;"|'''N50'''
|
+ |
| align="center" style="background:#f0f0f0;"|'''N50^'''
|
| |
| align="center" style="background:#f0f0f0;"|'''cor.Units''' |
|
| align="center" style="background:#f0f0f0;"|'''cor.Units''' |
| - |
| align="center" style="background:#f0f0f0;"|'''cor.N50'''
|
+ |
| align="center" style="background:#f0f0f0;"|'''cor.N50^'''
|
| |
| align="center" style="background:#f0f0f0;"|'''Errors,(Indel>=5,Inv,Rel)''' |
|
| align="center" style="background:#f0f0f0;"|'''Errors,(Indel>=5,Inv,Rel)''' |
| |
|- |
|
|- |
| Line 126: |
| |
|- |
|
|- |
| |
|} |
|
|} |
| |
|
+ |
[^] Please note that the ContigN50 calculated by Mauve Assembly Metrics is incorrect (off-by-one error). We have followed the definition of N50 (''A contig N50 is calculated by first ordering every contig by length from longest to shortest. Next, starting from the longest contig, the lengths of each contig are summed, until this running sum equals one-half of the total length of all contigs in the assembly. The contig N50 of the assembly is the length of the shortest contig in this list.'' [http://www.ncbi.nlm.nih.gov/pubmed/22510764 ref]) to calculate N50s. GAGE's N50 was calculated using the total reference genome length rather than the sum total of contig lengths. The GAGE's cor.N50 values were computed after correcting contigs by breaking them at each error.
|
| |
|
|
|
| |
[#] Please note that GAA and minimus2 were designed to merge two assemblies at a time, we thus performed [http://sb.nhri.org.tw/CISA/en/H_volc_all2 all runs] and took the average scores. |
|
[#] Please note that GAA and minimus2 were designed to merge two assemblies at a time, we thus performed [http://sb.nhri.org.tw/CISA/en/H_volc_all2 all runs] and took the average scores. |
