fix an error with the tag, and address limits for derepFasta that req…

…uire limiting the number of processed reads to 100k
h3abionet · Feb 8, 2024 · e29ba5d · e29ba5d
1 parent 90cdee8
commit e29ba5d
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/main.nf b/main.nf
@@ -732,7 +732,7 @@ if (params.reads != false || params.input != false ) { // TODO maybe we should c
     } else {
 
         process PerSampleInferDerepAndMerge {
-            tag { "PerSampleInferDerepAndMerge: ${meta.id}" }
+            tag { "PerSampleInferDerepAndMerge:${meta.id}" }
             publishDir "${params.outdir}/dada2-Derep-Single/Per-Sample", mode: "copy", overwrite: true
 
             input:

diff --git a/templates/PerSampleDadaInfer.R b/templates/PerSampleDadaInfer.R
@@ -11,13 +11,13 @@ if (!is.na(dadaOpt)) {
 cat("Processing:", "${meta.id}", "\\n")
 
 errF <- readRDS("errors.R1.RDS")
-derepF <- derepFastq("${reads[0]}")
+derepF <- derepFastq("${reads[0]}", n=100000)
 ddF <- dada(derepF, err=errF, multithread=${task.cpus}, pool=as.logical("${params.pool}"))
 saveRDS(ddF, "${meta.id}.dd.R1.RDS")
 
 if (file.exists ("errors.R2.RDS")) {
   errR <- readRDS("errors.R2.RDS")
-  derepR <- derepFastq("${reads[1]}")
+  derepR <- derepFastq("${reads[1]}", n=100000)
   ddR <- dada(derepR, err=errR, multithread=${task.cpus}, pool=as.logical("${params.pool}"))
   saveRDS(ddR, "${meta.id}.dd.R2.RDS")