Genetic data were obtained at the Center for Anchored Phylogenomics at Florida State University (http://www.anchoredphylogeny.com; Last accessed June 22, 2018) using the general methods outlined in Lemmon et al. (2012) (link). First, each genomic DNA sample was sonicated to a fragment size of 150–350 bp using a Covaris E220 focused-ultrasonicator. Subsequently, library preparation and indexing were performed following a protocol outlined in (Meyer and Kircher 2010 (link)). Indexed samples were then pooled at equal quantities determined using a Qubit fluorometer (16 samples per pool), and samples in each pool were enriched using an Agilent Custom SureSelect kit (Agilent Technologies) with custom designed probes. The general enrichment process uses streptavidin-coated magnetic beads to isolate targeted genomic fragments that hybridize with the oligonucleotide probes (Gnirke et al. 2009 (link)). Enriched fragments were pooled in groups of three (three pools total with 48 samples in each) before sequencing on the PE150 Illumina HiSeq2500 (three lanes, 48 samples per lane). Sequencing was performed in the Translational Science Laboratory in the College of Medicine at Florida State University.