The effect of spin-orbit coupling (SOC) in monolayer graphene has been thoroughly studied in the past. Its impact is commonly disregarded since the electronic structure of graphene is almost unaffected. Consequently, SOC has also been ignored in other graphitic systems. In this work, we calculate the intrinsic SOC in stacks of graphene layers from a model that includes all the valence orbitals of carbon as well as the intra-atomic SOC.The leading spin-orbit term in a graphene multilayer combines a Rashba-like spin-flip process with interlayer hopping. We find that this intrinsic term modifies the phase diagram and Fermi surfaces of finite rhombohedral stacks, and it leads to a distinct density-dependence of the magnetization in weak magnetic fields. This Letter shows that, unlike in the single-layer case, SOC has important consequences in multilayer graphitic systems and it should not be neglected when studying their physical properties.