Orthogonal Time-Frequency Space modulation stands out as a promising waveform for 6G and beyond wireless communication systems, offering superior performance over conventional methods, particularly in high-mobility scenarios and dispersive channel conditions. Error performance analysis remains crucial for accurately characterizing the reliability of wireless communication systems under practical constraints. In this paper, we systematically investigate the bit error rate performance of OTFS modulation over Nakagami-m fading channels in both single-user and multi-user scenarios. In analytical approaches, mathematical frameworks are employed for distinct receiver configurations: the Single-input Single-output scenario leverages Erlang probability density function of squared-Nakagami variables to derive closed-form BER expressions, while the Single-input Multiple-output case applies moment matching techniques with Gamma approximation to model multiple user interference, subsequently yielding Signal-to-interference-plus-noise Ratio characterizations through Meijer-G functions. This study examines single-path and multi-path channel conditions, evaluating the relationship between path multiplicity and error performance metrics while considering various fading intensities through Nakagami-m fading parameters. The derived closed-form BER expressions are validated through maximum likelihood detection based Monte Carlo simulations, demonstrating strong correlation between analytical and numerical results across various SNR regions. Furthermore, comparative benchmark evaluations against conventional orthogonal frequency division multiplexing with MLD reveal that OTFS consistently achieves superior error performance in high-mobility scenarios. In multipath fading environments, OTFS achieves superior diversity gain compared to conventional OFDM, which refers to enhanced error performance.