We study aspects of chaos and thermodynamics at strong coupling in a 1+1d QFT using numerical Hamiltonian truncation methods. We find that our eigenstate spectrum satisfies Wigner-Dyson statistics and that the coefficients describing eigenstates in our basis satisfy Random Matrix statistics both at weak and strong coupling. We also find scar states, but only at weak coupling. We then use these chaotic states to compute thermodynamic observables, obtaining results consistent with CFT expectations at temperatures above the scale of relevant interactions. Finally, we test the Eigenstate Thermalization Hypothesis in a new regime by being able to access the expectation value of local QFT operators in eigenstates with high energy densities. These expectation values show universality and are found to be consistent with analytic finite temperature CFT results.