We calculate the thermophoretic drift of a charged single colloidal particle with hydrodynamically slipping surface immersed in an electrolyte solution in response to a small temperature gradient. Here we rely on a linearized hydrodynamic approach for the fluid flow and the motion of the electrolyte ions while keeping the full nonlinearity of the Poisson-Boltzmann equation of the unperturbed system to account for possible large surface charging. The partial differential equations are transformed into a coupled set of ordinary differential equations in linear response. Numerical solutions are elaborated for parameter regimes of small and large Debye shielding and different hydrodynamic boundary conditions encoded in a varying slip length. Our results are in good agreement with predictions from recent theoretical work and successfully describe experimental observations on thermophoresis of DNA. We also compare our numerical results with experimental data on polystyrene beads.