Bayesian inference of mean velocity fields and turbulence models from flow MRI