The emission lineshape of the B850 band in the light harvesting complex 2 of purple bacteria is calculated by extending the approach of 2nd order time-nonlocal quantum master equation [S. Jang and R. J. Silbey, J. Chem. Phys. 118, 9312 (2003)10.1063/1.1569239]. The initial condition for the emission process corresponds to the stationary excited state density where exciton states are entangled with the bath modes in equilibrium. This exciton-bath coupling, which is not diagonal in either site excitation or exciton basis, results in a new inhomogeneous term that is absent in the expression for the absorption lineshape. Careful treatment of all the 2nd order terms are made, and explicit expressions are derived for both full 2nd order lineshape expression and the one based on secular approximation that neglects off-diagonal components in the exciton basis. Numerical results are presented for a few representative cases of disorder and temperature. Comparison of emission line shape with the absorption line shape is also made. It is shown that the inhomogeneous term coming from the entanglement of the system and bath degrees of freedom makes significant contributions to the lineshape. It is also found that the perturbative nature of the theory can result in negative portion of lineshape in some situations, which can be removed significantly by inclusion of the inhomogeneous term and completely by using the secular approximation. Comparison of the emission and absorption lineshapes at different temperatures demonstrates the role of thermal population of different exciton states and exciton-phonon couplings.