QUAST

Quality Assessment Tool for Genome Assemblies

{"date":"23 May 2025, Friday, 02:46:39","assembliesNames":["genome.nextpolish"],"referenceName":"","order":[0],"report":[["Alignment-based statistics",[]],["Reads mapping",[]],["Misassemblies",[]],["Unaligned",[]],["Per base quality",[{"metricName":"# N's per 100 kbp","quality":"Less is better","values":["0.00"],"isMain":true},{"metricName":"# N's","quality":"Less is better","values":[0],"isMain":false}]],["Statistics without reference",[{"metricName":"# contigs","quality":"Equal","values":[102],"isMain":true},{"metricName":"# contigs (>= 0 bp)","quality":"Equal","values":[102],"isMain":false},{"metricName":"# contigs (>= 1000 bp)","quality":"Equal","values":[102],"isMain":false},{"metricName":"# contigs (>= 5000 bp)","quality":"Equal","values":[44],"isMain":false},{"metricName":"# contigs (>= 10000 bp)","quality":"Equal","values":[13],"isMain":false},{"metricName":"# contigs (>= 25000 bp)","quality":"Equal","values":[13],"isMain":false},{"metricName":"# contigs (>= 50000 bp)","quality":"Equal","values":[12],"isMain":false},{"metricName":"Largest contig","quality":"More is better","values":[398012],"isMain":true},{"metricName":"Total length","quality":"More is better","values":[2490931],"isMain":true},{"metricName":"Total length (>= 0 bp)","quality":"More is better","values":[2490931],"isMain":false},{"metricName":"Total length (>= 1000 bp)","quality":"More is better","values":[2490931],"isMain":true},{"metricName":"Total length (>= 5000 bp)","quality":"More is better","values":[2264899],"isMain":false},{"metricName":"Total length (>= 10000 bp)","quality":"More is better","values":[2064296],"isMain":true},{"metricName":"Total length (>= 25000 bp)","quality":"More is better","values":[2064296],"isMain":false},{"metricName":"Total length (>= 50000 bp)","quality":"More is better","values":[2023342],"isMain":true},{"metricName":"N50","quality":"More is better","values":[205781],"isMain":false},{"metricName":"N90","quality":"More is better","values":[5305],"isMain":false},{"metricName":"auN","quality":"More is better","values":["196625.3"],"isMain":false},{"metricName":"L50","quality":"Less is better","values":[5],"isMain":false},{"metricName":"L90","quality":"Less is better","values":[40],"isMain":false},{"metricName":"GC (%)","quality":"Equal","values":["42.56"],"isMain":false}]],["K-mer-based statistics",[]],["BUSCO completeness",[{"metricName":"Complete BUSCO (%)","quality":"More is better","values":["0.00"],"isMain":true},{"metricName":"Partial BUSCO (%)","quality":"Less is better","values":["0.00"],"isMain":true}]],["Predicted genes",[{"metricName":"# predicted rRNA genes","quality":"More is better","values":["10 + 1 part"],"isMain":true}]],["Similarity statistics",[]],["Reference statistics",[]]],"subreferences":[],"subreports":[],"minContig":500}
{{ qualities }}
{{ mainMetrics }}
{"filenames":["genome.nextpolish"],"lists_of_lengths":[[398012,346038,244984,211011,205781,159536,85234,83521,80948,76551,70648,61078,40954,9326,9031,8519,8137,7435,7207,7160,7032,6936,6936,6932,6713,6595,6227,6221,6208,6169,6145,5994,5936,5787,5768,5610,5604,5602,5435,5305,5199,5171,5165,5098,4983,4868,4847,4579,4455,4452,4449,4424,4383,4377,4350,4330,4306,4295,4271,4233,4219,4215,4168,4151,4148,4116,4102,4073,4058,4052,4038,4014,3981,3901,3900,3891,3834,3797,3780,3675,3662,3632,3609,3581,3574,3554,3524,3492,3468,3452,3452,3415,3365,3357,3343,3341,3341,3302,3269,3268,3143,2173]]}
{{ assembliesLengths }}
{{ referenceLength }}
{"tickX":1}
{"coord_x":[[0.0,1e-10,15.978443401282492,15.978443401382492,29.870357709627445,29.870357709727443,39.70539529196112,39.70539529206112,48.17656530831243,48.176565308412435,56.437773667757156,56.43777366785716,62.842447261686495,62.842447261786496,66.26422008477955,66.26422008487955,69.61722343975003,69.61722343985004,72.86693208282365,72.86693208292365,75.940120380693,75.940120380793,78.77632901112074,78.77632901122074,81.22834394047848,81.22834394057848,82.87246816551723,82.87246816561724,83.24686633230708,83.24686633240708,83.60942153756969,83.60942153766969,83.9514221790969,83.9514221791969,84.27808718908713,84.27808718918713,84.57656996520578,84.57656996530578,84.865899537161,84.865899537261,85.15334226439833,85.15334226449833,85.4356463507018,85.4356463508018,85.71409645630489,85.71409645640489,86.27083608498188,86.27083608508188,86.54033371458303,86.54033371468303,86.80509415957327,86.80509415967327,87.05508101187868,87.05508101197869,87.30482699039034,87.30482699049034,87.55405107568214,87.55405107578214,87.80170948131442,87.80170948141442,88.04840439177158,88.04840439187159,88.28903731175211,88.28903731185211,88.52734178505948,88.52734178515948,88.75966455915479,88.75966455925479,88.99122456623648,88.99122456633648,89.21644156341544,89.21644156351545,89.44141768680063,89.44141768690064,89.66631351892124,89.66631351902124,89.88450503044845,89.88450503054845,90.09747760977723,90.09747760987723,90.30619475208266,90.30619475218266,90.5137878166838,90.5137878167838,90.72114000749117,90.72114000759117,90.92580244093473,90.92580244103473,91.12584812666428,91.12584812676428,91.32127706467983,91.32127706477984,91.51586294441717,91.51586294451717,91.69968979469925,91.69968979479926,91.87853858657667,91.87853858667667,92.05726694155719,92.05726694165719,92.23587485964083,92.23587485974083,92.41347913691708,92.41347913701708,92.5894374432692,92.5894374433692,92.76515487582755,92.76515487592755,92.93978837631391,92.93978837641392,93.11361896415437,93.11361896425437,93.28648605681973,93.28648605691973,93.45891154752982,93.45891154762982,93.63037354306482,93.63037354316482,93.8003100045726,93.8003100046726,93.96968442722822,93.96968442732822,94.13889826735466,94.13889826745466,94.3062252627632,94.30622526286321,94.47286978242272,94.47286978252272,94.63939386518534,94.63939386528534,94.80463328771451,94.80463328781451,94.96931067139154,94.96931067149154,95.132823831732,95.132823831832,95.295734807588,95.295734807688,95.45840490965024,95.45840490975024,95.62051297286035,95.62051297296036,95.78165754089535,95.78165754099535,95.9414773030646,95.9414773031646,96.09808541465019,96.09808541475019,96.25465338060347,96.25465338070347,96.41086003586611,96.41086003596611,96.56477839008788,96.56477839018788,96.71721135591471,96.71721135601472,96.86896184599252,96.86896184609252,97.01649704467928,97.01649704477929,97.1635103501462,97.1635103502462,97.30931928664423,97.30931928674423,97.45420487359947,97.45420487369947,97.59796638285043,97.59796638295043,97.74144687267531,97.74144687277531,97.8841244498543,97.8841244499543,98.0255976580644,98.0255976581644,98.16578620604103,98.16578620614104,98.30501125884258,98.30501125894259,98.44359398152739,98.44359398162739,98.71927403850206,98.71927403860207,98.85436409117716,98.85436409127716,98.98913297879388,98.98913297889388,99.12333982755845,99.12333982765846,99.25746638505844,99.25746638515844,99.52415382039888,99.52415382049888,99.65538989237358,99.65538989247358,99.78658581871598,99.78658581881598,99.9127635410214,99.91276354112141,100.0,100.0000000001]],"coord_y":[[398012,398012,398012,346038,346038,244984,244984,211011,211011,205781,205781,159536,159536,85234,85234,83521,83521,80948,80948,76551,76551,70648,70648,61078,61078,40954,40954,9326,9326,9031,9031,8519,8519,8137,8137,7435,7435,7207,7207,7160,7160,7032,7032,6936,6936,6932,6932,6713,6713,6595,6595,6227,6227,6221,6221,6208,6208,6169,6169,6145,6145,5994,5994,5936,5936,5787,5787,5768,5768,5610,5610,5604,5604,5602,5602,5435,5435,5305,5305,5199,5199,5171,5171,5165,5165,5098,5098,4983,4983,4868,4868,4847,4847,4579,4579,4455,4455,4452,4452,4449,4449,4424,4424,4383,4383,4377,4377,4350,4350,4330,4330,4306,4306,4295,4295,4271,4271,4233,4233,4219,4219,4215,4215,4168,4168,4151,4151,4148,4148,4116,4116,4102,4102,4073,4073,4058,4058,4052,4052,4038,4038,4014,4014,3981,3981,3901,3901,3900,3900,3891,3891,3834,3834,3797,3797,3780,3780,3675,3675,3662,3662,3632,3632,3609,3609,3581,3581,3574,3574,3554,3554,3524,3524,3492,3492,3468,3468,3452,3452,3415,3415,3365,3365,3357,3357,3343,3343,3341,3341,3302,3302,3269,3269,3268,3268,3143,3143,2173,2173,0.0]],"filenames":["genome.nextpolish"]}
{{ coordNGx }}
{{ coordNAx }}
{{ coordNGAx }}
{{ coordmisassemblies }}
{{ featuresInContigs }}
{{ 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 }}
{"filenames":["genome.nextpolish"],"reference_index":null,"list_of_GC_distributions":[[[0.0,1.0,2.0,3.0,4.0,5.0,6.0,7.0,8.0,9.0,10.0,11.0,12.0,13.0,14.0,15.0,16.0,17.0,18.0,19.0,20.0,21.0,22.0,23.0,24.0,25.0,26.0,27.0,28.0,29.0,30.0,31.0,32.0,33.0,34.0,35.0,36.0,37.0,38.0,39.0,40.0,41.0,42.0,43.0,44.0,45.0,46.0,47.0,48.0,49.0,50.0,51.0,52.0,53.0,54.0,55.0,56.0,57.0,58.0,59.0,60.0,61.0,62.0,63.0,64.0,65.0,66.0,67.0,68.0,69.0,70.0,71.0,72.0,73.0,74.0,75.0,76.0,77.0,78.0,79.0,80.0,81.0,82.0,83.0,84.0,85.0,86.0,87.0,88.0,89.0,90.0,91.0,92.0,93.0,94.0,95.0,96.0,97.0,98.0,99.0,100.0],[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,7,10,10,19,37,38,50,81,96,148,174,213,231,243,337,389,464,536,599,724,754,871,975,1029,1173,1339,1433,1412,1434,1379,1280,1214,1044,890,779,655,534,412,364,318,282,221,161,147,102,83,64,49,34,21,15,12,11,5,4,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]]],"list_of_GC_contigs_distributions":[[[0,5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95,100],[0,0,0,0,0,1,1,3,21,18,52,4,2,0,0,0,0,0,0,0,0]]],"lists_of_gc_info":null}
{"links":["icarus_viewers/contig_size_viewer.html"],"links_names":["View in Icarus contig browser"]}
{ "# 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 than or equal to 0 bp.", "Total length (>= 0 bp)" : "is the total number of bases in the contigs having size greater than or equal to 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. Similarly, Nx for any x between 0 and 100 could be computed (x% of assembly).", "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. Similarly, NGx for any x between 0 and 100 could be computed (x% of the reference).", "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. Similarly, Lx for any x between 0 and 100 could be computed.", "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. Similarly, LGx for any x between 0 and 100 could be computed.", "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. Similarly, NAx for any x between 0 and 100 could be computed.", "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. Similarly, NGAx for any x between 0 and 100 could be computed.", "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. Similarly, LAx for any x between 0 and 100 could be computed.", "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. Similarly, LGAx for any x between 0 and 100 could be computed.", "auN" : "is the area under the Nx curve.", "auNG" : "is the area under the NGx curve.", "auNA" : "is the area under the NAx curve.", "auNGA" : "is the area under the NGAx curve.", "Average %IDY" : "is the average of alignment identity percent (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).", "# large block misassemblies" : "is the number of misassemblies between alignments with length greater than or equal to 3 kbp and with the misassembly threshold equal to 5 kbp (instead of default 1 kbp for regular misassemblies).", "# 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 (based on alignments to the combined reference).", "# 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.", "# large relocations" : "is the number of relocation events among all large block misassemblies. Relocation is a misassembly where the left flanking sequence aligns over 5 kbp away from the right flanking sequence on the reference, or they overlap by more than 5 kbp and both flanking sequences align on the same chromosome.", "# large translocations" : "is the number of translocation events among all large block misassemblies. Translocation is a misassembly where the flanking sequences align on different chromosomes.", "# large i/s translocations" : "is the number of interspecies translocation events among all large block misassemblies. Interspecies translocation is a misassembly where the flanking sequences align on different references (based on alignments to the combined reference).", "# large inversions" : "is the number of inversion events among all large block misassemblies. 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 the misassembly threshold (1 kbp by default).
  3. The left and right flanking sequences both are on the same strand of the same chromosome of the reference genome.
", "# scaffold gap ext. mis." : "is the number of scaffold gap size extensive misassemblies. We define a scaffold gap size misassembly as a breakpoint where the flanking sequences are combined in a scaffold on the wrong distance. These misassemblies are not included in the total number of misassemblies. ", "# scaffold gap loc. mis." : "is the number of scaffold gap size local misassemblies. Such breakpoints satisfy the local misassembly conditions but they occur inside scaffolds. These misassemblies are not included in the total number of local misassemblies. ", "# possibly misassembled contigs": "is the number of contigs that contain large unaligned fragment (default min length is 500 bp) and thus could possibly contain interspecies translocation with unknown reference.", "# possible misassemblies" : "is the number of putative interspecies translocations in possibly misassembled contigs if each large unaligned fragment is supposed to be a fragment of unknown reference.", "# intergenomic misassemblies" : "is the number of all found and putative (possible) interspecies translocations.", "# structural variations" : "is the number of misassemblies matched with structural variations.", "# possible TEs" : "is the number of misassemblies possibly caused by transposable elements (TEs). We define a possible TE as an event that satisfies these conditions:
  1. There are two misassembly breakpoints of the same type around a short alignment (less than 7 kbp by default)
  2. The gap between two long flanking sequences on the sides of the short alignment is less than 7 kbp.
  3. The long flanking sequences both are on the same strand of the same chromosome of the reference genome.
", "# unaligned mis. contigs" : "is the number of contigs that have the number of unaligned bases more than 50% of the contig length and a misassembly event in their aligned fragment. Note that such misassemblies are not counted in # misassemblies and other misassemblies statistics.", "# 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. Uncalled bases (N's) are not counted.", "# partially unaligned contigs" : "is the number of contigs that have at least one alignment to the reference sequence but also have at least one unaligned fragment of length ≥ unaligned-part-size threshold (uncalled bases (N's) are not counted towards the fragment length).", "Partially unaligned length" : "is the total number of unaligned bases in all partially unaligned contigs. Uncalled bases (N's) are not counted.", "# 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 100,000 aligned bases in the assembly.", "# mismatches" : "is the number of mismatches in all aligned bases in the assembly.", "# indels per 100 kbp" : "is the average number of indels per 100,000 aligned bases in the assembly.", "# indels" : "is the number of indels in all aligned bases in the assembly.", "# indels (<= 5 bp)" : "is the number of indels of length less than 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.", "# genomic features" : "is the number of genomic features (genes, transcripts, CDS) in the assembly (complete and partial), based on a user-provided annotated list of gene positions in the reference genome. A feature counts as 'partially covered' if the assembly contains at least 100 bp of this feature but not the whole feature.", "# 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 than or equal to 0 bp.", "Complete BUSCO (%)" : "is the percent of BUSCO (Universal Single-Copy Ortholog) genes found in the assembly in a complete form.", "Partial BUSCO (%)" : "is the percent of BUSCO (Universal Single-Copy Ortholog) genes found in the assembly in a partial form.", "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, 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 100,000 assembly bases.", "# mapped" : "is the number of reads that mapped to the assembly.", "Mapped (%)" : "is the percent of reads that mapped to the assembly.", "# properly paired" : "is the number of reads that mapped within expected range and orientation of each other to the assembly.", "Properly paired (%)" : "is the percent of reads that mapped within expected range and orientation of each other to the assembly.", "# singletons" : "is the number of reads which do not overlap any other reads.", "Singletons (%)" : "is the percent of reads which do not overlap any other reads.", "# misjoint mates" : "is the number of reads with mate mapped to a different contig.", "Misjoint mates (%)" : "is the percent of reads with mate mapped to a different contig.", "Avg. coverage depth" : "is the average depth of coverage.", "Coverage >= 1x (%)" : "is the total number of bases with at least 1x coverage, divided by the total length of the assembly.", "# 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)." }