Table 2 Genetic diversity indices for each clade of Bombus morio (Main: main clade; TS: Teodoro Sampaio clade) and B. pauloensis (N: north clade; C: central clade; S: south clade) obtained from 1,575 bp of the following concatenated mitochondrial regions: Cytochrome C oxidase I (COI), Cytochrome B (CytB), the large ribosomal RNA subunit (16S), and cluster 4 of tRNA (covering a region of COII and ATPase 8 genes and tRNA^lys and tRNA^asp)

Ns

181

2

183

21

110

90

221

h

87

2

89

6

36

30

78

Hd

0.963

-

0.964

0.72381

0.91009

0.87491

0.95483

k

4.349

-

4.92374

1.20952

3.30525

2.89588

8.378

π

0.00278

-

0.00315

0.00077

0.00184

0.00170

0.00537

Tajima`s D (p-value)

-1.85387 (0.00600)

0.00000 (1.00000)

-1.97254 (0.00200)

-1.22665 (0.10600)

-1.54169 (0.03300)

-1.69378 (0.02000)

-0.61631 (0.30300)

Fu`s Fs (p-value)

-25.40396 (0.00000)

0.00000 (0.22800)

-25.18616 (0.00100)

-1.60993 (0.11000)

-25.66400 (0.00000)

-17.96060 (0.00000)

-24.21740 (0.00000)