sequences of bacteria
16S rRNA gene V3-V4 region was amplified using fusion primers containing Illumina adapters and barcode sequences. PCR amplicons were purified, pooled equimolarly, and sequenced on Illumina platform with paired-end 250 bp (PE250) strategy.