Phylogenetic tree of Gobiidae fish based on 12 protein-coding genes excluding ND6 gene. The tree was built based on the TIM3+1+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 analysis. 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.