The tidal stresses in a radially heterogeneous moon are calculated numerically as a function of time and location by using the latest moon model obtained from seismic data. Theoretical results are compared with the observed features of deep moonquakes and in particular with those of of the most active A1 hypocenter. The results show that (1) the magnitudes of tidal stresses are small, of the order of 1 bar and less, (2) they reach maxima in seismically active depth range of 600--1200 km, (3) the geographic distribution of moonquake epicenters does not appear to be controlled solely by the highest tidal stresses, and the presence of heterogeneities and weak zones may be important in determining moonquake locations, (4) the occurrence times of moonquakes correlate with one or more of the tidal stress components, and (5) the presence of an ambient tectonic stress of about the same order of magnitude as the tidal stress is necessary to explain the observed reversal in polarity of the A1 moonquakes.