Skip to content

Commit 9318bd2

Browse files
committed
Fixed storeDir issue
Given a default, single accession run followed by a run with all the accessions, the use of generic file output glob in conjunction with storeDir directive caused an issue with read files being incorrectly classed as present and reads for accession 1 output for every accession, causing pipeline to fail due tu duplicate file names at multiqc step where the files are collected.
1 parent 1c36afe commit 9318bd2

File tree

1 file changed

+2
-2
lines changed

1 file changed

+2
-2
lines changed

main.nf

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -64,14 +64,14 @@ process bgzip_chromosome_subregion {
6464

6565
process extract_reads {
6666
tag { accession }
67-
storeDir { "${params.outdir}/downloaded_reads" } //use with care, caching will not work as normal so changes to input may not take effect
67+
storeDir { "${workDir}/downloaded_reads" } //use with care, caching will not work as normal so changes to input may not take effect
6868

6969
input:
7070
val accession from accessionsChannel
7171
//e.g. ACBarrie
7272

7373
output:
74-
set val(accession), file('*.fastq.gz') into (extractedReadsChannelA, extractedReadsChannelB)
74+
set val(accession), file("${accession}_R?.fastq.gz") into (extractedReadsChannelA, extractedReadsChannelB)
7575
//e.g. ACBarrie, [ACBarrie_R1.fastq.gz, ACBarrie_R2.fastq.gz]
7676

7777
script:

0 commit comments

Comments
 (0)