@@ -55,32 +55,26 @@ Paired-end illumina reads and PacBio is optional <br/>
5555
5656
5757# Pipeline in 6 steps: <br />
58- ### 1. Downsampling of reads to a designated coverage using ``` Filtlong ```
58+ ### 1. Downsampling of reads to a designated coverage using [ Filtlong] ( https://github.com/rrwick/Filtlong )
5959#####   ;   ; -default is set to 100X (-x); which provided better assemblies compared to the typical 30-50X
60- ### 2. _ Optional_ : Polishing of downsampled reads with the paired-end illumina reads using ``` Ratatosk correct ```
60+ ### 2. _ Optional_ : Polishing of downsampled reads with the paired-end illumina reads using [ Ratatosk] ( https://github.com/DecodeGenetics/Ratatosk )
6161#####   ;   ; -uses a baseline quality score (-Q) of 90 and therefore assumes mildly recent ONT data (e.g. R10 or high-accuracy basecalling)
6262### 3. Genome Assembly
63- #### 3.a. Assembly with``` Flye ```
63+ #### 3.a. Assembly with [ Flye] ( https://github.com/mikolmogorov/Flye )
6464#####   ;   ; -removed the hard coded maximium value for the minimum overlap threshold (previously 10kb)
6565#####   ;   ; -by default the minimum overlap value is automatically provided as the read N90 after polishing
66- #### 3.b. Assembly with ``` Hifiasm ```
66+ #### 3.b. Assembly with [ Hifiasm] ( https://github.com/chhylp123/hifiasm )
6767#####   ;   ; -if Hifi reads are provided: uses the ``` --ul ``` option, with both polished ONT and Hifi reads
6868#####   ;   ; -without Hifi: uses the ``` --ont ``` option, with only the polished ONT reads
69- ### 4. 'Patch' the Flye assembly (target) using the the Hifiasm assembly (query) with ``` Ragtag patch ```
69+ ### 4. 'Patch' the Flye assembly (target) using the the Hifiasm assembly (query) with [ Ragtag] ( https://github.com/malonge/RagTag/wiki/ patch)
7070#####   ;   ; -uses a minimum unique alignment length (-f) of 25000 to be conservative during patching
71- ### 5. _ Optional_ : Polishing of assembly with PacBio Hifi and paired-end illumina reads using ``` NextPolish2 ```
72- ### 6. Filtering (minimum contig length 10kb), reordering and renaming using ``` Seqkit ``` and ``` awk ```
71+ ### 5. _ Optional_ : Polishing of assembly with PacBio Hifi and paired-end illumina reads using [ NextPolish2] ( https://github.com/Nextomics/NextPolish2 )
72+ ### 6. Filtering (minimum contig length 10kb), reordering and renaming using [ Seqkit] ( https://bioinf.shenwei.me/seqkit ) and ``` awk ```
73+ ### 7. Comprehensive evaluation of all assemblies using [ PAQman] ( https://github.com/SAMtoBAM/PAQman )
7374
7475## Schematic
7576
7677<p align =" center " >
77- <img src="https://github.com/SAMtoBAM/fusemblr/blob/main/figures/fusemblr_schematic.png" width=70%>
78+ <img src="https://github.com/SAMtoBAM/fusemblr/blob/main/figures/fusemblr_schematic.updated. png" width=70%>
7879</p >
7980
80-
81-
82- Following assembly it is recommended that you run [ PAQman] ( https://github.com/SAMtoBAM/PAQman ) on your resulting assembly to comprehensively check the quality <br />
83- It is recommended to feed your resulting assembly to PAQman alongside the 1.filtlong/* .fz.gz set of reads <br />
84- This can also help you compare any assemblies you have to check for the best.
85-
86-
0 commit comments