La forge institutionnelle d'INRAE étant en production depuis le 10 juin 2025, nous vous invitons à y créer vos nouveaux projets.
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]])