We give a counting formula in terms of modified Hall-Littlewood polynomials and the chromatic quasisymmetric function for the number of points on an arbitrary Hessenberg variety over a finite field. As a consequence, we express the Poincaré polynomials of complex Hessenberg varieties in terms of a Hall scalar product involving the symmetric functions above. We use these results to give a new proof of a combinatorial formula for the modified Hall-Littlewood polynomials.