RNA_Seq.nf: adding urqt block

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
+"""
+}
+