We present the first official release of the nCTEQ nuclear parton distribution functions (nPDFs) with errors. The main addition to the previous nCTEQ PDFs is the introduction of PDF uncertainties based on the Hessian method. Another important improvement is the inclusion of pion production data from RHIC giving us a handle to constrain gluon PDF. In this presentation we briefly discuss the framework of our analysis and concentrate on the comparison of our results with those of other groups providing nPDFs. Additionally we present predictions for selected results from the LHC heavy ion collisions.