map $_ /= ( $mutRate * 100), @$_ for @lnPAM;
 @PAM1 = map[ map exp( $_ ), @$_ ], @lnPAM;
 $mutRate = sum map{ $freq[ $_ ] * ( 1 - $PAM1[ $_ ][ $_ ] ) } 0..19;