We study the distribution of work induced by the two-point measurement protocol for a driven open quantum system. We first derive a general form for the generating function of work for the total system, bearing in mind that the Hamiltonian does not necessarily commute with its time derivative. Using this result, we then study the first few moments of work by using the master equation of the reduced system, invoking approximations similar to the ones made in the microscopic derivation of the reduced density matrix. Our results show that already in the third moment of work, correction terms appear that involve commutators between the Hamiltonian and its time derivative. To demonstrate the importance of these terms, we consider a sinusoidally, weakly driven and weakly coupled open two-level quantum system, and indeed find that already in the third moment of work, the correction terms are significant. We also compare our results to those obtained with the quantum jump method and find a good agreement.