We investigate the implications of the latest LHCb measurement of $R_K$ for NP explanations of the $B$ anomalies. The previous data could be explained if the $b to s mu^+ mu^-$ NP is in (I) $C_{9,{rm NP}}^{mumu}$ or (II) $C_{9,{rm NP}}^{mumu} = -C_{10,{rm NP}}^{mumu}$, with scenario (I) providing a better explanation than scenario (II). This continues to hold with the new measurement of $R_K$. However, for both scenarios, this measurement leads to a slight tension of $O(1sigma)$ between separate fits to the $b to s mu^+ mu^-$ and $R_{K^{(*)}}$ data. In this paper, we investigate whether this tension can be alleviated with the addition of NP in $b to s e^+ e^-$. In particular, we examine the effect of adding such NP to scenarios (I) and (II). We find several scenarios in which this leads to improvements in the fits. $Z$ and LQ models with contributions to both $b to s mu^+ mu^-$ and $b to s e^+ e^-$ can reproduce the data, but only within scenarios based on (II). If the tension persists in future measurements, it may be necessary to consider NP models with more than one particle contributing to $b to s ell^+ ell^-$.