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 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 } \mbox { $M _ { \hbox { $ \odot$ } } $ } and the 50 % Bayesian credible interval is ( 0.62 , 0.81 ) \times 10 ^ { 12 } \mbox { $M _ { \hbox { $ \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 \ > \mbox { kpc } ) = 0.26 ~ { } ( 0.20 , 0.36 ) \times 10 ^ { 12 } \mbox { $M _ { \hbox { $ \odot% $ } } $ } , M ( r = 50 \ > \mbox { kpc } ) = 0.37 ~ { } ( 0.29 , 0.51 ) \times 10 ^ { 12 } \mbox { $M _ { \hbox { $ \odot% $ } } $ } , M ( r = 100 \ > \mbox { kpc } ) = 0.53 ~ { } ( 0.41 , 0.74 ) \times 10 ^ { 12 } \mbox { $M _ { \hbox { $% \odot$ } } $ } , etc . ) , making it easy to compare our results directly to other studies .