Phase 3
Following amplification and size verification, the libraries were sequenced using a 2X300 bp paired-end run (MiSeq Reagent Kit v3) on an Illumina MiSeq platform. The Miseq paired-end reads were further processed to obtain trimmed paired-end reads. The quality control of raw sequence data (FastQ) was performed to remove low-quality reads (quality score, <30) by using FastQC v0.11.4.