In this work we show the strong convergence of propagation of chaos for the particle approximation of McKean-Vlasov SDEs with singular L p -interactions as well as the moderate interaction particle system in the level of particle trajectories. One of the main obstacles is to establish the strong well-posedness of particle system with singular interaction. In particular, we develop the theory of strong well-posedness of Krylov and Röckner [25] in the case of mixed L p -drifts, where the heat kernel estimates play a crucial role. Moreover, when the interaction kernel is bounded measurable, we also obtain the optimal rate of strong convergence, which is partially based on Jabin and Wang's entropy method [19] and Zvonkin's transformation.