To obtain an extended set of sequences that included not only SNPs but also information about insertions and deletions, we extracted FLM genomic sequences (including 2 kb upstream and 0.38 kb downstream sequence) from the GEBrowser 3.0 resource (http://signal.salk.edu/atg1001/3.0/gebrowser.php). A set of 850 sequences was manually curated and aligned using MEGA7.0.14 (Tamura et al., 2011). Some sequences were excluded since they possessed a high number of ambiguous bases. The core sequence set consisted of 840 sequences. This sequence set was used for association analysis (Figure 3—figure supplements 4–7) and all further analysis on FLM PRO2+ and INT6+ variation (Figure 8 and Figure 8—figure supplement 1). Accession identifiers, geographic data, and ADMIXTURE group membership (k = 9) were obtained from the 1001 genomes tool collection (http://tools.1001genomes.org/) (The 1001 Genomes Consortium, 2016).