We study the impact of including the baryonic decay $\Lambda_b\to \Lambda(\to p \pi^-)\mu^+\mu^-$ in a Bayesian analysis of $|\Delta B | = |\Delta S| = 1$ transitions. We perform fits of the Wilson coefficients $C_{9}$, $C_{9'}$, $C_{10}$ and $C_{10'}$, in addition to the relevant nuisance parameters. Our analysis combines data for the differential branching fraction and three angular observables of $\Lambda_b\to \Lambda(\to p \pi^-)\mu^+\mu^-$ with data for the branching ratios of $B_s \to \mu^+\mu^-$ and inclusive $b \to s\ell^+\ell^-$ decays. Newly available precise lattice QCD results for the full set of $\Lambda_b \to \Lambda$ form factors are used to evaluate the observables of the baryonic decay. Our fits prefer shifts to $C_{9}$ that are opposite in sign compared to those found in global fits of only mesonic decays, and the posterior odds show no evidence of physics beyond the Standard Model. We investigate a possible hadronic origin of the observed tensions between theory and experiment.