A new moving mesh scheme based on the Lagrange–Galerkin method for the approximation of the one-dimensional convection–diffusion equation is studied. The mesh movement is prescribed by a discretized dynamical system for the nodal points. This system is related to the velocity and diffusion coefficient in the convection–diffusion equation such that the nodal points follow the convective flow of the model. It is shown that under a restriction of the time step size the mesh movement cannot lead to an overlap of the elements and therefore an invalid mesh. Using a piecewise linear approximation, optimal error estimates in the ℓ∞(L2)∩ℓ2(H01)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\ell ^\\infty (L^2) \\cap \\ell ^2(H_0^1)$$\\end{document} norm are proved in case of both, a first-order backward Euler method and a second-order two-step method in time. These results are based on new estimates of the time dependent interpolation operator derived in this work. Preservation of the total mass is verified for both choices of the time discretization. Numerical experiments are presented that confirm the error estimates and demonstrate that the proposed moving mesh scheme can circumvent limitations that the Lagrange–Galerkin method on a fixed mesh exhibits.