  The LHCb collaboration has very recently released a new study of $B^+ \to
K^{+} \ell^+ \ell^-$ and $B \to K^{*0} \ell^+ \ell^-$ ($\ell = e,\mu$) decays,
testing lepton universality with unprecedented accuracy using the whole Run 1
and 2 dataset. In addition, the CMS collaboration has recently reported an
improved analysis of the branching ratios $B_{(d,s)}\to\mu^+\mu^-$. While these
measurements offer, per se, a powerful probe of New Physics, global analyses of
$b \to s \ell^+ \ell^-$ transitions also rely on the assumptions about
nonperturbative contributions to the decay matrix elements. In this work, we
perform a global Bayesian analysis of New Physics in (semi)leptonic rare $B$
decays, paying attention to the role of charming penguins which are difficult
to evaluate from first principles. We find data to be consistent with the
Standard Model once rescattering from intermediate hadronic states is included.
Consequently, we derive stringent bounds on lepton universality violation in
$|\Delta B | = | \Delta S| = 1$ (semi)leptonic processes.
