In standard treatments of open quantum systems, the reduced dynamics is described starting from the assumption that the system and the environment are initially uncorrelated. This assumption, however, is not always guaranteed in realistic scenarios and several theoretical approaches to characterize initially correlated dynamics have been introduced. For the uncorrelated scenario, stochastic unravelings are a powerful tool to simulate the dynamics. So far they have not been used in the most general case in which correlations are initially present since they cannot be applied to nonpositive operators or noncompletely positive maps. In our work, we employ the bath positive (B+) or one-sided positive decomposition (OPD) formalism as a starting point to generalize stochastic unraveling in the presence of initial correlations. Noticeably, our approach does not depend on the particular unraveling technique, but holds for both piecewise deterministic and diffusive unravelings. This generalization allows not only for more powerful simulations for the reduced dynamics, but also for a deeper theoretical understanding of open system dynamics.