Active learning for data-driven reduced models of parametric differential systems with Bayesian operator inference