dada2_post-processing.R 305 Bytes
Newer Older
Bastien Macé's avatar
Bastien Macé committed
1
2
3
4
5
6
7
8
tab <- read.table(Aquarium_2.fasta, header=T)
seqtab_1 <- makeSequenceTable(tab)
seqtab_2 <- removeBimeraDenovo(seqtab_1, verbose=T)
## processes the bimera removal

uniqueSeqs <- getUniques(seqtab_2)
uniquesToFasta(uniqueSeqs, paste0(sample.names, ".fasta")
## creates the new file without bimeras