Error Correction using Racon
The software package Racon was developed to complement the minimap2/miniasm pipeline but can be used for any long-read assembly. It provides a fasta consensus algorithm that uses either 2nd generation short reads or raw noisy long-reads to correct draft assemblies.
Change into the miniasm directory in practicals/error_correction_practical/racon. For convenience the trimmed nanofilt reads as well as the miniasm assembly were copied into this directory. Use minimap2 to map the trimmed nanofilt reads back to the miniasm assembly and then use the filtered reads and the mapping to build the consensus assembly.
minimap2 ./miniasm_unitigs.fasta nanofilt_result.fastq \ > minimap.racon.paf racon nanofilt_result.fastq \ minimap.racon.paf miniasm_unitigs.fasta \ > ./miniasm.racon.consensus.fasta
Analyse the assembly quality by comparing the consensus assembly to the published sequence using dnadiff and Assemblytics (for details see the Assembly using minimap/miniasm tutorial).
dnadiff ~/course_data/precompiled/chr17.fasta -p racon_first \ miniasm.racon.consensus.fasta
Create a directory round_2 in the racon/miniasm directory and change into it. Now repeat the process but this time use the consensus assembly you just created with racon.
mkdir round_2 cd round_2 minimap2 ../miniasm.racon.consensus.fasta \ ../nanofilt_result.fastq > round2.paf racon ../nanofilt_result.fastq round2.paf \ ../miniasm.racon.consensus.fasta \ > miniasm.racon.consensus.round2.fasta
Analyse the second consensus assembly as before.