A workflow for mapping and consensus generation of SARS-CoV2 whole genome amplicon nanopore data implemented in the Nextflow framework. Reads are mapped to a reference genome using Minimap2 after trimming the amplicon primers with a fixed length at both ends of the amplicons using Cutadapt. The consensus is called using Pysam based on a majority read support threshold per position of the Minimap2 alignment and positions with less than 30x coverage are masked using ‘N’.
Creator: David F. Nieuwenhuijse, Alexey Sokolov
Submitter: Ross Thorne
Galaxy version of pre-processing of reads from COVID-19 samples.
QC + human read cleaning
Submitter: ignacio eguinoa