nf-core/cutandrun
Analysis pipeline for CUT&RUN and CUT&TAG experiments that includes QC, support for spike-ins, IgG controls, peak calling and downstream analysis.
3.1
). The latest
stable release is
3.2.2
.
Define where the pipeline should find input data and save output data.
Path to comma-separated file containing information about the samples in the experiment.
string
^\S+\.csv$
You will need to create a design file with information about the samples in your experiment before running the pipeline. Use this parameter to specify its location. It has to be a comma-separated file with 3 columns, and a header row. See usage docs.
The output directory where the results will be saved. You have to use absolute paths to store on Cloud infrastructure.
string
MultiQC report title. Printed as page header, used for filename if not otherwise specified.
string
Save genome reference data to the output directory
boolean
Save any technical replicate FASTQ files that were merged to the output directory
boolean
Save trimmed FASTQ files to the output directory
boolean
Save BAM files aligned to the spike-in genome to the output directory
boolean
Save unaligned sequences to the output directory
boolean
Save alignment intermediates to the output directory (WARNING: can be very large)
boolean
Email address for completion summary.
string
^([a-zA-Z0-9_\-\.]+)@([a-zA-Z0-9_\-\.]+)\.([a-zA-Z]{2,5})$
Set this parameter to your e-mail address to get a summary e-mail with details of the run sent to you when the workflow exits. If set in your user config file (~/.nextflow/config
) then you don't need to specify this on the command line for every run.
Reference genome related files and options.
Name of iGenomes reference.
string
If using a reference genome configured in the pipeline using iGenomes, use this parameter to give the ID for the reference. This is then used to build the full paths for all required reference genome files e.g. --genome GRCh38
.
See the nf-core website docs for more details.
Path to bowtie2 index
string
Path to GTF annotation file
string
Path to gene BED file
string
Path to genome blacklist
string
Name of the igenome reference for the spike-in genome
string
K12-MG1655
Path to spike-in bowtie2 index
string
Path to spike-in fasta
string
Path to FASTA genome file.
string
^\S+\.fn?a(sta)?(\.gz)?$
This parameter is mandatory if --genome
is not specified. If you don't have a BWA index available this will be generated for you automatically. Combine with --save_reference
to save BWA index for future runs.
Directory / URL base for iGenomes references.
string
s3://ngi-igenomes/igenomes
Do not load the iGenomes reference config.
boolean
Do not load igenomes.config
when running the pipeline. You may choose this option if you observe clashes between custom parameters and those supplied in igenomes.config
.
Run pipeline up to input checking
boolean
Run pipeline up to reference preparation
boolean
Run pipeline up to pre-alignment
boolean
Run pipeline up to alignment
boolean
Run pipeline up to q-filtering
boolean
Run pipeline up to peak calling
boolean
Skips fastqc reporting
boolean
Skips trimming
boolean
Skips de-duplication
boolean
Skips reporting
boolean
Skips preseq reporting
boolean
Skips igv session generation
boolean
Skips deeptools QC repoting
boolean
Skips deeptools heatmap generation
boolean
Skips peak QC reporting
boolean
Skips multiqc
boolean
Instructs Trim Galore to remove bp from the 5' end of read 1 (or single-end reads).
integer
Instructs Trim Galore to remove bp from the 5' end of read 2 (paired-end reads only).
integer
Instructs Trim Galore to remove bp from the 3' end of read 1 AFTER adapter/quality trimming has been performed.
integer
Instructs Trim Galore to remove bp from the 3' end of read 2 AFTER adapter/quality trimming has been performed.
integer
Instructs Trim Galore to apply the --nextseq=X option, to trim based on quality after removing poly-G tails.
integer
This enables the option Cutadapt --nextseq-trim=3'CUTOFF
option via Trim Galore, which will set a quality cutoff (that is normally given with -q instead), but qualities of G bases are ignored. This trimming is in common for the NextSeq- and NovaSeq-platforms, where basecalls without any signal are called as high-quality G bases.
Select aligner
string
bowtie2
Normalisation constant for spike-in read normalisation
integer
10000
Filter reads below a q-score threshold
integer
De-duplicate target reads AND control reads (default is control only)
boolean
Sets the target read normalisation mode. Options are: ["Spikein", "RPKM", "CPM", "BPM", "None" ]
string
If normsalisation option is one of "RPKM", "CPM", "BPM" - then the binsize that the reads count is calculated on is used.
integer
1
Selects the peak caller for the pipeline. Options are: [seacr, macs2]. More than one peak caller can be chosen and the order specifies which is a primary peak called (the first) that will be used downstream. Any secondary peak callers will be run and outputed to the results folder.
string
seacr
Specifies whether to use a control to normalise peak calls against (e.g. IgG)
boolean
true
Specifies whether to extend paired-end fragments between the read mates when calculating coveage tracks
boolean
true
Specifies whether the background control is scaled prior to being used to normalise peaks.
number
0.5
SEACR specifies returns the top n fraction (between 0 and 1) of peaks based on total signal within peaks. This is only used if there are no controls included with the samples and if --use_control
is false
number
0.05
SEACR normalization.
string
SEACR stringency.
string
Q-value threshold for MACS2 peak caller.
number
0.01
P-value threshold for macs2 peak caller. If set it will overide the qvalue.
number
parameter required by MACS2. If using an iGenomes reference these have been provided when --genome
is set as GRCh37, GRCh38, GRCm38, WBcel235, BDGP6, R64-1-1, EF2, hg38, hg19 and mm10. Otherwise the gsize will default to GRCh38.
number
2700000000
Determines whether MACS2 broad or narrow peak mode is used for the peak caller
boolean
true
MACS2 broad cutoff parameter
number
0.1
Specifies what samples to group together for consensus peaks. Options are [group, all]
string
Minimum number of overlapping replicates needed for a consensus peak
number
1
Show gene names instead of symbols in IGV browser sessions
boolean
true
Deeptools multiBamSummary bam bin size
integer
500
Deeptools Correlation Plot statistical calculation method
string
pearson
Deeptools heatmap gene plot before length (bases)
integer
3000
Deeptools heatmap gene plot body length (bases)
integer
5000
Deeptools heatmap gene plot after length (bases)
integer
3000
Deeptools heatmap peak plot before length (bases)
integer
3000
Deeptools heatmap peak plot after length (bases)
integer
3000
Flag for whether to generate a heatmap for all samples together
boolean
true
Minimum fragment overlap for FriP score
number
0.2
Minimum peak overlap for peak reproducibility plot
number
0.2
Parameters used to describe centralised config profiles. These should not be edited.
Git commit id for Institutional configs.
string
master
Base directory for Institutional configs.
string
https://raw.githubusercontent.com/nf-core/configs/master
If you're running offline, Nextflow will not be able to fetch the institutional config files from the internet. If you don't need them, then this is not a problem. If you do need them, you should download the files from the repo and tell Nextflow where to find them with this parameter.
Institutional config name.
string
Institutional config description.
string
Institutional config contact information.
string
Institutional config URL link.
string
Set the top limit for requested resources for any single job.
Maximum number of CPUs that can be requested for any single job.
integer
16
Use to set an upper-limit for the CPU requirement for each process. Should be an integer e.g. --max_cpus 1
Maximum amount of memory that can be requested for any single job.
string
128.GB
^\d+(\.\d+)?\.?\s*(K|M|G|T)?B$
Use to set an upper-limit for the memory requirement for each process. Should be a string in the format integer-unit e.g. --max_memory '8.GB'
Maximum amount of time that can be requested for any single job.
string
240.h
^(\d+\.?\s*(s|m|h|day)\s*)+$
Use to set an upper-limit for the time requirement for each process. Should be a string in the format integer-unit e.g. --max_time '2.h'
Less common options for the pipeline, typically set in a config file.
Display help text.
boolean
Display version and exit.
boolean
Method used to save pipeline results to output directory.
string
The Nextflow publishDir
option specifies which intermediate files should be saved to the output directory. This option tells the pipeline what method should be used to move these files. See Nextflow docs for details.
Email address for completion summary, only when pipeline fails.
string
^([a-zA-Z0-9_\-\.]+)@([a-zA-Z0-9_\-\.]+)\.([a-zA-Z]{2,5})$
An email address to send a summary email to when the pipeline is completed - ONLY sent if the pipeline does not exit successfully.
Send plain-text email instead of HTML.
boolean
File size limit when attaching MultiQC reports to summary emails.
string
25.MB
^\d+(\.\d+)?\.?\s*(K|M|G|T)?B$
Do not use coloured log outputs.
boolean
Incoming hook URL for messaging service
string
Incoming hook URL for messaging service. Currently, MS Teams and Slack are supported.
Custom config file to supply to MultiQC.
string
Directory to keep pipeline Nextflow logs and reports.
string
${params.outdir}/pipeline_info
Boolean whether to validate parameters against the schema at runtime
boolean
true
Show all params when using --help
boolean