diff --git a/src/RNA_Seq.nf b/src/RNA_Seq.nf index 0eb73bbc91895a705c47b3199818d6a69af38f0a..112097780539633d49102ec6551c04f64b24c847 100644 --- a/src/RNA_Seq.nf +++ b/src/RNA_Seq.nf @@ -37,3 +37,43 @@ process adaptor_removal { """ } +/* +* urqt : +* Imputs : fastq files +* Output : fastq files +*/ +/* quality trimming */ + +/* +* for paired-end data +*/ + +params.fastq = "$baseDir/data/fastq/*_{1,2}.fastq" + +log.info "fastq files : ${params.fastq}" + +Channel + .fromFilePairs( params.fastq ) + .ifEmpty { error "Cannot find any fastq files matching: ${params.fastq}" } + .set { fastq_files } + +process trimming { + tag "${reads}" + cpus 4 + publishDir "results/fastq/trimming/", mode: 'copy' + + input: + file reads from fastq_files + + 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 +""" +} +