In this paper, we give the expression of the least square solution of the linear quaternion matrix equation AXB = C subject to a consistent system of quaternion matrix equations D1X = F1, XE2 =F2, and derive the maximal and minimal ranks and the leastnorm of the above mentioned solution. The finding of this paper extends some known results in the literature.