This is something many people struggle with. First, your premises are correct; the reason that a single application of the phase encode gradient is not sufficient is NOT because phase does not vary across this dimension each time we apply the gradient. It is because  we need a way to detect the phase difference. This is done with the Fourier transform, which requires multiple samples of the MR signal that differ in the way the signal is encoded along the dimension we are trying to discriminate.