For Illumina sequencing, genomic DNA was randomly interrupted into ~350 bp fragments using a Covaris ultrasonic crusher. Libraries were prepared with the NEBNext  UltraTM DNA Library Prep Kit for Illumina (NEB, Cat. No. E7370L) and sequenced on the Illumina Novaseq 6000 platform (Illumina, Inc., San Diego, CA, USA) with a paired-end (PE) 150 bp mode.