We present new mass estimates and cumulative mass profiles (CMPs) with Bayesian credible regions for the Milky Way (MW) Galaxy, given the kinematic data of globular clusters as provided by (1) the $\textit{Gaia}$ DR2 collaboration and the HSTPROMO team, and (2) the new catalog in Vasiliev (2019). We use globular clusters beyond 15kpc to estimate the CMP of the MW, assuming a total gravitational potential model $\Phi(r) = \Phi_{\circ}r^{-\gamma}$, which approximates an NFW-type potential at large distances when $\gamma=0.5$... We compare the resulting CMPs given data sets (1) and (2), and find the results to be nearly identical. The median estimate for the total mass is $M_{200}= 0.70 \times 10^{12} M_{\odot}$ and the $50\%$ Bayesian credible interval is $(0.62, 0.81)\times10^{12}M_{\odot}$. However, because the Vasiliev catalog contains more complete data at large $r$, the MW total mass is slightly more constrained by these data. In this work, we also supply instructions for how to create a CMP for the MW with Bayesian credible regions, given a model for $M(<r)$ and samples drawn from a posterior distribution. With the CMP, we can report median estimates and $50\%$ Bayesian credible regions for the MW mass within any distance (e.g., $M(r=25\text{kpc})= 0.26~(0.20, 0.36)\times10^{12}M_{\odot}$, $M(r=50\text{kpc})= 0.37~(0.29, 0.51) \times10^{12}M_{\odot}$, $M(r=100\text{kpc}) = 0.53~(0.41, 0.74) \times10^{12}M_{\odot}$, etc. ), making it easy to compare our results directly to other studies. (read more)