nf-core/hicar
Pipeline for HiCAR data, a robust and sensitive multi-omic co-assay for simultaneous measurement of transcriptome, chromatin accessibility and cis-regulatory chromatin contacts.
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.
Metho for the experiment.
string
HiCAR
It can be HiCAR, HiChIP, ChIA-PET, PLAC-Seq. If it is not HiCAR, user may want to define the anchor peaks. Otherwise, the anchor peaks will be called by R2 reads.
Path to anchor peaks
string
^\S+\.(narrowPeak|boradPeak)$
The anchor peaks can be a narrowPeak or broadPeak.
The output directory where the results will be saved. You have to use absolute paths to storage on Cloud infrastructure.
string
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.
MultiQC report title. Printed as page header, used for filename if not otherwise specified.
string
Reference genome related files and options required for the workflow.
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 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
.
Path to bwa index file.
string
If you have no genome reference available, the pipeline can build one using a FASTA file. This requires additional time and resources, so it's better to use a pre-build index if possible.
Path to annotation gtf file.
string
Annotation gtf or gff file is required for annotation and quality control. This parameter is mandatory if --genome
is not specified.
Path to annotation gff file.
string
Annotation gff or gff file is required for annotation and quality control. This parameter is mandatory if --genome
is not specified.
Path to annotation gene bed file.
string
Path to BED file containing gene intervals. This will be created from the GTF file if not specified.
Path to genome mappability file.
string
If you have no genome mappability available, the pipeline can build one using a FASTA file. This requires additional time and resources, so it's better to use a pre-build index if possible.
Effective genome size parameter required by MACS2.
string
Effective genome size 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. For other genomes, if this parameter is not specified then the MACS2 peak-calling and differential analysis will be skipped.
Path to blacklist regions in BED format, used for filtering alignments.
string
If provided, alignments that overlap with the regions in this file will be filtered out (see ENCODE blacklists). The file should be in BED format. Blacklisted regions for GRCh37, GRCh38, GRCm38, hg19, hg38, mm10 are bundled with the pipeline in the blacklists
directory, and as such will be automatically used if any of those genomes are specified with the --genome
parameter.
Parameters used to describe the experiment designs.
Specifies that the cutting position has to be using.
string
CviQI
Default CviQI digestion. Available enzymes are MboI, DpnII, BglII, HindIII, and CviQI.
Specifies that the cutting sequence has to be using.
string
^TAC
It will be automatically assigned by the enzyme inputed.
Specifies that the cutoff value used for mappability filter.
number
0.5
Default is 0.5. It is the average over just covered bases.
shift size for MACS2
integer
-75
How reads will be shift to increase the peak/noise ratios. See help from MACS2 documentation
extsize for MACS2
integer
150
How reads will be extend in 5'->3' direction to increase the fragment size. See help from MACS2 documentation
cutoff qvalue
number
0.01
The minimum FDR cutoff to call significant peak for R2 reads. See help from MACS2 documentation
read peak calling for fragment (R1) reads by MACS2
boolean
cutoff pvalue for fragment (R1)
number
0.1
The minimum FDR cutoff to call significant peak for R1 reads. R1 reads will be more noisy compare to R2 reads. Here we need a loose cutoff. And this will be fixed in the following steps. See help from MACS2 documentation
resolution bin size
string
5000_10000
Bin size for singal resampling.
output of restriction_cut_multipleenzyme.py.
string
None
See https://github.com/ijuric/MAPS/blob/master/bin/utils/genomic_features_generator/scripts/restriction_cut_multipleenzyme.py.
MAPS regression cutoff value
integer
12
Minimal cutoff value for coverage in the called connections. See help from MAPS documentation, default 12.
MAPS regression fold change cutoff value
number
2
Miniaml cutoff value for fold change from signal to expected values. See help from MAPS documentation, default 2.
MAPS regression -log10(fdr) cutoff value
number
2
Minimal cutoff -log10(FDR) value in the called connections. See help from MAPS documentation, default 2 (=fdr 0.05).
MAPS regression filter file name
string
None
file containing bins that need to be filtered out. Format: two columns 'chrom', 'bin'. 'chrom' contains 'chr1','chr2',.. 'bin' is bin label
MAPS regression type
string
pospoisson for positive poisson regression, negbinom for negative binomial. default is pospoisson.
source code path for merge_map.py
string
https://raw.githubusercontent.com/ijuric/MAPS/91c9c360092b25a217d91b9ea07eba5dd2ac72f4/bin/utils/genomic_features_generator/scripts/merge_map.py
The source code is available at https://github.com/ijuric/MAPS/blob/master/bin/utils/genomic_features_generator/scripts/merge_map.py
feature_frag2bin source path
string
https://raw.githubusercontent.com/ijuric/MAPS/91c9c360092b25a217d91b9ea07eba5dd2ac72f4/bin/utils/genomic_features_generator/scripts/feature_frag2bin.py
The source code is available at https://github.com/ijuric/MAPS/blob/master/bin/utils/genomic_features_generator/scripts/feature_frag2bin.py
make_maps_runfile source path
string
https://raw.githubusercontent.com/ijuric/MAPS/18c1a337f222130d7c5735d051614e2a253d5319/bin/MAPS/make_maps_runfile.py
The source code is available at https://github.com/ijuric/MAPS/blob/master/bin/MAPS/make_maps_runfile.py
remove duplicates for high resolution peaks or not
boolean
Type of snow cluster to use
string
SOCK
Possible values are 'SOCK' (default), 'MPI' and 'FORK'. This will be used by snow package for HiPeak calling
The block number of peak pair
number
1000000000
The block number of peak pair in each thread to count the reads. Smaller number will decrease the memory cost and increase the computation time for parepare the counts in hipeak calling.
Options for tracks viewed by igv, ucsc genome browser and virtual 4c plot tools
create track files for virtual 4c or not
boolean
max events to plot for virtual 4c
integer
25
JVM heap parameters for juicer
string
-Xms512m -Xmx4096m
For 64bit system, you may want to set the Xmx > 4g, eg. '-Xms1024m -Xmx32g'
The juicer_tools path
string
https://s3.amazonaws.com/hicfiles.tc4ga.com/public/juicer/juicer_tools_1.22.01.jar
The juicer_tools can be downloaded from juicer_tools download page
skip peak QC or not
boolean
skip plot profile or not
boolean
skip creat IGV files or not
boolean
skip create trackhub files or not
boolean
cutoff value for false discovery rate of enrichment analysis
number
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
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
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
Run this workflow with Conda. You can also use '-profile conda' instead of providing this parameter.
boolean
To skip some processes or workflows.
skip trim 5'end TAC
boolean
skip fastqc or not
boolean
skip peak annotation or not
boolean
skip enrichment or not
boolean
true
skip differential analysis or not
boolean
skip multiqc or not
boolean