Bioinformatics and statistical analysis
The quality control and sequence filtering of raw reads was performed according to the barcode matching and sequence overlapping using QIIME (version 1.9.1). The filtered high-quality reads were clustered into operational taxonomic units (OTUs) with a similarity level of 97%, then the obtained OTUs were used for taxonomical assignments by RDP classifier (