Population relationships were further assessed by modeling migration events with TREEMIX given a number of migration edges from m=1 to m=14. Using the Evanno method we identified that 10 migration events best explained the sample covariance of allele frequencies among louse populations (Figure 8). Gene flow with the highest migration weight across continents was between Hungary and Egypt. The analysis also showed a second gene flow event between Hungary and Turkey; however, the migration weight was lower. Hungary and Egypt still appear to be an outlier (not geographically clustered) in the ML tree generated through TREEMIX. Another intercontinental gene flow event was found between the Philippines and Colombia. All other admixture events were between populations within continents (Figure 8). Overall, the population relationships between the TREEMIX ML tree and SNP phylogeny are in agreement with each other.