added function to calculate a(x, P(i))

Now you can estimate a(x, P(i)) between two site patterns. this is
super messy right now, there has to be a faster way, but I just don’t
know what that it right now. It works…