A method of magnetic resonance (MR) imaging of an object includes: generating MR signals by subjecting the object to a number N of shots of a multi-echo imaging sequence comprising multi-slice RF pulses for simultaneously exciting two or more spatially separate image slices, with a phase offset in the slice direction being imparted to the MR signals; acquiring the MR signals that are received in parallel via a set of at least two RF coils having different spatial sensitivity profiles; and reconstructing a MR image for each image slice from the acquired MR signals using a parallel reconstruction algorithm, wherein the MR signal contributions from the different image slices are separated on the basis of the spatial encodings of the MR signals according to the spatial sensitivity profiles of the RF coils and of the phase offsets attributed to the respective image slices and shots.