diff --git a/src/RNASeq.nf b/src/RNASeq.nf index e8f6bb458804bce7ab9237df604a1455e8836db5..c130936426fd3e557cd15672b24faff21c3de5f9 100644 --- a/src/RNASeq.nf +++ b/src/RNASeq.nf @@ -38,66 +38,3 @@ process adaptor_removal { """ } -/* 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 "$pair_id" - publishDir "results/fastq/trimming/", mode: 'copy' - - input: - set pair_id, file(reads) from fastq_files - - output: - file "*_trim_R{1,2}.fastq.gz" into fastq_files_trim - - script: - """ - cutadapt -q 20,20 \ - -o ${pair_id}_trim_R1.fastq.gz -p ${pair_id}_trim_R2.fastq.gz \ - ${reads[0]} ${reads[1]} > ${pair_id}_report.txt - """ -} - -/* -* for single-end data -*/ - -params.fastq = "$baseDir/data/fastq/*.fastq" - -log.info "fastq files : ${params.fastq}" - -Channel - .fromPath( params.fastq ) - .ifEmpty { error "Cannot find any fastq files matching: ${params.fastq}" } - .set { fastq_files } - -process trimming { - tag "$reads.baseName" - publishDir "results/fastq/trimming/", mode: 'copy' - - input: - file reads from fastq_files - - output: - file "*_trim.fastq.gz" into fastq_files_trim - - script: - """ - cutadapt -q 20,20 \ - -o ${reads.baseName}_trim.fastq.gz \ - ${reads} > ${reads.baseName}_report.txt - """ -}