Table 2 Maximum-likelihood estimations (MLEs) and 95% of the highest posterior density (HPD) for population parameters of the populations from the SCS (TTN), the Indian Ocean (TPK) and the Indo-Pacific boundary (MKJ) under the estimated substitution rate of 1.616 * 10⁻⁹ substitutions per site per year (s/s/y; 95% HPD: 1.19–2.09 * 10⁻⁹ s/s/y) using the isolation with migration model. For each parameter, 95% HPD was shown in parentheses

t₁ (MYA)

0.067 (0.000–9.222)

0.049 (0.000–6.791)

0.038 (0.000–5.241)

Ne _TTN

60 (0–1,015)

44 (0–748)

34 (0–577)

Ne _MKJ

2,927 (896–9,497)

2,155 (660–6993)

1,663 (509–5397)

Ne _TPK

537 (179–1,971)

396 (132–1451)

305 (102–1120)

2NM _{TTN -> MKJ}

2NM _{MKJ -> TTN}

2NM _{MKJ -> TPK}

0.142 ^**(0.043–0.347)

2NM _{TPK -> MKJ}

0.002 (0.000–0.603)

2NM _{TTN -> TPK}

0.013 (0.000–0.197)

2NM _{TPK -> TTN}

6.375 *10⁻⁴ (0.000–0.037)