design: 1 g of genomic DNA was sheared using 96 microTUBE-50 AFA Fiber plates (Covaris Inc., Woburn, MA, USA) to a target size of 300-400 bp with a Covaris E220 system (Covaris Inc., Woburn, MA, USA).  gDNA samples were cleaned on a small volume Mosquito liquid handler (TTP labtech) with a 1:1 ratio of Kapa Pure beads (Roche, Basel, Switzerland) and eluted in Tris-Cl, pH 8.0. Genomic libraries were generated with 500 ng sheared DNA using Kapa Hyper library prep kit (Roche, Basel, Switzerland) and barcodes were added using the Illumina UD 96 index kit (Illumina, San Diego, CA, USA). Final libraries were PCR-amplified for 5 cycles with Kapa HIFI PCR kit (Roche, Basel, Switzerland). Samples were pair-end sequenced on either Hiseq 4000 or Novaseq at a read length of 2 x 150 bp.