code for finding fixed sites/empirical means