####Supplemental Data: Generalized Linear Model (GLM) used for comparing ecological attributes impacting haplotype diversity.

glm(formula = Haplotype ~ pH + Organic_Carbon + Nitrogen + Altitude + 
    Precipitation + Temperature, family = gaussian(link = "identity"), 
    data = normalized_data)

Coefficients:
               Estimate Std. Error t value Pr(>|t|)   
(Intercept)     11.3939     3.9710   2.869  0.00668 **
pH              -5.7919     5.3741  -1.078  0.28794   
Organic_Carbon   0.6554     2.4890   0.263  0.79374   
Nitrogen        -1.1276     2.7727  -0.407  0.68652   
Altitude        -0.6691     3.7730  -0.177  0.86018   
Precipitation   -1.7526     7.2940  -0.240  0.81140   
Temperature     -3.1233     4.3241  -0.722  0.47453   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for gaussian family taken to be 10.75421)

    Null deviance: 459.20  on 44  degrees of freedom
Residual deviance: 408.66  on 38  degrees of freedom
AIC: 242.98

Number of Fisher Scoring iterations: 2