We proposed a mathematical model for an incompressible, viscous, natural convection, and stagnation point slip flow of MHD Prandtl fluid over an infinite plate. The governing flow equations are constructed using the Prandtl rheological model. In account of physical relevance, we investigated the Soret and Dufour effects on the flow field. The complex natured flow equations are transformed to a set of PDEs using a suitable similarity variables. The non-dimensionalized ruling problem together with physical boundary conditions is numerically analyzed via Crank-Nicolson scheme. The velocity, temperature and concentration of the diffusing species distributions are enhanced for higher values of unsteadiness parameter. It is noted that velocity is slightly decreasing for higher values of Reynolds number while smaller values of providing more dominant effects on the velocity, temperature and concentration of the diffusing species profiles and enhanced heat and mass transfer rates is noticed. The physical behavior of reduced Nusselt and Sherwood numbers, friction factor, for distinct values of emerging parameters is examined and representative set of graphs are presented.