Previously, several levels of phylogenetic relationships in an insect order Odonata have been estimated using morphological and molecular markers. For the molecular phylogeny rRNA sequences were mainly, but other markers were not frequently employed. In this study, we sequenced both two mitochondrial genes (COI and 16S rRNA) and nuclear genes (28S rRNA and elongation factor-1α), composed of ~4,002 bp from 71 species of Odonata, occurring mostly in South Korea. These concatenated sequences were utilized to test the previous phylogenetic hypotheses of Odonata via Bayesian Inference (BI) and Maximum Likelihood (ML) algorithms, along with the data partition option available in BI method. Each families and superfamilies represented by multiple taxa consistently supported monophylies with the highest nodal supports in both Anisoptera and Zygoptera. A close relationship of Anisozygoptera to Anisoptera represented by a single species was obvious. On the other hand, familial relationships within each suborder of Anisoptera and Zygoptera have shown two compelling topologies. The topology obtained by BI method with partitioning of the four genes showed an unresolved relationship among Gomphidae, Aeshnidae, and the suborder Anisozygoptera in Anisoptera clade, presenting the relationships ((((Libellulidae + Corduliidae) + Macromiidae) + (Gomphidae + Aeshnidae + Anisozygoptera)) + (((Coenagrionidae + Platycnemdidae) + Calopterygidae) + Lestidae)). Another topology obtained by both BI and ML methods without partitioning, on the other hand, placed Anisozygoptera the basal lineage of Anisoptera, but Lestidae in Zygoptera was placed as the sister to Anisoptera + Anisozygoptera, presenting the relationships (((((((Libellulidae + Corduliidae) + Macromiidae) + Aeshnidae) + Gomphidae) + Anisozygoptera) + Lestidae) + ((Coenagrionidae + Platycnemdidae) + Calopterygidae)). Topological test to find out better supported tree turned out a slight higher support for the former topology, but the monophyly of Zygoptera with the inclusion of Lestidae was supported only poorly (BPP = 0.68) in the former topology.