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