A quark coalescence model is presented based on semi-relativistic molecular dynamics with color interactions among quarks taken into account and applied to pp collisions to study the effects of this model. A phenomenological potential with two tunable parameters is introduced to describe the color interactions between quarks and antiquarks. The interactions drive the process of hadronization and finally make them form different color neutral clusters, which can be identified as hadrons based on some criteria. A Monte Carlo generator, PYTHIA is used to generate the quarks in the initial state of hadronization, and different values of tunable parameters are used to study their effects on the final state distributions and correlations. Baryon to meson ratios, transverse momentum spectra, pseudorapidity distributions and forward-backward multiplicity correlations of hadrons produced in the hadronization process from this model with different parameters are compared with those from PYTHIA.