I tried to find out HWE for different HLA loci. When I repeated the analysis using same set of data, I found that HWE p values have changed.
Pypop user guide says that "HardyWeinbergGuoThompson overall p-value test was numerically unstable
because it attempted to check for equality in greater than or
equal to constructs ("