Data-driven discovery of model error in chaotic systems by integrating Bayesian sparse regression and data assimilation