This talk will present an update on results from the high statistics data from simulations of 2+1 flavor p4 (N_T=8) and Asqtad (N_T=6, 8) improved staggered fermion actions using the exact RHMC algorithm. These results are from ongoing calculations being done by the HotQCD collaboration on the BlueGene L supercomputers at LLNL and the NYBlue at BNL and represent over 100M node-hours. We will present results for the equation of state (EOS) and discuss the behavior of a variety of other observables including the Polyakov Loop, chiral and quark number susceptibility across the crossover. A case will be made that results for the EOS have reached a level of precision that they should be used in the phenomenological analysis of Relativistic Heavy Ion collisions.