Alpha and beta diversity analyses were performed using QIIME’s core_diversity_analyses.py, at rarefaction depths of 5721 (to retain as many samples as possible) and 10,000 to confirm that the results were similar with more sequences per sample.