Skip to content

Commit 848305b

Browse files
committed
get basic trimming in line again w/ cutadapt outputs, in case this trimming option is chosen
1 parent b41fc88 commit 848305b

1 file changed

Lines changed: 5 additions & 3 deletions

File tree

modules/local/illumina_filterandtrim.nf

Lines changed: 5 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -8,7 +8,7 @@ process ILLUMINA_DADA2_FILTER_AND_TRIM {
88
tuple val(meta), path(reads)
99

1010
output:
11-
tuple val(meta), path("${meta.id}.R[12].trim.fastq.gz"), optional: true, emit: trimmed
11+
tuple val(meta), path("${meta.id}*.trim.fastq.gz"), optional: true, emit: trimmed
1212
path("*.trimmed.txt"), emit: trimmed_report
1313

1414
when:
@@ -17,14 +17,16 @@ process ILLUMINA_DADA2_FILTER_AND_TRIM {
1717
script:
1818
def args = task.ext.args ?: ''
1919
def prefix = task.ext.prefix ?: "${meta.id}"
20+
def r1prefix = meta.single_end ? "" : "_1"
21+
def r2prefix = meta.single_end ? "" : "_2"
2022
"""
2123
#!/usr/bin/env Rscript
2224
suppressPackageStartupMessages(library(dada2))
2325
2426
out <- filterAndTrim(fwd = "${reads[0]}",
25-
filt = "${meta.id}.R1.trim.fastq.gz",
27+
filt = "${meta.id}${r1prefix}.trim.fastq.gz",
2628
rev = if("${reads[1]}" == "null") NULL else "${reads[1]}",
27-
filt.rev = if("${reads[1]}" == "null") NULL else "${meta.id}.R2.trim.fastq.gz",
29+
filt.rev = if("${reads[1]}" == "null") NULL else "${meta.id}${r2prefix}.trim.fastq.gz",
2830
trimLeft = if("${reads[1]}" == "null") ${params.trim_for} else c(${params.trim_for}, ${params.trim_rev}),
2931
truncLen = if("${reads[1]}" == "null") ${params.trunc_for} else c(${params.trunc_for}, ${params.trunc_rev}),
3032
maxEE = if("${reads[1]}" == "null") ${params.maxEE_for} else c(${params.maxEE_for}, ${params.maxEE_rev}),

0 commit comments

Comments
 (0)