Genealogical relationships of DNA sequences were assessed
with phylogenetic trees (for both RAG1 and ND5) and a statistical
parsimony network (for ND5 only). Phylogenetic trees were constructed
with the program PAUP v4.0b10 (Swofford, 2003) (Neighbour
Joining, Maximum Parsimony and Maximum Likelihood
methods) and MrBayes v3.1 (Ronquist and Huelsenbeck, 2003) for
Bayesian method. All trees were rooted with outgroup sequences
of three Mobula species. Both transitions and transversions at all
codon positions were included in the analyses.