# Report significance from Tukey HSD

i run a Tukey HSD test inside the R package agricolae and got significant differences between my groups. Now I would like to report that in a figure with the basic significance asterisks. How can I obtain the p-values from Tukey HSD test? Or would you rather commend using a pairwise wilcoxon test since I directly obtain the p-values of each pair?

    > print(tukey_both_rarefied_Observed)

$statistics MSerror Df Mean CV 45.63004 37 12.45 54.25704$parameters
test  name.t ntr StudentizedRange alpha
Tukey Staging   3         3.452781  0.05

$means Observed std r Min Max Q25 Q50 Q75 nonlesional 14.608696 5.781940 23 4 26 10.00 14.0 19.0 Patch 11.000000 8.694826 11 1 31 5.50 9.0 15.5 Plaque 6.833333 6.274286 6 1 17 1.75 5.5 10.0$comparison
NULL

\$groups
Observed groups
nonlesional 14.608696      a
Patch       11.000000     ab
Plaque       6.833333      b

attr(,"class")
[1] "group"
$$$$
`