Scalable statistical methods for microbial metagenomics