WGS of Yellow mutant B. terrestris
We selected a subset of samples for the whole-genome resequencing approach (both wildtype (n=7) and yellow mutant (n=7)) from the specimens used for RAD-Seq. To optimize the potential to find fixed sites within the association peak, we identified unique RAD-Seq haplotypes within the peak of association and selected individuals which maximally represented the range of differences in RAD-Seq haplotypes. Genomic DNA extraction was conducted from thoracic tissues of haploid males using either Qiagen DNeasy or E.Z.N.A. (Omega Bio-tek) kit followed up by RNaseA treatment. 150 bp paired-end sequencing libraries were prepared using Illumina TruSeq DNA Nano kits for whole-genome re-sequencing approach implemented in Illumina HiSeq 2500 sequencer (5-10 individuals per lane aimed at ~30X per sample coverage) in Pennsylvania State University Genomics Core Facility (University Park, Pennsylvania, USA).