Skip to content
Snippets Groups Projects
Commit 322c97ac authored by nservant's avatar nservant
Browse files

update test-dataset

parent 4ad259c9
No related branches found
No related tags found
No related merge requests found
...@@ -13,15 +13,15 @@ params { ...@@ -13,15 +13,15 @@ params {
splitFastq = false splitFastq = false
bwt2_opts_end2end = '--very-sensitive -L 30 --score-min L,-0.6,-0.2 --end-to-end --reorder' bwt2_opts_end2end = '--very-sensitive -L 30 --score-min L,-0.6,-0.2 --end-to-end --reorder'
bwt2_opts_trimmed = '--very-sensitive -L 20 --score-min L,-0.6,-0.2 --end-to-end --reorder' bwt2_opts_trimmed = '--very-sensitive -L 20 --score-min L,-0.6,-0.2 --end-to-end --reorder'
min_mapq = 0 min_mapq = 0
// Digestion Hi-C // Digestion Hi-C
restriction_site = 'A^AGCTT' restriction_site = 'A^AGCTT'
ligation_site = 'AAGCTAGCTT' ligation_site = 'AAGCTAGCTT'
min_restriction_fragment_size = 100 min_restriction_fragment_size =
max_restriction_fragment_size = 100000 max_restriction_fragment_size =
min_insert_size = 100 min_insert_size =
max_insert_size = 600 max_insert_size =
// Hi-C Processing // Hi-C Processing
min_cis_dist = min_cis_dist =
......
...@@ -9,7 +9,7 @@ ...@@ -9,7 +9,7 @@
params { params {
config_profile_name = 'Hi-C test data from Dixon et al. (2012)' config_profile_name = 'Hi-C test data from Schalbetter et al. (2017)'
config_profile_description = 'Minimal test dataset to check pipeline function' config_profile_description = 'Minimal test dataset to check pipeline function'
// Limit resources so that this can run on Travis // Limit resources so that this can run on Travis
...@@ -19,10 +19,13 @@ params { ...@@ -19,10 +19,13 @@ params {
// Input data // Input data
readPaths = [ readPaths = [
['SRR400264_00', ['https://github.com/nf-core/test-datasets/raw/hic/SRR400264_00_R1.fastq.gz', 'https://github.com/nf-core/test-datasets/raw/hic/SRR400264_00_R2.fastq.gz']], ['SRR4292758_00', ['https://github.com/nf-core/test-datasets/raw/hic/data/SRR4292758_00_R1.fastq.gz', 'https://github.com/nf-core/test-datasets/raw/hic/data/SRR4292758_00_R2.fastq.gz']]
['SRR400264_01', ['https://github.com/nf-core/test-datasets/raw/hic/SRR400264_01_R1.fastq.gz', 'https://github.com/nf-core/test-datasets/raw/hic/SRR400264_01_R2.fastq.gz']] ]
]
// Annotations // Annotations
genome = 'GRCh37' fasta = 'https://github.com/nf-core/test-datasets/raw/hic/reference/W303_SGD_2015_JRIU00000000.fsa'
// Options
skip_cool = true
} }
...@@ -33,11 +33,11 @@ def helpMessage() { ...@@ -33,11 +33,11 @@ def helpMessage() {
Mandatory arguments: Mandatory arguments:
--reads Path to input data (must be surrounded with quotes) --reads Path to input data (must be surrounded with quotes)
--genome Name of iGenomes reference
-profile Configuration profile to use. Can use multiple (comma separated) -profile Configuration profile to use. Can use multiple (comma separated)
Available: conda, docker, singularity, awsbatch, test and more. Available: conda, docker, singularity, awsbatch, test and more.
References If not specified in the configuration file or you wish to overwrite any of the references. References If not specified in the configuration file or you wish to overwrite any of the references.
--genome Name of iGenomes reference
--bwt2_index Path to Bowtie2 index --bwt2_index Path to Bowtie2 index
--fasta Path to Fasta reference --fasta Path to Fasta reference
--chromosome_size Path to chromosome size file --chromosome_size Path to chromosome size file
...@@ -72,6 +72,10 @@ def helpMessage() { ...@@ -72,6 +72,10 @@ def helpMessage() {
--email 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 --email 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
-name Name for the pipeline run. If not specified, Nextflow will automatically generate a random mnemonic. -name Name for the pipeline run. If not specified, Nextflow will automatically generate a random mnemonic.
Step options:
--skip_cool Skip generation of cool files
--skip_multiQC Skip MultiQC
AWSBatch options: AWSBatch options:
--awsqueue The AWSBatch JobQueue that needs to be set when running on AWSBatch --awsqueue The AWSBatch JobQueue that needs to be set when running on AWSBatch
--awsregion The AWS Region for your AWS Batch job to run on --awsregion The AWS Region for your AWS Batch job to run on
...@@ -138,8 +142,7 @@ if (params.readPaths){ ...@@ -138,8 +142,7 @@ if (params.readPaths){
.from( params.readPaths ) .from( params.readPaths )
.map { row -> [ row[0], [file(row[1][0]), file(row[1][1])]] } .map { row -> [ row[0], [file(row[1][0]), file(row[1][1])]] }
.separate( raw_reads, raw_reads_2 ) { a -> [tuple(a[0], a[1][0]), tuple(a[0], a[1][1])] } .separate( raw_reads, raw_reads_2 ) { a -> [tuple(a[0], a[1][0]), tuple(a[0], a[1][1])] }
.println() }else{
}else{
raw_reads = Channel.create() raw_reads = Channel.create()
raw_reads_2 = Channel.create() raw_reads_2 = Channel.create()
...@@ -674,6 +677,9 @@ process generate_cool{ ...@@ -674,6 +677,9 @@ process generate_cool{
tag "$sample" tag "$sample"
publishDir "${params.outdir}/export/cool", mode: 'copy' publishDir "${params.outdir}/export/cool", mode: 'copy'
when:
!params.skip_cool
input: input:
set val(sample), file(vpairs) from all_valid_pairs_4cool set val(sample), file(vpairs) from all_valid_pairs_4cool
file chrsize from chromosome_size_cool.collect() file chrsize from chromosome_size_cool.collect()
...@@ -694,15 +700,18 @@ process generate_cool{ ...@@ -694,15 +700,18 @@ process generate_cool{
process multiqc { process multiqc {
publishDir "${params.outdir}/MultiQC", mode: 'copy' publishDir "${params.outdir}/MultiQC", mode: 'copy'
when:
!params.skip_multiqc
input: input:
file multiqc_config from ch_multiqc_config file multiqc_config from ch_multiqc_config
file ('input_*/*') from all_mstats.concat(all_mergestat).collect() file ('input_*/*') from all_mstats.concat(all_mergestat).collect()
file ('software_versions/*') from software_versions_yaml file ('software_versions/*') from software_versions_yaml
file workflow_summary from create_workflow_summary(summary) file workflow_summary from create_workflow_summary(summary)
output: output:
file "*multiqc_report.html" into multiqc_report file "*multiqc_report.html" into multiqc_report
file "*_data" file "*_data"
script: script:
rtitle = custom_runName ? "--title \"$custom_runName\"" : '' rtitle = custom_runName ? "--title \"$custom_runName\"" : ''
......
...@@ -24,7 +24,8 @@ params { ...@@ -24,7 +24,8 @@ params {
readPaths = false readPaths = false
chromosome_size = false chromosome_size = false
restriction_fragments = false restriction_fragments = false
skip_cool = false
skip_multiqc = false
// Boilerplate options // Boilerplate options
name = false name = false
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment