import simuPOP as sim
from simuPOP.utils import Exporter

Naatamo=sim.Population(
   size=[50,50],
   ploidy=2,
   loci=[17,2,11,5,7,5,3,3,14,15,5,8,16,13,6,6,2,8,4,6,6,5,3,9,5,8,4,4],
   lociPos=[
      4.9,21.4,23.9,26.0,34.1,34.4,48.9,50.8,78.6,82.4,118.2,119.9,120.1,122.3,131.1,136.1,151.3,
      24.7,43.3,
      28.1,28.6,32.1,37.9,40.3,59.6,63.6,71.9,74.8,77.4,91.8,
      9.0,20.0,47.5,47.7,55.8,
      22.3,22.7,23.7,34.7,34.8,39.9,54.5,
      35.9,38.9,55.9,59.8,68.0,
      5.5,18.0,47.0,
      18.5,19.8,22.0,
      18.2,18.3,31.5,47.8,70.2,84.4,94.4,96.3,103.9,110.5,111.5,120.1,128.7,139.0,
      9.3,17.8,19.3,20.8,24.2,29.1,31.7,37.0,64.1,67.9,74.6,74.9,75.0,75.4,113.0,
      5.7,14.5,19.4,23.2,57.6,
      26.4,34.2,34.3,36.9,61.5,62.5,66.6,78.4,
      9.3,12.5,20.1,22.2,23.2,27.6,44.8,53.1,53.9,55.8,71.0,75.0,75.1,75.2,76.5,85.5,
      14.8,15.2,23.9,34.1,37.0,54.3,54.6,55.4,60.7,61.1,64.8,80.1,82.3,
      5.5,40.9,41.9,48.7,52.8,60.0,
      16.9,42.0,68.1,69.0,71.0,81.9,
      27.0,50.8,
      13.9,23.3,24.2,27.1,28.2,53.6,56.3,66.4,
      11.4,51.2,51.4,72.5,
      33.0,33.8,47.4,50.5,63.1,69.1,
      14.4,20.0,20.1,38.8,43.6,51.6,
      6.5,7.5,18.4,20.1,44.2,
      30.8,33.6,40.2,
      1.3,1.4,11.1,12.1,12.4,17.0,19.8,34.4,47.0,
      11.9,14.9,34.1,34.5,47.1,
      13.5,19.8,21.0,23.9,25.8,26.4,27.0,28.9,
      16.1,21.4,22.6,26.2,
      8.6,13.2,28.1,35.8
   ],
 
   alleleNames=['1','2'],
   subPopNames=['escapee','wild']

)
sim.initSex(Naatamo)

sim.initGenotype(Naatamo,subPops=[0] , loci=[0], freq=[0.515,0.485])
sim.initGenotype(Naatamo,subPops=[0] , loci=[1], freq=[0.535,0.465])
sim.initGenotype(Naatamo,subPops=[0] , loci=[2], freq=[0.431,0.569])
sim.initGenotype(Naatamo,subPops=[0] , loci=[3], freq=[0.372,0.628])
sim.initGenotype(Naatamo,subPops=[0] , loci=[4], freq=[0.423,0.577])
sim.initGenotype(Naatamo,subPops=[0] , loci=[5], freq=[0.479,0.521])
sim.initGenotype(Naatamo,subPops=[0] , loci=[6], freq=[0.36,0.64])
sim.initGenotype(Naatamo,subPops=[0] , loci=[7], freq=[0.489,0.511])
sim.initGenotype(Naatamo,subPops=[0] , loci=[8], freq=[0.424,0.576])
sim.initGenotype(Naatamo,subPops=[0] , loci=[9], freq=[0.601,0.399])
sim.initGenotype(Naatamo,subPops=[0] , loci=[10], freq=[0.444,0.556])
sim.initGenotype(Naatamo,subPops=[0] , loci=[11], freq=[0.179,0.821])
sim.initGenotype(Naatamo,subPops=[0] , loci=[12], freq=[0.475,0.525])
sim.initGenotype(Naatamo,subPops=[0] , loci=[13], freq=[0.515,0.485])
sim.initGenotype(Naatamo,subPops=[0] , loci=[14], freq=[0.617,0.383])
sim.initGenotype(Naatamo,subPops=[0] , loci=[15], freq=[0.515,0.485])
sim.initGenotype(Naatamo,subPops=[0] , loci=[16], freq=[0.467,0.533])
sim.initGenotype(Naatamo,subPops=[0] , loci=[17], freq=[0.425,0.575])
sim.initGenotype(Naatamo,subPops=[0] , loci=[18], freq=[0.263,0.737])
sim.initGenotype(Naatamo,subPops=[0] , loci=[19], freq=[0.376,0.624])
sim.initGenotype(Naatamo,subPops=[0] , loci=[20], freq=[0.455,0.545])
sim.initGenotype(Naatamo,subPops=[0] , loci=[21], freq=[0.599,0.401])
sim.initGenotype(Naatamo,subPops=[0] , loci=[22], freq=[0.552,0.448])
sim.initGenotype(Naatamo,subPops=[0] , loci=[23], freq=[0.618,0.382])
sim.initGenotype(Naatamo,subPops=[0] , loci=[24], freq=[0.626,0.374])
sim.initGenotype(Naatamo,subPops=[0] , loci=[25], freq=[0.468,0.532])
sim.initGenotype(Naatamo,subPops=[0] , loci=[26], freq=[0.5,0.5])
sim.initGenotype(Naatamo,subPops=[0] , loci=[27], freq=[0.454,0.546])
sim.initGenotype(Naatamo,subPops=[0] , loci=[28], freq=[0.535,0.465])
sim.initGenotype(Naatamo,subPops=[0] , loci=[29], freq=[0.347,0.653])
sim.initGenotype(Naatamo,subPops=[0] , loci=[30], freq=[0.449,0.551])
sim.initGenotype(Naatamo,subPops=[0] , loci=[31], freq=[0.572,0.428])
sim.initGenotype(Naatamo,subPops=[0] , loci=[32], freq=[0.417,0.583])
sim.initGenotype(Naatamo,subPops=[0] , loci=[33], freq=[0.444,0.556])
sim.initGenotype(Naatamo,subPops=[0] , loci=[34], freq=[0.369,0.631])
sim.initGenotype(Naatamo,subPops=[0] , loci=[35], freq=[0.393,0.607])
sim.initGenotype(Naatamo,subPops=[0] , loci=[36], freq=[0.667,0.333])
sim.initGenotype(Naatamo,subPops=[0] , loci=[37], freq=[0.309,0.691])
sim.initGenotype(Naatamo,subPops=[0] , loci=[38], freq=[0.505,0.495])
sim.initGenotype(Naatamo,subPops=[0] , loci=[39], freq=[0.49,0.51])
sim.initGenotype(Naatamo,subPops=[0] , loci=[40], freq=[0.606,0.394])
sim.initGenotype(Naatamo,subPops=[0] , loci=[41], freq=[0.616,0.384])
sim.initGenotype(Naatamo,subPops=[0] , loci=[42], freq=[0.46,0.54])
sim.initGenotype(Naatamo,subPops=[0] , loci=[43], freq=[0.396,0.604])
sim.initGenotype(Naatamo,subPops=[0] , loci=[44], freq=[0.434,0.566])
sim.initGenotype(Naatamo,subPops=[0] , loci=[45], freq=[0.495,0.505])
sim.initGenotype(Naatamo,subPops=[0] , loci=[46], freq=[0.719,0.281])
sim.initGenotype(Naatamo,subPops=[0] , loci=[47], freq=[0.525,0.475])
sim.initGenotype(Naatamo,subPops=[0] , loci=[48], freq=[0.551,0.449])
sim.initGenotype(Naatamo,subPops=[0] , loci=[49], freq=[0.571,0.429])
sim.initGenotype(Naatamo,subPops=[0] , loci=[50], freq=[0.424,0.576])
sim.initGenotype(Naatamo,subPops=[0] , loci=[51], freq=[0.455,0.545])
sim.initGenotype(Naatamo,subPops=[0] , loci=[52], freq=[0.081,0.919])
sim.initGenotype(Naatamo,subPops=[0] , loci=[53], freq=[0.479,0.521])
sim.initGenotype(Naatamo,subPops=[0] , loci=[54], freq=[0.413,0.587])
sim.initGenotype(Naatamo,subPops=[0] , loci=[55], freq=[0.48,0.52])
sim.initGenotype(Naatamo,subPops=[0] , loci=[56], freq=[0.418,0.582])
sim.initGenotype(Naatamo,subPops=[0] , loci=[57], freq=[0.434,0.566])
sim.initGenotype(Naatamo,subPops=[0] , loci=[58], freq=[0.429,0.571])
sim.initGenotype(Naatamo,subPops=[0] , loci=[59], freq=[0.404,0.596])
sim.initGenotype(Naatamo,subPops=[0] , loci=[60], freq=[0.546,0.454])
sim.initGenotype(Naatamo,subPops=[0] , loci=[61], freq=[0.561,0.439])
sim.initGenotype(Naatamo,subPops=[0] , loci=[62], freq=[0.485,0.515])
sim.initGenotype(Naatamo,subPops=[0] , loci=[63], freq=[0.449,0.551])
sim.initGenotype(Naatamo,subPops=[0] , loci=[64], freq=[0.722,0.278])
sim.initGenotype(Naatamo,subPops=[0] , loci=[65], freq=[0.74,0.26])
sim.initGenotype(Naatamo,subPops=[0] , loci=[66], freq=[0.5,0.5])
sim.initGenotype(Naatamo,subPops=[0] , loci=[67], freq=[0.5,0.5])
sim.initGenotype(Naatamo,subPops=[0] , loci=[68], freq=[0.455,0.545])
sim.initGenotype(Naatamo,subPops=[0] , loci=[69], freq=[0.51,0.49])
sim.initGenotype(Naatamo,subPops=[0] , loci=[70], freq=[0.366,0.634])
sim.initGenotype(Naatamo,subPops=[0] , loci=[71], freq=[0.587,0.413])
sim.initGenotype(Naatamo,subPops=[0] , loci=[72], freq=[0.212,0.788])
sim.initGenotype(Naatamo,subPops=[0] , loci=[73], freq=[0.379,0.621])
sim.initGenotype(Naatamo,subPops=[0] , loci=[74], freq=[0.403,0.597])
sim.initGenotype(Naatamo,subPops=[0] , loci=[75], freq=[0.412,0.588])
sim.initGenotype(Naatamo,subPops=[0] , loci=[76], freq=[0.521,0.479])
sim.initGenotype(Naatamo,subPops=[0] , loci=[77], freq=[0.357,0.643])
sim.initGenotype(Naatamo,subPops=[0] , loci=[78], freq=[0.546,0.454])
sim.initGenotype(Naatamo,subPops=[0] , loci=[79], freq=[0.51,0.49])
sim.initGenotype(Naatamo,subPops=[0] , loci=[80], freq=[0.48,0.52])
sim.initGenotype(Naatamo,subPops=[0] , loci=[81], freq=[0.54,0.46])
sim.initGenotype(Naatamo,subPops=[0] , loci=[82], freq=[0.591,0.409])
sim.initGenotype(Naatamo,subPops=[0] , loci=[83], freq=[0.419,0.581])
sim.initGenotype(Naatamo,subPops=[0] , loci=[84], freq=[0.485,0.515])
sim.initGenotype(Naatamo,subPops=[0] , loci=[85], freq=[0.611,0.389])
sim.initGenotype(Naatamo,subPops=[0] , loci=[86], freq=[0.556,0.444])
sim.initGenotype(Naatamo,subPops=[0] , loci=[87], freq=[0.444,0.556])
sim.initGenotype(Naatamo,subPops=[0] , loci=[88], freq=[0.5,0.5])
sim.initGenotype(Naatamo,subPops=[0] , loci=[89], freq=[0.516,0.484])
sim.initGenotype(Naatamo,subPops=[0] , loci=[90], freq=[0.444,0.556])
sim.initGenotype(Naatamo,subPops=[0] , loci=[91], freq=[0.546,0.454])
sim.initGenotype(Naatamo,subPops=[0] , loci=[92], freq=[0.612,0.388])
sim.initGenotype(Naatamo,subPops=[0] , loci=[93], freq=[0.576,0.424])
sim.initGenotype(Naatamo,subPops=[0] , loci=[94], freq=[0.696,0.304])
sim.initGenotype(Naatamo,subPops=[0] , loci=[95], freq=[0.561,0.439])
sim.initGenotype(Naatamo,subPops=[0] , loci=[96], freq=[0.288,0.712])
sim.initGenotype(Naatamo,subPops=[0] , loci=[97], freq=[0.47,0.53])
sim.initGenotype(Naatamo,subPops=[0] , loci=[98], freq=[0.273,0.727])
sim.initGenotype(Naatamo,subPops=[0] , loci=[99], freq=[0.429,0.571])
sim.initGenotype(Naatamo,subPops=[0] , loci=[100], freq=[0.466,0.534])
sim.initGenotype(Naatamo,subPops=[0] , loci=[101], freq=[0.51,0.49])
sim.initGenotype(Naatamo,subPops=[0] , loci=[102], freq=[0.419,0.581])
sim.initGenotype(Naatamo,subPops=[0] , loci=[103], freq=[0.389,0.611])
sim.initGenotype(Naatamo,subPops=[0] , loci=[104], freq=[0.409,0.591])
sim.initGenotype(Naatamo,subPops=[0] , loci=[105], freq=[0.571,0.429])
sim.initGenotype(Naatamo,subPops=[0] , loci=[106], freq=[0.49,0.51])
sim.initGenotype(Naatamo,subPops=[0] , loci=[107], freq=[0.52,0.48])
sim.initGenotype(Naatamo,subPops=[0] , loci=[108], freq=[0.719,0.281])
sim.initGenotype(Naatamo,subPops=[0] , loci=[109], freq=[0.561,0.439])
sim.initGenotype(Naatamo,subPops=[0] , loci=[110], freq=[0.526,0.474])
sim.initGenotype(Naatamo,subPops=[0] , loci=[111], freq=[0.609,0.391])
sim.initGenotype(Naatamo,subPops=[0] , loci=[112], freq=[0.229,0.771])
sim.initGenotype(Naatamo,subPops=[0] , loci=[113], freq=[0.657,0.343])
sim.initGenotype(Naatamo,subPops=[0] , loci=[114], freq=[0.612,0.388])
sim.initGenotype(Naatamo,subPops=[0] , loci=[115], freq=[0.934,0.066])
sim.initGenotype(Naatamo,subPops=[0] , loci=[116], freq=[0.526,0.474])
sim.initGenotype(Naatamo,subPops=[0] , loci=[117], freq=[0.46,0.54])
sim.initGenotype(Naatamo,subPops=[0] , loci=[118], freq=[0.475,0.525])
sim.initGenotype(Naatamo,subPops=[0] , loci=[119], freq=[0.516,0.484])
sim.initGenotype(Naatamo,subPops=[0] , loci=[120], freq=[0.566,0.434])
sim.initGenotype(Naatamo,subPops=[0] , loci=[121], freq=[0.542,0.458])
sim.initGenotype(Naatamo,subPops=[0] , loci=[122], freq=[0.268,0.732])
sim.initGenotype(Naatamo,subPops=[0] , loci=[123], freq=[0.571,0.429])
sim.initGenotype(Naatamo,subPops=[0] , loci=[124], freq=[0.515,0.485])
sim.initGenotype(Naatamo,subPops=[0] , loci=[125], freq=[0.52,0.48])
sim.initGenotype(Naatamo,subPops=[0] , loci=[126], freq=[0.459,0.541])
sim.initGenotype(Naatamo,subPops=[0] , loci=[127], freq=[0.495,0.505])
sim.initGenotype(Naatamo,subPops=[0] , loci=[128], freq=[0.495,0.505])
sim.initGenotype(Naatamo,subPops=[0] , loci=[129], freq=[0.771,0.229])
sim.initGenotype(Naatamo,subPops=[0] , loci=[130], freq=[0.475,0.525])
sim.initGenotype(Naatamo,subPops=[0] , loci=[131], freq=[0.427,0.573])
sim.initGenotype(Naatamo,subPops=[0] , loci=[132], freq=[0.613,0.387])
sim.initGenotype(Naatamo,subPops=[0] , loci=[133], freq=[0.479,0.521])
sim.initGenotype(Naatamo,subPops=[0] , loci=[134], freq=[0.429,0.571])
sim.initGenotype(Naatamo,subPops=[0] , loci=[135], freq=[0.46,0.54])
sim.initGenotype(Naatamo,subPops=[0] , loci=[136], freq=[0.576,0.424])
sim.initGenotype(Naatamo,subPops=[0] , loci=[137], freq=[0.49,0.51])
sim.initGenotype(Naatamo,subPops=[0] , loci=[138], freq=[0.304,0.696])
sim.initGenotype(Naatamo,subPops=[0] , loci=[139], freq=[0.647,0.353])
sim.initGenotype(Naatamo,subPops=[0] , loci=[140], freq=[0.939,0.061])
sim.initGenotype(Naatamo,subPops=[0] , loci=[141], freq=[0.52,0.48])
sim.initGenotype(Naatamo,subPops=[0] , loci=[142], freq=[0.583,0.417])
sim.initGenotype(Naatamo,subPops=[0] , loci=[143], freq=[0.4,0.6])
sim.initGenotype(Naatamo,subPops=[0] , loci=[144], freq=[0.571,0.429])
sim.initGenotype(Naatamo,subPops=[0] , loci=[145], freq=[0.394,0.606])
sim.initGenotype(Naatamo,subPops=[0] , loci=[146], freq=[0.369,0.631])
sim.initGenotype(Naatamo,subPops=[0] , loci=[147], freq=[0.51,0.49])
sim.initGenotype(Naatamo,subPops=[0] , loci=[148], freq=[0.552,0.448])
sim.initGenotype(Naatamo,subPops=[0] , loci=[149], freq=[0.546,0.454])
sim.initGenotype(Naatamo,subPops=[0] , loci=[150], freq=[0.469,0.531])
sim.initGenotype(Naatamo,subPops=[0] , loci=[151], freq=[0.248,0.752])
sim.initGenotype(Naatamo,subPops=[0] , loci=[152], freq=[0.429,0.571])
sim.initGenotype(Naatamo,subPops=[0] , loci=[153], freq=[0.571,0.429])
sim.initGenotype(Naatamo,subPops=[0] , loci=[154], freq=[0.432,0.568])
sim.initGenotype(Naatamo,subPops=[0] , loci=[155], freq=[0.53,0.47])
sim.initGenotype(Naatamo,subPops=[0] , loci=[156], freq=[0.59,0.41])
sim.initGenotype(Naatamo,subPops=[0] , loci=[157], freq=[0.588,0.412])
sim.initGenotype(Naatamo,subPops=[0] , loci=[158], freq=[0.313,0.687])
sim.initGenotype(Naatamo,subPops=[0] , loci=[159], freq=[0.35,0.65])
sim.initGenotype(Naatamo,subPops=[0] , loci=[160], freq=[0.404,0.596])
sim.initGenotype(Naatamo,subPops=[0] , loci=[161], freq=[0.712,0.288])
sim.initGenotype(Naatamo,subPops=[0] , loci=[162], freq=[0.438,0.562])
sim.initGenotype(Naatamo,subPops=[0] , loci=[163], freq=[0.51,0.49])
sim.initGenotype(Naatamo,subPops=[0] , loci=[164], freq=[0.51,0.49])
sim.initGenotype(Naatamo,subPops=[0] , loci=[165], freq=[0.68,0.32])
sim.initGenotype(Naatamo,subPops=[0] , loci=[166], freq=[0.364,0.636])
sim.initGenotype(Naatamo,subPops=[0] , loci=[167], freq=[0.415,0.585])
sim.initGenotype(Naatamo,subPops=[0] , loci=[168], freq=[0.97,0.03])
sim.initGenotype(Naatamo,subPops=[0] , loci=[169], freq=[0.516,0.484])
sim.initGenotype(Naatamo,subPops=[0] , loci=[170], freq=[0.755,0.245])
sim.initGenotype(Naatamo,subPops=[0] , loci=[171], freq=[0.443,0.557])
sim.initGenotype(Naatamo,subPops=[0] , loci=[172], freq=[0.397,0.603])
sim.initGenotype(Naatamo,subPops=[0] , loci=[173], freq=[0.434,0.566])
sim.initGenotype(Naatamo,subPops=[0] , loci=[174], freq=[0.525,0.475])
sim.initGenotype(Naatamo,subPops=[0] , loci=[175], freq=[0.495,0.505])
sim.initGenotype(Naatamo,subPops=[0] , loci=[176], freq=[0.495,0.505])
sim.initGenotype(Naatamo,subPops=[0] , loci=[177], freq=[0.51,0.49])
sim.initGenotype(Naatamo,subPops=[0] , loci=[178], freq=[0.546,0.454])
sim.initGenotype(Naatamo,subPops=[0] , loci=[179], freq=[0.48,0.52])
sim.initGenotype(Naatamo,subPops=[0] , loci=[180], freq=[0.621,0.379])
sim.initGenotype(Naatamo,subPops=[0] , loci=[181], freq=[0.369,0.631])
sim.initGenotype(Naatamo,subPops=[0] , loci=[182], freq=[0.516,0.484])
sim.initGenotype(Naatamo,subPops=[0] , loci=[183], freq=[0.592,0.408])
sim.initGenotype(Naatamo,subPops=[0] , loci=[184], freq=[0.641,0.359])
sim.initGenotype(Naatamo,subPops=[0] , loci=[185], freq=[0.755,0.245])
sim.initGenotype(Naatamo,subPops=[0] , loci=[186], freq=[0.459,0.541])
sim.initGenotype(Naatamo,subPops=[0] , loci=[187], freq=[0.454,0.546])
sim.initGenotype(Naatamo,subPops=[0] , loci=[188], freq=[0.374,0.626])
sim.initGenotype(Naatamo,subPops=[0] , loci=[189], freq=[0.505,0.495])
sim.initGenotype(Naatamo,subPops=[0] , loci=[190], freq=[0.341,0.659])
sim.initGenotype(Naatamo,subPops=[0] , loci=[191], freq=[0.561,0.439])
sim.initGenotype(Naatamo,subPops=[0] , loci=[192], freq=[0.418,0.582])
sim.initGenotype(Naatamo,subPops=[0] , loci=[193], freq=[0.505,0.495])
sim.initGenotype(Naatamo,subPops=[0] , loci=[194], freq=[0.542,0.458])
sim.initGenotype(Naatamo,subPops=[0] , loci=[195], freq=[0.663,0.337])
sim.initGenotype(Naatamo,subPops=[0] , loci=[196], freq=[0.628,0.372])
sim.initGenotype(Naatamo,subPops=[0] , loci=[197], freq=[0.369,0.631])
sim.initGenotype(Naatamo,subPops=[0] , loci=[198], freq=[0.429,0.571])
sim.initGenotype(Naatamo,subPops=[0] , loci=[199], freq=[0.537,0.463])
sim.initGenotype(Naatamo,subPops=[1] , loci=[0], freq=[0.089,0.911])
sim.initGenotype(Naatamo,subPops=[1] , loci=[1], freq=[0.996,0.004])
sim.initGenotype(Naatamo,subPops=[1] , loci=[2], freq=[0.034,0.966])
sim.initGenotype(Naatamo,subPops=[1] , loci=[3], freq=[0.206,0.794])
sim.initGenotype(Naatamo,subPops=[1] , loci=[4], freq=[0.688,0.312])
sim.initGenotype(Naatamo,subPops=[1] , loci=[5], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[6], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[7], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[8], freq=[0.835,0.165])
sim.initGenotype(Naatamo,subPops=[1] , loci=[9], freq=[0.991,0.009])
sim.initGenotype(Naatamo,subPops=[1] , loci=[10], freq=[0.839,0.161])
sim.initGenotype(Naatamo,subPops=[1] , loci=[11], freq=[0.758,0.242])
sim.initGenotype(Naatamo,subPops=[1] , loci=[12], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[13], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[14], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[15], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[16], freq=[0.902,0.098])
sim.initGenotype(Naatamo,subPops=[1] , loci=[17], freq=[0.803,0.197])
sim.initGenotype(Naatamo,subPops=[1] , loci=[18], freq=[0.216,0.784])
sim.initGenotype(Naatamo,subPops=[1] , loci=[19], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[20], freq=[0.97,0.03])
sim.initGenotype(Naatamo,subPops=[1] , loci=[21], freq=[0.078,0.922])
sim.initGenotype(Naatamo,subPops=[1] , loci=[22], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[23], freq=[0.013,0.987])
sim.initGenotype(Naatamo,subPops=[1] , loci=[24], freq=[0.553,0.447])
sim.initGenotype(Naatamo,subPops=[1] , loci=[25], freq=[0.629,0.371])
sim.initGenotype(Naatamo,subPops=[1] , loci=[26], freq=[0.961,0.039])
sim.initGenotype(Naatamo,subPops=[1] , loci=[27], freq=[0.844,0.156])
sim.initGenotype(Naatamo,subPops=[1] , loci=[28], freq=[0.014,0.986])
sim.initGenotype(Naatamo,subPops=[1] , loci=[29], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[30], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[31], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[32], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[33], freq=[0.685,0.315])
sim.initGenotype(Naatamo,subPops=[1] , loci=[34], freq=[0.923,0.077])
sim.initGenotype(Naatamo,subPops=[1] , loci=[35], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[36], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[37], freq=[0.741,0.259])
sim.initGenotype(Naatamo,subPops=[1] , loci=[38], freq=[0.076,0.924])
sim.initGenotype(Naatamo,subPops=[1] , loci=[39], freq=[0.201,0.799])
sim.initGenotype(Naatamo,subPops=[1] , loci=[40], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[41], freq=[0.858,0.142])
sim.initGenotype(Naatamo,subPops=[1] , loci=[42], freq=[0.885,0.115])
sim.initGenotype(Naatamo,subPops=[1] , loci=[43], freq=[0.029,0.971])
sim.initGenotype(Naatamo,subPops=[1] , loci=[44], freq=[0.816,0.184])
sim.initGenotype(Naatamo,subPops=[1] , loci=[45], freq=[0.924,0.076])
sim.initGenotype(Naatamo,subPops=[1] , loci=[46], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[47], freq=[0.994,0.006])
sim.initGenotype(Naatamo,subPops=[1] , loci=[48], freq=[0.873,0.127])
sim.initGenotype(Naatamo,subPops=[1] , loci=[49], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[50], freq=[0.85,0.15])
sim.initGenotype(Naatamo,subPops=[1] , loci=[51], freq=[0.918,0.082])
sim.initGenotype(Naatamo,subPops=[1] , loci=[52], freq=[0.432,0.568])
sim.initGenotype(Naatamo,subPops=[1] , loci=[53], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[54], freq=[0.858,0.142])
sim.initGenotype(Naatamo,subPops=[1] , loci=[55], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[56], freq=[0.014,0.986])
sim.initGenotype(Naatamo,subPops=[1] , loci=[57], freq=[0.946,0.054])
sim.initGenotype(Naatamo,subPops=[1] , loci=[58], freq=[0.13,0.87])
sim.initGenotype(Naatamo,subPops=[1] , loci=[59], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[60], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[61], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[62], freq=[0.013,0.987])
sim.initGenotype(Naatamo,subPops=[1] , loci=[63], freq=[0.862,0.138])
sim.initGenotype(Naatamo,subPops=[1] , loci=[64], freq=[0.312,0.688])
sim.initGenotype(Naatamo,subPops=[1] , loci=[65], freq=[0.508,0.492])
sim.initGenotype(Naatamo,subPops=[1] , loci=[66], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[67], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[68], freq=[0.945,0.055])
sim.initGenotype(Naatamo,subPops=[1] , loci=[69], freq=[0.092,0.908])
sim.initGenotype(Naatamo,subPops=[1] , loci=[70], freq=[0.658,0.342])
sim.initGenotype(Naatamo,subPops=[1] , loci=[71], freq=[0.077,0.923])
sim.initGenotype(Naatamo,subPops=[1] , loci=[72], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[73], freq=[0.886,0.114])
sim.initGenotype(Naatamo,subPops=[1] , loci=[74], freq=[0.979,0.021])
sim.initGenotype(Naatamo,subPops=[1] , loci=[75], freq=[0.991,0.009])
sim.initGenotype(Naatamo,subPops=[1] , loci=[76], freq=[0.934,0.066])
sim.initGenotype(Naatamo,subPops=[1] , loci=[77], freq=[0.943,0.057])
sim.initGenotype(Naatamo,subPops=[1] , loci=[78], freq=[0.997,0.003])
sim.initGenotype(Naatamo,subPops=[1] , loci=[79], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[80], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[81], freq=[0.323,0.677])
sim.initGenotype(Naatamo,subPops=[1] , loci=[82], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[83], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[84], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[85], freq=[0.982,0.018])
sim.initGenotype(Naatamo,subPops=[1] , loci=[86], freq=[0.903,0.097])
sim.initGenotype(Naatamo,subPops=[1] , loci=[87], freq=[0.967,0.033])
sim.initGenotype(Naatamo,subPops=[1] , loci=[88], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[89], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[90], freq=[0.931,0.069])
sim.initGenotype(Naatamo,subPops=[1] , loci=[91], freq=[0.955,0.045])
sim.initGenotype(Naatamo,subPops=[1] , loci=[92], freq=[0.965,0.035])
sim.initGenotype(Naatamo,subPops=[1] , loci=[93], freq=[0.054,0.946])
sim.initGenotype(Naatamo,subPops=[1] , loci=[94], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[95], freq=[0.013,0.987])
sim.initGenotype(Naatamo,subPops=[1] , loci=[96], freq=[0.817,0.183])
sim.initGenotype(Naatamo,subPops=[1] , loci=[97], freq=[0.839,0.161])
sim.initGenotype(Naatamo,subPops=[1] , loci=[98], freq=[0.881,0.119])
sim.initGenotype(Naatamo,subPops=[1] , loci=[99], freq=[0.887,0.113])
sim.initGenotype(Naatamo,subPops=[1] , loci=[100], freq=[0.727,0.273])
sim.initGenotype(Naatamo,subPops=[1] , loci=[101], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[102], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[103], freq=[0.847,0.153])
sim.initGenotype(Naatamo,subPops=[1] , loci=[104], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[105], freq=[0.088,0.912])
sim.initGenotype(Naatamo,subPops=[1] , loci=[106], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[107], freq=[0.322,0.678])
sim.initGenotype(Naatamo,subPops=[1] , loci=[108], freq=[0.326,0.674])
sim.initGenotype(Naatamo,subPops=[1] , loci=[109], freq=[0.087,0.913])
sim.initGenotype(Naatamo,subPops=[1] , loci=[110], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[111], freq=[0.983,0.017])
sim.initGenotype(Naatamo,subPops=[1] , loci=[112], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[113], freq=[0.155,0.845])
sim.initGenotype(Naatamo,subPops=[1] , loci=[114], freq=[0.995,0.005])
sim.initGenotype(Naatamo,subPops=[1] , loci=[115], freq=[0.446,0.554])
sim.initGenotype(Naatamo,subPops=[1] , loci=[116], freq=[0.955,0.045])
sim.initGenotype(Naatamo,subPops=[1] , loci=[117], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[118], freq=[0.944,0.056])
sim.initGenotype(Naatamo,subPops=[1] , loci=[119], freq=[0.003,0.997])
sim.initGenotype(Naatamo,subPops=[1] , loci=[120], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[121], freq=[0.198,0.802])
sim.initGenotype(Naatamo,subPops=[1] , loci=[122], freq=[0.698,0.302])
sim.initGenotype(Naatamo,subPops=[1] , loci=[123], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[124], freq=[0.02,0.98])
sim.initGenotype(Naatamo,subPops=[1] , loci=[125], freq=[0.039,0.961])
sim.initGenotype(Naatamo,subPops=[1] , loci=[126], freq=[0.964,0.036])
sim.initGenotype(Naatamo,subPops=[1] , loci=[127], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[128], freq=[0.92,0.08])
sim.initGenotype(Naatamo,subPops=[1] , loci=[129], freq=[0.437,0.563])
sim.initGenotype(Naatamo,subPops=[1] , loci=[130], freq=[0.999,0.001])
sim.initGenotype(Naatamo,subPops=[1] , loci=[131], freq=[0.857,0.143])
sim.initGenotype(Naatamo,subPops=[1] , loci=[132], freq=[0.977,0.023])
sim.initGenotype(Naatamo,subPops=[1] , loci=[133], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[134], freq=[0.873,0.127])
sim.initGenotype(Naatamo,subPops=[1] , loci=[135], freq=[0.085,0.915])
sim.initGenotype(Naatamo,subPops=[1] , loci=[136], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[137], freq=[0.128,0.872])
sim.initGenotype(Naatamo,subPops=[1] , loci=[138], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[139], freq=[0.12,0.88])
sim.initGenotype(Naatamo,subPops=[1] , loci=[140], freq=[0.789,0.211])
sim.initGenotype(Naatamo,subPops=[1] , loci=[141], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[142], freq=[0.916,0.084])
sim.initGenotype(Naatamo,subPops=[1] , loci=[143], freq=[0.87,0.13])
sim.initGenotype(Naatamo,subPops=[1] , loci=[144], freq=[0.08,0.92])
sim.initGenotype(Naatamo,subPops=[1] , loci=[145], freq=[0.788,0.212])
sim.initGenotype(Naatamo,subPops=[1] , loci=[146], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[147], freq=[0.178,0.822])
sim.initGenotype(Naatamo,subPops=[1] , loci=[148], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[149], freq=[0.923,0.077])
sim.initGenotype(Naatamo,subPops=[1] , loci=[150], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[151], freq=[0.87,0.13])
sim.initGenotype(Naatamo,subPops=[1] , loci=[152], freq=[0.921,0.079])
sim.initGenotype(Naatamo,subPops=[1] , loci=[153], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[154], freq=[0.956,0.044])
sim.initGenotype(Naatamo,subPops=[1] , loci=[155], freq=[0.859,0.141])
sim.initGenotype(Naatamo,subPops=[1] , loci=[156], freq=[0.092,0.908])
sim.initGenotype(Naatamo,subPops=[1] , loci=[157], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[158], freq=[0.645,0.355])
sim.initGenotype(Naatamo,subPops=[1] , loci=[159], freq=[0.935,0.065])
sim.initGenotype(Naatamo,subPops=[1] , loci=[160], freq=[0.832,0.168])
sim.initGenotype(Naatamo,subPops=[1] , loci=[161], freq=[0.012,0.988])
sim.initGenotype(Naatamo,subPops=[1] , loci=[162], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[163], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[164], freq=[0.981,0.019])
sim.initGenotype(Naatamo,subPops=[1] , loci=[165], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[166], freq=[0.961,0.039])
sim.initGenotype(Naatamo,subPops=[1] , loci=[167], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[168], freq=[0.541,0.459])
sim.initGenotype(Naatamo,subPops=[1] , loci=[169], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[170], freq=[0.617,0.383])
sim.initGenotype(Naatamo,subPops=[1] , loci=[171], freq=[0.931,0.069])
sim.initGenotype(Naatamo,subPops=[1] , loci=[172], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[173], freq=[0.94,0.06])
sim.initGenotype(Naatamo,subPops=[1] , loci=[174], freq=[0.135,0.865])
sim.initGenotype(Naatamo,subPops=[1] , loci=[175], freq=[0.075,0.925])
sim.initGenotype(Naatamo,subPops=[1] , loci=[176], freq=[0.821,0.179])
sim.initGenotype(Naatamo,subPops=[1] , loci=[177], freq=[0.058,0.942])
sim.initGenotype(Naatamo,subPops=[1] , loci=[178], freq=[0.929,0.071])
sim.initGenotype(Naatamo,subPops=[1] , loci=[179], freq=[0.833,0.167])
sim.initGenotype(Naatamo,subPops=[1] , loci=[180], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[181], freq=[0.936,0.064])
sim.initGenotype(Naatamo,subPops=[1] , loci=[182], freq=[0.912,0.088])
sim.initGenotype(Naatamo,subPops=[1] , loci=[183], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[184], freq=[0.237,0.763])
sim.initGenotype(Naatamo,subPops=[1] , loci=[185], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[186], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[187], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[188], freq=[0.059,0.941])
sim.initGenotype(Naatamo,subPops=[1] , loci=[189], freq=[0.113,0.887])
sim.initGenotype(Naatamo,subPops=[1] , loci=[190], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[191], freq=[0.078,0.922])
sim.initGenotype(Naatamo,subPops=[1] , loci=[192], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[193], freq=[0.982,0.018])
sim.initGenotype(Naatamo,subPops=[1] , loci=[194], freq=[0.878,0.122])
sim.initGenotype(Naatamo,subPops=[1] , loci=[195], freq=[0.107,0.893])
sim.initGenotype(Naatamo,subPops=[1] , loci=[196], freq=[0,1])
sim.initGenotype(Naatamo,subPops=[1] , loci=[197], freq=[0.833,0.167])
sim.initGenotype(Naatamo,subPops=[1] , loci=[198], freq=[1,0])
sim.initGenotype(Naatamo,subPops=[1] , loci=[199], freq=[0.924,0.076])



sim.utils.export(Naatamo, format='PED', output='Naatamo_ref.ped')

