design: DNA extractions were performed on bee legs and/or muscles using Omega Biotek EZNA DNA extraction or Qiagen DNeasy kits with standard protocols with RNase removal. All genomes were sequenced with 150 bp paired end Illumina sequencing at mean ~17 times coverage (ranging from 11.7 to 44.9x).