In this study, performance analysis of variable-gain amplify-and-forward relaying orthogonal frequency division multiplexing-based systems over independently but not necessarily identically distributed Nakagami-m fading channels using non-linear power amplifier (PA) at the relay and maximal ratio combining scheme at destination. Specifically, novel closed-form expressions of outage probability and asymptotic outage probability are investigated for the considered system model. Further, an average symbol error rate (ASER) expression is derived for general-order rectangular quadrature amplitude modulation (RQAM) scheme using cumulative distribution function-based approach, which is also the novel contribution of the work. Since the transfer function of PA used at relay is nonlinear, this introduces non-linear distortion (NLD) in the considered system model. Thus, the impact of NLD and fading parameters on outage probability, diversity gain of the system and ASER performances are analysed. Further, the impact of NLD on constellation order for RQAM scheme is also discussed. The theoretical results are compared with computer simulations to verify the accuracy of the derivations. The derived expressions are generalised over variety of fading environments for integer-valued fading parameter.