In this paper, a single-server queue with negative customers is considered. The arrival of a negative customer will remove one positive customer that is being served, if any is present. An alternative approach will be introduced to derive a set of equations which will be solved to obtain the stationary queue length distribution. We assume that the service time distribution tends to a constant asymptotic rate when time t goes to infinity. This assumption will allow for finding the stationary queue length of queueing systems with non-exponential service time distributions. Numerical examples for gamma distributed service time with fractional value of shape parameter will be presented in which the steady-state distribution of queue length with such service time distributions may not be easily computed by most of the existing analytical methods.