For nuclear sequences, two alignments were created: first, heterozygous sites were coded using the IUPAC ambiguity codes (unphased dataset); and second, gametic phases of nuclear sequences were inferred using the Bayesian algorithm implemented in PHASE 2.1 [34] (phased dataset).