Genomic DNA extraction was performed according to TaKaRa MiniBEST Universal Genomic DNA Extraction Kit protocol. The concentration and quality of extracted DNA was determined using NanoDrop2000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA), Quant-iT PicoGreen dsDNA assay kit (Thermo Fisher Scientific, Wilmington, DE, USA) and agarose gel electrophoresis. Genomic DNA was then fragmented using Covaris to a targeted size of 350 bp and upon adapter ligation, a library containing fragments of 400 bp was generated. The library size was determined using Agilent Bioanalyzer 2100 (Agilent, Santa Clara, CA, USA). Paired-end DNA sequencing was performed with Illumina NvoSeq 6000.