For singularly perturbed two-point boundary value problems, bounds for the solution and its derivatives are established, and a B-spline collocation method on Shishkin meshes is proposed for numerical approximation. Stability and convergence estimates in the maximum norm are derived, showing convergence of order 2 up to a logarithmic factor. Numerical examples demonstrate both the accuracy and the convergence order of the proposed numerical method.