We present results on the spontaneous crystallization of freely-jointed polymers of hard spheres obtained in an unprecedentedly long Monte Carlo (MC) simulation on a system of 54 chains of 1000 monomers. Starting from a purely amorphous configuration and after a transitory dominance of the hexagonal closed packed (HCP) polymorph, the system crystallizes in a final, stable, face centered cubic (FCC) crystal of very high perfection. Through refined metrics we gauge the degree of ordering and identify the regions of the phase transition and the corresponding morphologies. An analysis of chain conformational characteristics, of the spatial distribution of monomers and of the volume accessible to them shows that the phase transition is caused by an increase in translational entropy that is larger than the loss of conformational entropy of the chains in the crystal compared to the amorphous state. Polymer chains in the crystal adopt ideal random walk statistics as their great length renders local conformational details, imposed by the geometry of the FCC crystal, irrelevant.