A method for stochastic unraveling of general time-local quantum master equations (QMEs) is proposed. The present kind of jump algorithm allows a numerically efficient treatment of QMEs which are not in Lindblad form, i.e. are not positive semidefinite by definition. The unraveling can be achieved by allowing for trajectories with negative weights. Such a property is necessary, e.g. to unravel the Redfield QME and to treat various related problems with high numerical efficiency. The method is successfully tested on the damped harmonic oscillator and on electron transfer models including one and two reaction coordinates. The obtained results are compared to those from a direct propagation of the reduced density matrix (RDM) as well as from the standard quantum jump method. Comparison of the numerical efficiency is performed considering both the population dynamics and the RDM in the Wigner phase space representation.