This paper presents a hybrid stress approach for the analysis of laminated composite plates. The plate mechanical model is based on the so called First-order Shear Deformation Theory, rationally deduced from the parent three-dimensional theory. Within this framework, a new quadrilateral four-node finite element is developed from a hybrid stress formulation involving, as primary variables, compatible displacements and elementwise equilibrated stress resultants. The element is designed to be simple, stable and locking-free. The displacement interpolation is enhanced by linking the transverse displacement to the nodal rotations and a suitable approximation for stress resultants is selected, ruled by the minimum number of parameters. The transverse stresses through the laminate thickness are reconstructed a posteriori by simply using three-dimensional equilibrium. To improve the results, the stress resultants entering the reconstruction process are first recovered using a superconvergent patch-based procedure called Recovery by Compatibility in Patches, that is properly extended here for laminated plates. This preliminary recovery is very efficient from the computational point of view and generally useful either to accurately evaluate the stress resultants or to estimate the discretization error. Indeed, in the present context, it plays also a key role in effectively predicting the shear stress profiles, since it guarantees the global convergence of the whole reconstruction strategy, that does not need any correction to accommodate equilibrium defects. Actually, this strategy can be adopted together with any plate finite element. Numerical testing demonstrates the excellent performance of both the finite element and the reconstruction strategy.