diff --git a/makeSequenceTable.R b/makeSequenceTable.R index 364397cdfe2cde364d826e45b0239d6eef9e5a77..cdbfbde918b9165018d579645359dd5cb83e1872 100644 --- a/makeSequenceTable.R +++ b/makeSequenceTable.R @@ -9,7 +9,7 @@ seqtab.nochim <- removeBimeraDenovo(seqtab, method="consensus", multithread = sn saveRDS(seqtab.nochim, snakemake@output$rds) -uniquesToFasta(seqtab.nochim, snakemake@output$fasta) +uniquesToFasta(seqtab.nochim, snakemake@output$fasta, ids = getSequences(seqtab.nochim)) write.table(t(seqtab.nochim), snakemake@output$tsv, sep = "\t", quote = FALSE)