Gaussian processes are powerful tools for modeling and predicting various numerical data. Hence, checking their quality of fit becomes a vital issue. In this article, we introduce a testing methodology for general Gaussian processes based on a quadratic form statistic. We illustrate the methodology on three statistical tests recently introduced in the literature, which are based on the sample autocovariance function, time average mean-squared displacement, and detrended moving average statistics. We compare the usefulness of the tests by taking into consideration three very important Gaussian processes: the fractional Brownian motion, which is self-similar with stationary increments (SSSIs), scaled Brownian motion, which is self-similar with independent increments (SSIIs), and the Ornstein-Uhlenbeck (OU) process, which is stationary. We show that the considered statistics' ability to distinguish between these Gaussian processes is high, and we identify the best performing tests for different scenarios. We also find that there is no omnibus quadratic form test; however, the detrended moving average test seems to be the first choice in distinguishing between same processes with different parameters. We also show that the detrended moving average method outperforms the Cholesky method. Based on the previous findings, we introduce a novel procedure of discriminating between Gaussian SSSI, SSII, and stationary processes. Finally, we illustrate the proposed procedure by applying it to real-world data, namely, the daily EURUSD currency exchange rates, and show that the data can be modeled by the OU process.