Phylogenetic tree of Gobiidae fishes based on 12 protein-coding genes excluding ND6 gene. The tree was built based on the TIM3+I+G model (-ln L= 119,585.39, A:C:G:T= 0.3172, 0.3187, 0.1016, 0.2625, I= 0.2830, G= 0.4140). The topology generated by Bayesian inference is shown, and is identical to the result of maximum likelihood analyses. Bayesian posterior probabilities and bootstrap support values above 50% replicates are shown adjacent to each internode. M. swinhonis and P. glenii were used as outgroups.