Illumina HiSeq 2500 paired end sequencing
DNA was sheared with an S2 focused-ultrasonicator (Covaris, Woburn, MA, USA) to achieve a target fragment size of around 100-900 bp according to the manufacturer’s instructions. The KAPA High Throughput Library Preparation Kit with Standard PCR Library Amplification (KAPA Biosystems, Wilmington, MA, USA) was used to generate Illumina sequencing libraries according to the manufacturer’s instructions. All libraries were amplified with seven PCR cycles. Illumina libraries were subjected to a dual size selection using AMPure XP beads (Beckman Coulter, Brea, CA, USA) with a target size window of 300 to 1000 bp.