Detailed knowledge about equilibration processes is of interest for various fields of physics, including heavy ion collision experiments and quantum quenched condensed matter systems. We study the approach to equilibrium at late times within two types of strongly coupled thermal systems in 3 + 1 dimensions: systems in the presence of (i) a nonzero charge density, or (ii) a magnetic field at vanishing charge density. Utilizing the gauge/gravity correspondence, we map the aforementioned problem to the computation of quasinormal frequencies around two particular classes of black branes within the Einstein-Maxwell theory. We compute (i) the tensor and vector quasinormal modes of Reissner-Nordstrom black branes and (ii) the scalar, as well as tensor quasinormal modes of magnetic black branes. Some of these quasinormal modes correspond to the late-time relaxation of the above systems after starting with initial pressure anisotropy. We provide benchmarks which need to be matched at late-times by all holographic thermalization codes with the appropriate symmetries.