This study presents a fitted operator numerical method for solving singularly perturbed boundary value problems with integral boundary condition. The stability and parameter uniform convergence of the proposed method are proved. To validate the applicability of the scheme, a model problem is considered for numerical experimentation and solved for different values of the perturbation parameter, $\varepsilon$ and mesh size, $h$. The numerical results are tabulated in terms of maximum absolute errors and rate of convergence and it is observed that the present method is more accurate and $\varepsilon$-uniformly convergent for $h\geq \varepsilon$ where the classical numerical methods fails to give good result and it also improves the results of the methods existing in the literature.