In this paper, a stabilized space-time finite element method for solving linear parabolic evolution problems is analyzed. The proposed method is developed on a base of a space-time variational setting, that helps on the simultaneous and unified discretization in space and in time by finite element techniques. Stabilization terms are constructed by means of classical bubble spaces. Stability of the discrete problem with respect to an associated mesh dependent norm is proved, and a priori discretization error estimates are presented. Numerical examples confirm the theoretical estimates.