We consider the J/psi photo-production data collected at HERA in the light of next-to-leading order predictions for the color-singlet yield and polarization. We find that, while the shapes of inclusive distributions in the transverse momentum and inelasticity are well reproduced, the experimental rates are larger than those given by the color-singlet contribution alone. Furthermore, the next-to-leading order calculation predicts the J/psis to be mostly longitudinally polarized at high transverse momentum in contrast with the trend of the preliminary data from the ZEUS collaboration.