Refactoring mendelian error function.

Now all as numpy array operations

TO DO: remove old map function
New tests, fixed bug in homref vs homalt

Added new multivariant test.

Moved benchmark code to src

Fixed bug where had forgotton to index in the homref vs homalt case