Skip to content
Snippets Groups Projects
Verified Commit c82cceba authored by Laurent Modolo's avatar Laurent Modolo
Browse files

UrQt: nf add file_id var for single-end like for paired-end

parent 7ed402cc
No related branches found
No related tags found
No related merge requests found
...@@ -5,24 +5,25 @@ log.info "fastq files : ${params.fastq}" ...@@ -5,24 +5,25 @@ log.info "fastq files : ${params.fastq}"
Channel Channel
.fromPath( params.fastq ) .fromPath( params.fastq )
.ifEmpty { error "Cannot find any fastq files matching: ${params.fastq}" } .ifEmpty { error "Cannot find any fastq files matching: ${params.fastq}" }
.map { it -> [(it.baseName =~ /([^\.]*)/)[0][1], it]}
.set { fastq_files } .set { fastq_files }
process trimming { process trimming {
tag "$reads.baseName" tag "$file_id"
cpus 4 cpus 4
input: input:
file reads from fastq_files set file_id, file(reads) from fastq_files
output: output:
file "*_trim.fastq.gz" into fastq_files_trim set file_id, "*_trim.fastq.gz" into fastq_files_trim
script: script:
""" """
UrQt --t 20 --m ${task.cpus} --gz \ UrQt --t 20 --m ${task.cpus} --gz \
--in ${reads} \ --in ${reads} \
--out ${reads.baseName}_trim.fastq.gz \ --out ${file_id}_trim.fastq.gz \
> ${reads.baseName}_trimming_report.txt > ${file_id}_trimming_report.txt
""" """
} }
...@@ -49,26 +49,24 @@ log.info "fastq files : ${params.fastq}" ...@@ -49,26 +49,24 @@ log.info "fastq files : ${params.fastq}"
Channel Channel
.fromPath( params.fastq ) .fromPath( params.fastq )
.ifEmpty { error "Cannot find any fastq files matching: ${params.fastq}" } .ifEmpty { error "Cannot find any fastq files matching: ${params.fastq}" }
.map { it -> [(it.baseName =~ /([^\.]*)/)[0][1], it]}
.set { fastq_files } .set { fastq_files }
process trimming { process trimming {
tag "$reads.baseName" tag "$file_id"
cpus 4 cpus 4
publishDir "results/fastq/trimming/", mode: 'copy'
input: input:
file reads from fastq_files set file_id, file(reads) from fastq_files
output: output:
file "*_trim.fastq.gz" into fastq_files_trim set file_id, "*_trim.fastq.gz" into fastq_files_trim
script: script:
""" """
UrQt --t 20 --m ${task.cpus} --gz \
UrQt --t 20 --m ${task.cpus} --gz \ --in ${reads} \
--in ${reads} \ --out ${file_id}_trim.fastq.gz \
--out ${reads.baseName}_trim.fastq.gz \ > ${file_id}_trimming_report.txt
> ${reads.baseName}_trimming_report.txt """
"""
} }
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment