Purpose: The purpose of this paper is to discuss about evaluating the integrals involving B-spline wavelet on the interval (BSWI), in wavelet finite element formulations, using Gauss Quadrature. Design/methodology/approach: In the proposed scheme, background cells are placed over each BSWI element and Gauss quadrature rule is defined for each of these cells. The nodal discretization used for BSWI WFEM element is independent to the selection of number of background cells used for the integration process. During the analysis, background cells of various lengths are used for evaluating the integrals for various combination of order and resolution of BSWI scaling functions. Numerical examples based on one-dimensional (1D) and two-dimensional (2D) plane elasto-statics are solved. Problems on beams based on Euler Bernoulli and Timoshenko beam theory under different boundary conditions are also examined. The condition number and sparseness of the formulated stiffness matrices are analyzed. Findings: It is found that to form a well-conditioned stiffness matrix, the support domain of every wavelet scaling function should possess sufficient number of integration points. The results are analyzed and validated against the existing analytical solutions. Numerical examples demonstrate that the accuracy of displacements and stresses is dependent on the size of the background cell and number of Gauss points considered per background cell during the analysis. Originality/value: The current paper gives the details on implementation of Gauss Quadrature scheme, using a background cell-based approach, for evaluating the integrals involved in BSWI-based wavelet finite element method, which is missing in the existing literature. © 2019, Emerald Publishing Limited.