In this paper, we develop and characterize the fully dissipative Lattice Boltzmann method for ultra-relativistic fluids in two dimensions using three equilibrium distribution functions: Maxwell-Juttner, Fermi-Dirac and Bose-Einstein. Our results stem from the expansion of these distribution functions up to fifth order in relativistic polynomials. We also obtain new Gaussian quadratures for square lattices that preserve the spatial resolution. Our models are validated with the Riemann problem and the limitations of lower order expansions to calculate higher order moments are shown. The kinematic viscosity and the thermal conductivity are numerically obtained using the Taylor-Green vortex and the Fourier flow respectively and these transport coefficients are compared with the theoretical prediction from Grads theory. In order to compare different expansion orders, we analyze the temperature and heat flux fields on the time evolution of a hot spot.