diff --git a/src/RNAseq.nf b/src/RNAseq.nf index e1d9f8729e10a68a67850ac2ddde4992df2cbf7f..ebe2ecdf33d82b64275882220671bd471dbc541d 100644 --- a/src/RNAseq.nf +++ b/src/RNAseq.nf @@ -332,35 +332,35 @@ process fastqc_genome { /* Deduplication of reads */ //////////////////////////// -process dedup_genome { - tag "$file_id" - publishDir "${params.output}/03_hisat2/dedup/", mode: 'copy' - - input: - set file_id, file(bam) from HISAT_ALIGNED_DEDUP - - output: - set file_id, "*.{bam, bai}" into DEDUP_GENOME - file "*.log" into DEDUP_LOG - - when: - params.do_dedup - - """ - umi_tools dedup -I ${bam[0]} \ - -S ${file_id}.dedup.bam \ - --paired \ - 2> ${file_id}_dedup.log - """ -} - -if (! params.do_dedup) { - HISAT_ALIGNED_DEDUP.into{DEDUP_GENOME} +if (params.do_dedup) { + process dedup_genome { + tag "$file_id" + publishDir "${params.output}/03_hisat2/dedup/", mode: 'copy' + + input: + set file_id, file(bam) from HISAT_ALIGNED_DEDUP + + output: + set file_id, "*.{bam, bai}" into DEDUP_GENOME + file "*.log" into DEDUP_LOG + + when: + params.do_dedup + + """ + umi_tools dedup -I ${bam[0]} \ + -S ${file_id}.dedup.bam \ + --paired \ + 2> ${file_id}_dedup.log + """ + } +} else { + HISAT_ALIGNED_DEDUP.set{DEDUP_GENOME} } DEDUP_GENOME.into{DEDUP_GENOME_HTSEQ; DEDUP_GENOME_RNASEQC - } + } /////////// /* HTseq */ @@ -474,6 +474,9 @@ hisat2 -x ${index2_id} \ -2 ${fastq_unaligned[1]} \ --rna-strandness ${params.strand} \ --dta\ + --no-softclip\ + --trim3 1\ + --trim5 1\ 2> ${file_id}_postgenome.txt \ | samtools view -bS -F 4 - \ | samtools sort -@ ${task.cpus} -o ${file_id}.bam \ @@ -491,27 +494,30 @@ POSTGENOME_ALIGNED.into{POSTGENOME_ALIGNED_FASTQC; //////////////////////////// /* Deduplication of reads */ //////////////////////////// +if (params.do_dedup) { + process dedup_postgenome { + tag "$file_id" + publishDir "${params.output}/05_post_genome_hisat2/dedup/", mode: 'copy' -process dedup_postgenome { - tag "$file_id" - publishDir "${params.output}/05_post_genome_hisat2/dedup/", mode: 'copy' + input: + set file_id, file(bam) from POSTGENOME_ALIGNED_DEDUP - input: - set file_id, file(bam) from POSTGENOME_ALIGNED_DEDUP + output: + set file_id, "*.{bam, bai}" into DEDUP_POSTGENOME + file "*.log" into DEDUP_POSTGENOME_LOG - output: - set file_id, "*.{bam, bai}" into DEDUP_POSTGENOME - file "*.log" into DEDUP_POSTGENOME_LOG + when: + params.do_dedup - when: - params.do_dedup + """ + umi_tools dedup -I ${bam[0]} \ + -S ${file_id}.dedup.bam \ + --paired \ + 2> ${file_id}_dedup.log + """ + } +} else { - """ - umi_tools dedup -I ${bam[0]} \ - -S ${file_id}.dedup.bam \ - --paired \ - 2> ${file_id}_dedup.log - """ } process fastqc_postgenome { diff --git a/src/RibosomeProfiling.nf b/src/RibosomeProfiling.nf index a145b80c27b7378975cee891ee9df9154f84ac14..48ce6e55c45388e684f6a5efb984f197dd5e3600 100644 --- a/src/RibosomeProfiling.nf +++ b/src/RibosomeProfiling.nf @@ -458,6 +458,9 @@ hisat2 -x ${index2_id} \ -U ${fastq_unaligned[0]} \ --rna-strandness ${params.strand} \ --dta\ + --no-softclip\ + --trim3 1\ + --trim5 1\ 2> ${file_id}_postgenome.txt \ | samtools view -bS -F 4 -F 256 - \ | samtools sort -@ ${task.cpus} -o ${file_id}.bam \