kkonganti@1: <tool id="cfsan_flye" name="CFSAN_Flye" version="2.8.1+galaxy0">
kkonganti@0:     <description>de novo assembler for single molecule sequencing reads</description>
kkonganti@0:     <requirements>
kkonganti@0:     	<requirement type="package" version="2.8.1">flye</requirement>
kkonganti@0:     </requirements>    
kkonganti@0:     <version_command>flye --version</version_command>
kkonganti@0:     <command detect_errors="exit_code"><![CDATA[
kkonganti@0:         #if $inputs.is_of_type('fastqsanger', 'fastq'):
kkonganti@0:             #set $ext = 'fastq'
kkonganti@0:         #elif $inputs.is_of_type('fastqsanger.gz', 'fastq.gz'):
kkonganti@0:             #set $ext = 'fastq.gz'
kkonganti@0:         #elif $inputs.is_of_type('fasta.gz'):
kkonganti@0:             #set $ext = 'fasta.gz'
kkonganti@0:         #elif $inputs.is_of_type('fasta'):
kkonganti@0:             #set $ext = 'fasta'
kkonganti@0: 	#end if
kkonganti@0: 	infile=\$(basename $inputs '.dat').${ext} &&
kkonganti@0: 	ln -s $inputs ./\${infile} &&   
kkonganti@0: 	flye
kkonganti@0:         $mode_conditional.mode
kkonganti@0: 	\${infile}
kkonganti@0:         -o out_dir
kkonganti@0:         -t \${GALAXY_SLOTS:-4}
kkonganti@0:         -i $iterations
kkonganti@0:         #if $mode_conditional.mode == '--pacbio-hifi' and $mode_conditional.hifi_error:
kkonganti@0:             --hifi-error $mode_conditional.hifi_error
kkonganti@0:         #end if
kkonganti@0:         #if $min_overlap:
kkonganti@0:             -m $min_overlap
kkonganti@0:         #end if
kkonganti@0:         #if $asm.asm_select == 'true':
kkonganti@0:             --asm-coverage $asm.asm_coverage
kkonganti@0:             -g '${asm.genome_size}'
kkonganti@0:         #end if
kkonganti@0:         $meta
kkonganti@0:         $scaffold
kkonganti@0:     ]]></command>
kkonganti@0:     <inputs>
kkonganti@0:         <param name="inputs" type="data" format="fasta,fasta.gz,fastq,fastq.gz,fastqsanger.gz,fastqsanger" label="Input reads" />
kkonganti@0:         <conditional name="mode_conditional">
kkonganti@0:             <param name="mode" type="select" label="Mode">
kkonganti@0:                 <option value="--nano-raw">Nanopore raw (--nano-raw)</option>
kkonganti@0:                 <option value="--nano-corr">Nanopore corrected (--nano-corr)</option>
kkonganti@0:                 <option value="--nano-hq">Nanopore HQ (--nano-hq)</option>
kkonganti@0:                 <option value="--pacbio-raw">PacBio raw (--pacbio-raw)</option>
kkonganti@0:                 <option value="--pacbio-corr">PacBio corrected (--pacbio-corr)</option>
kkonganti@0:                 <option value="--pacbio-hifi">PacBio HiFi (--pacbio-hifi)</option>
kkonganti@0:             </param>
kkonganti@0:             <when value="--nano-raw"/>
kkonganti@0:             <when value="--nano-corr"/>
kkonganti@0:             <when value="--nano-hq"/>
kkonganti@0:             <when value="--pacbio-raw"/>
kkonganti@0:             <when value="--pacbio-corr"/>
kkonganti@0:             <when value="--pacbio-hifi">
kkonganti@0:                 <param argument="--hifi-error" type="float" min="0" max="1" optional="true" label="Expected HiFi reads error rate" help="Default: 0.01"/>
kkonganti@0:             </when>
kkonganti@0:         </conditional>
kkonganti@0:         <param argument="--iterations" type="integer" value="1" label="Number of polishing iterations" 
kkonganti@0:             help="Polishing is performed as the final assembly stage. By default, Flye runs one polishing iteration. Additional iterations 
kkonganti@0:                 might correct a small number of extra errors (due to improvements on how reads may align to the corrected assembly). If the 
kkonganti@0:                 parameter is set to 0, the polishing is not performed"/>
kkonganti@0:         <param argument="--min-overlap" type="integer" min="1000" max="10000" optional="true" label="Minimum overlap between reads" 
kkonganti@0:             help="This sets a minimum overlap length for two reads to be considered overlapping. By default it is chosen 
kkonganti@0:                 automatically based on the read length distribution (reads N90) and does not require manual setting. Typical 
kkonganti@0:                 value is 3k-5k (and down to 1k for datasets with shorter read length). Intuitively, we want to set this 
kkonganti@0:                 parameter as high as possible, so the repeat graph is less tangled. However, higher values might lead to assembly gaps.
kkonganti@0:                 In some rare cases it makes sense to manually increase minimum overlap for assemblies of big genomes with long reads and high coverage." />
kkonganti@0:         <param argument="--keep-haplotypes" type="boolean" truevalue="--keep-haplotypes" falsevalue="" checked="False" label="Keep haplotypes"
kkonganti@0:             help="By default, Flye collapses graph structures caused by alternative haplotypes (bubbles, superbubbles, roundabouts) to produce longer 
kkonganti@0:                 consensus contigs. This option retains the alternative paths on the graph, producing less contigouos, but more detailed assembly."/>
kkonganti@0:         <param argument="--scaffold" type="boolean" truevalue="--scaffold" falsevalue="" label="Enable scaffolding using graph" 
kkonganti@0:             help="Starting from the version 2.9 Flye does not perform scaffolding by default, which guarantees that all assembled sequences do not have any gaps" />
kkonganti@0:         <param argument="--meta" type="boolean" truevalue="--meta" falsevalue="" checked="False" label="Perform metagenomic assembly" 
kkonganti@0:             help="It is designed for highly non-uniform coverage and is sensitive to underrepresented sequence at low coverage (as low as 2x).
kkonganti@0:                 In some examples of simple metagenomes, we observed that the normal mode assembled more contigious bacterial 
kkonganti@0:                 consensus sequence, while the metagenome mode was slightly more fragmented, but revealed strain mixtures"/>
kkonganti@0:         <conditional name="asm">
kkonganti@0:             <param name="asm_select" type="select" label="Reduced contig assembly coverage" help="Typically, assemblies of large genomes at high coverage require a hundreds of RAM. For high coverage assemblies, 
kkonganti@0:                         you can reduce memory usage by using only a subset of longest reads for initial contig extension stage (usually, the memory bottleneck)">
kkonganti@0:                 <option value="true">Enable reduced coverage for initial disjointing assembly</option>
kkonganti@0:                 <option value="false" selected="true">Disable reduced coverage for initial disjointing assembly</option>
kkonganti@0:             </param>
kkonganti@0:             <when value="true">
kkonganti@0:                 <param argument="--asm-coverage" type="integer" min="0" value="30"
kkonganti@0:                     label="Reduced coverage for initial disjointing assembly" 
kkonganti@0:                     help="This parameter specifies the target coverage of the longest reads. For a typical assembly, 30x is enough to produce good 
kkonganti@0:                         initial contigs. Regardless of this parameter, all reads will be used at the later pipeline stages."/>
kkonganti@0:                 <param argument="--genome-size" type="text" optional="true" label="Estimated genome size"
kkonganti@0:                     help="For example, 5m or 2.6g. No longer required as input. However, it must be used in conjunction with --asm-coverage option.">
kkonganti@0:                     <validator type="regex" message="Genome size must be a float  or integer, optionally followed by the a unit prefix (kmg)">^([0-9]*[.])?[0-9]+[kmg]?$</validator>
kkonganti@0:                 </param>
kkonganti@0:             </when>
kkonganti@0:             <when value="false" />
kkonganti@0:         </conditional>
kkonganti@0:         <param name="generate_log" type="boolean" truevalue="true" falsevalue="false" checked="false" label="Generate a log file"/>
kkonganti@0:     </inputs>
kkonganti@0:     <outputs>
kkonganti@0:         <data name="consensus" format="fasta" from_work_dir="out_dir/assembly.fasta" label="${tool.name} on ${on_string}: consensus"/>
kkonganti@0:         <data name="assembly_graph" format="graph_dot" from_work_dir="out_dir/assembly_graph.gv" label="${tool.name} on ${on_string}: assembly graph"/>
kkonganti@0:         <data name="assembly_gfa" format="txt" from_work_dir="out_dir/assembly_graph.gfa" label="${tool.name} on ${on_string}: graphical fragment assembly"/>
kkonganti@0:         <data name="assembly_info" format="tabular" from_work_dir="out_dir/assembly_info.txt" label="${tool.name} on ${on_string}: assembly info"/>
kkonganti@0:         <data name="flye_log" format="txt" from_work_dir="out_dir/flye.log" label="${tool.name} on ${on_string}: log">
kkonganti@0:             <filter>generate_log</filter>
kkonganti@0:         </data>
kkonganti@0:     </outputs>
kkonganti@0:     <tests>
kkonganti@0:         <!--Test 01: pacbio-raw-->
kkonganti@0:         <test expect_num_outputs="5">
kkonganti@0:             <param name="inputs" ftype="fastq.gz" value="ecoli.fastq.gz"/>
kkonganti@0:             <param name="mode" value="--pacbio-raw"/>
kkonganti@0:             <param name="iterations" value="0"/>
kkonganti@0:             <param name="generate_log" value="true"/>
kkonganti@0:             <output name="assembly_info" file="result1_assembly_info.txt" ftype="tabular" compare="sim_size"/>
kkonganti@0:             <output name="assembly_graph" file="result1_assembly_graph.dot" ftype="graph_dot" compare="sim_size"/>
kkonganti@0:             <output name="assembly_gfa" file="result1_assembly_graph.gfa" ftype="txt" compare="sim_size"/>
kkonganti@0:             <output name="consensus" file="result1_assembly.fasta" ftype="fasta" compare="sim_size"/>
kkonganti@0:             <output name="flye_log" file="result1.log" ftype="txt" compare="sim_size"/>
kkonganti@0:         </test>
kkonganti@0:         <!--Test 02: nano raw-->
kkonganti@0:         <test expect_num_outputs="4">
kkonganti@0:             <param name="inputs" ftype="fasta.gz" value="nanopore.fasta.gz"/>
kkonganti@0:             <param name="mode" value="--nano-raw"/>
kkonganti@0:             <param name="iterations" value="0"/>
kkonganti@0:             <output name="assembly_info" ftype="tabular">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="95" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_graph" ftype="graph_dot">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="803" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_gfa" ftype="txt">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="35047" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="consensus" ftype="fasta">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="35573" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:         </test>
kkonganti@0:         <!--Test 03: reduce coverage-->
kkonganti@0:         <test expect_num_outputs="4">
kkonganti@0:             <param name="inputs" ftype="fastq.gz" value="ecoli_hifi.fastq.gz"/>
kkonganti@0:             <conditional name="mode_conditional">
kkonganti@0:                 <param name="mode" value="--nano-raw"/>
kkonganti@0:             </conditional>
kkonganti@0:             <conditional name="asm">
kkonganti@0:                 <param name="asm_select" value="true" />
kkonganti@0:                 <param name="asm" value="30"/>
kkonganti@0:                 <param name="genome_size" value="3980000"/>
kkonganti@0:             </conditional>
kkonganti@0:             <output name="assembly_info" ftype="tabular">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="286" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_graph" ftype="graph_dot">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="1840" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_gfa" ftype="txt">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="420752" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="consensus" ftype="fasta">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="427580" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:         </test>
kkonganti@0:         <!--Test 04: metagenomic mode-->
kkonganti@0:         <test expect_num_outputs="4">
kkonganti@0:             <param name="inputs" ftype="fastq.gz" value="ecoli.fastq.gz"/>
kkonganti@0:             <conditional name="mode_conditional">
kkonganti@0:                 <param name="mode" value="--pacbio-raw"/>
kkonganti@0:             </conditional>
kkonganti@0:             <param name="meta" value="true"/>
kkonganti@0:             <output name="assembly_info" ftype="tabular">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="95" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_graph" ftype="graph_dot">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="367" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_gfa" ftype="txt">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="418729" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="consensus" ftype="fasta">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="425667" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:         </test>
kkonganti@0:         <!--Test 05: nanopore HQ mode-->
kkonganti@0:         <test expect_num_outputs="4">
kkonganti@0:             <param name="inputs" ftype="fastq.gz" value="ecoli_hifi.fastq.gz"/>
kkonganti@0:             <conditional name="mode_conditional">
kkonganti@0:                 <param name="mode" value="--nano-hq"/>
kkonganti@0:             </conditional>
kkonganti@0:             <param name="min_overlap" value="1000"/>
kkonganti@0:             <output name="assembly_info" ftype="tabular">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="286" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_graph" ftype="graph_dot">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="1248" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_gfa" ftype="txt">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="420252" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="consensus" ftype="fasta">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="427129" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:         </test>
kkonganti@0:         <!--Test 06: hifi error option-->
kkonganti@0:         <test expect_num_outputs="4">
kkonganti@0:             <param name="inputs" ftype="fastq.gz" value="ecoli_hifi.fastq.gz"/>
kkonganti@0:             <conditional name="mode_conditional">
kkonganti@0:                 <param name="mode" value="--pacbio-hifi"/>
kkonganti@0:                 <param name="hifi_error" value="0.21"/> 
kkonganti@0:             </conditional>
kkonganti@0:             <param name="min_overlap" value="1000"/>
kkonganti@0:             <output name="assembly_info" ftype="tabular">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="286" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_graph" ftype="graph_dot">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="1273" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_gfa" ftype="txt">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="420252" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="consensus" ftype="fasta">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="427129" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:         </test>
kkonganti@0:         <!--Test 07: keep haplotypes-->
kkonganti@0:         <test expect_num_outputs="4">
kkonganti@0:             <param name="inputs" ftype="fastq.gz" value="ecoli_hifi.fastq.gz"/>
kkonganti@0:             <conditional name="mode_conditional">
kkonganti@0:                 <param name="mode" value="--pacbio-corr"/>
kkonganti@0:                 <param name="hifi_error" value="0.21"/>
kkonganti@0:             </conditional>
kkonganti@0:             <param name="min_overlap" value="1000"/>
kkonganti@0:             <param name="keep-haplotypes" value="true"/>
kkonganti@0:             <output name="assembly_info" ftype="tabular">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="286" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_graph" ftype="graph_dot">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="1273" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_gfa" ftype="txt">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="420252" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="consensus" ftype="fasta">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="427129" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:         </test>
kkonganti@0:         <!--Test 08: scaffolding mode-->
kkonganti@0:         <test expect_num_outputs="4">
kkonganti@0:             <param name="inputs" ftype="fastq.gz" value="ecoli_hifi.fastq.gz"/>
kkonganti@0:             <param name="mode" value="--nano-hq"/>
kkonganti@0:             <param name="min_overlap" value="1000"/>
kkonganti@0:             <param name="scaffolding" value="true"/>
kkonganti@0:             <output name="assembly_info" ftype="tabular">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="286" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_graph" ftype="graph_dot">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="1248" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="assembly_gfa" ftype="txt">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="420252" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:             <output name="consensus" ftype="fasta">
kkonganti@0:                 <assert_contents>
kkonganti@0:                     <has_size value="427129" delta="100"/>
kkonganti@0:                </assert_contents>
kkonganti@0:             </output>
kkonganti@0:         </test>
kkonganti@0:     </tests>
kkonganti@0:     <help><![CDATA[
kkonganti@0: 
kkonganti@0: .. class:: infomark
kkonganti@0: 
kkonganti@0: **Purpose**
kkonganti@0: 
kkonganti@0: Flye is a de novo assembler for single molecule sequencing reads, such as those produced by PacBio and Oxford Nanopore Technologies. 
kkonganti@0: It is designed for a wide range of datasets, from small bacterial projects to large mammalian-scale assemblies. The package represents 
kkonganti@0: a complete pipeline: it takes raw PacBio/ONT reads as input and outputs polished contigs. Flye also has a special mode for metagenome 
kkonganti@0: assembly.
kkonganti@0: 
kkonganti@0: ----
kkonganti@0: 
kkonganti@0: .. class:: infomark
kkonganti@0: 
kkonganti@0: **Quick usage**
kkonganti@0: 
kkonganti@0: Input reads can be in FASTA or FASTQ format, uncompressed or compressed with gz. Currently, PacBio (raw, corrected, HiFi) and ONT reads 
kkonganti@0: (raw, corrected) are supported. Expected error rates are <30% for raw, <3% for corrected, and <1% for HiFi. Note that Flye was primarily 
kkonganti@0: developed to run on raw reads. You may specify multiple files with reads (separated by spaces). Mixing different read types is not yet supported. The *--meta* o
kkonganti@0: ption enables the mode for metagenome/uneven coverage assembly.
kkonganti@0: 
kkonganti@0: Genome size estimate is no longer a required option. You need to provide an estimate if using *--asm-coverage* option.
kkonganti@0: 
kkonganti@0: To reduce memory consumption for large genome assemblies, you can use a subset of the longest reads for initial disjointig assembly by 
kkonganti@0: specifying *--asm-coverage* and *--genome-size* options. Typically, 40x coverage is enough to produce good disjointigs.
kkonganti@0: 
kkonganti@0: ----
kkonganti@0: 
kkonganti@0: .. class:: infomark
kkonganti@0: 
kkonganti@0: **Outputs**
kkonganti@0: 
kkonganti@0: The main output files are:
kkonganti@0: 
kkonganti@0:     ::
kkonganti@0: 
kkonganti@0:         - Final assembly: contains contigs and possibly scaffolds (see below).
kkonganti@0:         - Final repeat graph: note that the edge sequences might be different (shorter) than contig sequences, because contigs might include multiple graph edges.
kkonganti@0:         - Extra information about contigs (such as length or coverage).
kkonganti@0: 
kkonganti@0: Each contig is formed by a single unique graph edge. If possible, unique contigs are extended with the sequence from flanking unresolved repeats on the graph. Thus, 
kkonganti@0: a contig fully contains the corresponding graph edge (with the same id), but might be longer then this edge. This is somewhat similar to unitig-contig relation in 
kkonganti@0: OLC assemblers. In a rare case when a repetitive graph edge is not covered by the set of "extended" contigs, it will be also output in the assembly file.
kkonganti@0: 
kkonganti@0: Sometimes it is possible to further order contigs into scaffolds based on the repeat graph structure. These ordered contigs will be output as a part of scaffold in 
kkonganti@0: the assembly file (with a scaffold prefix). Since it is hard to give a reliable estimate of the gap size, those gaps are represented with the default 100 Ns. 
kkonganti@0: assembly_info.txt file (below) contains additional information about how scaffolds were formed.
kkonganti@0: 
kkonganti@0: Extra information about contigs/scaffolds is output into the assembly_info.txt file. It is a tab-delimited table with the columns as follows:
kkonganti@0: 
kkonganti@0:     ::
kkonganti@0: 
kkonganti@0:         - Contig/scaffold id
kkonganti@0:         - Length
kkonganti@0:         - Coverage
kkonganti@0:         - Is circular, (Y)es or (N)o
kkonganti@0:         - Is repetitive, (Y)es or (N)o
kkonganti@0:         - Multiplicity (based on coverage)
kkonganti@0:         - Alternative group
kkonganti@0:         - Graph path (graph path corresponding to this contig/scaffold).
kkonganti@0: 
kkonganti@0: Scaffold gaps are marked with ?? symbols, and * symbol denotes a terminal graph node. Alternative contigs (representing alternative haplotypes) will have the same alt. 
kkonganti@0: group ID. Primary contigs are marked by *.
kkonganti@0: 
kkonganti@0: ----
kkonganti@0: 
kkonganti@0: .. class:: infomark
kkonganti@0: 
kkonganti@0: **Algorithm Description**
kkonganti@0: 
kkonganti@0: This is a brief description of the Flye algorithm. Please refer to the manuscript for more detailed information. The draft contig extension is organized as follows:
kkonganti@0: 
kkonganti@0:     ::
kkonganti@0: 
kkonganti@0:         - K-mer counting / erroneous k-mer pre-filtering
kkonganti@0:         - Solid k-mer selection (k-mers with sufficient frequency, which are unlikely to be erroneous)
kkonganti@0:         - Contig extension. The algorithm starts from a single read and extends it with a next overlapping read (overlaps are dynamically detected using the selected solid k-mers).
kkonganti@0: 
kkonganti@0: Note that we do not attempt to resolve repeats at this stage, thus the reconstructed contigs might contain misassemblies. Flye then aligns the reads on these draft 
kkonganti@0: contigs using minimap2 and calls a consensus. Afterwards, Flye performs repeat analysis as follows:
kkonganti@0: 
kkonganti@0:     ::
kkonganti@0: 
kkonganti@0:         - Repeat graph is constructed from the (possibly misassembled) contigs
kkonganti@0:         - In this graph all repeats longer than minimum overlap are collapsed
kkonganti@0:         - The algorithm resolves repeats using the read information and graph structure
kkonganti@0:         - The unbranching paths in the graph are output as contigs
kkonganti@0: 
kkonganti@0: If enabled, after resolving bridged repeats, Trestle module attempts to resolve simple unbridged repeats (of multiplicity 2) using the heterogeneities between repeat copies. 
kkonganti@0: Finally, Flye performs polishing of the resulting assembly to correct the remaining errors:
kkonganti@0: 
kkonganti@0:     ::
kkonganti@0: 
kkonganti@0:         - Alignment of all reads to the current assembly using minimap2
kkonganti@0:         - Partition the alignment into mini-alignments (bubbles)
kkonganti@0:         - Error correction of each bubble using a maximum likelihood approach
kkonganti@0: 
kkonganti@0: 
kkonganti@0: The polishing steps could be repeated, which might slightly increase quality for some datasets.
kkonganti@0: 
kkonganti@0: 
kkonganti@0:   ]]></help>
kkonganti@0:     <citations>
kkonganti@0:             <citation type="doi">10.1073/pnas.1604560113</citation>
kkonganti@0:             <citation type="bibtex">
kkonganti@0:                 @misc{githubFlye,
kkonganti@0:                 author = {Kolmogorov, Mijhail},
kkonganti@0:                 year = {2021},
kkonganti@0:                 title = {Flye},
kkonganti@0:                 publisher = {GitHub},
kkonganti@0:                 journal = {GitHub repository},
kkonganti@0:                 url = {https://github.com/fenderglass/Flye}}
kkonganti@0:             </citation>
kkonganti@0:         </citations>
kkonganti@0: </tool>