A maximum likelihood phylogeny of the derived from 13 coding protein of mtDNA sequences for Anura. Bayesian inference obtained the same topology. Numbers above the lines or beside the nodes were rapid bootstrap proportions calculated with 1000 replicates and Bayesian posterior probabilities, respectively. The branch lengths are to scale. Mt genomic features that are useful markers for inferring relationships, and genomic features of each species are labeled on tree as follows: (A) LTPF cluster; (B) trnS (AGY) pseudogene; (C) translocation of trnH; (D) translocation of nad5; (E) LTPF changed into TLPF cluster; (F) duplication of CR and trnM; (G) duplication of trnM; (H) translocation of trnE and loss one copy of trnM; (I) trnA, trnN, trnC and trnE gene loss; (J) trnQ, trnA, trnN, trnC, trnY, and trnP gene loss; (K) translocation of nad5; (L) translocation of trnP; (M) translocation of trnL(CUN); (N) WANCY changed into WACYN; (O) translocation of nad5, trnE, and trnP.