“…Regions with ambiguous positions that could not be aligned were difficult to identify manually and were removed using Gblocks (Castresana, 2000), with the 37, 95, 102-104, 165-176, 182, 206, 211-262, 268-84, 298, 325-328, 366, 382, 397, 421, 530, 551, 741-1008, 1022-1023, 1045, 1121-1124, 1128, 1170-1171, 1189-1191, 1293, 1395, 1719-1923, 2254-2264, 2272-2273, 2284-2286. Substitution saturation was measured using Xia's test (Xia et al, 2003;Xia and Lemey, 2009) implemented in DAMBE (Xia and Xie, 2001). No saturation was detected in the 16S alignment (321 sites) from which gaps and ambiguous regions had been removed (Iss significantly < Iss.c).…”