DNA of honey
All DNA was extracted from the honey sample. Extracted genomic DNA was fragmented to 200 bp using Covaris S2 focused-ultrasonicator (Covaris, MA, USA). Whole-genome libraries were prepared using the FOCUS (Fragmented DNA Compact Sequencing Assay) NIPT protocol at CCHT with 12 cycles for the final PCR enrichment step. Following library quantification, 36 samples were pooled equimolarly, and the quality and quantity of the pool were assessed on Agilent 2200 TapeStation (Agilent Technologies, USA). WGS was performed on the NextSeq 550 or NovaSeq 6000 SP instrument (Illumina Inc.) 85 bp single-end.