Massive multiple-input multiple-output (MIMO) techniques with a large number of antenna elements at base station (BS) have been proved as an alternative to provide potential opportunity to increase the spectrum and energy efficiency. However, in the system, there generally exists a spatial correlation effect due to insufficient antenna elements spacing and/or the lack of rich scattering at BS. The minimum mean square error (MMSE) method performs signal detection at the expense of large-scale matrix inversion operation. Thus, the conjugate gradient (CG) method has received a lot of attentions to realize the MMSE detection efficiently. Unfortunately, this efficiency can be compromised due to the ill-conditioned equalization matrix of MMSE method over the correlated channel environments. Moreover, the hard output signal detection exhibits a sharply degradation in performance for higher-order quadrature amplitude modulation (QAM). Therefore, the modern communication systems use the soft-output information, i.e., log-likelihood ratio (LLR) along with the forward error-correcting code (FEC) to achieve satisfactory performance. The LLR computation along with a higher-order QAM remains challenging due to the exhaustive search of symbol in the modulation constellation. In this paper, a low-complexity soft-output signal detector based on approximate inverse symmetric successive over-relaxation preconditioned conjugate gradient (AI-SSOR-CG-SOD) method is proposed to realize MMSE method detection for uplink multiuser massive MIMO correlated channel. In the proposed method, a new preconditioner, an AI-SSOR, which is based on the Neumann series approximation of the inverse of the conventional SSOR preconditioner is firstly developed to handle ill-conditioned matrix, and then incorporated with CG method to improve the convergence rate and performance. According to the characteristic of the Gray-coding that adjacent symbols in the constellation set have only one different bit, the constellation set is divided multiple times based on the bits of the inphase and the quadrature components of the symbol, which reduces the complexity of the LLR computation of the transmitted bits by avoiding the exhaustive search process. Simulation results show that the AI-SSOR preconditioner is robust against spatial correlation effect, and the proposed detector converges at 3 iterations. Simulation results also show that the proposed detector achieves a better trade-off between the complexity and the performance compared to other existing detectors.