Assembly and intrahost/low-frequency variant calling for viral samples
nfcore/viralrecon is a bioinformatics analysis pipeline used to perform assembly and intra-host/low-frequency variant calling for viral samples. The pipeline supports short-read Illumina sequencing data from both shotgun (e.g. sequencing directly from clinical samples) and enrichment-based library preparation methods (e.g. amplicon-based: ARTIC SARS-CoV-2 enrichment protocol; or probe-capture-based).
The pipeline is built using Nextflow, a workflow tool to run tasks across multiple compute infrastructures in a very portable manner. It comes with Docker containers making installation trivial and results highly reproducible. Furthermore, automated continuous integration tests that run the pipeline on a full-sized dataset using AWS cloud ensure that the code is stable.
- Download samples via SRA, ENA or GEO ids (
parallel-fastq-dump; if required)
- Merge re-sequenced FastQ files (
cat; if required)
- Read QC (
- Adapter trimming (
- Variant calling
- Read alignment (
- Sort and index alignments (
- Primer sequence removal (
iVar; amplicon data only)
- Duplicate read marking (
picard; removal optional)
- Alignment-level QC (
- Genome-wide and amplicon coverage QC plots (
- Choice of multiple variant calling and consensus sequence generation routes (
iVar variants and consensus||
- Intersect variants across callers (
- Read alignment (
- De novo assembly
- Primer trimming (
Cutadapt; amplicon data only)
- Removal of host reads (
- Choice of multiple assembly tools (
- Primer trimming (
- Present QC and visualisation for raw read, alignment, assembly and variant calling results (
NB: The pipeline has a number of options to allow you to run only specific aspects of the workflow if you so wish. For example, you can skip all of the assembly steps with the
--skip_assemblyparameter. See the usage docs for all of the available options when running the pipeline.
Numerous QC and reporting steps are included in the pipeline in order to collate a full summary of the analysis within a single MultiQC report. You can see an example MultiQC report here, generated using the parameters defined in this configuration file. The pipeline was run with these samples, prepared from the ncov-2019 ARTIC Network V1 amplicon set and sequenced on the Illumina MiSeq platform in 301bp paired-end format.
Download the pipeline and test it on a minimal dataset with a single command:
Please check nf-core/configs to see if a custom config file to run nf-core pipelines already exists for your Institute. If so, you can simply use
-profile <institute>in your command. This will enable either
singularityand set the appropriate execution settings for your local compute environment.
Start running your own analysis!
Typical command for shotgun analysis:
Typical command for amplicon analysis:
See the usage documentation for all of the available options when running the pipeline.
The nf-core/viralrecon pipeline comes with documentation about the pipeline, found in the
- Pipeline configuration
- Running the pipeline
- Output and how to interpret the results
These scripts were originally written by Sarai Varona, Miguel Juliá and Sara Monzon from BU-ISCIII and co-ordinated by Isabel Cuesta for the Institute of Health Carlos III, Spain. Through collaboration with the nf-core community the pipeline has now been updated substantially to include additional processing steps, to standardise inputs/outputs and to improve pipeline reporting; implemented primarily by Harshil Patel from The Bioinformatics & Biostatistics Group at The Francis Crick Institute, London.
Many thanks to others who have helped out and contributed along the way too, including (but not limited to)*:
* Listed in alphabetical order
Contributions and Support
If you would like to contribute to this pipeline, please see the contributing guidelines.
If you use nf-core/viralrecon for your analysis, please cite it using the following doi: 10.5281/zenodo.3901628
An extensive list of references for the tools used by the pipeline can be found in the
You can cite the
nf-core publication as follows:
The nf-core framework for community-curated bioinformatics pipelines.
Philip Ewels, Alexander Peltzer, Sven Fillinger, Harshil Patel, Johannes Alneberg, Andreas Wilm, Maxime Ulysse Garcia, Paolo Di Tommaso & Sven Nahnsen.