Table 3 Posterior probabilities for all models based on multinomial logistic regression; models specifically tested for nuclear loci only

Equal gene flow between both chilensis and pallatangae, and chilensis and albiceps

Model D

Higher rate of gene flow between chilensis and pallatangae compared to chilensis and albiceps (asymmetric gene flow model)

Model F

Asymmetric gene flow model with gene flow until the last 20,000 years and isolation thereafter (isolation-migration model).

Model G