@@ -19,13 +19,15 @@ Original pipelining and code forked from the CCBR Exome-seek Pipeline [Exome-see
19
19
[ singularity] ( https://singularity.lbl.gov/all-releases ) must be installed on the target system. Snakemake orchestrates the execution of each step in the pipeline. To guarantee the highest level of reproducibility, each step relies on versioned images from [ DockerHub] ( https://hub.docker.com/orgs/nciccbr/repositories ) . Nextflow uses singularity to pull these images onto the local filesystem prior to job execution, and as so, nextflow and singularity are the only two dependencies.
20
20
21
21
## Setup
22
- LOGAN is installed on the Biowulf in the ccbrpipeliner module.
23
- Please clone this repository to your local filesystem using the following command:
22
+ LOGAN can be used with the Nextflow pipelining software
23
+ Please clone this repository to your local filesystem using the following command on Biowulf :
24
24
``` bash
25
25
# start an interactive node
26
26
sinteractive --mem=2g --cpus-per-task=2 --gres=lscratch:200
27
- # load the ccbrpipeliener module
28
- module load ccbrpipeliner
27
+ git clone https://github.com/CCBR/LOGAN
28
+ module load nextflow
29
+ # #Example run
30
+ nextflow run /data/LOGAN//main.nf
29
31
```
30
32
31
33
## Usage
@@ -89,17 +91,17 @@ c130889189_TUMOR c130889189_PBMC
89
91
90
92
#### 2. Tumor only mode
91
93
92
- No flags are required
94
+ No addtional flags for sample sheet are required as all samples will be used to call variants
93
95
94
96
#### Calling Mode
95
97
96
98
Adding flags determines SNV (germline and/or somatic), SV, and/or CNV calling modes
97
99
98
100
` --vc ` - Enables somatic SNV calling using mutect2, vardict, varscan, octopus, strelka (TN only), MUSE (TN only), and lofreq (TN only)
99
101
100
- ` --germline ` - Enables germline using DV
102
+ ` --germline ` - Enables germline using Deepvariant
101
103
102
- ` --sv ` - Enables somatic SV calling using Manta and SVABA
104
+ ` --sv ` - Enables somatic SV calling using Manta, SVABA, and GRIDSS (coming soon)
103
105
104
106
` --cnv ` - Enables somatic CNV calling using FREEC, Sequenza, and Purple (hg19/hg38 only)
105
107
@@ -111,33 +113,32 @@ Adding flags determines SNV (germline and/or somatic), SV, and/or CNV calling mo
111
113
` --callers ` - Comma separated argument for callers, the default is to use all available.
112
114
Example: ` --callers mutect2,octopus `
113
115
114
- ` --cnvcallers ` - - Comma separated argument for cnvcallers . Adding flag allows only certain callers to run.
116
+ ` --cnvcallers ` - - Comma separated argument for CNV callers to use . Adding flag allows only certain callers to run.
115
117
Example: ` --cnvcallers purple `
116
118
119
+ ` --svcallers ` - - Comma separated argument for SV callers. Adding flag allows only certain callers to run.
120
+ Example: ` --cnvcallers manta `
121
+
117
122
118
123
## Running LOGAN
119
124
Example of Tumor_Normal calling mode
120
125
``` bash
121
- # copy the logan config files to your current directory
122
- logan init
123
126
# preview the logan jobs that will run
124
- logan run --mode local -profile ci_stub --genome hg38 --sample_sheet samplesheet.tsv --outdir out --fastq_input " *R{1,2}.fastq.gz" -preview --vc --sv --cnv
127
+ nextflow run /data/LOGAN/main.nf --mode local -profile ci_stub --genome hg38 --sample_sheet samplesheet.tsv --outdir out --fastq_input " *R{1,2}.fastq.gz" -preview --vc --sv --cnv
125
128
# run a stub/dryrun of the logan jobs
126
- logan run --mode local -profile ci_stub --genome hg38 --sample_sheet samplesheet.tsv --outdir out --fastq_input " *R{1,2}.fastq.gz" -stub --vc --sv --cnv
129
+ nextflow run /data/LOGAN/main.nf --mode local -profile ci_stub --genome hg38 --sample_sheet samplesheet.tsv --outdir out --fastq_input " *R{1,2}.fastq.gz" -stub --vc --sv --cnv
127
130
# launch a logan run on slurm with the test dataset
128
- logan run --mode slurm -profile biowulf,slurm --genome hg38 --sample_sheet samplesheet.tsv --outdir out --fastq_input " *R{1,2}.fastq.gz" --vc --sv --cnv
131
+ nextflow run /data/LOGAN/main.nf --mode slurm -profile biowulf,slurm --genome hg38 --sample_sheet samplesheet.tsv --outdir out --fastq_input " *R{1,2}.fastq.gz" --vc --sv --cnv
129
132
```
130
133
131
134
Example of Tumor only calling mode
132
135
``` bash
133
- # copy the logan config files to your current directory
134
- logan init
135
136
# preview the logan jobs that will run
136
- logan run --mode local -profile ci_stub --genome hg38 --outdir out --fastq_input " *R{1,2}.fastq.gz" --callers octopus,mutect2 -preview --vc --sv --cnv
137
+ nextflow run /data/LOGAN/main.nf --mode local -profile ci_stub --genome hg38 --outdir out --fastq_input " *R{1,2}.fastq.gz" --callers octopus,mutect2 -preview --vc --sv --cnv
137
138
# run a stub/dryrun of the logan jobs
138
- logan run --mode local -profile ci_stub --genome hg38 --outdir out --fastq_input " *R{1,2}.fastq.gz" --callers octopus,mutect2 -stub --vc --sv --cnv
139
+ nextflow run /data/LOGAN/main.nf --mode local -profile ci_stub --genome hg38 --outdir out --fastq_input " *R{1,2}.fastq.gz" --callers octopus,mutect2 -stub --vc --sv --cnv
139
140
# launch a logan run on slurm with the test dataset
140
- logan run --mode slurm -profile biowulf,slurm --genome hg38 --outdir out --fastq_input " *R{1,2}.fastq.gz" --callers octopus,mutect2 --vc --sv --cnv
141
+ nextflow run /data/LOGAN/main.nf --mode slurm -profile biowulf,slurm --genome hg38 --outdir out --fastq_input " *R{1,2}.fastq.gz" --callers octopus,mutect2 --vc --sv --cnv
141
142
```
142
143
143
144
0 commit comments