Quantitative MRI methods and learning-based algorithms require exact forward simulations. One critical factor to correctly describe magnetization dynamics is the effect of slice-selective RF pulses. While contemporary simulation techniques correctly capture their influence, they only provide final magnetization distributions, require to be run for each parameter set separately, and make it hard to derive general theoretical conclusions and to generate a fundamental understanding of echo formation in the presence of slice-profile effects. This work aims to provide a mathematically exact framework, which is equally intuitive as extended phase graphs (EPGs), but also considers slice-profiles through their natural spatial representation. We show, through an analytical, hybrid Bloch-EPG formalism, that the spatially-resolved EPG approach allows to exactly predict the signal dependency on off-resonance, spoiling moment, microscopic dephasing, and echo time. We also demonstrate that our formalism allows to use the same phase graph to simulate both gradient-spoiled and balanced SSFP-based MR sequences. We present a derivation of the formalism and identify the connection to existing methods, i.e. slice-selective Bloch, slice-selective EPG, and the partitioned EPG. As a use case, the proposed hybrid Bloch-EPG framework is applied to MR Fingerprinting.