dada2.R 281 bytes
library(dada2)
derep <- derepFastq(snakemake@input[[1]], verbose = TRUE)
err <- learnErrors(snakemake@input[[1]], multithread = snakemake@threads, verbose = TRUE)
dada <- dada(derep, err = err, multithread = snakemake@threads, verbose = TRUE)
saveRDS(dada, snakemake@output[[1]])