Practical and Rigorous Uncertainty Bounds for Gaussian Process Regression