The optimal Jastrow-type ground-state wave function of spin-aligned atomic hydrogen is calculated using the pair potential of Kolos and Wolniewicz. The optimization is performed by solving the Euler equation in the hypernetted chain approximation. Accurate energies as well as pair-distribution functions are obtained. The Bose-Einstein condensate fraction is evaluated from the one-particle momentum distribution. The pair distribution function is also used to obtain stability criteria for the system and minimal values for the aligning magnetic field are calculated at low densities. The resulting values of the minimal aligning fields are considerably higher than those obtained previously.