QUAST

Quality Assessment Tool for Genome Assemblies by CAB

Contigs are ordered from largest (contig #1) to smallest.

Contigs are broken into nonoverlapping 100 bp windows. Plot shows numbers of windows for each GC percentage.

{"assembliesWithNs":null,"minContig":500,"report":[["Genome statistics",[{"values":["96.042","96.000"],"quality":"More is better","isMain":true,"metricName":"Genome fraction (%)"},{"values":["1.030","1.002"],"quality":"Less is better","isMain":true,"metricName":"Duplication ratio"},{"values":[55106,55421],"quality":"More is better","isMain":true,"metricName":"Largest alignment"},{"values":[178143,172255],"quality":"More is better","isMain":true,"metricName":"Total aligned length"},{"values":[49658,49658],"quality":"More is better","isMain":false,"metricName":"NA50"},{"values":[17534,27260],"quality":"More is better","isMain":false,"metricName":"NA75"},{"values":[2,2],"quality":"Less is better","isMain":false,"metricName":"LA50"},{"values":[4,3],"quality":"Less is better","isMain":false,"metricName":"LA75"}]],["Misassemblies",[{"values":[0,0],"quality":"Less is better","isMain":true,"metricName":"# misassemblies"},{"values":[0,0],"quality":"Less is better","isMain":false,"metricName":" # relocations"},{"values":[0,0],"quality":"Less is better","isMain":false,"metricName":" # translocations"},{"values":[0,0],"quality":"Less is better","isMain":false,"metricName":" # inversions"},{"values":[0,0],"quality":"Less is better","isMain":false,"metricName":" # interspecies translocations"},{"values":[0,0],"quality":"Less is better","isMain":false,"metricName":"# possibly misassembled contigs"},{"values":[0,0],"quality":"Less is better","isMain":false,"metricName":"# misassembled contigs"},{"values":[0,0],"quality":"Less is better","isMain":true,"metricName":"Misassembled contigs length"},{"values":[0,0],"quality":"Less is better","isMain":false,"metricName":"# local misassemblies"}]],["Unaligned",[{"values":[2,1],"quality":"Less is better","isMain":false,"metricName":"# fully unaligned contigs"},{"values":[2730,885],"quality":"Less is better","isMain":false,"metricName":"Fully unaligned length"},{"values":[0,0],"quality":"Less is better","isMain":false,"metricName":"# partially unaligned contigs"},{"values":[0,0],"quality":"Less is better","isMain":false,"metricName":" # with misassembly"},{"values":[0,0],"quality":"Less is better","isMain":false,"metricName":" # both parts are significant"},{"values":[0,0],"quality":"Less is better","isMain":false,"metricName":"Partially unaligned length"}]],["Mismatches",[{"values":[57,47],"quality":"Less is better","isMain":false,"metricName":"# mismatches"},{"values":[2,0],"quality":"Less is better","isMain":false,"metricName":"# indels"},{"values":[2,0],"quality":"Less is better","isMain":false,"metricName":"Indels length"},{"values":["32.97","27.20"],"quality":"Less is better","isMain":true,"metricName":"# mismatches per 100 kbp"},{"values":["1.16","0.00"],"quality":"Less is better","isMain":true,"metricName":"# indels per 100 kbp"},{"values":[2,0],"quality":"Less is better","isMain":false,"metricName":" # indels (<= 5 bp)"},{"values":[0,0],"quality":"Less is better","isMain":false,"metricName":" # indels (> 5 bp)"},{"values":[14,0],"quality":"Less is better","isMain":false,"metricName":"# N's"},{"values":["7.74","0.00"],"quality":"Less is better","isMain":true,"metricName":"# N's per 100 kbp"}]],["Statistics without reference",[{"values":[20,12],"quality":"Less is better","isMain":true,"metricName":"# contigs"},{"values":[47,35],"quality":"Less is better","isMain":false,"metricName":"# contigs (>= 0 bp)"},{"values":[11,7],"quality":"Less is better","isMain":false,"metricName":"# contigs (>= 1000 bp)"},{"values":[5,6],"quality":"Less is better","isMain":false,"metricName":"# contigs (>= 5000 bp)"},{"values":[5,5],"quality":"Less is better","isMain":false,"metricName":"# contigs (>= 10000 bp)"},{"values":[3,3],"quality":"Less is better","isMain":false,"metricName":"# contigs (>= 25000 bp)"},{"values":[1,1],"quality":"Less is better","isMain":false,"metricName":"# contigs (>= 50000 bp)"},{"values":[55106,55421],"quality":"More is better","isMain":true,"metricName":"Largest contig"},{"values":[180875,173140],"quality":"More is better","isMain":true,"metricName":"Total length"},{"values":[187773,179095],"quality":"More is better","isMain":false,"metricName":"Total length (>= 0 bp)"},{"values":[173623,169301],"quality":"More is better","isMain":true,"metricName":"Total length (>= 1000 bp)"},{"values":[161370,165854],"quality":"More is better","isMain":false,"metricName":"Total length (>= 5000 bp)"},{"values":[161370,158866],"quality":"More is better","isMain":true,"metricName":"Total length (>= 10000 bp)"},{"values":[132733,132339],"quality":"More is better","isMain":false,"metricName":"Total length (>= 25000 bp)"},{"values":[55106,55421],"quality":"More is better","isMain":true,"metricName":"Total length (>= 50000 bp)"},{"values":[49658,49658],"quality":"More is better","isMain":false,"metricName":"N50"},{"values":[17535,27260],"quality":"More is better","isMain":false,"metricName":"N75"},{"values":[2,2],"quality":"Less is better","isMain":false,"metricName":"L50"},{"values":[4,3],"quality":"Less is better","isMain":false,"metricName":"L75"}]],["Predicted genes",[]],["Similarity statistics",[{"values":[4,4],"quality":"Equal","isMain":false,"metricName":"# similar correct contigs"},{"values":[0,0],"quality":"Equal","isMain":false,"metricName":"# similar misassembled blocks"}]],["Reference statistics",[{"values":[179997,179997],"quality":"Equal","isMain":false,"metricName":"Reference length"}]]],"referenceName":"combined_reference","date":"07 November 2016, Monday, 15:05:48","order":[0,1],"assembliesNames":["meta_contigs_1","meta_contigs_2"]}
{{ qualities }}
{{ mainMetrics }}
{"lists_of_lengths":[[55106,49658,27969,17535,11102,3898,3427,1610,1141,1120,1057,987,929,862,850,832,747,714,694,637],[55421,49658,27260,16099,10428,6988,3447,885,817,729,714,694]],"filenames":["meta_contigs_1","meta_contigs_2"]}
{"assemblies_lengths":[180875,173140],"filenames":["meta_contigs_1","meta_contigs_2"]}
{"reflen":[54999,64999,59999]}
{"tickX":1}
{"coord_y":[[55106,55106,55106,49658,49658,27969,27969,17535,17535,11102,11102,3898,3898,3427,3427,1610,1610,1141,1141,1120,1120,1057,1057,987,987,929,929,862,862,850,850,832,832,747,747,714,714,694,694,637,637,0.0],[55421,55421,55421,49658,49658,27260,27260,16099,16099,10428,10428,6988,6988,3447,3447,885,885,817,817,729,729,714,714,694,694,0.0]],"coord_x":[[0.0,1e-10,30.466344160331722,30.46634416043172,57.92066344160332,57.92066344170332,73.38382861091914,73.38382861101914,83.07836903939184,83.07836903949185,89.21630960608155,89.21630960618155,91.37138908085694,91.37138908095694,93.26606772633033,93.26606772643034,94.15618521078093,94.15618521088093,94.78700760193504,94.78700760203505,95.4062197650311,95.4062197651311,95.99060124395301,95.99060124405301,96.53628196268141,96.53628196278142,97.04989633724948,97.04989633734948,97.52646855563235,97.52646855573235,97.99640635798202,97.99640635808203,98.45639253628197,98.45639253638197,98.86938493434693,98.86938493444693,99.26413268832066,99.26413268842066,99.64782308223911,99.64782308233912,100.0,100.0000000001],[0.0,1e-10,32.00935659004274,32.009356590142744,60.69019290747372,60.690192907573724,76.4346771398868,76.4346771399868,85.7329328866813,85.7329328867813,91.75580455123021,91.75580455133021,95.79184474991337,95.79184475001337,97.78271918678526,97.78271918688526,98.29386623541643,98.29386623551643,98.76573870855955,98.76573870865955,99.18678526048285,99.18678526058285,99.5991683031073,99.5991683032073,100.0,100.0000000001]],"filenames":["meta_contigs_1","meta_contigs_2"]}
{{ coordNGx }}
{"coord_y":[[55106,55106,55106,49658,49658,27969,27969,17534,17534,11101,11101,3898,3898,3427,3427,1141,1141,1057,1057,987,987,929,929,862,862,850,850,832,832,747,747,714,714,694,694,637,637,0.0],[55421,55421,55421,49658,49658,27260,27260,16099,16099,10428,10428,6988,6988,3447,3447,817,817,729,729,714,714,694,694,0.0]],"coord_x":[[0.0,1e-10,30.466344160331722,30.46634416043172,57.92066344160332,57.92066344170332,73.38382861091914,73.38382861101914,83.07781617138907,83.07781617148908,89.21520387007602,89.21520387017603,91.37028334485142,91.37028334495142,93.2649619903248,93.26496199042481,93.89578438147892,93.89578438157892,94.48016586040083,94.48016586050083,95.02584657912924,95.02584657922924,95.5394609536973,95.5394609537973,96.01603317208017,96.01603317218017,96.48597097442986,96.48597097452986,96.94595715272979,96.94595715282979,97.35894955079475,97.35894955089475,97.75369730476848,97.75369730486848,98.13738769868694,98.13738769878694,98.48956461644782,98.48956461654782],[0.0,1e-10,32.00935659004274,32.009356590142744,60.69019290747372,60.690192907573724,76.4346771398868,76.4346771399868,85.7329328866813,85.7329328867813,91.75580455123021,91.75580455133021,95.79184474991337,95.79184475001337,97.78271918678526,97.78271918688526,98.25459165992838,98.25459166002838,98.67563821185168,98.67563821195168,99.08802125447615,99.08802125457615,99.48885295136884,99.48885295146884]],"filenames":["meta_contigs_1","meta_contigs_2"]}
{{ coordNGAx }}
{{ genesInContigs }}
{{ operonsInContigs }}
[{{ num_contigs }}, {{ Largest_alignment }}, {{ Total_aligned_length }}, {{ num_misassemblies }}, {{ Misassembled_contigs_length }}, {{ num_mismatches_per_100_kbp }}, {{ num_indels_per_100_kbp }}, {{ num_N's_per_100_kbp }}, {{ Genome_fraction }}, {{ Duplication_ratio }}, {{ NGA50 }}]
{{ allMisassemblies }}
{{ krona }}
{{ gcInfos }}
{"links_names":["View in Icarus contig browser"],"links":["icarus.html"]}
{ "# contigs" : "is the total number of contigs in the assembly.", "Largest contig" : "is the length of the longest contig in the assembly.", "Total length" : "is the total number of bases in the assembly.", "Reference length" : "is the total number of bases in the reference.", "# contigs (>= 0 bp)" : "is the total number of contigs in the assembly that have size greater or equal than 0 bp.", "Total length (>= 0 bp)" : "is the total number of bases in the contigs having size greater or equal than 0 bp.", "N50" : "is the contig length such that using longer or equal length contigs produces half (50%) of the bases of the assembly. Usually there is no value that produces exactly 50%, so the technical definition is the maximum length x such that using contigs of length at least x accounts for at least 50% of the total assembly length.", "NG50" : "is the contig length such that using longer or equal length contigs produces half (50%) of the bases of the reference genome. This metric is computed only if a reference genome is provided.", "N75" : "is the contig length such that using longer or equal length contigs produces 75% of the bases of the assembly. Usually there is no value that produces exactly 75%, so the technical definition is the maximum length x such that using contigs of length at least x accounts for at least 75% of the total assembly length.", "NG75" : "is the contig length such that using longer or equal length contigs produces 75% of the bases of the reference genome. This metric is computed only if a reference genome is provided.", "L50" : "is the minimum number of contigs that produce half (50%) of the bases of the assembly. In other words, it's the number of contigs of length at least N50.", "LG50" : "is the minimum number of contigs that produce half (50%) of the bases of the reference genome. In other words, it's the number of contigs of length at least NG50. This metric is computed only if a reference genome is provided.", "L75" : "is the minimum number of contigs that produce 75% of the bases of the assembly. In other words, it's the number of contigs of length at least N75.", "LG75" : "is the minimum number of contigs that produce 75% of the bases of the reference genome. In other words, it's the number of contigs of length at least NG75. This metric is computed only if a reference genome is provided.", "NA50" : "is N50 where the lengths of aligned blocks are counted instead of contig lengths. I.e., if a contig has a misassembly with respect to the reference, the contig is broken into smaller pieces. This metric is computed only if a reference genome is provided.", "NGA50" : "is NG50 where the lengths of aligned blocks are counted instead of contig lengths. I.e., if a contig has a misassembly with respect to the reference, the contig is broken into smaller pieces. This metric is computed only if a reference genome is provided.", "NA75" : "is N75 where the lengths of aligned blocks are counted instead of contig lengths. I.e., if a contig has a misassembly with respect to the reference, the contig is broken into smaller pieces. This metric is computed only if a reference genome is provided.", "NGA75" : "is NG75 where the lengths of aligned blocks are counted instead of contig lengths. I.e., if a contig has a misassembly with respect to the reference, the contig is broken into smaller pieces. This metric is computed only if a reference genome is provided.", "LA50" : "is L50 where aligned blocks are counted instead of contigs. I.e., if a contig has a misassembly with respect to the reference, the contig is broken into smaller pieces.", "LGA50" : "is LG50 where aligned blocks are counted instead of contigs. I.e., if a contig has a misassembly with respect to the reference, the contig is broken into smaller pieces.", "LA75" : "is L75 where aligned blocks are counted instead of contigs. I.e., if a contig has a misassembly with respect to the reference, the contig is broken into smaller pieces.", "LGA75" : "is LG75 where aligned blocks are counted instead of contigs. I.e., if a contig has a misassembly with respect to the reference, the contig is broken into smaller pieces.", "Average %IDY" : "is the average of alignment identity percent (Nucmer measure of alignment accuracy) among all contigs.", "# misassemblies" : "is the number of positions in the assembled contigs where the left flanking sequence aligns over 1 kbp away from the right flanking sequence on the reference (relocation) or they overlap on more than 1 kbp (relocation) or flanking sequences align on different strands (inversion) or different chromosomes (translocation).", "# misassembled contigs" : "is the number of contigs that contain misassembly events.", "Misassembled contigs length" : "is the number of total bases contained in all contigs that have one or more misassemblies.", "# relocations" : "is the number of relocation events among all misassembly events. Relocation is a misassembly where the left flanking sequence aligns over 1 kbp away from the right flanking sequence on the reference, or they overlap by more than 1 kbp and both flanking sequences align on the same chromosome.", "# translocations" : "is the number of translocation events among all misassembly events. Translocation is a misassembly where the flanking sequences align on different chromosomes.", "# interspecies translocations" : "is the number of interspecies translocation events among all misassembly events. Interspecies translocation is a misassembly where the flanking sequences align on different references (--meta only).", "# inversions" : "is the number of inversion events among all misassembly events. Inversion is a misassembly where it is not a relocation and the flanking sequences align on opposite strands of the same chromosome.", "# local misassemblies" : "is the number of local misassemblies. We define a local misassembly breakpoint as a breakpoint that satisfies these conditions:
  1. Two or more distinct alignments cover the breakpoint.
  2. The gap between left and right flanking sequences is less than 1 kbp.
  3. The left and right flanking sequences both are on the same strand of the same chromosome of the reference genome.
", "# scaffold gap size misassemblies" : "is the number of scaffold gap size misassemblies. We define scaffold gap size misassembly as a breakpoint where the flanking sequences combined in scaffold on the wrong distance. These misassemblies are not included in the total number of misassemblies. ", "# possibly misassembled contigs": "is the number of contigs that contain large unaligned fragment and thus could possibly contain interspecies translocation with unknown reference.", "# structural variations" : "is the number of misassemblies matched with structural variations.", "# fully unaligned contigs" : "is the number of contigs that have no alignment to the reference sequence.", "Fully unaligned length" : "is the total number of bases contained in all fully unaligned contigs.", "# partially unaligned contigs" : "is the number of contigs that are not fully unaligned ones but have fragments with no alignment to the reference.", "# with misassembly" : "is the number of partially unaligned contigs that contain misassembly events in their aligned fragment. Note that such misassemblies are not counted in # misassemblies and other misassemblies statistics.", "# both parts are significant" : "is the number of partially unaligned contigs that contain both aligned and unaligned fragments of length ≥ min-contig threshold.", "Partially unaligned length" : "is the total number of unaligned bases in all partially unaligned contigs.", "# ambiguous contigs" : "is the number of contigs that have reference alignments of equal quality in multiple locations on the reference.", "Ambiguous contigs length" : "is the total number of bases contained in all ambiguous contigs.", "Genome fraction (%)" : "is the total number of aligned bases in the reference, divided by the genome size. A base in the reference genome is counted as aligned if there is at least one contig with at least one alignment to this base. Contigs from repeat regions may map to multiple places, and thus may be counted multiple times in this quantity.", "GC (%)" : "is the total number of G and C nucleotides in the assembly, divided by the total length of the assembly.", "Reference GC (%)" : "is the total number of G and C nucleotides in the reference, divided by the total length of the reference.", "# mismatches per 100 kbp" : "is the average number of mismatches per 100000 aligned bases.", "# mismatches" : "is the number of mismatches in all aligned bases.", "# indels per 100 kbp" : "is the average number of indels per 100000 aligned bases.", "# indels" : "is the number of indels in all aligned bases", "# indels (<= 5 bp)" : "is the number of indels of length less or equal to 5 bp", "# indels (> 5 bp)" : "is the number of indels of length greater than 5 bp", "Indels length" : "is the number of total bases contained in all indels", "# genes" : "is the number of genes in the assembly (complete and partial), based on a user-provided annotated list of gene positions in the reference genome. A gene counts as 'partially covered' if the assembly contains at least 100 bp of this gene but not the whole gene.", "# operons" : "is the number of operons in the assembly (complete and partial), based on a user-provided annotated list of operon positions in the reference genome. An operon counts as 'partially covered' if the assembly contains at least 100 bp of this operon but not the whole operon.", "# predicted genes (unique)" : "is the number of unique genes in the assembly found by a gene prediction tool.", "# predicted genes (>= 0 bp)" : "is the number of found genes having length greater or equal than 0 bp.", "Cumulative length" : "plot shows the growth of assembly contig lengths. On the x-axis, contigs are ordered from largest (contig #1) to smallest. The y-axis gives the size of the x largest contigs in the assembly.", "Nx" : "plot shows the Nx metric value as x varies from 0 to 100. Nx is the minimum contig length y such that using contigs of length at least y accounts for at least x% of the total assembly length.", "NGx" : "plot shows the NGx metric value as x varies from 0 to 100. NGx is the minimum contig length y such that using contigs of length at least y accounts for at least x% of the bases of the reference genome. This metric is computed only if a reference genome is provided.", "NAx" : "plot shows the NAx metric value as x varies from 0 to 100. NAx is computed similarly to Nx, but based on lengths of aligned blocks instead of contig lengths. Contigs are broken into aligned blocks at misassembly breakpoints. NAx is the minimum block length y such that using blocks of length at least y accounts for at least x% of the bases of the assembly. This metric is computed only if a reference genome is provided.", "NGAx" : "plot shows the NGAx metric value as x varies from 0 to 100.NGAx is computed similarly to NGx, but based on lengths of aligned blocks instead of contig lengths. Contigs are broken at misassembly breakpoints. NGAx is the minimum block length y such that using blocks of length at least y accounts for at least x% of the bases of the reference genome. This metric is computed only if a reference genome is provided.", "GC content" : "plot shows the distribution of GC percentage among the contigs, i.e., the total number of bases in contigs with such GC content. Typically, the distribution is approximately Gaussian. However, for some genomes it is not Gaussian. For assembly projects with contaminants, the GC distribution of the contaminants often differs from the reference genome and may give a superposition of multiple curves with different peaks.", "Duplication ratio" : "is the total number of aligned bases in the assembly (i.e. Total length - Fully unaligned length - Partially unaligned length), divided by the total number of aligned bases in the reference (see the Genome fraction (%) metric). If the assembly contains many contigs that cover the same regions of the reference, its Duplication ratio may be much larger than 1. This may occur due to overestimating repeat multiplicities and due to small overlaps between contigs, among other reasons.", "Largest alignment" : "is the length of the largest continuous alignment in the assembly. This metric is always equal to the Largest contig metric but it can be smaller if the largest contig of the assembly contains a misassembly event.", "Total aligned length" : "is the total number of aligned bases in the assembly.", "Avg contig read support" : "is the average coverage of contigs that have large unique alignments to the reference. Read coverage is extracted from contig names (SPAdes/Velvet naming scheme only).", "# N's" : "is the total number of uncalled bases (N's) in the assembly.", "# N's per 100 kbp" : "is the average number of uncalled bases (N's) per 100000 assembly bases.", "# similar correct contigs" : "is the number of correct contigs similar among > 50% assemblies (see Icarus for visualization).", "# similar misassembled blocks" : "is the number of misassembled blocks similar among > 50% assemblies (see Icarus for visualization)." }