diff --git a/src/Rnaseq.nf b/src/Rnaseq.nf index 303cdb6cb036c2dfc41d439697a597ef2dc8ffad..851cc91c0b25c5078e1d412b3cdcb8dee5eb3012 100644 --- a/src/Rnaseq.nf +++ b/src/Rnaseq.nf @@ -24,3 +24,24 @@ process adaptor_removal { ${reads[0]} ${reads[1]} > ${pair_id}_report.txt """ } + +process trimming { + tag "${reads}" + cpus 4 + publishDir "results/fastq/trimming/", mode: 'copy' + + input: + set pair_id, file(reads) from fastq_files_cut + + output: + file "*_trim_R{1,2}.fastq.gz" into fastq_files_trim + + script: +""" +UrQt --t 20 --m ${task.cpus} --gz \ +--in ${reads[0]} --inpair ${reads[1]} \ +--out ${reads[0].baseName}_trim_R1.fastq.gz --outpair ${reads[1].baseName}_trim_R2.fastq.gz \ +> ${reads[0].baseName}_trimming_report.txt +""" +} +