The massive MIMO communication system for cellular networks has shown excellent performance in terms of high spectral and energy efficiencies in recent years. This achievement has been due to the assiduous research undertaken in terms of formulating mathematical expressions to compute these efficiencies. The purpose of this work is to find a generalized approach to compute the error rates of a user terminal connected to a linear receiver of a Massive MIMO system. In line with the objectives set for present and future cellular networks to achieve ultra-reliable communication systems the issue of finding bit error rates required per user terminal is evaluated for a range of BS antennas is addressed in this paper. The approach employed to tackle this issue is to consider three practical linear receivers such as Maximum Ratio Combining (MRC), Zero-Forcing (ZF), Minimum Mean Square Error (MMSE) altogether with an uncorrelated Rayleigh fading environment in the case of imperfect Channel State Information (CSI). Based on solid mathematical foundations, the approach to find the Bit Error Rate (BER) closed-form expressions which are simple to evaluate and have not been reported in literatures yet, are derived for MPSK and MQAM modulation techniques. We verified the validity of our analysis by simulating all mathematical formulations using Matlab R2018 and we showed interesting results such as the number of BS antennas required to achieve the desired BER for a specific uplink receiver and modulation scheme. The mathematical analysis presented in this work can be applied to future cellular systems to compute the distribution of error rates for a user with a given number of BS antennas and modulation scheme. Moreover, Machine Learning algorithms could be applied to the analysis presented to predict the bit error rates for any user terminal in the cellular network.