We derive the triply differential spectrum for the inclusive rare decay B → X s ℓ + ℓ − in the shape function region, in which X s is jet-like with m 2 X < ∼ m b Λ QCD . Experimental cuts make this a relevant region. The perturbative and non-perturbative parts of the matrix elements can be defined with the Soft-Collinear Effective Theory, which is used to incorporate α s corrections consistently. We prove that, with a suitable power counting for the dilepton invariant mass, the same universal jet and shape functions appear as in B → X s γ and B → X u ℓν decays. Parts of the usual α s (m b ) corrections go into the jet function at a lower scale, and parts go into the non-perturbative shape function. For B → X s ℓ + ℓ − , the perturbative series in α s are of a different character above and below µ = m b . We introduce a "split matching" method that allows the series in these regions to be treated independently.