bayes_abundance(db, bayes, reads, rank = ['family','order','class','phylum','species']);
[as DatabaseReader]
[as AbundanceEstimate]
all reads data in one sample
[as Environment]