Applying Nyquist's sampling theorem to a real signal
I do not see specifically what you think is incorrect. I follow your thinking and do not see a flaw but let's review the key characteristics below in case that fills in a hole in the thought process: 

For periodic systems we can consider a time span that extends over one period (0 to T) or from $-\infty$ to $+\infty$. Both are equivalent as long as we agree that the $\infty$ case must be periodic. 

Similarly for sampled systems we can consider a frequency span that extends over one frequency range, 0 to $F_s$, where $F_s$ represents the sampling frequency, or from $-\infty$ to $+\infty$. Both are equivalent as long as we agree that the $\infty$ case must be similarly periodic in frequency.  (Note the base frequency that repeats can also be given as $-F_s/2$ to $+F_s/2$ but the range 0 to $F_s$ is most common when referring to the DFT.)

So the DFT itself is both periodic in time and in frequency, representing a sampled system that is periodic in time. 

For any other systems; aperiodic or not sampled, we cannot use the DFT and call that equivalent, we can only approximate those systems with the DFT. 

Below is a summary of the various Fourier Transforms that we can use in those other cases:

CASE 1: Continuous and Aperiodic in time: The Fourier Transform

First case we have is the Fourier Transform, where the time function is by definition continuous and extending from $-\infty$ to $+\infty$. Since it is continuous in time (not sampled) the frequency result will be aperiodic (not repeating). Since it is aperiodic in time, the frequency result will be continuous. 



CASE 2: Continuous and periodic in time: The Fourier Series Expansion

Note that x(t) can still be periodic in the case above; there is nothing that excludes taking the Fourier Transform as written in the graphic above on a periodic function. This will result in the continuous Fourier Transform result as being 0 except for the frequencies at multiples of 1/T, including zero, which may have non-zero values at only these locations.  Knowing this, we can solve for just those frequencies in this case, resulting in the Fourier Series Expansion:



CASE 3: Discrete and aperiodic in time: The Discrete Time Fourier Transform

The third case is to now consider the sampled version of our first case, where we have an aperiodic time domain function again, but this time it is sampled in time. This is the Discrete Time Fourier Transform. We have a similar result to the first case, the Fourier Transform, with the only change that it is now repeating in frequency. As previously described the repetition boundaries occur at multiples of $F_s =  1/T$, so if our frequency spectrum in our first case extended beyond $F_s$, we would get frequency folding (aliasing) as a result and therefore would not be able to represent our time domain waveform exactly (the sampling rate was not high enough!).

(Side note: Zero padding prior to taking the DFT approximates the DTFT, the more zeros you add in, the more the frequency result approximates a continuous function, and that function is samples of the DTFT!).



CASE 4: Discrete and periodic in time: The Discrete Fourier Transform

Finally, we come to the DFT (the FFT is an algorithm to compute the DFT exactly). In this case, when considering the same time axis as above that extends from $-\infty$ to $+\infty$, the time domain waveform must be repeating (periodic).  The periodicity in time results in discrete frequencies that can only exist at integer multiples of 1/T Hz where T represents one period in time in seconds. Note that like the Fourier Series Expansion, the periodicity is implied; the actual transform considers N samples in time and provides N samples in frequency, but the periodicity described is evident when both axis are extended to infinity.



I hope this answered your question. If not please comment on the remaining points of confusion and I will try and update accordingly.