Extraction of parton distributions and $\alpha_{\mathrm{s}}$ from DIS data within a Bayesian treatment of systematic errors