In this article, the authors present finite element analysis and approximation of Burgers’‐Fisher equation. Existence and uniqueness of weak solution is proved by Galerkin's finite element method for non‐smooth initial data. Next, a priori error estimates of semi‐discrete solution in
L
∞
(
0
,
T
;
L
2
(
Ω
)
)
norm, are derived and the convergence of semi‐discrete solution is established. Then, fully discretization of the problem is done with the help of Euler's backward method. The nonlinearity is removed by lagging it to previous known level. The scheme is found to be convergent. Positivity of fully discrete solution is discussed, and bounds on time step are discovered for which the solution preserves its positivity. Finally, numerical experiments are performed on some examples to demonstrate the effectiveness of the scheme. The proposed scheme found to be fast, easy and accurate.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1652–1677, 2017