initialLambdaVec is now scaled with mutation rate
put ???? when triophasing non-mendelian sites