The evolutionary model used
was Hasegawa–Kishino–Yano (HKY) (Hasegawa et al., 1985) for
the RAG1 data set and HKY + U (a = 0.2044) for the ND5 data set
as inferred by hierarchical likelihood ratio test implemented in
Modeltest v3.7 (Posada and Crandall, 1998). A statistical parsimony
haplotype network of ND5 data was constructed with the software
TCS v1.21 (Clement et al., 2000) with a 95% parsimony connection
limit.