Arviz is the companion package to pymc that handles these sorts of analysis-of-posterior jobs.
You could try az.summary for a tabular description of the mean, variance, and HDI of your posterior, or az.plot_posterior to see the actual plot with the mean and HDI labeled automatically.