Fitting the observational data of the twin peak kHz quasiperiodic oscillations (QPO) from low mass X-ray binaries (LMXBs) by the relativistic precession model gives a substantially higher neutron star mass estimate, M similar to 2 M-circle dot, than the"canonical value", M similar to 1.4 M-circle dot. Using a fully general relativistic approach we discuss the non-geodesic corrections to the orbital and epicyclic frequencies of slightly charged circularly orbiting test particles caused by the presence of a neutron star magnetic field. We show that consideration of such non-geodesic corrections can bring down the neutron star mass estimate and improve the quality of twin peak QPO data fits based on relativistic precession frequency relations.