Fecal microbiota of NZBWF1
Data analysis was carried out using the modified protocol, according to Chi-Cheng Huang. Amplicon sequencing was performed by using 300-bp paired-end raw reads and assembled using FLASH v.1.2.7. De-multiplexing was carried out based on barcode identification. For quality control, Q < 20 reads were discarded from the pipeline of QIIME (v1.9.1). If three consecutive bases were < Q20, the read was classified as truncated. The read was preserved in the data set if it had more than 75% of the original length using split_libraries_fastq.py script in QIIME. Using UCHIME algorithms, sequences were chimera-checked to acquire the effective tags and filtered from the data set by using the UPARSE function in the peptide of USEARCH (v.7). The information of operational taxonomic unit (OUT) abundance was normalized to the variation and rarefied to the minimum sequence depth using the QIIME script (single_rarefaction.py).