The pyrosequencing data was initially processed as follows: (i) the terminal adaptor sequences and the reverse primer sequences were removed, as were reads containing multiple Ns, poly A/T tails, and low quality sequences; (ii) the barcode index sequences and forward primer sequences were removed; (iii) sequences with lengths