The uncertainty in the properties estimated using atomistic simulations are estimated by propagating the uncertainties for the free energy terms involved in the calculation of the property. The method of propagation of the uncertainties from reweighted observables using the covariance matrix has been already described in our previous work. Only the uncertainties of the properties estimated from the simulations contribute to the uncertainty calculation give the uncertainties of IAPWS95 derived properties are very small.