All nucleotide sequences determined by Miseq (referred to as “reads”) were converted to FASTAQ format on MiSeq reporter V2.3 and subsequently analyzed using CLC Genomics Workbench 6.0 (CLC bio, Cambridge, MA, USA). Briefly, the ends of all reads were trimmed to remove adaptor sequences located at both ends of each read. Trimmed reads were assembled into contigs using a de novo assembly algorithm. Contigs generated by de novo assembly algorithm were analyzed using BlastN.