Bayesian Spline Learning for Equation Discovery of Nonlinear Dynamics with Quantified Uncertainty