diff --git "a/src/RNAseq_sen1\342\210\206.nf" "b/src/RNAseq_sen1\342\210\206.nf" index eda4e8a0afb8fcb5336f4acc2db0e1ec8adfe754..7d3fb5ab13f7db7218f690a174acd004a6e814df 100644 --- "a/src/RNAseq_sen1\342\210\206.nf" +++ "b/src/RNAseq_sen1\342\210\206.nf" @@ -36,3 +36,43 @@ process adaptor_removal { ${reads[0]} ${reads[1]} > ${pair_id}_report.txt """ } + +/* +* 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: + set pair_id, file(reads) from fastq_files_cut + + output: + set pair_id, "*_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 ${pair_id}_trim_R1.fastq.gz --outpair ${pair_id}_trim_R2.fastq.gz \ +> ${pair_id}_trimming_report.txt +""" +}