This repository contains a Nextflow variant calling pipeline for analyzing Next-Generation Sequencing (NGS) data using deepvariant.
%%{init: {'theme':'base'}}%%
flowchart TB
subgraph "input"
v0["trimmed reads"]
v1["reference"]
v30["indices"]
v31["vep_DB"]
end
subgraph "output"
v18["vcf"]
v25["html"]
end
v8([ALIGN])
v9([FLAGSTAT])
v10([QUALIMAP])
v13([BAMINDEX])
v14([VARCALL])
v17([ANNOTATE])
v24([REPORT])
v21(( ))
v22(( ))
v23(( ))
v0 --> v8
v1 --> v8
v30 --> v8
v8 --> v9
v8 --> v10
v8 --> v13
v8 --> v14
v9 --> v21
v10 --> v22
v13 --> v14
v1 --> v14
v14 --> v17
v17 --> v18
v17 --> v23
v31 --> v14
v30 --> v14
v21 --> v24
v22 --> v24
v23 --> v24
v24 --> v25
The pipeline is implemented in Nextflow and includes several stages for NGS data analysis:
- ALIGN: Sequence alignment using BWA mem.
- FLAGSTAT: Alignment quality control using Samtools flagstat.
- QUALIMAP: Alignment quality control using Qualimap.
- FAINDEX: Fai index creation using Samtools faidx.
- BAMINDEX: Bai index creation using Samtools index.
- PREPARE: File processing and preparation using Samtools.
- VARCALL: Variant calling using deepvariant.
- ANNOTATE: Annotation using VEP (Variant Effect Predictor).
- REPORT: Compiling report using MultiQC.
To quickly run the pipeline, use the following command:
nextflow run glebus-sasha/deepvariant_nf \
-profile <docker/singularity> \
--reference <path-to-reference> \
--reads <path-to-reads-folder> \
--outdir results- Nextflow (https://www.nextflow.io/docs/latest/install.html)
- Docker (https://docs.docker.com/engine/install/) or
- Singularity (https://github.com/sylabs/singularity/blob/main/INSTALL.md)
- Install all the necessary dependencies such as Nextflow, Singularity.
- Clone this repository:
git clone https://github.com/glebus-sasha/deepvariant_nf.git - Navigate to the pipeline directory:
cd deepvariant_nf - Edit the
nextflow.configfile to set the required parameters, if necessary. - Run the pipeline, setting the required parameters, for example:
nextflow run main.nfThis project is licensed under the MIT License.