-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathCGE_auton_part1.html
233 lines (165 loc) · 103 KB
/
CGE_auton_part1.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<!-- saved from url=(0014)about:internet -->
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<title>CGE Final Exam: Adam Auton</title>
<base target="_blank"/>
<style type="text/css">
body, td {
font-family: sans-serif;
background-color: white;
font-size: 12px;
margin: 8px;
}
tt, code, pre {
font-family: 'DejaVu Sans Mono', 'Droid Sans Mono', 'Lucida Console', Consolas, Monaco, monospace;
}
h1 {
font-size:2.2em;
}
h2 {
font-size:1.8em;
}
h3 {
font-size:1.4em;
}
h4 {
font-size:1.0em;
}
h5 {
font-size:0.9em;
}
h6 {
font-size:0.8em;
}
a:visited {
color: rgb(50%, 0%, 50%);
}
pre {
margin-top: 0;
max-width: 95%;
border: 1px solid #ccc;
}
pre code {
display: block; padding: 0.5em;
}
code.r {
background-color: #F8F8F8;
}
table, td, th {
border: none;
}
blockquote {
color:#666666;
margin:0;
padding-left: 1em;
border-left: 0.5em #EEE solid;
}
hr {
height: 0px;
border-bottom: none;
border-top-width: thin;
border-top-style: dotted;
border-top-color: #999999;
}
@media print {
* {
background: transparent !important;
color: black !important;
filter:none !important;
-ms-filter: none !important;
}
body {
font-size:12pt;
max-width:100%;
}
a, a:visited {
text-decoration: underline;
}
hr {
visibility: hidden;
page-break-before: always;
}
pre, blockquote {
padding-right: 1em;
page-break-inside: avoid;
}
tr, img {
page-break-inside: avoid;
}
img {
max-width: 100% !important;
}
@page :left {
margin: 15mm 20mm 15mm 10mm;
}
@page :right {
margin: 15mm 10mm 15mm 20mm;
}
p, h2, h3 {
orphans: 3; widows: 3;
}
h2, h3 {
page-break-after: avoid;
}
}
</style>
</head>
<body>
<h1>CGE Final Exam: Adam Auton</h1>
<h2>Plotting SNP density</h2>
<p>SNP densities calculated using VCFtools. I used <em>region3</em> for this question. SNP density was calculated for 10kb regions for all populations, CEU, YRI, and CHB populations.</p>
<p>Density information was first loaded into R.</p>
<h3>Density Plot</h3>
<p>Plotting density of rare SNPs (maf <= 0.001) for each population.</p>
<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAGwCAMAAABy/YdMAAAAnFBMVEX9/v0AAAAAADkAAGUAAP8AOY8AZo8AZrUA/wA5AAA5ADk5AGU5OY85ZmU5ZrU5j485j9plAABlADllAGVlOY9lZgBlZjllZmVltbVltf2POQCPOTmPOWWPZgCPjzmPj9qPtY+P27WP29qP2/21ZgC1tWW124+1/rW1/tq1/v3ajzna/rXa/tra/v39tWX924/9/rX9/tr9/v3/AABaVr5AAAAANHRSTlP//////////////////////////////////////////////////////////////////wD/9XEwmgAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3QmDpLh1AOBolkkz3o3Ts+s47raTOJ6Ok+6kr8f//28poes9XUigoqjSeztbB4cQ+kAIFdD/ABxdxj9cOgMclwmG7zQYvtNg+E6D4TsNhu80GL7TYPhOg+E7DYbvNBi+02D4ToPhOw2G7zQYvtNg+E6D4TsNhu80GL7TYPhOg+E7DYbvNBi+02D4ToPhO4194T++C/H1Vb49AHw+/vQsB5ziAU9ymiCTwE/PAE8P+kMs3r+JLz/yKaRnTi3RRHbJSwnRFaMZKSiZtrErvFqZO/l+Wkm0ekgqDz/HU7Y48mO3wS+kvZBQBr6kZNrGrvDv376+qvUV4l6vniyNp9M3GZ+PQvxpHvKiNvbTJP8lhCyc034s309z/8+jLCJVZHoyM1YnIctOjrmXX7/8URfemypGXN4ueTe9HiIzN1vN06p5P70l30cyCGhRbqRbsXhGsiWDU24XO8OrPea0iv/y03MAPxfszPY0f7g3Q+70HjEXDoI3k5mxYOGf9HxmOwCcAoKnC7xzQwi8nvf//CXfhRkEtCg70i0nkZFcyeCUG8au8HMBqP3279/ugwrt/dtp3MtpHeX2P2/2pxnuT3uIHgK64pUVrvxgJzNj55Bj37/Nu9S8k9655ZuyRvDzVvOAp1dD6B4fLhlPjzOIF4VG6hVLZSRXMn7KjWLnVv3TvH5yfZ++/M2unj6QvUkkWdhvQq/0qab+oYvfFowrfjuZGasX8aBSOhX1w5yADjmVD69q7Xs0vR7iwat50ZJx+jiDeFF2pF2xVEZyJUNTbha7n87JSk2t8z/aCk1HDF4X/3ygm4lj8GbsnEoCfp5kJbyZNwbvZRAvyo4k8LGMZEsGr1y72BX+RRbAk1o92Tzy4GlVL8OVqwwzJ6nqbZgWt1fVmxJ7mZtJFVX9KZkXfSg288aqei+DeFF2JKnqYxlZKplNpxPx2LtxZzd2Were6vltLdTEInM+0SaWGzunMpeRa3wZeFU9BPCpxp1N21Ut4ZLvIMggXpQdSRp3sYzkSoasXLvYt6qXK6HPZGQJ+Nu1LKB/1q1ZYU5oXU1q5jx9/vq/qidHF6ceO4faOdzpnKkjT9Pev5jKFJ/OzTkgp2dyiDy0/skc4/W8aMlueppBvCg30q1YPCPZksEpt4uuu2zVET0/5FaD4ReG3Gow/MKQW42u4XsOhu80GL7TYPhOg+E7DYbvNBi+02D4ToPhOw2G7zQYvtNg+E6D4TsNhu80GL7TYPhOg+E7DYbvNBi+02D4TmMLvOA4cpwRfsO8HOcOhu80GL7TYPhOYyv8x2/P5pFMlUlzXDQawEt7eP99bdIcF40G8O+/vKo9vy5pjovGZvjvX/72F7nH/xLU9Qx/5NjeuJufHPUWuYuU4Y8c3KrvNM4AX9gpSOPj1x+n48XPz6d/YB7scr9i4XN8PsqnTj3d2efDqFQ5UBxlj3/713sML1/nbWFdvP/uh0vs5esrwwdxEPjPP//nH149+M9H+fCw+bmFv/175fN9nu7g6R5sYgwfxPZWva7Xw9adS3qIB5r2dFLw9ODBy932pPcmHy319ZU0H8d4oGz9+tdfXoH3+HRs3uM/H1OH4po9/uUEfOcf4x9Uv+CvPz6+P5zqhKpdXj59TCfz9ZWP8WFsr+qTh+IKeP2QXrrH6/eT+LyIp6qHN6o64/Ty9uUHMHwYxzjGz90/Tw8BvNnjT/CVe7yFhxdu1cfiGPAvcm9+uwvg7TH+Tr7VLNzBy3O764EflidpE4eA/5SdvrJBFsDrVv2vf6x9aqeDPx3p7xk+iEPAL8WGM/prC4bHwfDt4yrgO4phL3mGP1YwfKfRHfyb7mKb+9q+/PAa9/Xxhrvs5Gmi7gs8evQG/zL/QQnXqb4VnqZn+gdk7//BozN41Wz/fHxoBO+lZ+Are/8uEbcFP8XDTWq75Qrhl24H89Kze3x4YeDR4rbgF+PN/L1H83vaxj3eS0+Yvw13lj/K2DTq4DdsJceAt7tio6reS89W9Y+Hb90NVZhXDx8c49WA1VVz4hhf+dPuJaIz+KBVL6+dylzjUZmegb+Crt/e4Ol5vPljkKvdr/g8vjd4DhUM32kwfKfB8J0Gw3cacfgUMMPfTFTBb+nfZfhDxRBHZvhbD4bvNBi+0xggqszwtx4M32kwfKfB8J0Gw3caDN9pMHynEYdPATP8zcQAMXy/A998Y/ibiSL4geFvKAb7ugQ/8B5/S8HwnQbDdxoE3j+op+A3XILD8AeJEH5A4wYyIcPfUDB8p7ECPnHVRlkw/EEigB8Yvotw8Iqd4TuJOLwZWgxfvCUw/EEihHewDH/DMThIta8zfB9RCo/69Rj+KiPslaXw4Bp4DH9LwfCdRgA/MHwXkYFHw8wHMlsOvlSe4S8VZ4Ev771n+EtFGh4PMx/IbAx/zVENP7i+PNvbE4Uvsr8c/Ppu5tuIIUCWZiOdhMADw99CBPDzyxgMtMAMf9VhV7sKXo8FDI+O9n4yBcHwO8f54cvO6Bh+51iCH+m0Rpfhrz2S8BDCA8PfTqyCHxn+6oPhGT4YvAZ+IBMz/GEjC38i3gpve30WguF3jiV4GIPJGf4GYtgIP6CpIQ5fVrIMv2+0h0fJuCs4loPh9w0E7+2ssBl+uBL4LuU3wdt5EvDA8IeNAH4w36AS3pwAMvxVRAiP9+EW8MDwRwwffrCCkII//RtHhr/ySMAPZmQ9PHgpAsMfMdwz6mAf+FQpM/y+4cEbQQN/Uo/An7aGMnj/AzD8USILDxbe6Q+F8MMSvFfcDL9vROCRGcO3jENtYBgencotwsOwBE+a9xB+bg2v/+L7T8/BGIaPxHr4AdbBD5GBsB3+8/Fhfn/7+uqPYvhIUHj78frgP357Ju/lSV8dfMVJcjoPS/BCw494hiPCd7THt4d3HxH8LD/iXT6ER+1C8ilM2y7Qz/bmY/zH916O8TXw/kSOez94Org5fDpuE74s2+eHH9z7AN4CLwQvTOQnY3g09gzwA+/xTdPZCh+DcGlqeK/H/trh6bJTJdh8e+gc3n0gsfl07ruu18PW3WHgG6V4Rnj76XrgT+dz94kxftLedb/XCV948fJmeGK/Bt5d2XUmePj49Ud8BMOH3+PwLhD8WARv0CPw9u1M8MnwkvaztwzfCqwtfFG2/GUObnBTeLg5+GBl1scF4IN6wa1UMfxyVT9CAl7nkoy7HLy/9L7gXeEXwkM5fJipED5Y3GHgB/qhHXzhUXk5nUvAK3w55Ibgw+MS+sDw3oeO4BtU0aSwtqZ1XnhNnIRX62LmUvDjLcIPLcSOBe/qMlgJD7cAP3YBb5jXwY8p+PkDgve3MwIfzfRx4GlRHRDeFOR54UVT+AF9pXEUeFRUg37dIIaL/Tzw+TR3ho/1FzF8U3hcey5Ni3PC8C6uCN7Vp1XwgxOJwHuPt9L9dTF4dUKPFjuWwdtVZ/gtyR0HXoOPkIEnx/4rgh+gBfwQh69O1qVj/Bb6GHaCd2V4FfA6Dwvw0dyWxk7wyYQGo47hqVEAP/8cR+DNL3Su71bV8HCt8OZAFYXHjefDwKvMYHhX6SeWXwk/AoIf3UD1huEj3TrxDGP4cLL94U0R+PAkk5XwAS069F0WHlxO9Pa+Gl53zWfgXV1hGkqHhB8X4YuJ8vBhRVgZ+8GLSYC7tjoGT3+zy8OT/HtRBv/5KIS4y09akvQh4Kt//DHVj0lpWII3i83DE/kEvHW+DLy6oPIlvD0uG7Gk3SWAi/AA1w/v1qoUHjA8UPhxd3h1K2zkhthsrIQft8DTY5qy8ltA1aeJq+GhCl6+bIUf28LD07zHV9b1NfCuGNfAD3vDLxyJMvDjNnjYEz5zz0Q28vDD6As3gMet2GFH+CGYoTW815ujPyXh6fDLt+rxvhjA65dDwsMQhR/QFMEMaKVUO09GGh4OCP/xH/L18y8N9ngrejD43AIL4P0lBPC22yoLP26EH20fQBv4uVX/8iVxx0wqLgs/nBEeV+U0A3iG5vD4gA8L8ItrXHoe/0/fUrfIJSNIelyGH9bBu+K17/Znnjr42BLDPd4tJJrkgIZF4Ydt8Pbl7PC6XV8XZfCmYObXwYMvlPd/jDgzvEsmC0/WrAqeDjf7dxSehIYPE46s1p6teh9+Xv/zwJteswj8EE11EX7w5zwevFFvBb82QvhxCd6v6svgB4jDD6aEd4Anaabh1cNJrwT+tMc/5KeKRz18cIzfHz41zrYV0dCN8D6whQc6vAV8uFKFrfpTRV/ZU39peGgGbxmX4QdvDiAZAgw/Bnt2Al5PNY4QwvvGreHneNt+jB8J8IjKjg4/N3xK9xzwZnkWXn8phLdN/ADebQE0De+vlJLiRLH3MV5/PAu8JfXhg3r4DPAo3cFrx4JrYPrw3pl5BN7b4+0p/D7wL19fX0TtkT7YdqEQfmgA74bvCO9O3QvhxwJ4HRXwtCrZAv/x64/Tv/eft1X1MtshPIXIwqe2giz8AE3gkavL6lp42dS7Dvjfnk/7fAt48/lo8MPZ4eeYi8DB+52w0B5efVlf1YsvP962V/UUXu7/CXj03b2uhMdHGDuVn5TfBeRGuAsWF+DdAAUPi/BwXnj7dQP8qsjBz2t+bvi5lKPwQVKRvj9vmbb6GNEsJGMIfiyD9xt3OXgohDf1CIGPxDXADxXwNLkIPNmU0NA4/ICmUPDQEbz8szMr+u4K4Un1qv/iUhU83dHL4fGUgzdzITzZjAaXaDG8qZnL4EktnoQHCp+QL4CXf4SitmEXS3oNvClv8zkWAbybbiu8d7zPw0Mcfp7sAvDjZnh5de3nnyuvwogknYIHsoen4aP96LAJHm1TBN5cIlcETzZMNSc9hxqANOWuC/6xvq5vBx9MRoLuvWXw82TkKN0A3i2fdp4Z+LEIfoquJMOjEfjTsAO8HIxWYqiFL9rj0+dz4OCdrj8FHjXCRvhGF2JcEN6fWiVGjtI+/LAMbycmk9F6QQ/BLfcMPJTDowYcnsJ92A6/NuLwOhtF8GMOnjqvgscHCAqvNVfCe6W9DG+nFvouqkh0AG/P3xrAa4UoPETh1Yhy+MHrB6QT6cQbwJupyYc4PNpIDlDVm2yofGTh9fcxC4/BR2gK7zUA3Lz18CO0gDdDKuBH2Ayv42VjX30RPLkiScFbnXr4cQHefT87/Hit8B/fK/f3EnhdZDH4eeoIPKob9ME4By9HUPiBzA67wMuVRfCq6ZiGj/zBvgAePQHPxkg+NoJ/EZuvq0/D4wIuhbcV9Rng7RW1AbxZCTvnQBfopsrCQwN48uZ/bgW/YncPk7bZMGXhPFLwMIbwbhYHP2bgoRreJbMM71VDbqoaeHTafw74uHwJ/NuK3T1Muhpeizn4kfS4GPgBCDz2Oi68ewSOqdpU+PCCvJnCm4rhR38iP+VUnKFVrxs8aXhTFghelxWFdwfYZXh8iNAfyJAAHmrg/aiGt3MWw+Pkw89t4NdGFN7kpB5+gAL4sQR+JMtcCz82hxfN4O3FfEeAV1kohLc9mQTedo97Lat94XVBjomJIAHvhuwFj4o7jBL4j+938mKM2gZeGh68xvY+8Ja2Cbw9O4ms+Wr4KQcvdod/up9/nWvw8KPW8Jq/CN71we4IPxbB69DwEMALCj8tw0MLePmzrLwCp8XjzuwNQQG8HuHg1a7gwcMSvJrDHk5dE8d1vo8xeHVNt9889yICH6vps/BiEX7C8AK9RuDx7AvwHn8hvHy44dbr6uelY3h3fwWE8JCCdy8Gyp3vm08xeEAbBw6VUi28hVyA1xsdgTfr6UcIL0wlL6ASHi09vtsXVfUP80Nwnho85y4D7woEw9te3Sy86fbfAA9hH+31w/vPSyILy4Zu3Imvr7KFVxcF8OGYDPy4CX6eabBDXRD4YIyff/OyFT524FXwkIQXtfBoxCr4lVEFbwpmoFSmGTDDm0nMSxJ+9OAtUhZ+KIV36cYDLbEGXmh4yMDPc/UFj65ncqdmKXi7p1N4YHiUdDZ2hFffEvDqx04P3pxLYfj5r+8m4fXkhfCx/LsUKuChDt7V9foSPAsPGfii/KClJSef42zwfm5q4O2uPrh5W8FHi88v3TL4cayEN+06Au+uzVDwgjQC0vC49GO5PRC8qhVHwDTr4ccUfFAKQzV8fB+iE42Hgo8tLRu7wNuPo4G3+RyN50CyP5TCh1vYgeGnIvjTRDcC7365SsCPlfC6wMcK+FH/9BeFH6OluwSPFh+Bj8wdh5/6hR/DPZ5U9fMc+gqs88ADyg9EPgbTQxJeZOB112wSHkQEPpY1PWtsBaLjI9EcnmRltB9dafjwTiqWfTuNgwcL75K6DLzZ+LbAT+DD019xbN78BA8K7x77E4OnBZ2FH+Pw8w5XDe/f9GbnkPe0+VVHY3hB4CczcJLfsvAQT/Dw8COmcMZo+nXwUAsPlfBp+RAeCuCnEF7W9CH85JVsb/DkNuBRX2xbBk8Tct3B/t2Odg4JP/ozpg49m+FBwQsFD0vwal2vDR7Lqk/B5hyHdz85am77VJk4vL8tuXDdwXH4UTeoI/Bxebut2QZmCTwgeNk3J6rgg0zQpv+IL9vRy8vGAeBtecfgR/vuwdtF+PBBDXBOeLvZroGf+3ChCD6Sj9uAt8j+iBBe/wxgZyLwwaFkjqEYnh5rUpV9CO8mLYaXjfl5A6Dw87bgPTXDz4Y9TpikjwuPqM4Cjy5JyMJHK00FH/bPB1tCdCUD+FgO1AgfHjS8ML/WCH349x+X4iUWwM8THBx+xLOJhvCkbvHgTWoE3u6ZM/y4M7xQ64/hp1J4Uz2475eHt1EP7xWzhdfj6d1o+tNm+NHA05QBjYusG51iXAE/l8KExmTh/aTg+PC6EHBRFcPDWA8/BsU+evB0brUMXYqx80lbL9Cs4UzqyqIGHjx4dS3e5MGnCvhq4HEXbSm8mz0NbyfMwJsN42zwMFJ4/StNE/hECUfh/UricvBTPTxY22p4dLAFQhOHx0sQYhmebJVxeLNKcfiJws+0kxk1n8EJ3Ymnmm56MoCQNAEvwokycWF4AeXwbk+tgLc5oPDjaA/Qej1EJGUCH/3tdiW8PJEDBD9VwYsQHtrDy7st5H3Ukb9RdTB4QG2sBXg7pa7oYREe0a6DF4DhyR5vyNfDz8eN9vDzrVXvvw9GLcOjppcdjI9cCfhxdJJ2xjy8GeGgzVsOXn1NwruE/K2JrmkZvHuVsQVetRHOD//+yyvgG+uEidx8JfDozy7SfTwP75U7+uQqZCvswwsxkn3ehxdeyng7cpsemWzETYDN8BOeehadvHR2gv/+5W/y78rP+hVJj+qXzovC+3U1zA8R9uHHZXiUqH6nkwXw9sc0AagN4cFbUAFReNP2z8ObI4ivsb1x9/ko7uAtcvP8Evzk4PFgyMHbs8BieLRIVONicvR1XvqYhRcoPZoSfsvB43wXwU8KHlbCCzgPfDKW4MUh4B06lMLr41gUXvcU5OHdZ9uiWwM/dQQPrh4+Hzxkq3oMRVPSxjn4mQTBe+nF4SEFr+6xJBMKQeHFAeEnEcIDhdcMFN68EOJ6eKRVBi80PF25cRxRUnou3LRqCA/+Hi/lo/ACNQAPCa+/oOE+vJ3CTmvhIQoPWfggvRBe58/AWxKBG3d45UYMP38TdOV9eJ1EBTwUwgv7n20IMDyZPAmvF2mtEDzuXqKw9r0KXj3qCFrB6wmE6uabHPx0WHg8PAVPC3c/eIjDx9fOHIDy8CYFDZKAp/0ySXi3hVwFvLyIaAt82DaogycJRuEhBo/W6czwXmNd9eV78EDgRRbe77+5FnjcxwNReDgCvMlpBD5YJWFoozukDw+l8OqfPgzMXxh+G7y4LPy0CG/VwZz76ep/EpH0d4c35ZiAd2U2ZeDhDPDzIYTA4+zIIj81691FT2Tt/KZeDfwU3SEDeGWIr71Jw8/HhclciK3qfj/9y8ALdaFwbIpFePXVn60EPj49hnf73lZ4MigOD3Xw8zxpeKeO4NV0h4DXrdM4PGrsVMOP6XGxiMDrmpzAm9PsOni/7l4FH8v1hODFZBKZXGEZeP2iprk2eFJ8FfBF4bUZ9KfgaEvhoQgeKuGhOfw8bjJ34kym3y/IYybOAD//C+FtU0TPp8rDTnI+eLfzR+EhAQ9+5snXiTyEOg0/NYBXu7f6guB1VeJ+4vHzmI2d4UUpfBiV8FAEL9CnLDxdV3yZzJzwCvhoxOCFhZ+zqA749joODR9eib8zvC5ZCi/smAPA22cM0o2gAbxZdzzRWvjJwqtzdl2G9sZ6Hz5+KW4mzgGvDpwl8BNcCF6I4LIV0QIep6ngp1XwIAL4CWcDMetN4CjwICLw+IdkcUl4+6NZCE8uiHIfPfgpBT+5OmwzvFmymCvGq4AX+nCE4PHPyPoIRuGXllYLj+c08KNQNY3d/jbBu4EI3q20Zp9gC7y6HD8Hr76GTcVjwVt7U5e1gc+vZABv8oNmnGrghep9QZijqUb6hgcNb75l4KEYvjYnLgj8NLkzdzfjNOGbVYNEQ3iYCuGhCt7258XhU6lcB7x+GMQmeEHe81f4F8Dr1lMA7xaDWysO3t5cEIMHDF/oHoGf7GHxCuABcNtJfZwovOlzWAkvMPkCvOu8JfBoP9bw8+kwbtYvwdsho7nMyp55UfjSHd4u3fxM58En5+oGHtbAqyfMTAGn2plCeLEA74ZYeAjgYSM86Noz2nS3c3UGb97OD0+XN5f0FMKLqSG8+o7hU+4HgHeD0P1cRfDVixb2LdJRTcODd4kswEMcXnfJTagXR8FDAt48tbogLPyE4eHQ8JO5DEyN1PCTbsjpHxKEedKXeURzWZGI8CuGX1jL0aRQC49q+HXwOv1zw8d+3N9zj8/Bm9Pn9vAiHJ9KIQ6vDqbmt44pmMstIQevqrXm8Kbf9sDwIoRXNbruA9ed5JPpNG0GP+0Fbw8wEXjIwMNGeL0X5eDDQReGF3F4MbWCVwx7wgvA8CZhtR6b4d2S0WEIdNoHhQfceJ3siZbw4HX3qJl9uUjMbGgA2v1AnZqXZXYU+JyIwGuzJXiz8zE8+R7Cq6M66tJaBU/b7fRcDlrC+6UryIs91gfw5udcD95fUkWYGa4H3pzBWOYIPGyAn4TtJ9FLm0qOF2XwfvcYw8cihHfVm/uDOhMYeDOTLt1FeF2AglaVDl6Yo20iERF82QBvDy1mjgnD22voE/BLXUzpzKPcpN1zCSSi6emcsPCuQ1zBk8vDSekuw4Mg8ObSBJN6KbzaXdPwEINH3KTWQf/78CYzeOl1ZGT264DXe0IE3pZoDh63hyAG7xp1eilxeFdVuDxG4PXy3Pct8G7WDuEhAQ8OflqAnxLwttLVv/ZOSD3ISgX8tAgvgk9F8OFlr+viauBtCVj4CbLwKBkx4RY6OmvTT4MxP+Jj+CmEF/gnAkETC+DxpkbhVQUV/kIHCXizDfUIb9u6CB4svKmUJ7o6GH6CLDyYfjrTrLbLI02pyc4Yg8cH+QmS8MFtzLRtF8Dbo4a9k/GycSR4mKYcvO44w/BCPS/M9P7MP/RMGXjT0hJb4NWQOnhws0xmg71sHAreDojC69NBAq8+efCQhdedwiZBgdKOwQfHeGO6AG/XzRVF5/D6Xw38pNNz8CICL1xnug/vKgdAF+8m4S1QHN61VcjaLsJD9/B+e8s74E05eHxod3eKof5+27Iy+vMkbktxnTQhvK7ESV2NM2bhp2Cc17aznwQdx/Cb4IX+BU8lp1vyDh7B2aubDHzkGRLgw5tju2+7Bh78vxDYL7y5egBNuQRvqmsg8JO77xDtyhjeWk1A4MElGIGHJfhYbcDwuaQnZ5KDp0d823azh3Z1gY6+yUo7CbQVIHj0Zpwx/IThJwQfO4wz/OqkyWUOySkX4Gcw97QHrWQbiN6eivz1xX12qHA3MepHgjnYOHzyC8Th9fYmfPiSH4rPHPvCB6UVHeqVrzl2O/h5iAashLeHANB7PIJHt8mYjSqXz/i6QAAv/C1ZLfjCcR3wQpj2GehbR8zeOdXAhwO2wSf0luAnhi+J+elN6CzO1fA6QXNrahweLTmEn8B1nQsMD9AOfvLgp8VS3yOuAF4d0/UJMT606wQL4c0jCZLw/gV1i/Ap+QA+qIIYviCmKQcPBh6JZndEXJvn4cN0VsGntsQLx1XAz52xFn7y4KEcfgrhTddaeO10AXx8Ud7+zfArw/bQ6USDsvRaYtGHgUIMHvRBohTeS68M/oDoMq4GfnKJ+j2mFfBTBj6oupfgU5MwfLugvbiRIi+CD44QFD45Zyq9MBsQDKtLc7/oDX6i3wH0QyUUfGWeQFcjQV7S3w4TVwsfTEHq8Hgi/omVbjboW6NW1cpTpDY/KDWNK4EPej0zE+wM7zcYr8O9O3g6QL2JDfD0uMPwBUnXxWJpFjSiUpOIaQu8fTGVPsNvmDeIJvCJ4QwfRj/wmzrYJvt6zG66SNwO/IY0xFZ4t98z/N4/QW0pbncb75al2/P6LQntFgwPDB+JPuBhM7xtGTJ8b/DkmvsriJuB3xQNrY76o4wfDC+jJXy7pM4aDN9pMHynwfCdBsN3GgzfaTB8p8HwnQbDdxoM32kwfKfB8J0Gw3caDN9pMHynwfCdBsN3GgzfaWyGf/8mZPz0XJs0x0VjK/zn48P8/vb1tTJpjovGVviP357Ju0zSxMascZwzeI/vNDYf4z++8zH+GoNb9Z0Gw3caDN9pMHynwfCdBsN3GgzfaTB8p3FOeI4jx/ng/e2gXVK3kvaBs83wDL81DryWl0r7wNlmeIbfGgdey0ulfeBsMzzDb40Dr+Wl0j5wthvCc1xTMHynwfCdBsN3GgzfaTB8p8HwnQbDdxoM32kwfKdRB//+8zPAm7qf6kXM7/rt47uQN9fl3irS/nwUX37kU9uadqt8o6Rf5gtfHs6SbZN2s+Kugn+Ty5S5ebkDeFI3U6o3eWvlaWDurTJQqWNqAAACY0lEQVTtt6+vBYmuTrtVvknSAOfLtkq7WXFXwT99+evP6tbJU24+//xjzoJ6kzdRnwbm3irS1rdkFyS6Ou1G+SZJy/R+/XGebOu0WxX3qqr+FKctar6J9sG8vf/yKjOWe6tI+/2Xf5PVcUGiq9Nulm+UtH49S7b1a7viXgP//k2q/G7e/vSbrIVOS8u91aT9bV63gkRXp90s3yhptVPCWbKt025X3Ov2eJOwPc433LzLdpataTfLNymS+ekRZ8k2fjJFk+JeB09bdqe3BgcdlPbHH+a8tzpYxtJulm9SJE/30K5pEku7WbZXwKMa5fMvz/rt8/FetSjTbxVpyxU9bbgFia5Ou1m+cdKq6XWebKu02xX3ij3+dCqpT4TRW6OTVpeaPOVue0JM026Vb5K0qknOk22VdrPi5p67ToPhOw2G7zQYvtNg+E6D4TsNhu80GL7TYPhOg+E7DYbvNBi+02D4ToPhOw2G7zQYvtNg+E6D4TsNhu80GL7T6BFe/bXce3/oz8/q30XytHt0CS9xY/ccdKMOHcPLm0r1ldbqRuT3n//79PXvPz/rwfImO/Fw6byeLbqFl3ehPd3LOxTkrSfybkpb1avB79/uY39Q91aiS/hv+gEGJ3F5p5aq9C28Hiw3gRuu+7uER7edyzuTThuCvI3WwOvBDH9zYW9CVbu2/Cirdt7jbz3cTajzUVweyDG8OcYz/K2F5dSt+ifdqn/+fMSteobnuMFg+E6D4TsNhu80GL7TYPhOg+E7DYbvNBi+02D4ToPhOw2G7zQYvtNg+E6D4TsNhu80GL7TYPhO4/8BCfBG3S4qGb0AAAAASUVORK5CYII=" alt="plot of chunk unnamed-chunk-3"/> </p>
<p>Average SNP densities for each population are shown in Table 1.</p>
<p>Table 1: Average Density per population</p>
<table><thead>
<tr>
<th>Population</th>
<th>Average Density</th>
</tr>
</thead><tbody>
<tr>
<td>All_pops</td>
<td>8.63857142857143</td>
</tr>
<tr>
<td>CEU</td>
<td>3.17238095238095</td>
</tr>
<tr>
<td>YRI</td>
<td>6.27619047619048</td>
</tr>
<tr>
<td>CHB</td>
<td>3.79714285714286</td>
</tr>
</tbody></table>
<h3>Interpretation of density plot</h3>
<p>The average SNP density is different for each population. The YRI population contains the most variation among the three populations while the CEU and CHB populations have similar SNP densities. The mean densities are shown in <em>Table 1</em>.</p>
<p>This may be due to the reduction of genomic diversity as populations migrated out of Africa.</p>
<h2>Frequency spectrum</h2>
<p>Frequency for each population was calculated using VCFtools.</p>
<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAtAAAALQCAMAAACOibeuAAAAqFBMVEX9/v0AAAAAADkAAGUAOTkAOY8AZrU5AAA5ADk5AGU5OQA5OTk5OWU5OY85ZmU5ZrU5j485j7U5j9plAABlADllAGVlOY9lZgBlZjllZmVlZrVltf2POQCPOTmPOWWPZgCPjzmPj7WPtY+P29qP2/21ZgC1Zjm1ZmW1jzm1tWW124+1/rW1/v3ajznaj2XatWXa/rXa/tra/v39tWX924/9/rX9/tr9/v15RaA0AAAAOHRSTlP/////////////////////////////////////////////////////////////////////////ADtcEcoAAAAJcEhZcwAACxIAAAsSAdLdfvwAAB1ASURBVHic7d0JY+LItYbhqBnjSTqeTmB6kg7uyQZObtqaTMz2///ZVVVJYhdScUo6PrzPze1hE8vHRyGVQfxmCxjym6HvACCJQsMUCg1TKDRModAwhULDFAoNUyg0TKHQMIVCwxQKDVMoNEyh0DCFQsMUCg1TKDRModAwhULDFAoNUyg0TOm10KvHh7ftdpnNiv8VR9fTrDA5f9lwkbDYp/1jZy64O/fwcovs8mKWFY+7eOCrx+I/H+arxyLh/MN8d3bCxI6eYXcw3JXi5OI+ND7jMnotdP7d9/ODQo+LCPej3rML2l/ssstPz3rq0r07m+fQHt9k17H2hb45saNn2N345nn0unosbt8XuuEZl9FnoTfPf3weHxfa5V2MJcWDXGY/PBav3mXx/4vRq7tI7l7Om+cse/hvGHGKS60e/aUcf3Z1ddV1VIfczWXZ6Nv0u0f33PqTigW+uEv7Gzi6PjeeTVwX1tOkI0h6VUXPFDpZYuVtHj/Du/vwcVwX2p2WZ6nePPss9OpxlhexnBTa/bd4tMvs4S2v0ysusipe7cUhd3ZxrLxU8e71v+I1764unB2uzqVULhUObcN4s576J8KftHocF++p1Q3UJ/rrK28kD6+kd60ajU8LnS6xur+Hz7Bbq5z40//w+XWv0P6a0gzUfRY6d2HMTlc5cv+yn7jYyv+UI7RbAwspuOTCpXbpl2eHq8vDYLM7tK2ennF4Ot05WXnj/hrKE3fPq7vG8hl516p0/Dr0w9vBKkfixI6fYX8nwusjn+ytcqwek4XcY6HDu9P4eKNw5tMrznAn+v/sRo7R6+Kg0MXBXaHLs8snJfObG7tD2/2nJ5xUP4Xh6QknlteXV++Z40XjGvs70DBCJ07s5Bn2dydsmK6e/l1tFIaRKkvT6R4L7UeMbO8dqtrauzRCu9N3qxwnI3R59m6U2W63u0Pbo/HGnXn49IQTD8eb4qK/e+er0Mfr0H6bIO8lsZNnuLo/fpGP5QhdytOs2vVYaP8I/EryUaHdmlq1Du1fzW41Lrzhraf1CF1eaq/Q4ezdGqG/jurQdu/pKU8qrqHY7JmVN1CfOKleNe7gMku6Dd6HMMsRgvBHR6/+pJSJhSf05Bn2p5SLLLNdoZdh0zTFo++v0Bu/KVdEdFLow1mOYitk9CdfaLfRXWxjuxP+b2+Woyp0eXZ5dWHCc3dou/f0VCct6+v7k990CVOk9RZ8NrYx1eeCqeah/Wpc9faeLjH/hJ4+w2FK3L/xFufujdDLZNPRiv5SuEw7417dSPO4cGXS+w69s8Qo9IE80abKO/bOElNUaOB2FBqmUGiYQqFhCoWGKRQaplBomEKhYQqFhikUGqZQaJhCoWEKhYYpFBqmUGiYQqFhCoWGKRQaplBomEKhYQqFhikUGqZQaJhCoWEKhYYpFBqmUGj9sgND3xvlKLR+2cUjOEGh9aPQHVBo/Sh0BxRaPwrdAYXWj0J3QKH1o9AdUGj9KHQHFFo/Ct0BhdaPQndAofWj0B1QaP0odAcUWj8K3cEtheZDM/2g0B3cVOiLRyCJnDug0PqRcwcUWj9y7oBC60fOHVBo/ci5AwqtHzl3QKH1I+cOGgu99BPMs0tnE7QQcpbTVOjN13nx7+q3rxfOJ2gZ5CyoqdDrn96Kfzc/E3Ra5CyosdA/MnL0gZwFsQ6tADnLYZZDP3LugBFaAXKWwyzH8MhZELMcwyNnQcxyDI+cBbEOrQA5y2GWQz9y7oARWgFylsMsx/DIWVCXWY7jr3cTtAxyFtR9loOgpZGzoO7r0AQtjpzldJ/lIOh+kHMURmgFyFlO91kOgpZGzoK6f5aDoKWRsyBmOYZHzoJYh1aAnOUwy6EVOUe5NkLPWLdLj5zlXJvlWMwIOjVyFnR1luPlnwSdGDkLujoPvfn6e4JOi5wFNa5D+5kkppOSI2c5zHJoRc5RKLRW5ByFQmtFzlEotFbkHIVCa0XOUSi0VuQchUJrRc5RKLRW5ByFQmtFzlEotFbkHIVCa0XOUSi0VuQchUJrRc5RKLRW5ByFQmtFzlEotFbkHIVCa0XOUSi0VuQchUJrRc5RKLRW5ByFXYEpQM5yQqHdniFOf8mU3bxKI+fkfKEXfoR4eDs6j928CiPn9HYj9Cl28yqNnJNjHVoBcpZTFtq9GY4u/fDjIYK+ATmnVq5yfD4bMiOHMHJOLhT6dMvbn8rWtzByTq4coafn3grZ+pZGzsnxo0HDI2dBzHIoQM5ymlY5ziPoeOSc3G6EXv/l5ExGjgTIOam9Qp9MKbH1nQI5J7W3yjE+Pu946zsL6vMJuityTo5ZjuGRsyBmORQgZzlloZengV5C0Dcg59R2n+U43Vgp1u2Klb7jWSaCjkfOyV0p9Mv85DO8BB2PnJNrWuUIQfMZA0HknFrjLMd0/AsjR3rkLKgsdD7ZPJ/ZWFlm4zCntIegb0DOqTWtQ59H0PHIObndB/wJOj1yTq5xo3DqJ/xHbKyIIefUGv9SuHo68617ghZHznKa922Xn9mAIWh55CyGnTVqRc5RKLRW5ByFQmtFzlHqWY7xgk+BpUfOqVV/WPnPJ+ZH0yPn5Ch0n8g5OT7g3ytyTo2NQq3IOQq70+0VOafWuDvdswg6HjknV47Q/yLoXpBzaqxy9IqcU2OVo0/knFw5Qr/Mmy+2h6BvQM6psTvdPpFzcsxDa0XOUSi0VuQchT9994qcU2M3Bn0i5+QodJ/IOTlWOXpFzqmxw3MFyFlOU6H5MZt+kLOgemeNi2xyfB4/2SuOnFOrNgq/Ff93siNufsxGGDknV3+n8OnX061v1u2EkXNy9cdHZ4uTt8LzCPoG5JwasxwKkLOcpnlotr7FkXNqV34Fa8vWtyRyTq6x0Gx9CyPn5Br/9M26nTRyTo3PQ2tFzlEav/XNyCGNnFNr+tY3W9/SyDm5ph3NHG99Z0F9PkF3Rs6pNa1ysPUtjpxTa9zRDOt2wsg5OXY00ytyTq1xRzOMHMLIOTm+sTI8chbUVGg+Y9APchbU+OP1bH1LI+fUGvfLwbqdMHJOjh3N9Imck2NHM70i59QaNwqn/q1wxLpdWuQsqFzl8NvZJ1ZPZ04m6HjknFzzHvzzM2+PBB2PnJPjA/5akXMUCq0VOUfhdwp7Rc6p8TuFfSLn5Php5F6Rc2qu0OfnQS8h6Fjk3AM2CrUi5ygUWityjuILvfp+vs1bvhMSdDxyTs8VOnxjYvlw9s+yJwg6Fjn3wG8Unv3GxCUEHYuce+BXOfxO5XNGjtTIOb2wUej+gDVuuQRBxyPn5Jjl0Iqco1Borcg5CoXWipyjUGityDkKhdaKnKNQaK3IOQqF1oqco1Borcg5CoXWipyjUGityDkKhdaKnKM0Fpq9YvaDnOWwB//hkbMg9uA/PHIW1Fho9izfC3IWxDq0AuQsh1kOrcg5CiO0AuQsh1mO4ZGzIGY5hkfOgpjlGB45C2IdWgFyltNlliMLdscPzpS6RyDnGzBCK0DOcpjlGB45C2KWY3jkLIhZjuGRsyDWoRUgZzl8lkMrco5yZR16PT35jRuClkbOgq4U+mV+8oPrBC2NnAVdLzRb36mRs6DGQk/HvzBypEfOgpo3CpfZOMwp7SFoeeQshlkOrcg5CoXWipyjUGityDkKhdaKnKNQaK3IOQqF1oqco1Borcg5CoXWipyjUGityDkKhdaKnKNQaK3IOQqF1oqco1Borcg5CoXWipyjUGityDkKhdaKnKNQaK3IOQqF1oqco1Borcg5CoXWipyjUGityDkKhdaKnKNQaK3IOYpcoQ/dfs/uHYWOkmiEJvabUegoFForCh2FQmtFoaNQaK3YVolCobUi5ygUWityjkKhtSLnKKkKzYrercg5Sh8jNEHHIOcovReaMaUlCh2ll0Jf7HDWoPMdM4ZCO507keiXZBtG6JiF9BF98Unl3HbguPHu9qlzJ5oKvfka/RvUAoW+Xdura4jgIskXX5qcLz/EgVYB2z45Dfe1xa1c+Z3C4t/d7+cd32BEx4xqETQ5C2gRZmOhfzw3ckAaOQvqvg4NceQsp/ssB6AYhYYpFBqmUGiYQqFhCoWGKRQaptxS6CH/ZqSL2NNBzo1ahHVToVnohoUSX/29LkShh1oo8dXf60IUeqiFEl/9vS5EoYdaKPHV3+tCFHqohRJf/b0uRKGHWijx1d/rQsxDwxQKDVMoNEyh0DCFQsMUCg1TKDRModAwJarQi+zD/PBAh4U2z1n28Nb1lrbbfNL57uVZNu68UJfHtF1/fu2+UGvkXGmfc0yhlw9vSx9VfaDTQsVjX7TcB8XeDSyztkHfdktdHlNxp0avR/dTEDnXS7XPOabQxWs4vGLqA10W8kdaPvzdQqun/7YNenf3vrQfOeqFVk9vy7YLbf5e7r+rSxDtkXOpS85RhZ6VO2KrD3RZaBte0t0WKh7GsnXQ1UKLcYentL57y/Zv1PUO6boE0R4519rnPMQIvWi9wrW7JfcFnJZJ7xYqQm779OwWGnfZy5zqEfouc+5/HXrz3H4nbvs30HrkOFi3azty1At162Y12ihch77TnONnOdxOYDtvfRcLLToMArtb6hD0bqG8/Q3VC7nJgfZNcEF3DaI9cq60z5l5aJhCoWEKhYYpFBqmUGiYQqFhCoWGKRQaplBomEKhYQqFhikUGqZQaJhCoWEKhYYpFBqmUGiYorjQ/itn68+v9Td1Dr6ys0jxFZG7ZCtnxYVefJlcDjrBd1Lvla2c9RZ6/dOvRZhl0Av3rbX6wNZ9NW30t4+j13C0OPLn+erTdvsyDxdYfXx0X5LPs+zhr7OtOweXGMtZb6HzmftfCDp8r7g+sPUjRz4uv29c/Lv+sQw6nLJ6ett8nRf/bvMfPrX+RvJ9Mpaz2kKvp5nbN1sIOuwuoj6wDUFPyt1IuCBfyqDDKWXo7oKbn7+9t3fNXlnLWW2h/V5/FrNy5PBvf/u79QlBh6PFRcM4UYwf4ZQQtBs5luP8S+uv2N8jazlrLXSR3NbtWOTXet0um9UHtlXQ5dFF9uHjfPPs1vDCKeV6XjGMjF6L9Id+MIqZy1lrobt6uZimGz8gRX3O5gudv7N5VOXU52yl0IBHoWEKhYYpFBqmUGiYQqFhCoWGKRQaplBomEKhYQqFhikUGqZQaJhCoWEKhYYpFBqmUGiYQqFhCoWGKRQaplBomEKhYQqFhikUGqZQaJhCoWEKhYYpFBqmUGiYQqFhCoWGKRQaplBomEKhYUqfhQ6/PLN6nPgfMFg9Vj8dVllms2U22x05XPSmH8FbPT68HdzAsrwv7mcjP8z9b5tlKn7FScDm+cN8u55++EeZ8EGU1kPur9Cb5/C4doWeHP0yx+Ws19OH236RJv/u+/n+DRQ3vnkeva4ei9v3WY/f369aX7bM3E8LTqqEWxbaRsj9FboK76DQ5Yl5NZIs/Ss6Kx717pB7LWTZ6Nv0u0c/rruTigW+uEtPtovRa3m51eMPj24AcEP/xL181tNJef2b5z8+j4+yDv/5OK6zdqfl2Y2jlApFjb5NXZMm+4/5PkLur9DVYHxQaJdUcVLxyi4OhSjc6eFIeWgbBo/11KfqT1o9jos3yCrr+sSH/xUDgsutWDivri/cpju+l3Xx9jfxp//h8+te1v6aDAzUy93K3aIaoe8k5P4KHVILr+0sO1qHXrpTQhR5GDl2h7ZV1uPwqnDnFKcu66zLE3dPUrgZH2+Qu2ditv92W9x6eOryyd674epxpOHnfW9WjLbjKulZ/d54FyEPPEIvy9fw6HVRZ535ou8ObfezDifVz0fIOpxYZp1Xb4Djxbi86fDWOD5cf8z9i2qyevp3tb0yC0ObhU77h1k8Oj/iliP0fYQ88Dq0D9B3vX43LIu/O7Q9GjzcmYdZhxMPB4/ior+rtqjDULX/9ljegl/kYzl4lHIDK9F1oYsD4zL5Owm571mO8PiORuji1PW0Gjz8ZoM/Uh7a7mVdnlSs3hXvqrNyCKpP9Fm7y4Wnsn4z9fHl/trrDaRZtYi7XJX18niS672qC+2m8OoQ7iHkHueh3Wb0yTz0uDrnSx3Fsl7xW1abw3XW1UlLv0CxzTL6k39Pra643Bx3V+tnoXxwbu7IPU3jvcFj4d9nq6e8HjyWRqaj60K7ITo85jsJ+f3+pfDai3z/HQ6R3l/IZgudm9i4G9r7C/n9Fho4g0LDFAoNUyg0TKHQMIVCwxQKDVMoNEyh0DCFQsMUCg1TKDRModAwhULDFAoNUyg0TKHQMIVCwxQKDVMoNEyh0DCFQsMUCg1TKDRModAwhULDFAoNUyg0TKHQMIVCwxQKDVMoNEyh0DCFQsOUWwqdHRC7S0C8mwp98QgwEAoNUyg0TKHQMIVCwxQKDVMoNEyh0DCFQsMUCg1TKDRModAwpbHQS/+po9mlsym0kCs5o4OmQm++zot/V799vXA+hZZxLWd00FTo9U9vxb+bnyl0WtdyRgeNhf6REboP13JGB6xDK8A6tBxmOWAKI7QCjNBymOUYHrMcgpjlGN7VnPkycnvMcgyPnAWxDq0AOcthlkM/cu6AEVoBcpbDLMfwyFlQl1mO481sgpZBzoK6z3IQtDRyFtR9HZqgxZGznO6zHATdD3KOwgitADnL6T7LQdDSyFlQ989yELQ0chbELMfwyFkQ69AKkLMcZjm0Iuco10boGet26ZGznGuzHIsZQadGzoKuznK8/JOgEyNnQVfnoTdff0/QaZGzoMZ1aD+TxHRScuQsh1kOrcg5CoXWipyjUGityDkKhdaKnKNQaK3IOQqF1oqco1Borcg5CoXWipyjUGityDkKhdaKnKNQaK3IOQqF1oqco1Borcg5CoXWipyjUGityDkKhdaKnKNQaK3IOQqF1oqco1Borcg5CoXWipyjsCswBchZTii02zPE6S+ZsptXaeScnC/0wo8QD29H57GbV2HknN5uhD7Fbl6lkXNyrEMrQM5yykK7N8PRpR9+PETQNyDn1MpVjs9nQ2bkEEbOyYVCn255+1PZ+hZGzsmVI/T03FshW9/SyDk5fjRoeOQsiFkOBchZTtMqx3kEHY+ck9uN0Ou/nJzJyJEAOSe1V+iTKSW2vlMg56T2VjnGx+cdb31nQX0+QXdFzskxyzE8chbELIcC5CynLPTyNNBLCPoG5Jza7rMcpxsrxbpdsdJ3PMtE0PHIObkrhX6Zn3yGl6DjkXNyTascIWg+YyCInFNrnOWYjn9h5EiPnAWVhc4nm+czGyvLbBzmlPYQ9A3IObWmdejzCDoeOSe3+4A/QadHzsk1bhRO/YT/iI0VMeScWuNfCldPZ751T9DiyFlO877t8jMbMAQtj5zFsLNGrcg5CoXWipyjUGityDlKPcsxXvApsPTIObXqDyv/+cT8aHrknByF7hM5J8cH/HtFzqmxUagVOUdhd7q9IufUGnenexZBxyPn5MoR+l8E3QtyTo1Vjl6Rc2qscvSJnJMrR+iXefPF9hD0Dcg5NXan2ydyTo55aK3IOQqF1oqco/Cn716Rc2rsxqBP5Jwche4TOSfHKkevyDk1dniuADnLaSo0P2bTD3IWVO+scZFNjs/jJ3vFkXNq1Ubht+L/TnbEzY/ZCCPn5OrvFD79err1zbqdMHJOrv746Gxx8lZ4HkHfgJxTY5ZDAXKW0zQPzda3OHJO7cqvYG3Z+pZEzsk1Fpqtb2HknFzjn75Zt5NGzqnxeWityDlK47e+GTmkkXNqTd/6ZutbGjkn17SjmeOt7yyozyfozsg5taZVDra+xZFzao07mmHdThg5J8eOZnpFzqk17miGkUMYOSfHN1aGR86CmgrNZwz6Qc6CGn+8nq1vaeScWuN+OVi3E0bOybGjmT6Rc3LsaKZX5Jxa40bh1L8Vjli3S4ucBZWrHH47+8Tq6czJBB2PnJNr3oN/fubtkaDjkXNyfMBfK3KOQqG1Iuco/E5hr8g5NX6nsE/knBw/jdwrck7NFfr8POglBB2LnHvARqFW5ByFQmtFzlF8oVffz7d5y3dCgo5Hzum5QodvTCwfzv5Z9gRBxyLnHviNwrPfmLiEoGORcw/8KoffqXzOyJEaOacXNgrdH7DGLZcg6HjknByzHFqRcxQKrRU5R6HQWpFzFAqtFTlHodBakXMUCq0VOUeh0FqRcxQKrRU5R6HQWpFzFAqtFTlHodBakXOUxkKzV8x+kLMc9uA/PHIWxB78h0fOghoLzZ7le0HOgliHVoCc5TDLoRU5R2GEVoCc5TDLMTxyFsQsx/DIWRCzHMMjZ0GsQytAznK6zHJkwe74wZlS9wjkfANGaAXIWQ6zHMMjZ0HMcgyPnAUxyzE8chbEOrQC5CyHz3JoRc5RrqxDr6cnv3FD0NLIWdCVQr/MT35wnaClkbOg64Vm6zs1chbUWOjp+BdGjvTIWVDzRuEyG4c5pT0ELY+cxTDLoRU5R6HQWpFzFAqtFTlHodBakXMUCq0VOUeh0FqRcxQKrRU5R6HQWpFzFAqtFTlHodBakXMUCq0VOUeh0FqRcxQKrRU5R6HQWpFzFAqtFTlHodBakXMUCq3VpZwP9X63tKPQWrXKmdSPUWitKHQUCq0VhY5CobWi0FEotFYUOgqF1opCR6HQWlHoKBRaKwodhUJrRaGjUGitKHQUCq0VhY5CobWi0FEotFbtCs0nlY7IFbrB7XfzDnUfocl5m2yEJuibReTMKNJToQk6wo0Dx71+iLr3Edp0mpJuzflyhU0/BcMW+v2OHOnvuGjOlxeyJtEvyd4+cjRtZCp4EYi+6yTJuWEhdWle1vnuNRV68zX6N6gFgo65hpZXF6Xhqb/xynvP+fKRKFHX0Hahyw/jvCu/U1j8u/v9vOObiHv4FrUImpwFtAizsdA/nhs5II2cBXVfh4Y4cpbTfZYDUIxCwxQKDVMoNEyh0DCFQsMUCg1Tbin0kH8z0kXs6SDnRi3CuqnQLHTDQomv/l4XotBDLZT46u91IQo91EKJr/5eF6LQQy2U+OrvdSEKPdRCia/+Xhei0EMtlPjq73Uh5qFhCoWGKRQaplBomEKhYQqFhikUGqZQaJgSVehF9mF+eKDDQpvnLHt463pL220+6Xz38iwbd16oy2Parj+/dl+oNXKutM85ptDLh7elj6o+0Gmh4rEvWu6DYu8GllnboG+7pS6PqbhTo9ej+ymInOul2uccU+jiNRxeMfWBLgv5Iy0f/m6h1dN/2wa9u3tf2o8c9UKrp7dl24U2fy/339UliPbIudQl56hCz8odsdUHuiy0DS/pbgsVD2PZOuhqocW4w1Na371l+zfqeod0XYJoj5xr7XMeYoRetF7h2t2S+wJOy6R3CxUht316dguNu+xlTvUIfZc5978OvXluvxO3/RtoPXIcrNu1HTnqhbp1sxptFK5D32nO8bMcbiewnbe+i4UWHQaB3S11CHq3UN7+huqF3ORA+ya4oLsG0R45V9rnzDw0TKHQMIVCwxQKDVMoNEyh0DCFQsMUCg1TKDRModAwhULDFAoNUyg0TKHQMIVCwxQKDVMoNExRXGj/lbP159f6mzoHX9lZpPiKyF2ylbPiQi++TC4HneA7qffKVs56C73+6dcizDLohfvWWn1g676aNvrbx9FrOFoc+fN89Wm7fZmHC6w+ProvyedZ9vDX2dadg0uM5ay30PnM/S8EHb5XXB/Y+pEjH5ffNy7+Xf9YBh1OWT29bb7Oi3+3+Q+fWn8j+T4Zy1ltodfTzO2bLQQddhdRH9iGoCflbiRckC9l0OGUMnR3wc3P397bu2avrOWsttB+rz+LWTly+Le//d36hKDD0eKiYZwoxo9wSgjajRzLcf6l9Vfs75G1nLUWukhu63Ys8mu9bpfN6gPbKujy6CL78HG+eXZreOGUcj2vGEZGr0X6Qz8YxczlrLXQXb1cTNONH5CiPmfzhc7f2TyqcupztlJowKPQMIVCwxQKDVMoNEyh0DCFQsMUCg1TKDRModAwhULDFAoNUyg0TPl/D3nVuTeBD2EAAAAASUVORK5CYII=" alt="plot of chunk unnamed-chunk-7"/> </p>
<p>Upon examination of the SNP frequencies for each population, there are many SNPs that are not polymorphic in some populations. The following plots show allele frequencies after removing SNPs that have an allele frequency of 0 in the given population. </p>
<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAtAAAALQCAMAAACOibeuAAAAq1BMVEX9/v0AAAAAADkAAGUAOTkAOY8AZrU5AAA5ADk5AGU5OQA5OTk5OWU5OY85ZmU5ZrU5j485j7U5j9plAABlADllAGVlOY9lZgBlZjllZmVlZrVltdpltf2POQCPOTmPOWWPZgCPjzmPj7WPtY+P29qP2/21ZgC1Zjm1ZmW1jzm1tWW124+1/rW1/v3ajznaj2XatWXa/rXa/tra/v39tWX924/9/rX9/tr9/v2Jv+W/AAAAOXRSTlP//////////////////////////////////////////////////////////////////////////wA6rGVAAAAACXBIWXMAAAsSAAALEgHS3X78AAAeDklEQVR4nO3dC2PiyJWG4agZ40k6ns7C9CSd4J4kG3CyaWsyMbf//8tWpRISV6E6VEnHh/fZbA/GiMvH56JUNuI3W8CQ3wx9B4CYKDRModAwhULDFAoNUyg0TKHQMIVCwxQKDVMoNEyh0DCFQsMUCg1TKDRModAwhULDFAoNUyg0TKHQMIVCw5ReC716fHjbbpfZrPhf8eV6mhUm5y/rL+I3+7T/1ZkLNt89vNwiu7yZLcUjLR7q6rH4z4f56rHINP8wb76dLKOlewZn1dUXz9Oluzd6rS8+8zUo7l/r8y/Va6Hz776fHxR6XDzW/eD3NLGXF7vs8pO1nrqfnzuwefYNKZvsxo3uhb4xo/q6ihMtz1NV6OKObZ5Hr6vH4r6VhW55/qX6LPTm+X+ex8eFdukXI0vxsJbZD4/Fz+uy+P8iAHeR3P0Ab56z7OE/fvwpLrV6LC/llN/eXd3uOnan3M1l2ejb9LtH90yXZxUbfHGXLm/g6Prc6DZxzVhPI48Z6e1qdabQyTLyt1n/u8z+4p6nN7/B+vAaR3+sC+3/83FcF9qdl2fxXkr7LPTqcZYXj+yk0O6/xeNbFnnkdZbFRVbFeF6cct+uBoByFHr473OZT/Vtf3Uul2orf2rrR5/1tHxayrNWj+PiFXZ3A/WZ5fVVN5L7n6R3ZjcanxY6XUYnhS6foGqDw2uclF+WN1YORavHP3x+3St0ecFYA3Wfhc7do5udTjnychCYuBCr/1QjtJtz+cftcvSXap6L6tv+6nI/9DSntrsna+yfXPedrLrx8hqqM5tn2V1j9TP3zuzyKOfQbpTcm3KkzWhZ3mBT6GqDg2v8MN+/g77p+WRvyrF6jBh5j4X2r1Xj453CWfnIi2+4M8v/NOPI6HVxUOjiZFPo6tvVU5SVOxjNqe3+k+XPqp9Q/2T5M6vry/2TtXkeL1pn7Cq1jNCJMzoeoasNjq+xGSRyv9O6evrXbqfQj1tZrE73WOhy/Mj2Xq92exGXRujqZ/vSCF19uxlzttttc2p7NPq4bx4+Wf7Mw9GnuOjv3t0U+ngOXe4F5L1kdFJov8H+Ne6P0P6y5dV9rEboSh5rotdjocv7nNeLPE2h3bxtN4cuf37dnMu//K2n9QhdXWqv0P7bzfywvI7dqe3ek1WdVVxDsRM0q26gPnOy+6lxJ5dZ5L3uPvhVDv/Qyy9Hr+VZKTO6PIc+2N2rn8+tvze7q1tmTaGXfrc1Thb9FXpT7soVgZ0U+nCVo9gncXvFxUXcLnixx+3O+L+9VY5doatvV1fnlzibU9u9J2t31rK+vj+WO6d+UbTen8/G73Wpz0WxW4cuJ267l/B0GZ0W2m3zuqzuwN41Nqsc5XJ5+TJc3LG9EXoZcTla0W8Kl7HX2M/fSPtIcGXR+y6864wo9IE82s7JO/auM1JUaOB2FBqmUGiYQqFhCoWGKRQaplBomEKhYQqFhikUGqZQaJhCoWEKhYYpFBqmUGiYQqFhCoWGKRQaplBomEKhYQqFhikUGqZQaJhCoWEKhYYpFBqmUGiYQqFhCoWGKRQaplBomEKhYQqFhikUGqZQaJhCoWEKhYYpFBqmUGiYQqFhyi2Fzg5Eu0uA3E2FvvgFMBAKDVMoNEyh0DCFQsMUCg1TWgu9LNfjZpe+TaEjuZIzArQVevN1Xvy7+u3rhe9T6Diu5YwAbYVe//RW/Lv5mUKndS1nBGgt9I+M0H24ljMCMIdWgDl0PKxywBRGaAUYoeNhlWN4rHJExCrH8FjliIhVjuGxyhERc2gFmEPHwyoHTGGEVoAROh5WOYbHKkdEIascx++GpdBxXMsZAcJXOSh0bNdyRoDwOTSFju5KzggQvspBoftBoCKM0AowQscTvspBoWO7ljMChP8tB4WO7VrOCMAqx/BY5YiIObQCzKHjYZVDKwIVuTZCz5hDp3clZwS4tsqxmFHo1K7ljABXVzle/kGhE7uWMwJcXYfefP09hU7rWs4I0DqHLleSWLZL7krOCMAqh1YEKkKhtSJQEQqtFYGKUGitCFSEQmtFoCIUWisCFaHQWhGoCIXWikBFKLRWBCpCobUiUBEKrRWBilBorQhUhEJrRaAiFForAhWh0FoRqAiF1opARSi0VgQqQqG1IlARCq0VgYpQaK0IVIRDgSnAocDi8YV2R4Y4/SRTDqcbmyxnBCgLvShHiIe3o+9xON3IhDkjQDNCn+JwurHJckYA5tAKMIeOpyq0ezEcdfvgRwp9A1HOCFBNOT6fDZkROjJZzgjgC326512eyypHZLKcEaAaoafnXgpZ5YhNljMC8KFBw2OVIyJWORRgDh1P25TjPAotJ8sZAZoRev3nk28yQicQnjMC7BX6ZEmJVY4UwnNGgL0px/j4e8d735lXf59Ch5LljACscgyPVY6IWOVQgDl0PFWhl6eBXkKhbyDKGQGav+U43Vkp5nbFpO94lYlCy8lyRoArhX6Zn/wNL4WWk+WMAG1TDh80f8sRkShnBGhd5ZiOf2GETu9azghQFTqfbJ7P7Kwss7FfU9pDoW8gyhkB2ubQ51FoOVnOCND8gT+FTk+WMwK07hROywX/ETuF0YhyRoDW3xSuns68655CR3clZwRoP7ZdfmYHhkLH154zAnCwRq0IVIRCa0WgIhRaKwIVqVc5xgv+2i49Uc4IsPvFyr8/sQ6dnixnBKDQfaLQyfEH/r3iD/xTY6dQKwIV4XC6veJwuqm1Hk73LAotJ8sZAaoR+p8UuheinBGAKUevmHKkxpSjT0w5kqtG6Jd5+8X2UOgbiHJGAA6n2ycOp5sc69BaEagIhdaKQEX41Xev+NV3ahzGoE8cxiA5Ct0nCp0cU45eMeVIjQOeK8ABz+NpKzQfGtQPPjQoovpgjYtscvw9Pho5OlHOCLDbKfxW/N/Jgbj50KDIZDkjQP2ewqdfT/e+mUNHJssZAeo/H50tTl4Kz6PQNxDljACscijACB1P2zo0qxzRiXJGgCufgrVllSMmWc4I0FpoVjkik+WMAK2/+mYOHZsoZwTg76G1IlCR1nd9M0LHJsoZAdre9c0qR2yynBGg7UAzx3vfmVd/n0IHE+WMAG1TDlY5ohPljACtB5phDh2ZLGcE4EAzveJAM6m1HmiGEToyWc4IwDtWhscqR0RtheZvOfrB33JE1Prh9axyxCbKGQFaj8vBHDoyWc4IwIFm+sSBZpLjQDO94kAzqbXuFE7Ll8IRc+i0ruWMANWUo9zPPrF6OnM2hZaT5YwA7Ufwz8+8PFJoOVnOCMAf+GtFoCIUWisCFeFzCnvF5xSmxucU9onPKUyOj0buFR+NnJor9Pl10EsotJQ0ZwRgp1ArAhWh0FoRqEhZ6NX3823e8ZWQQssJc0YAV2j/jonlw9lfy56g0FLSnBGg3Ck8+46JSyi0lDRnBCinHOVB5XNG6NSEOSOA3yl0v8Aad9yCQsvJckYAVjm0IlARCq0VgYpQaK0IVIRCa0WgIhRaKwIVodBaEagIhdaKQEUotFYEKkKhtSJQEQqtFYGKUGitCFSktdAcfbQfHH00Ho7gPzyO4B8RR/AfHkfwj6i10BzBvxccwT8i5tAKMIeOh1UOrQhUhBFaAUboeFjlGB6rHBGxyjE8VjkiYpVjeKxyRMQcWgHm0PGErHJkXvP1wTdj3SMc54wAjNAKMELHwyrH8FjliIhVjuGxyhERqxzDY5UjIubQCjCHjoe/5dCKQEWuzKHX05PPuKHQsV3LGQGuFPplfvKB6xQ6tms5I8D1QrPKkdq1nBGgtdDT8S+M0OldyxkB2ncKl9nYryntodDxteeMAKxyaEWgIhRaKwIVodBaEagIhdaKQEUotFYEKkKhtdrL+cCAd+k9oNBakbMIhdaKnEUotFbkLEKhtSJnEQqtFTmLUGityFmEQmtFziIUWityFqHQWpGzCIXWipxFKLRW5CxCobUiZxEKrRU5i1BorchZhEJrRc4iFForchah0FqRs0i8QvO+irgotEiiEZrYb0ahRSi0VhRahEJrRaFFKLRW7KuIUGityFmEQmtFziIUWityFklVaCZ6tyJnkT5GaIKWIGeR3gvNmNIRhRbppdAXO5y1CL5jxlBoJ7gTiT5JtmWElmykT9Qfvlg5dx04bry7fQruRFuhN1/Fn/UdodC363p1LRFcFPOHL03Olx/iQFPArk9Oy33tcCtXPqew+Lf5/LzjGxR0zKgOQZNzBB3CbC30j+dGDsRGzhGFz6ERHTnHE77KAShGoWEKhYYpFBqmUGiYQqFhCoWGKbcUesjfGekS7ekg51Ydwrqp0Gx0w0aJr/5eN6LQQ22U+OrvdSMKPdRGia/+Xjei0ENtlPjq73UjCj3URomv/l43otBDbZT46u91I9ahYQqFhikUGqZQaJhCoWEKhYYpFBqmUGiYIir0IvswPzwRsNHmOcse3kJvabvNJ8F3L8+ycfBGIY9pu/78Gr5RZ+S80z1nSaGXD2/LMqr6RNBGxWNfdDwGxd4NLLOuQd92SyGPqbhTo9ej+xkROddbdc9ZUujiZ9j/xNQnQjYqv+j48JuNVk//6Rp0c/e+dB856o1WT2/Lrhtt/l4dvyskiO7IuRKSs6jQs+pAbPWJkI22/kc6bKPiYSw7B73baDEOeErru7fs/kJdH5AuJIjuyLnWPechRuhF5wlXc0vuDTgdk242KkLu+vQ0G41DjjKneoS+y5z7n0NvnrsfxG3/BjqPHAdzu64jR71RWDd3o43COfSd5ixf5XAHgQ3e+y42WgQMAs0tBQTdbJR3v6F6I7c40L0JLujQILoj553uObMODVMoNEyh0DCFQsMUCg1TKDRModAwhULDFAoNUyg0TKHQMIVCwxQKDVMoNEyh0DCFQsMUCg1TFBe6fMvZ+vNr/U6dg7fsLFK8ReQu2cpZcaEXXyaXg07wntR7ZStnvYVe//RrEWYV9MK9a60+sXVvTRv97ePo1X9ZfPGn+erTdvsy9xdYfXx0b5LPs+zhr7Ot+w4uMZaz3kLnM/c/H7R/X3F9YluOHPm4er9x8e/6xypof87q6W3zdV78u81/+NT5Hcn3yVjOagu9nmbu2Gw+aH+4iPrE1gc9qQ4j4YJ8qYL251Shuwtufv723l41e2UtZ7WFLo/6s5hVI0f58rd/WB8ftP+yuKgfJ4rxw5/jg3Yjx3Kcf+n8Fvt7ZC1nrYUuktu6A4v8Ws/tsll9YrsLuvpykX34ON88uxmeP6ea5xXDyOi1SH/oB6OYuZy1FjrUy8U03fiBWNTnbL7Q+TtbR1VOfc5WCg2UKDRModAwhULDFAoNUyg0TKHQMIVCwxQKDVMoNEyh0DCFQsMUCg1TKDRModAwhULDFAoNUyg0TKHQMIVCwxQKDVMoNEyh0DCFQsMUCg1TKDRModAwhULDFAoNUyg0TKHQMIVCwxQKDVMoNExJW2j/WTOrx0n5kQWrx92Hhe0ss9kymzVfHG56w8feLd0tzaqrv/yBkIvRa31xdz/dh0h+mJefdJap+EynjjbPH+bb9fTD/1YJH0R5XyGnLPTm2d/3ptCTo8/iuJz1evpwy2fQ1NdVnFhPx5cuVmVd3LHN8+h19VjctzLr8Xv7jOtl5j5McLJLuGOhLYacstC7B3xQ6OrMfDeSlD/dj1nxyJpT7mchy0bfpt89luO6O6vY4Iu79KQMyJ+5evzh0f2Qu6F/4n581tOJv/7632X2l2f3wZJ+g/XhNY7+WGft//NxXGftzsuzm0awHhVV+TZ1bZnsF/geQ05Z6N1gfFBo//O6+n7uTvms3fn+i+rU1g8e62mZannW6nFcvEDusq7PfPhv8UPvsik2znfXtz3Iuhw8qg0Or3FSflneWPnat3r8w+fXvazLC76TgXrZTO4WuxH6LkNOWejd5MnPnY/m0Et3jg8j9yNHc2q7y3rsfyrcd4pzl3XW1ZnNk+RvxsVdX3sxYjRZVxscXOOH+f499E9CPtl7NVw9jjR82G8nxWg73iU9q18b7zDk3kfoZfVzOnpd1FlnZdGbU9v9rP1Z9fPhs/ZnVlnnPuvN83ixm8gdDx7VBsfXWD83xZnl/Vs9/Wu3vzLzz9h76XT5iHcDYjVC32PIvc+h/d5D9ZPbDBn+4s0e48Hg4b55mLU/83DwKC76u8nBTR8OHtvqxqvT+4OHv2x5dR+rwaOSv5NJdF3o4sS4evh3GXL6VQ7/GI5G6NxFsBs8yl2D8ovq1HYv6+qsYnpXvKrOqiGoPrOatI2rp7J+MT03vTvYEyk33k3vXKDV1bnr2GW9PF4A06wutFvCqwp9jyEnXYd2u9En69Dj3Xe+1GuYy3rit9zt8tZZ785alhsU+yzlXnOzD1TtjrurLVehTrN227z6DQ6usdkBL9fLJ3Ud6sFj+Y6Wo+tCuyHaP/y7DPk9/abw2g9yy1oounrvIRsqdP5uduA0e+8hv6dCA1dRaJhCoWEKhYYpFBqmUGiYQqFhCoWGKRQaplBomEKhYQqFhikUGqZQaJhCoWEKhYYpFBqmUGiYQqFhCoWGKRQaplBomEKhYQqFhikUGqZQaJhCoWEKhYYpFBqmUGiYQqFhCoWGKRQaplBomEKhYQqFhikUGqZQaJhCoWEKhYYpFBqmUGiYQqFhCoWGKRQaplBomEKhYUproZeZM+vrvtwtco6nrdCbr/Pi39VvX/u6M3eKnCNqK/T6p7fi383PBJ0WOUfUWugfGTn6QM4RMYdWgJzjYZUDpjBCK0DO8bDKMTxyjohVjuGRc0SscgyPnCNiDq0AOcfDKgdMYYRWgJzjYZVjeOQcUcgqR+b1cK/uCzlHFL7KQdCxkXNE4XNogo6OnOMJX+Ug6H6QswgjtALkHE/4KgdBx0bOEYX/LQdBx0bOEbHKMTxyjog5tALkHA+rHFqRs8i1EXrG3C49co7n2irHYkbQqZFzRFdXOV7+QdCJkXNEV9ehN19/T9BpkXNErXPociWJ5aTkyDkeVjm0ImcRCq0VOYtQaK3IWYRCa0XOIhRaK3IWodBakbMIhdaKnEUotFbkLEKhtSJnEQqtFTmLUGityFmEQmtFziIUWityFqHQWpGzCIXWipxFKLRW5CxCobUiZxEKrRU5i1BorchZhEJrRc4iHApMAXKOxxfaHRni9JNMOcxrbOScXFnoRTlCPLwdfY/DvEZGzuk1I/QpDvMaGzknxxxaAXKOpyq0ezEcdfvgR4K+ATmnVk05Pp8NmZEjMnJOzhf6dM+7PJe978jIOblqhJ6eeylk7zs2ck6ODw0aHjlHxCqHAuQcT9uU4zyCliPn5JoRev3nk28yciRAzkntFfpkSYm97xTIOam9Kcf4+HvHe9+Z1+O9s4ack2OVY3jkHBGrHAqQczxVoZengV5C0Dcg59Sav+U43Vkp5nbFpO94lYmg5cg5uSuFfpmf/A0vQcuRc3JtUw4fNH9jEBE5p9a6yjEd/8LIkR45R1QVOp9sns/srCyzsV9T2kPQNyDn1Nrm0OcRtBw5J9f8gT9Bp0fOybXuFE7LBf8ROyvRkHNqrb8pXD2dedc9QUdHzvG0H9suP7MDQ9DxkXM0HKxRK3IWodBakbMIhdaKnEXqVY7xgr8CS4+cU9v9YuXfn1gfTY+ck6PQfSLn5PgD/16Rc2rsFGpFziIcTrdX5Jxa6+F0zyJoOXJOrhqh/0nQvSDn1Jhy9IqcU2PK0SdyTq4aoV/m7RfbQ9A3IOfUOJxun8g5OdahtSJnEQqtFTmL8KvvXpFzahzGoE/knByF7hM5J8eUo1fknBoHPFeAnONpKzQfZtMPco6oPljjIpscf4+P7I2OnFPb7RR+K/7v5EDcfJhNZOScXP2ewqdfT/e+mdtFRs7J1X8+OlucvBSeR9A3IOfUWOVQgJzjaVuHZu87OnJO7cqnYG3Z+46JnJNrLTR735GRc3Ktv/pmbhcbOafG30NrRc4ire/6ZuSIjZxTa3vXN3vfsZFzcm0Hmjne+868nu6ZSeScWtuUg73v6Mg5tdYDzTC3i4yck+NAM70i59RaDzTDyBEZOSfHO1aGR84RtRWavzHoBzlH1Prh9ex9x0bOqbUel4O5XWTknBwHmukTOSfHgWZ6Rc6pte4UTsuXwhFzu7TIOaJqylHuZ59YPZ05m6DlyDm59iP452deHglajpyT4w/8tSJnEQqtFTmL8DmFvSLn1Picwj6Rc3J8NHKvyDk1V+jz66CXELQUOfeAnUKtyFmEQmtFziJloVffz7d5x1dCgpYj5/Rcof07JpYPZ38te4Kgpci5B+VO4dl3TFxC0FLk3INyylEeVD5n5EiNnNPzO4XuF1jjjlsQtBw5J8cqh1bkLEKhtSJnEQqtFTmLUGityFmEQmtFziIUWqsm5+zAgHfpPaDQWmVnTh1/gRMUWisKLUKhtaLQIhRaKwotQqG1otAirYXmqJj9uJIzhQ7AEfyHdy1nCh2AI/gP71rOFDpAa6E5snwvruVMoQMwh1aAOXQ8rHJoRaFFGKEVYISOh1WO4bHKERGrHMNjlSMiVjmGxypHRMyhFWAOHU/IKgd/Yd6P45wpdABGaAUYoeNhlWN4rHJExCrH8FjliIhVjuGxyhERc2gFmEPHw99yaEWhRa7ModfTk8+4IdHYruVMoQNcKfTL/OQD10k0tms5U+gA1wvNKkdq13Km0AFaCz0d/8IInd61nCl0gPadwmU29mtKe0g0vvacKXSAW1Y5OIhgShRa5KZCnz8bUZCzCIXWipxFKLRW5CxCobUiZxEKrRU5i1BorchZhEJrRc4iFForchah0FqRswiF1oqcRSi0VuQsQqG1ImcRCq0VOYtQaK3IWYRCa0XOIhRaq0s581forSi0Vp1yJvVjFForCi1CobWi0CIUWisKLUKhtaLQIhRaKwotQqG1otAiFForCi1CobWi0CIUWisKLUKhtaLQIhRaKwotQqG16lZo/lLpSLxCt7j9bt6h8BGanLfJRmiCvpkgZ0aRngpN0AI3Dhz3+kfUvY/QptOM6dacL1fY9FMwbKHf78iR/o5HzfnyRtbc8kmySUeOtp1MBT8EUV91kuTcspG6NC8Lvntthd58FX8GdYSgJdfQ8epEWp76G6+895wvfyEiuoauG11+GOdd+ZzC4t/m8/OOb0L28C3qEDQ5R9AhzNZC/3hu5EBs5BxR+Bwa0ZFzPOGrHIBiFBqmUGiYQqFhCoWGKRQaplBomHJLoYf8nZEu0Z4Ocm7VIaybCs1GN2yU+OrvdSMKPdRGia/+Xjei0ENtlPjq73UjCj3URomv/l43otBDbZT46u91Iwo91EaJr/5eN2IdGqZQaJhCoWEKhYYpFBqmUGiYQqFhCoWGKaJCL7IP88MTARttnrPs4S30lrbbfBJ89/IsGwdvFPKYtuvPr+EbdUbOO91zlhR6+fC2LKOqTwRtVDz2RcdjUOzdwDLrGvRttxTymIo7NXo9up8RkXO9VfecJYUufob9T0x9ImSj8ouOD7/ZaPX0n65BN3fvS/eRo95o9fS27LrR5u/V8btCguiOnCshOYsKPasOxFafCNlo63+kwzYqHsayc9C7jRbjgKe0vnvL7i/U9QHpQoLojpxr3XMeYoRedJ5wNbfk3oDTMelmoyLkrk9Ps9E45Chzqkfou8y5/zn05rn7Qdz2b6DzyHEwt+s6ctQbhXVzN9oonEPfac7yVQ53ENjgve9io0XAINDcUkDQzUZ59xuqN3KLA92b4IIODaI7ct7pnjPr0DCFQsMUCg1TKDRModAwhULDFAoNUyg0TKHQMIVCwxQKDVMoNEyh0DCFQsMUCg1TKDRModAwRXGhy7ecrT+/1u/UOXjLziLFW0Tukq2cFRd68WVyOegE70m9V7Zy1lvo9U+/FmFWQS/cu9bqE1v31rTR3z6OXv2XxRd/mq8+bbcvc3+B1cdH9yb5PMse/jrbuu/gEmM56y10PnP/80H79xXXJ7blyJGPq/cbF/+uf6yC9uesnt42X+fFv9v8h0+d35F8n4zlrLbQ62nmjs3mg/aHi6hPbH3Qk+owEi7Ilypof04Vurvg5udv7+1Vs1fWclZb6PKoP4tZNXKUL3/7h/XxQfsvi4v6caIYP/w5Pmg3cizH+ZfOb7G/R9Zy1lroIrmtO7DIr/XcLpvVJ7a7oKsvF9mHj/PNs5vh+XOqeV4xjIxei/SHfjCKmctZa6FDvVxM040fiEV9zuYLnb+zdVTl1OdspdBAiULDFAoNUyg0TKHQMIVCwxQKDVMoNEyh0DCFQsMUCg1TKDRModAw5f8BNDQ5oI70S1AAAAAASUVORK5CYII=" alt="plot of chunk unnamed-chunk-8"/> </p>
<h3>Interpretation of Frequency plots</h3>
<p>The frequency plots show a large number of rare variants within the human genome. When looking at all populations together, there is a large spike in the number of subjects that cary a rare allele. The cause of this is partly due to the presence of variants that are only seen in a subset of populations. Each population has more rare variants than common and some variants are unique to individual populations. The presence of variants that are population private causes the enrichment of low allele frequency when the populations are combined.</p>
<h2>The result of combined low and high coverage</h2>
<p>One pitfall of combining low and high sequencing coverage is that the low coverage areas are more suceptible to errors in DNA amplification and sequencing. In low coverage areas, sequencing errors may result in some reads being called with an alternative allele and there is a higher chance that the resulting SNP will be called heterzygous when it is in reality, homozygous. </p>
<p>Another consequence of low coverage is that there is reduced power to detect rare variants. </p>
<h2>SNPs common to multiple populations</h2>
<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAtAAAALQCAMAAACOibeuAAAAsVBMVEX9/v0AAAAAADkAAGUAOTkAOWUAOY8AZrU5AAA5ADk5AGU5OQA5OTk5OWU5OY85ZmU5ZrU5j485j7U5j9plAABlADllAGVlOQBlOY9lZgBlZjllZmVlZrVltdpltf2POQCPOTmPOWWPZgCPjzmPj2WPj7WPtY+P29qP2/21ZgC1Zjm1ZmW1jzm124+1/rW1/v3ajznaj2XatWXa/rXa/tra/v39tWX924/9/rX9/tr9/v0+LOa4AAAAO3RSTlP/////////////////////////////////////////////////////////////////////////////AKHEOx4AAAAJcEhZcwAACxIAAAsSAdLdfvwAAB3WSURBVHic7d0LY+NG2YZhlJQ4hSVdcLoUFpy2wIedlo2gjWzr//+wT3PQwbIsa0an0ev7orSObNnKqyfj0dga/SYFBPnN3BsADIlAQxQCDVEINEQh0BCFQEMUAg1RCDREIdAQhUBDFAINUQg0RCHQEIVAQxQCDVEINEQh0BCFQEMUAg1RCDRE6RLo/ePDe5om0Sb7xy4qbym76ORHH72eIonU2noz47vt4Tn7MVqf3F/9Kd1/7LpJ+lmz/9xt94/ZU2RP3vzI+7c0Pb5k9x6e7/7v0bx80rsm3Uy5d7qU5Oz3rpX70s4yy3vrEuj4q6+3bSU7PKuS9tLrKbLKHl/u3/aPWXF1jVZZ6auFPg20vv+q44upuN5t+jWuBFq/Shyt80dOFejp9o5LSaorn5T70s6yy3tuaadAH1/+9LIqSpb9eWYvX95S90fR/Zfnrx7Vr6oXxVH0WT16rXe1Wbh//OZR5Ur9ea9VaQ7PtizZvfop3lJ7b1q8SLbKw8+muStu1p4vtZVV//mwKmqklsX537xee115rd/+Wb3+/ZvaDrvN9WfNQ9Gw92LdpuQPjaP7v+htz/bHl2e1q9a1fI1qwr1zuSS2wAf1Mv+M/p5t0DH7v66TWrkocsvOsst76xDo/eMmzn4fUzL1ulkdylup+QNW8Ujtov3jKntzyktWLHz4NfsLVL9AtnKsnyI1f7/699JtgL3XLkyySsTRyrxecdPeaZ8v1WvpDO4f//jprVIj/cK23mrtfIuLJ9BtTbF51a1My6bnfO/tsxaxfGh2h/7d0/LNdK3bnIkCPeHeuVwSW2D9Mmb9w/PG1kk9R7k1l3ZWvry3DoGO1e+7sSWLTQNQ3ipKtjK/nbonso/Wv4VdWNa6shvMKjoTumT2XrswsT1Rne3ipr2zUiPVjpgdFq8r72L7x/wdLNEN6tpuianwJvntnzeJebrytymftfLcpompvL8meRfS/nr2sVljtLIrmAoMsYMC2juXS2ILrF9GPdvX2+Rua+tk/ygq7XvTzsqX93Y90OYNaJUXKtKd+PLWScnMoqKspmRmYb7/za+VvSPtdNcq1kdSq7wNKAOdLbwc6OzOaqDViqa5ffo5P87YmBbT9m43ahW7JWpTjy/reJVVVDUsZptPtjJtbaGzP5TdSaBjux16n+VN9jSBnnTvtLTQG/M4E2hVXRV4XScd6GJrLu4su7x3Qa4H2rQ4+q/X/uWb36A8HDhpA9SdpyUzC0/bgOyhv1/nJevdQps19ZIP9o/eMnuo0kKn9ihlp9rn3/55VfSLT7YyPesw6j5l2WyVb9aVFroIdHZjNVGgZ9g7DSWptdBpvNqt8zrlLXTFhZ01SMmuB1pnItb1shtiDuHtrbRSMrso66Vl774b21QVC9fVHlUS5ZWs96F1WFe1FFdu2jsrgY7zO8zT5jXK17FdPPOA/F0u0a2E7fKW95XPag7pza+pfzSHkOblDs/FQyt96DLQaghvmkBPuXcul8QWuAi0GcMwdcqfsziibtxZ+fLeBbkaaDOWkv1m9jg6KTqIxcBhUbJ8UaKPo7NDD338Xx4r2aNq1dMsBpPM8W/+o703P1hvCnRxXF620Ds77GCSVNQoyd/m8lEOu3nZdr2pnWmCWG5zZStTMzyQD7rqd1Tbf9F3fC6728UoRyXQqomeJNDT7p2LJSlGOWygdc5tnXS5K4PMF3aWXd7bWJ8UXtud3QaDMY6h904yzAjFEGYKdDzIES08Db13biDQwCwINEQh0BCFQEMUAg1RCDREIdAQhUBDFAINUQg0RCHQEIVAQxQCDVEINEQh0BClNdBJVJ7BCCxBW6CP3+uTnH/HV/GxGG2BPnynT438gUBjMVoD/S0tNBaGPjREYZQDotBCQxRGOSAKoxwQhVEOiEIfGqK4jHJExmjbgmbRibm3JnDuLTQVnVp08QeccR/loKJTI9AO3Ec5qOjUCLQD91EOKjo1Au2APnT4CLQD9+9yUNGpEWgHBDp8BNpBax9aX0OufnkCKjo1Au2gtYXePzVc9pyKTo1AO2jvcsQNH3tT0akRaAf0ocNHoB0Q6PARaAcEOnwE2gGBDsCVr+kSaAcEen7XTnUj0A4I9PyunepGoB0Q6PldO9WNQDsg0AGgDz0cAh0+Au2AQAeAFno4BHp+jHIMiEDPj1GOAQ0W6IiT7X0xyjGg4QLd5UFoRB96OAQ6PPU3OQLtgEDPL+tDH57rJwalBNoLgZ5fFujXrT00LBFoLwR6fibQlyf0IdAOCPT8Ds+r/9JCD4RAhyCJVmbsroJAeyHQoSLQXgh0qAi0FwI9v2sT+hBoBwQ6AFcm9CHQDgh0CNon9CHQDgh0qAi0FwIdKgLtZbAJzwn0wAi0l8EuGkSgB0agvQx20SACPTAC7WWwiwYR6IERaC/0oUNFoL0wyhEqAu2FFjpUBNoLoxyhItBeXEY5WmcoINADI9BeGOUIFYH2Qh86VATaC6McoSLQXq600Ks4urt4rlvLYsreG4H20jrK8cNbvGo5G7llMWXvjUB7uTLKkaz5LsdMCLQXWuhQEWgv9KFDRaC9MMoRKgLthUCHikB7IdChItBeCHSoCLQXAh0qAu2FQIeKQHsh0KEi0F4IdKgItBcCHSoC7YVAh4pAe+kT6JZrx1L23gi0l16Bbl58/hM8EGgvBDpUBNoLgQ4VgfZCoENFoL0Q6FARaC8m0OqslPqpVpcQaH+D1xl1OtA7Pe720HAppgYE2tvwdUZd2UJ3RaD9DV5n1NGHDhWB9mIDvTu/kuklBLqHoeuMOtvl+NStyAqB9jd4nVFnAt31yFsh0P4GrzPqbAv9TJdjCoPXGXUcFIaKQHsh0KEi0F7ockyJLsfoyhb68LduaxDofgatM+oqge44pESg+2mo85VLfxBoB5Uux+rsTv9CU/ZGzXW+dvk8Au2gz3UKCfQw6pfPMwi0F5frFBoEemjXLp9HoB3YQCfnXYtehabszRrrTB96QOV3OQY9WKHsjS7UuRGB9tIa6EYE2t/gDQfq2roctNCDa6wzoxwDYpRjfoxyDMgGOl4fX84PCmuFrs/5RaCdNdeZUY7htPWhGeUYGn3o0ZVf8GeUY3wX6lzj8E6IutaDwkYEugfGocfWfiXZh19fWq4kS6CHwSjHgK6McsTrlmt9E+hhMMoxoCvf5Yg3nt/luDwXOuoY5RhQa5fjp7f9x3T/1LeFZh9cQR96OK2BbjxjiEBPg0B7KUY5VrvRRjnYB4Wh64y6/IOV/3wc78tJ7INcc531G+HZWyGB9kKgp3ShzvXDFI1Ae5nigxX2QeFCneOGZQTayxQTzbAPfBBoL1NMp8s+KDCd7timmE6XfZBjOt3R2Rb6JwI9iaHrjDq6HJOiyzE2uhxTossxOttCv27bH1ZBoHsYus6om2I6XfZBjul0R8c4dKgItBcCHSoC7YWPvic1dJ1RN8VUYOyDHHPbjY5AT4lAj44ux6TocoyNg8JQEWgvBDpUBNpLMVnjLlp3W4NA9zB0nVGXHxR+yf7HQeHYBq8z6opzCp9+IdCjG7zOqCu+PrrZ0eUY39B1Rh0HhaEi0F4Yh54U49Bj45PCKfFJ4egI9JQI9OjockyKLsfYOCgMFYH20nrW90AXDWIfFDjre2xtZ30PdeFN9kGOs75H1zbRTI9rf7APmjHRzNjauhxDXXiTfVCgyzG21olm6EMPjC7H6JhoZlJMNDO21olmaKEHxkQzo7ty4c2UPvRcCLSXKxfeTKujHA4XVWcf9EagvbR99M0ox+D46HtsrV9Oog89ML6cNLqr37Y7u04kgfZHoEfX2uXwvyAk+6AZXY6xtX7bLl4N30JHLRy2Wz4C7cV2Ob5ruJRpJlkNH+iOK4l0qc5NCLSXKzP4Jw8/Eujh8MHK6K59wX//SKDnQaC9TH7GCoHuiEB7mfw6hbcdaL4+OrbJr1N404Hm66Ojm+TSyBfH5m4t0JyxMjoV6OYPUC6hhfY1Qp1Rx0FhqAi0FwIdKgLtRQd6//U2jTu+ExJof5fq7P+tRtSpQJszU5KHbh/LEmhfl+rc48wg1OmDwsb5Ny4ZMdDCv6l0qc495j9Bne5y6Enl47BaaIk77kKde5wZhDpzUKg+wFp1XINA+7tQZ/rQwwl2lOPmdxyB9kKgA0ALPRwCPT9GOQZEoOfXY/4T1BHo+THKMSACHQD60MMh0KEi0F4I9PyyPvTh+exbpQTaC4GeXxbo1+3ZFAcE2guBnp8JNN/lGASBnt/hefVfWuiBEOgQJNHKjN1VEGgvBDpUBNoLgQ4VgfZCoENFoL0Q6FARaC8EOlQE2ktroMe4xgqB7ohAe5n8OoUEuiMC7cXlOoUGgZ4GgfbSGugxrlNIoDsi0F7oQ/sZfw4RAu2FUQ4/428rgfZCC+2HQAeKUQ4/BDpQLqMcDmcjE+gBX2FJdZkdoxx+CHSg6EP7IdCBYpTDD4EO1LUWesMnhY0IdKCujXLsNgS6CYEO1NVRjtd/BxDo8Ob2J9CBujoOffz+DwEE+vJKMyHQgWrtQ+sRuxCG7cLbpQQ6UMsY5QhvlxLoQBFoPwQ6UOEGuuMV72dCoAMVbqA7rjQTAh0oAu2HQAeKQPsh0IEi0H4IdKAItB8CHaibDHTU/6N0Ah0oWYHuGNSWl+2KQAdKWKC7PY5Ay0WgCbQoBJpAi0KgCbQoBJpAi0KgCbQotxLolgE9Ai3J0gPd8hmJ18t2RaADtcRAd2xsCfQtWmKgx3zZrgh0oAh014h07dwMhUB7IdCdA93xlYZCoL0Q6AG+0USgw0Ggu+aFQC8CgSbQokgOdMeDuIUFOrx5/oIiOdC9V7r8sFBaaAJdR6AJtCgE+hYCfUP9FALd9jifD9mHMmigfVZaJgI96pFkDwTaC4EetZ/SA4H2QqAJtCgEmkCLQqCXHejLWlYafFsDsowryQoP9Ch17jhC019YQ4JXr4IV+EWDJAR61jr3N/6blYur1ymsXEm2/jfY8n53Y/rtA+rcVYditgb626aWA0OjzgNy70NjcNR5OO6jHEDACDREIdAQhUBDFAINUQg0RCHQEKVPoOf8zCgsg+0O6tyqQ7F6BZqVeqw08tPf6koEeq6VRn76W12JQM+10shPf6srEei5Vhr56W91JQI910ojP/2trkSg51pp5Ke/1ZUYh4YoBBqiEGiIQqAhCoGGKAQaohBoiEKgIYpXoHfR3fb0hsNKx5coenh3faU0jdfOmxdH0cp5JZffKT18enNfqTPqnOteZ59AJw/viS5VccNppex333Wcg6LyAknUtdD9Xsnld8o26v6ttp0Dos7FWt3r7BPo7G/Y/MUUN1xW0j90/PXLlfZP/+ta6HLzPndvOYqV9k/vSdeVjv+y83e5FKI76my51Nkr0Bs7EVtxw2Wl1PxJu62U/RpJ50LnK+1WDru02Lyk+xt1MSGdSyG6o86F7nWeo4Xede5wla+kTsDpWOlypazIXXdPudLKZZa5oFvom6zz9H3o40v3SdyqL9C55Tjp23VtOYqV3LKZtzYB9qFvtM7+oxxqEljno+9spZ1DI1C+kkOhy5Xi7i9UrKQGB7onQRXatRDdUedc9zozDg1RCDREIdAQhUBDFAINUQg0RCHQEIVAQxQCDVEINEQh0BCFQEMUAg1RCDREIdAQhUBDFAINUQIOtD7l7PDprThT5+SUnd0Yp4jcJFl1DjjQu8/ry4Ue4ZzUWyWrzuEG+vDdL1kxbaF36qy14kaqTk27/8eH+zfzY/bDX7f7j2n6ujUP2H94VCfJx1H08OMmVffgEmF1DjfQ8Ub9YwptzisubqS65YhX9nzj7N+Hb22hzZL90/vx+2327zT+5mPnM5Jvk7A6Bxvow3Ok5mYzhTbTRRQ3UlPotZ1GQhXy1RbaLLFFVw88/vBlae+ak5JW52ADrWf92W1sy6Hf/qrT+phCmx+zh5p2Ims/zBJTaNVyJKv4c+dT7G+RtDqHGuiscqmaWOSXom8XbYobaV5o++MuuvuwPb6oHp5ZYvt5WTNy/5ZVf+5fJmDi6hxqoF29Xqymaj8wlODrLD7Q8cLGUQMXfJ2lBBrQCDREIdAQhUBDFAINUQg0RCHQEIVAQxQCDVEINEQh0BCFQEMUAg1RCDREIdAQhUBDFAINUQg0RCHQEIVAQxQCDVEINEQh0BCFQEMUAg1RCDREIdAQhUBDFAINUQg0RCHQEIVAQxQCDVEINEQZN9D7x4f3NE2iTRLl17Arbym7KOp7cbuzp0jOnvJsSe0SkYm+JI7e2Phuqy90pq9plkT9Nw/TGjfQ8Vdfb9sCfXh+6HulGZ+nODyvqj/uH9fHl/u3/ePd1gR6ZS5xbZf33EBMatRAH1/+9LIqAr1/jLKUlLfU/VF0/+X5qyxJdlEcRZ/Vo9fpLkuSWbh//OZRNZjZT1GWsIf3w7O5SJ66Vz/FW2oXH9Rz/VP/ySSRXssu+Xu2Hcfs/+p6kWu1zsN7vhE60OY/H1ZFoNUyuxxLMmqg94+bOEubCbSKR5bS8lZqmtfDs46uXrR/XGVdiDzQxcKHX7OGUuUsWznWT5GaZlbHT7fQZrF+rsQG+uE9LpfE6gk2++wNY2eeqtwIdTHVtd7aP356qwbaLseSjBroWKVxYwMdm+a5vJXmgV6Zi4KpeyL7aB02u7D8S1DPaf9IzCrZXTbQZrF+rjzQa3O3WZJFObnT3Z/IRLbcCN3Wmz+qeF3tctjlWJIxA236Aas8xpE+1ipvpdVAm0VF6E2gzUIb6NikL+s47HQPWD02W9cG2iw+CfTG3G2WHF/Wsepm3L/tbKCLjdBPFqlez3r/9HN+ULgplo9YIQxuzECrTq/q4ZbtcqrDVV6j8aSFVneeBtosPG2hs4f+fm3+W22hzeLLLXQar3ZrvY7tctQuFJm/E3ywLXRl+YgVwuDGDLRuU2OdZpuXLEXlrbQSaLso60NnR2wbHUTTh9YL19WObxLlOa/2oc3ik0BnfWibTnMgqsYwVKfaLCw3Qm+nfZ0kqgTaLh+xQhjciIE2Q15Z7uwohxnVLW+llUDnixI9ypEdGN7/RR/MmfHhYswjWlXG6cwwRf6j/u9JoO0oh12i+j96WMU+/1tlkHmnOx/qdY4v1RZ6V+mUYBnC+6Tw2pt8bRD50uKEKN6kxQU6bh54OFtMoG9TeIEGeiDQEIVAQxQCDVEINEQh0BCFQEMUAg1RCDREIdAQhUBDFAINUQg0RCHQEIVAQxQCDVFaA51EEZNhYVHaAn38Xp0guv8dU1NgMdoCffhOz3jxA4HGYrQG+ltaaCwMfWiIwigHRKGFhiiMckAURjkgCqMcEKVPHzo6Mfy2Ac5cRjnqyT3JMIFGCNxbaAKNgLmPchBoBMx9lINAI2DuoxwEGgGjDw1R3L/LQaARsPYW+uHXl7ML9RFoBOzKKEe8toeGJQKNgF0Z5Yg3jHJgSVq7HD+97T+m+ydaaCxGa6D1Va/rV1Ej0AgYoxwQhUBDFAINUQg0RCHQEIVAQxQCDVEINEQh0BCFQEMUAg1RCDREIdAQhUBDFAINUTjrG6IwcxJEYeYkiMLMSRCFPjREYZQDotBCQxRGOSAKoxwQhVEOiEIfGqJwFSyIQgsNURjlgCiMckAURjkgCn1oiMJ3OSDKlT704ZkZ/LEkVwL9uuUqWFiS64FmlAML0hro59V/aaGxKO0HhUm0MmN3FQQaAWOUA6IQaIhCoCEKgYYoBBqiEGiIQqAhCoGGKAQaohBoiEKgIQqBhigEGqIQaIhCoCEKZ31DFGZOgijMnARRmDkJotCHhigm0KpzUe9aXEKgETAd6J1uiR/e63fSQmNpyhb6HKMcWBxGOSCKDbTqdNRnZWSUA8tjuxyfGg8I6UNjaUygu41wcBUsBM+20M9NXQ5aaCwO3+WAKIxyQJS2LgejHFicsoU+/O3sTvrQWJpKoJuH7na1PBNohKzS5VjV79NLz7oiBBoBa/36aLyihcayXLnGyopAY1FsoJPzgz+z/OFHAo0FKb/L0XxQuH8k0FiQa4E+R6ARsCtdjgYEGgFj5iSIYgMdr48vtNBYPvrQEKX8gj+BhgAcFEIUDgohCoGGKAQaohBoiEKgIUoxyrE6+57oBQQaAcs/WPnPR8ahIQCBhih8sAJRuCQFRGmbTpepwLA4bdPpMhUYFse20D81BpqpwLA0bV0O+tBYnNYZ/BsRaATMttCv26Y7aaGxNG3T6TLKgcVhwnOIwsXrIUrrR9/0obE0LtMYcFk3BK810LTQWJq2LgejHFgcRjkgCqMcEKWYrHEXrc/upA+NpckPCr9k/+MULCxecU7h0y9NoxwPP0Z3ta95EGgErPj66GZ31uXIDgfjh3dzaFgi0AjYlVGOZM0oB5bk+jg0LTQWpPWTQjViZ8buKgg0AsYlKSAKE81AFGYfhSitZ303ItAIGGd9Q5S2iWaaEWgEjC4HRKHLAVFaJ5ppRKARsLaJZpoRaASMcWiIQqAhCh99QxS+nARRCDREYW47iNJ2UMjMSVgc2+WonWaVVpZyTiEWpO2DFWZOwuJwJVmIwgcrEIXrFEKUtq+PMsqBxWm9NDKjHFgaFWg9xsEoByRglAOiuIxycBUsBE8Hev/1No0bTlihhcbSqECb0Yzkof7xN6McWBx9UNg4msEoBxZIdzn05P3xWQvNKAcWxxwUqg8KV+d30ofG0vBdDohybRx6Qx8aS3LtjJXdhkBjQa5e6/v13wQay3H1nMLj938g0FiM1j60HrFj2A4LwigHRCHQEIVAQxQCDVEINEQh0BCFQEMUAg1RCDREIdAQhUBDlOECfar/lgEeRmqhCTTmQaAhCoGGKAQaoow0WSOBxjxGuqwbgcY8rp4k6zeNAYHGPFoD7T8VGIHGPOhDQxRGOSAKLTREYZQDojDKAVEY5YAo9KEhykiXdSPQmActNERhlAOiMMoBURjlgCj0oSEK3+WAKAQaorT2oZ91l+OePjQWo/2iQU/1y3+nBBpBa+9yxJvzZQQaAaMPDVEINEQh0BCFQEMUAg1RCDREIdAQhUBDFAINUQg0RCHQEIVAQxQCDVEINEQh0BBlrEBzWVnMYooWmkBjMgQaohBoiDJvoCP62hjWzIFu+QnwMEmgLzbDBBoDm7yFJtAYE4GGKNMH+vJhIIFGb1NMp0sLjclMcUkKr0AzoAcfLpekOLsKFqwx9xCcuF+SAgiYex8aCJj7KAcQMAINUQg0RCHQEIVAQxQCDVEINETpE+iZP54LyGC7A331CjQr9VgJoyDQc62EURDouVbCKAj0XCthFAR6rpUwCgI910oYBePQEIVAQxQCDVEINEQh0BCFQEMUAg1RCDRE8Qr0Lrrbnt5wWOn4EkUP766vlKbx2nnz4ihaOa/k8julh09v7ithRD6BTh7eEx3J4obTSlnGdh3n+qi8QBJ1DXS/V3L5nbKNun+rbSfm5RPorK00LVNxw2Ul/UPHmJUr7Z/+1zXQ5eZ97t5CFyvtn96Trisd/2XnSXMpBEblFeiNnfCuuOGyUmqaTreVsrgknQOdr7RbOfzpFJuXdO8QFRP/uRQCo5qjhd517tiWr6ROdOqY6HKlLMxd/wzKlVYus/nRQodm+j708aX7ZHnVF+jcQp/0obu20MVKbtnMW3X60KHwH+VQk+06j3JkK+0cGtvylRwCXa4Ud3+hYiU1CNP9L04F2rUQGBPj0BCFQEMUAg1RCDREIdAQhUBDFAINUQg0RCHQEIVAQxQCDVEINEQh0BCFQEMUAg1RCDREIdAQJeBA61P7Dp/eijOiTk6N2nGKCBoEHOjd5/XlQHNOKhqFG+jDd79kobWB3qmzA4sbqToF8P4fH+7fzI/ZD3/d7j+m6evWPGD/4VFNRhBH0cOPm1TdgxsRbqDjjfrHBNqcv13cSHULHa/sed3Zvw/f2kCbJfun9+P32+zfafzNx85nfkOAYAN9eI7UHHgm0GZajuJGagK9ttN1qMC+2kCbJTbc6oHHH77QO7khwQZaz66029gWWnczqtMnmUCbH7OHmvY4a6fNEhNo1UInq/hz56kMsHyhBjpLaKomcPml6ENHm+JGmgfa/riL7j5sjy+qJ22W2P501lzfv2Upn/uXwXRCDbSr14upVe00bob4QMeMV98UKYEGNAINUQg0RCHQEIVAQxQCDVEINEQh0BCFQEMUAg1RCDREIdAQhUBDlP8HtXkw70Jb6gIAAAAASUVORK5CYII=" alt="plot of chunk unnamed-chunk-10"/> </p>
<h3>Interpretation of Frequency Plots</h3>
<p>SNPs that are common to the two populations have a larger proportion of SNPs with higher allele frequencies. The median allele frequency for CEU-YRI snps is: 0.1923 while the median allele frequencies for CEU and YRI are: 0.0118 and 0.0227 respectively. </p>
<p>The SNPs that are unique to individual populations are overwhelmingly SNPs with allele frequencies of <0.1. It is possible that more common SNPs point to regions of shared genetic material between populations while rare SNPs differentiate the differen populations.</p>
<h2>Squared Correlation Coefficient Analysis</h2>
<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAIAAAApSmgoAAAABnRSTlMA/QD+AP2iVEMGAAAACXBIWXMAAAsSAAALEgHS3X78AAAgAElEQVR4nOy9e5hc1XXo+avT1Q+9hbolgxFPoRZIcmyLZ1syRnaMUYQdYgjo5s44c78ZiJzri8OdcE0I0F/LNraD7xA0Hlsxk2/GztxEIuCQGAUFOzbm4caYh4wlgVpuLBAgIakl9OxXVfX8sffatc45Vd0tqVvVXVq/rz/pnLXX3medU1W7dq299tqZfCHPmGfDyprl34WW+zuevXVOUdy5enHzbe1wy/r8mmtG5kqPr6xZzvrCSDV3MvAPB4CW+zue+eKcQdWPq/Hkkx9DDY5Njvk242+8Mf6Uxrh5Rpqo0gaMITofWFKz/MFKW3GMPL4y9PLAwmb73I0/xuUbzxhXZCttgDEitNy/7ZnRGFxdsyafXzOmG6xK7CkZI0v1jeg7Vy+uqYn835IHOpPFDywJpUqhc/Ximubb2gF4cHlUU7N4dWdoauUGHl+ZrKIkKx8fziWADStFsmFlUFiyOmljEqXsDFNyPxJsv21uTU20ckO66jCuWN5gubS/aGhtg37IcavkoZWRxBtM3V3qYQ75gpZ4WmXuZfCHWerulqzuLCk8HsOO5Y1X4ikNarnS71y9JOik3w6xduJvA6kYbmSQd0WyvaRC+j0w9OMa6m1gnBj5Qn7s/62/BYCW+zti8o77WwC4Zb2XbPOCGKpWR5nyYlOxWiXVAVpakgU3rx/yEvl8Yf3NZV6FUD31V7JKy/3b1GPRzQyr+vANTjz5VGst96/31eXSxZfAGxlquTbjDZZ+wsWnMdQLmvwrqR/eHoM+zJJ311FaOLRhiXfsMb7x0m/4wS1PvxOE4r2rP2ksvMrKQt/mcbwrVGvFxjqG9ToO9TawvxP+G1cdfTn8u7ns2yvRz8obKPk90ZHos9T7z1dR71ffZpD4Rga/RPHjKgrqMzfYjac6o6AfPmXbSj66Ia84xDMp09Gry20LT03bk5bEv5zi3WUZ24Z8QRN/Sf3E0xvqYZa6u9LCoQ0r+dyG/cZLdvRDvg3CByQ1pEh/95eyv8wb9ZjeFeU7+qEe1xBvA/sbgb8q6uiTb5d8PnyE0u8YPUJKDrtSHX3681D2EzL4JUSklQfrqdMfoZT+MDv64Vyx1DMp3dHHbrb4jMMzvPmWmxOS0g3qgV76AR7TC1pSf/3N6tUc8mGWvLtSwmEYVuY36DDfeEP2pGW+PlUL6Tb1X2L8XvKOjvVdUbajH/JxDf42sL+R+Ks+Hz0PLi96+rz3s33LNlcWHOvHFOTwgVQsS1oSGMYlWi6aWzxpnt/iTOxI6XV2bAISsTSD6JdnsCse+zOJtcac5Te2AO0Pre+kc/1D7XDzdWuuuzkmoeXG5SWemK8KwHeXl/Q+M/gLOjjL1uRDpOywH2b87gYTHoNhx/fGO0bLaZlfwspSzPniXTcDtD/0WCfQ+Zh7ke6/fdlIGFyGso9reG8D40QYVx19OR+9o2NL++DVO1cvkXftzetT1UeEk3CJkWUkDJ5zrevpt2xj25Z2193Mnd+iJLTceG3Jb8Y5tz6bcge339bsPuRDvqCV4lgNG4vvimuuuwXKfRmPrMFDP65B3wbGSDCuOvrB8WOcknM4a64hvOFuXl/Ir1mmKg4yPD9WhneJ9le3paq0zG9OtTaneSEAmzo6h6VfnrJXHJFn4kdkDz76QMcmX9eJNnU88OiDlBvPS+0vPiMvk7gJ2h9a3zmMF7QMxcfVuXpJiC0ZuYd5zIad4EMeQcsV19x+v/96Dv18+DIe2U/K8B5X2beBMRJUUUcvPwAfXB7CyiTSLhZn9uCjLnLr8ftuG60R41CX+O5yiR7bsNINnUp3hTLsuu2PZWwzuH55hrjiCT4T6elvu60dbv79a/C9U/tttw3az4t/QILtfOfjvRTDfUETVtB+232uuPOxh9ppv22ui9UbuYd5rIZ5jvshj5zlGnnRvvK5h9qBW+5KrcMYpsFz5yce+wN/HNMf8nEN/jYwRoSKzxIM5++EwitT03FJkhEsQZieaB1SMvglRiq8Us2wDTvqpswVh3gmQ88KJp98cn4vYdtwwiuL7Q/xgqb+hgjjG/xhlry7Mrc8lGFDhKUmDEtNeA4vvLJslE5+qMnYfOq68XfgMb0rysdKDB5eWSq8LVG74j1P1fxV0YgemHPrM4k3Tcv9HcWfh9esScRHrvczUv4X4rLbi8XDnO5LMtQlRN6htbbFfyAnGkx8kGJ3NFzKX3FYBg+NuBeKo7DiDNsg47ISztmb1+eL6XqGeEFLtPdMopO6ZX2+mI9lZB7msRt24m+8kbNcs+y6EMh4e+wdeGzvimvWxGy7eX2q4x7icQ31NjBOmEx+PCQ1qxY2rIyWPwgjnnpsLF3RMIwxR3WN6A3DMIwU1tEbhmFUOdbRG4ZhVDnmozcMw6hybERvGIZR5VhHbxiGUeVYR28YhlHlWEdvGIZR5VhHbxiGUeVYR28YhlHlWEdvGIZR5VhHbxiGUeVYR28YhlHlWEdvGIZR5VhHbxiGUeVYR28YhkH0ZBQ9WbX9YdXemGEYhuGwjt4wjFOd6MmIp+CpStsxalhHbxjGKU20qvq7wWylDTAMw6gMp0IX7zhV7tMwDOOUxTp6wzBORU6d4Ty2laBhGKcgX8/EevmceLHvGChUxJ7Rxnz0hmGcQrgu/s62uPRKH3JTuKcCJp0ETqEfL4ZhGKcm5roxDONUIToavTwJ4MNHYvL9kzjtCEBhorluDMMwxifR0QjgmzzeCvDhb8ZKTxPXDea6MQzDMMYj1tEbhlHlnFKRlCUxH71hGNWMC7NZBsDjEm9zb2tMp+rDK62jNwyjOolFUv45AN/0XXz58Mrq7OhP9V80hmEYVY919IZhVCHml9eY68YwjKrCR1JCLGTexVP+OadmHL196RmGYVQ51tEbhlE9mMemJLYy1jCMKiE6GvHn4qUhvghWom5OzZWx1tEbhjHuifZFwJZGgPlf8MJ+998XALZM8kVaWOTplKS6sJ85hmEYVY6N6A3DGN9E+yK+BTDf+We+5eW1rcXTUKSFRa4UibluDMMwxhrRExEFuALgvUkA06/wRT5yPF6khUV+kZJUF9bRG4YxLomeiACegy/48fj0JXgJADWfKp6GIi0sEiZjrx5NiyuH+egNwzCqHOvoDcMYf1i8/DFhKRAMwxhPuJyUiAu+BvqhFlC5iB06I/EpnqbYfPSGYRhjgkF+ppxg/mTr6A3DGDd8PRMVU8m7OBk1GeuWv4ZVr7HU81I0VD760bK8spjrxjCM8cE7mWgjfEhODwIwFd6EswGY/QLAW5d4hY2A6IciLQwchKkAvL+irhsb0RuGceryr5kIeAXubCu610v46C+Bcj56KRrCRz9adzA8MqPVsM1cG4ZhVDk2ojcMY0zz9Uz0OQCubePeVqrZRz9qI3rr6A3DGKNE34iA6A663Xk9EVAvxRO8kAleWJgEEIlCJLV0kRYWmZCSVARz3RiGYRjHh43oDcMYi0TfiOgFKMjYnV4K4ISAH+f3QrcXRkdEgq8YTkORFhbpTkkqgrluDMM4dYi+EdHgjwdCHsp6BlA+lpwXckCEWZFIxeJpVrVAylGTN9eNYRjGySL6RuTH8j3QC71koAZq8MdOSC9kIQu9ME0kOcgVFTIofSmKCcNfjRyMXzpXL4lWbihRsGFlVGMdvWEYxtggU/5vcB5fWTP3tvYSBZ2rFy9nfd5cN4ZhjBWiNeKxmQYTvBc+A4dc8SwywCzRLnghR72wZxJAgyhkpJYu0sIiAylJRTgu182GlTXLv9ty//r7H1q+JVW4bUv7zdctMx+9YRhjg2hNxG7xlTu3ey/AAExxGrsZAHZLhfO9kIle2HBEJPiK4TQUaWGRC1OS8cM1a/L5NdC5+qF0WWfHppb5t5uP3jAMY6xw3K6bcnRscf4cG9EbhlFhXIr5zOcZkGF4BiLxzdzZxmOtAN9vLZXr5lP0yz7gPum8KJTIR99aPtfNUoA7Bkb85o6FEY+6aZ7fAlhHbxhGZflKJrrnpwAMQEblKhMf/apWbgDg2p+yqhWvDPwCgMvhAEwD2L8U4E5RWNUKeP1QpIVFnoIrAQpXjezNVZo5zQvbH91mHb1hGJUiejLCDaXdxtyzKProe4s++uxSn0OYp5QykuvmqWKum9MuEwm+YjgNRVpYJGwOftUI3dvxMfJx9HPntzz4lQduNx+9YRjG2GCkfPTFmPo5tz67fuFtzTaiNwyjAkSrIuctMYqcyIh+zq3PFEodc82agsXRG4Zx0om+EVEPfYDKauBSSErmA50CIeQtiKVA6PNCciKMRCIVi6eRaoFUwoO+sZECYdSwjt4wjJNH9EQEcDng3eL5T8npAJxTYjI2HxZMXU4+1EVNxrZ74f4JAKeJQl5q6SItLCKTsRXGkpoZhmFUOZbUzDCM8c4gm18bo4qN6A3DGHVcJCVXir/F5Ym8ElwfZOGVDnPdGIZhVDnW0RuGMU6Jnoh4WU4WAFAH+NF9fqmMzQtwvmwnko1Nxh50da8gjygDz3khz3vhe5MApotCXmrpIi0s8ouUpLqwjt4wjFEkeiLiORXTknLd1CBddhnXTc2nxHXzHDWfEmWka36u6LqZvkQk+IrhNBRpYZHgurn6+O90BLARvWEY4ws/9VrVI+URxqJuDMMwjOPDRvSGYYw8Fkl5PIzaiD6TL+RHq23DME49XHL5QAGiMC8qPgTXoTnh3a9w7+/4ooLKR18jCne+wv/9OwC74U5RRrLVu8ZrALgIgFflWne+Anj9UKSFgdDCHQMFKkf07bLfjoU/PSHDbERvGIYxNrDJWMMwxj5fz0RuC6ciVzB41M29S/2uT+Wibu5t5XMAzG7jXtkiyrdMLOqGPwf4g2/6creNlNeXopgwIFE3hXuO46bHAeZHMwxjZDC//Iky4nvGCjaiNwzjRLEufmSw8ErDMAzj+LARvWEYJ4QN50cMm4w1DGMMEj0Z2xEwT2oHj1+oBAPIDKrLdfMUQG6pVBmAC0tsPJILuW6uJEdqXvdKeNoL908COE0UclJLF2lhkWrfeMTi6A3DOB6ivgjYHN+B7w04J652EbwKC3r96bZ6kE0AnbCjnn5AIu5dPHsejsJEAGqheRvA5rnU4pVRofFnwtuuwX8B2PwZr1ALiH4o0sLAgl5/IxdVNo7+wfJx9DdbHL1hGEYVYCN6wzDGDgcz0fRHAfhoTP7lRu7uiqs+D5fB03LaAsgI82mAtutodVXcGP4AANOKrpu2Rj4JwEe6aGsUZdcycBl8C74AsKURYL4otDUCXj8UaWGRp/2NFGZUdET/t+VH9P/rCRlmHb1hGMfAVzIR8EfwD0DKJ/AX2/ja3JjEZRfIyanzzLiRqxPefR9fu90LB1TFkALhL+7jR7cDvAB/IcooZ36/eGk+BsDPgjH3AV4/FGlhICc38qXKum5GraM3141hGMbYwFw3hmFUnI2ZyMXLLHyYvTcANPXEFH7QwGfjkq4GGnvY2+BPm9w0aB7wwjfg4h6AIw1M6pHBZ47+Bmp7AF5s4OKfAtDCiw1e2bUMNPawqYGFTvjfAfjfvcKLDSCNh6KYUNjb4G+kUF/REf3/U35E/5/MdWMYxujzd5loO/xXAH4ELv/jm3GdZ2BJXPJ+eAfOllMXHuM6HSf8H+LnH4CMSnUTXDdPw+kAzC360n3LwDtwKfwSgLcAmC0Kbl7go/EiLQycLTdyVWVdN6PW0ZvrxjCMIfi7TAT8cRsvtfIGAJ+VEf0F8aHxew18PDWiv1CN6C9IjeivHWpEP02N6KelRvQX9rApXDQ+op9WakQ/rcyI/gI/oh/mIxkdbMGUYRhGlWM+esMwKsLGTPRDAD4Ni9pY1QoWdTM6RN8v77r5nPnoDcMYBTZmImBRG99rBXHd2GTs6DF6Hb25bgzDMMYGlqbYMIyTieWkrACjtvGIvZaGYSSxXr7KMNeNYRhFrIuvJBZeaRiGUeWYj94wjNHGhvPVio3oDcMA2JeJaI1J1gNwdiuLAdjbynq422Un/jRNbUByY6bNS/lse0zS+Au4nKZSO0w54bdbubgdYJLLgSDh8bWSpnj9UiYsBZj/U9Yv9cq+ZeByFh6AaQD7/xLgtJ+K/UsBrx+KtDDQFHaYuqrM0zk5mOvGMAyjyrGVsYZhjBJfzUTA8pR8URvAS63YytiTQ7S2/IKpFbYy1jCM42VjJnJZJL0rRvFSK9jK2JNItK58R3+TrYw1DOPY8VOvrUPpGScN89EbhmFUOdbRG4YxgmyRGJvH4J5HATalhvaPAVDXymfxCo/BPW0A77XSVGpz8M0NfPZITNL4U7iMpg1yrjYHd8L/6zouPgIwaU9sc/Dago+6eayB3qUAH+nisQav7FsGLmPhN2Rz8L8EtTn4Y85LcyRWpIWBpg1yI5XNRz9qmI/eME45tmSiheKRf6kV89GPFR/9w+V99DeYj94wjOFhfvkxja2MNQzDMI4PG9EbxqmCZTgY69hkrGEYx020JgKY5U/71fF+eL87mpWstV/k3XKwH2a44zNkDVS8b6pPSfJAprhgypdmQBZM7ZfTAir3ekZOnUKjFwZlZMEUGfqD8H0xk/brK75PtZAyMpeSVAZz3RiGYRjHh43oDaPKiX4eMTMmOQ+C5Bw4vReUJHAOXj5LDs6BJnf8Y06/BGBzPB4xk5JcBJvrWdDrT7fVA/RlAS88WO+rFCCq9zkS8jUchYn1AGcDE/y1zlbtX4QXnifCBbNiJp1N8TQUaWFgQa+XXDRAJTHXjWEYx0H0SMQe2B0TvnUD0zf744M3sLceoGlzsu7BGwDYzCE5OHgDLi5x+gJfa0E8TnFrQ1LS1cACFV45NxVeeWgY4ZW8BbDgrBLhlQt0eOWvARZc7RVikZRSVC68coEPr0w+hJOKJTUzDOOY+NtMBNRDLrizhQ/Cr+T4HVgCwOupFt4B4P3xg2UAvAZnArAvXmULzI9LpsN7wbMP7wHOHS/CZ6RKA/RIerJ+iERtC5wu+rr96dLgOfCGany6MgbRD0VaGJghN3J9ZZOa/XP5OPrftzh6wzCMKsBG9IZhDJ+/zUQ3u5WuC0i7bja1YitjNWNlZewPy4/oP20jesMwBFv7aqSxjt4wDGNsYFE3hmEMSfRIxAJwC5cWALAZfhd2xtRmh1KYCk0u9nEBCabi5VPkYGoIr9xM06cBNjfEqjSkJBfB5oZUeGUtSHjlvnpfpQBRg4RXZjkKExsApniL2dzAFNW+D69sYLYIF3wgZtIUiqehSAsDC3q9xMIrDcMY00SPREiIZO8N+OMF8IqFV57q4ZXW0RvGuCd6JALYXByV24ie8TiiHzUsBYJhGMbYIFP+bzA2rIxqaqKamsWrO9OFnauXRDXW0RvG+MZyUlYPx9PRd65evJz1+Xwh33HjQ80rN8RLN6yc+9CN2/IWR28Y45goFwF7agFm9vsD4OfwEbxwZ23SRTvzCfaII/s7cM/XAfbckWz8OwB8Hg4DMNkpPwHAUolMj28l+OXF3P1svJUX4BJ4Wk4vB/DR+E8DrLqDe3QV1yHVwAS/leCqxbg5gvnPsmqxUn4BgEvgCEwC6FoM0CgKqxYDXj8UaWGRp/2NFD5S0Tj6DeXj6K8pZ9iGldGj1xXWXAN0rl4yd8td7tjRuXrJ5/jes7eaj94wxiVdmQiY+c8AM2XedaZMqO5uZaYsmDpj0AVT17ay9w5A9BXXtnr5v8qCqWtb2XQ1uAVTd0BqwdQC4OK4nYtp7GHvYn8aWzC1GGC5VCk3Gbsc5suCqaBM6LvVZGzjV4Higqnl7r+LY0VaGNi7OHkjleE4JmM7Oza1zL/dHc9pXshDHZ1cM0dK5zQvbP/K+k7r6A3DMMYIx9HRd2xpT2bu0VyzpsDKqMY6esMYf0Tropk3ArAVgHlyPM8r1GvhEjgtVv0MipqToSmnGlFMxssnycFk2QKErb5WRzx8pS4lae6jo4Fm2Xlkey1I1I0TdmV9lQHIlIq6mewtpqOByar95j4vPEOEzRfFTHL2J4q0sGhkzksuGHdRN83zWwYp7Vy9ZO5DN27LW0dvGOOJaF0EsJU9EcBM1zuHkPkQR498ByyALUnXzU5o3OqPD8PeLEDTVhI41zxbOSIHh6EPgMZ5vlZz3OOxqSEp6WqgWcXRn5uKoz88nDj6nQDN55aIo2/WcfSvxkyKhcxLUbk4+uZxGkc/p3lh+6PbYA7Q2bGJhdfNUaUdW9pbbvzeHAuvNAzDGCMcT9TN3PktD37lgU6g87GH2m+57hpd2Dy/pf2h9Z2WvdIwxg8bM5FzLWRg0S6Al08HCMLgeNgGc0XYK0EugbmwTY7Xw93bAV46N3m59QAshwMATIP18GcA/AY+vB1IOoW+Mo27DsRbcSEx++XUuU5cdNB+gFXnco+r4jLQh5vp92qrpvmM+R8/wKppouxaBiZxaBpTDgBsnAbwIVFYNQ3w+qFIC4vs9zdSmFrRqJsfl4+6+d1BDNuwMlr+INByf8ezt84hHn7z+Mqa5Q+a68Ywxge/yUSL2uEDAPyad1sAPuw60NcAuFAO4MkWbmoX4a+lYxW2tLCo3R+/2MLec4GiJPBii5c/1AKwtJ0XW9iBv66r1dQdq9JM8lvl7Wmc2e2VgSbn+nGum3PBbSpSB3BwAlO7izuP9E+jttsrfPwQAHVFZdcycGY3O2B+HcCHXgSKu4q47UqIF2lhYO+5yRupDMeZAuGaNYX8Gi2Yc+sz4Xth2Zp8YY119IYx1omei4DDLfCE7KW0lcjNwf07IK75t4s++ilO2QlT4ZWFUAoXh7m8J0jgi56gWw4ubvG7PvHv+OD2eN+UTUkaEpJMstbFclqj5UrtYiRkvqWo7FsGMhSC0K0WklwOF+umpCgm1IxanpljYNRsMB+9YRhGlWMjesMY00T3RnwK3AaqWfFsZGVkLaf+WD7QA1pYm/ygF0uhIFuzpjuDIB+Qg+BCp7Z0rYGUpKCtBR8+mSm2H6rk3UGNV8xDbVYU8soSZXlSWBczKZiti2LCuJGVx7YSNIxTjR9nIiADmwA4ApfDRACOws8B+KQUAZPkAGiHFhG+rfbLdoQECcAnn6FtCSgJqgj40RJ/rY8o5U/Cj4Bit+y5FVbHJQWI1O7kzofgOjQnbH2Cr1xd1ERpuq78rif4+6sBXoe7RBkpdY07M1rk3pGK4PVDkRYGQgt3VnZz8CfLT8ZeZVsJGoZhVAE2ojeMU4pvZ6IzAfXZ/xV80K2EgnqJrzkfELdEn4olccpOOEVqBV6DC+X4CD6Q5kKShB8KrwDwO0r5fHg9bp7jEzI9HHDj9BD36YasA+rfc2WCWY/9ayAn49Bz8LudvAvniDJqRD8gZlwod4dURPRDkRYGQgu3V3ZE/7PyI/qP2YjeMKqLvZloovJmh9WaNRK23iulMwCf5JEZchCUnXAAcsSoVS6X0+Btpa9x18pLj+yuvkua3Q5Af7zKlNS1spBT4Z2udEDMAOpVlYJ0uDklrJcvsLfjyrKElgb5JuhXl0CeWy5eVB/XQSxJ3EhlsB2mDONU4NFMBOyHGfCZRwDoh98DmDyVj6wtjkgvXAEw18Ud/gaAC+QAulbw6bUiPJLcYeqFFVyy1h//bAXXAzBnLQl+tgLgY2vZvgLg02v52Qo+BsAlL3DeJRA2nBJ+AF+IS7qgEfbKqftmcgNXJ3xLquyDGdLh1sJhWVy1ET72KsCBi9io2u8CoFGtDpv2CMCB673CRkD0Q5EWBvambqQyWHilYRiGcXzYiN4wxgrRc9HkFoB+mBUrgODAccua5jHX+Wz/DXAZreDf5ECP4P4Nn0BAcZ46PgPm5JIKoSh90CjNzvkQQEe8C5kJ78Ybae6jo05lr8yCZEZzwlzWVxmAd8WDlIejMklwOn6o+y6crtr32SvrOE+E00SNUDFVpIVFI3P+RqaUfBAnDZuMNYzq5jeZaKd0bfvhPz7DC0sAInCr93fBcvFmHIbHAPgjAA4CMFUOgB/Cp0X4RjIhDY/BtXK8qI0vt0LYpkOxqA3gpVZ+CMCnYVEbq1r9df8BSA0V/2IbX5sbk+RlZtWhwui98O77+NrtXqiTBNfIxOxf3MePbgd4Af5ClFHTtpIUx7uVfhaMuQ/w+qFICwNh4vdLlZ2Mfbb8ZOziEzLMOnrDqDAbMxGwaC3AD1cA5GESfNJtipSDqwFeamFRWzHL/B63/dPDQIk0xQ+2crPsMJVOgaB3mHqplbMBaErtMPVSK8CiNr4nO0y91OrnRRc+zN4bILXD1A8a+GwqTXGjSlPclEpT/MYw0hRfLDtMpdMU6x2m/O6GssNULCOxFJVLU9zk0xRXZ0dvrhvDMIyxgbluDKMq+WYmct6MGqhV0SnTxJGNCpqcoELmu0UN0axTVXok51cdHFIBmo68CqasgaNASgflLu+Wq9fIJSaASyiZGGcm3C/IAtcwUi3Ei4DJyuNUkqly0Z64eyqSBsONh2XDoSKiH4q0MGEkcHdlXTft5Uf0LTaiN4xxyDczEfDf2sRDYq6bMeG6ST6E6sA6esMwjLGBuW4Mo5rYmIncIHWpBLd0tGJRN6d61M0vyrtuLreoG8MYV2zMRIva+KtWMNfNmHPdVLSjf758R3+Z+egNY5wQrYuAyStgniQgm+cH6WyFq/2M6H6YiSTi6vfbBw4A82CrrzWzILUo9v7+wM2sBuGS5JD+DIqak6EppxpRTBYLJylTG8VaV6ujIValLiVp7qOjQS2YqgXoqwVZMNWV9VUGINMgM8BZjsLEBjHjDH+tyap9v2CqgTNE2HxRzCRnf6JIC4tG5rzkgsRUcrVgHb1hGMbYwJKaGcZ4J1oXuQH4YWCrhEJu5bA7mAdv0NsCUO9iBF2+3ZxPZJZpUWP2reyJAGa605TrphcZ7C+ALUnXzU5o3OqPD8PeLEDTVhIcFguPKFOd2Y3zfK3muA9kU0NS0tVAs3LdnJty3Rwehr8TVh0AACAASURBVOvGJWVrPreE66ZZu25ejZkU89JIUTnXTfNYiLqxjt4wxi/RqghK+EYM4+RgHb1hGMbYwEb0hjFOidZFicnSydhkLNhkbBrLR28Y4xHvtDGMimIjesMYFVwkZXEMro5tMtYmY0tjrhvDMIwqxzp6wxhHhEhKUL5v89Gbj75CWEdvGCOM7+V1/+4w1425bgbHRvSGYRhVjmWvNIyxz75MBDTl2JsVf4iMl4Mc+A58Hr4DqIO7c3RmeQSAAvTC5wEVGNe0nb3n+mb3Zmn6OgC3AeytB2jq9QfAr+CDeOGu+uSArulR9l7nj78Dd7cB7G1N3k6w8D0ApjvlRwFYBs6AK2NVvryUu38ab+UXcDk8JadXALJ/ylMAq1q5x1UZgIxKSjnB73iyaik3ADD/p6xaKsquZeByOOB3YNm/FOA0UVi1FPD6oUgLizzlb6RwVUWTmr1SPqnZ71j2SsMYA+zLRD4B5DwGd9281Iplr8SyV6aIfl2+o/+AZa80jIriIimbbohPitpkLGCTsceE+egNwzCqHPPRG8bY5DeZyEWnrJdNmibD4TAcDrErIgcWtfFSK7bDlO0wlSDaXN51s8B89IZRIV7IRJd93zlKOHwTk9cB0Av1LsIR8KVFOTx/E5et44GbAL64judvArjs+3AW/2MpwBGYDb/3AwD6fZ/9q4l8cJ00W8/BmwCmuu7W9Q+F4tTt393E/7xOhGf68MTAb2/ivHX++PmbmA8gxiu8Yev4+5sA/mgdz9/ETADO+wcO/weAQ/EqP4GPxyXT4ABMkdM+eRKIcDN+kuAoTCy66P03BPAmXN4BsLOZN0UZP//KAegB54k5458Adv6BV3gTEP1QpIWBKXIj76tsR7+lfEc/33z0hnHSif4+Aqb/R+iBOQDdMHk2AJthTnFe1JUW5TAdmC198mymu4MeOJMmAHphKnAmADkOTAQX+jJbmp3jAlKY6rpMN4jeJgeuw58twjOlyxROC6Uw3ce2iPGK6WJhVpnqfxz0+lqd8SpTU5L50Kn8/68D4GZk58lVXJVa6IeJAByFbpgAuOmK0/y1Zqr254sBC2ETAGfMjpk0U1kYirQwME8k70s+g5OLJTUzDMMwjg8b0RvGMRM9GBHCNna64Tczc7BDSYK3ZGpSE2CHdzH7Kq60yw/pJrmx3VsA9DPtMMDkybBDmp3K+9xP+bcB8Vy8TbAqH1p+G95Kum6KpYjlcUmMHTIe3CF1gbd9rf3xLuQ9mBWvPaOPWXVyCTiSBXHgOOG+rK/iwujdb4/J4slxcroAZs1kQLU/ow9gVh35INwBMOtSr+A8/omimFA9gf1joS+0qBvDGAv4nJQLfC98mGLI496sZBEYKo7eUiBYCoTSmOvGMAzDOD6sozeM4WK7iBijS6b834lhrhvDGBbRIxELxB/S4wMlpwKdsACgqdcfsJWiJiIMcldrgayMXSB++074lI8z2eecyAsB6PdLZDOoZhcw0/lAgtPGHcu16rXwd5M++tnKqqnOcm0nxSInn6JMdXFBbKbp0wCb40tMG1KSi2BzAwsk2HRbPcjKWCfcV++rFCAqtTJ2ireYzQ1MUe1fhBfOFuGCD8RMmkLxNBRpYWBBr5dcZCtjDePUJLo3AomSdJOKff7HcMFJnDCsC8opIfGDnKolwkIo7ffHA24iMbRWm7pWTmYtQ9i5Ow7NJoRa7jzkIikgTcV1fFHcwkKYjM36WrWpKglJ1EdtXTG+s1ZuELnugBKGydgIarVy3tcdUO1HfQC1deSDMBczaUBdMRTFhIGalKS6sI7eMAxjbGAjesOoCNG6CB3Z4o7Psqgbi7oZBayjN4yTjJ96TeV0NIxxh3X0hmEYYwPLXmkYJ5ONGR9JOVlikA/KqcNtunRZGx2tWPbKgGWvPBGi35ZPanaeJTUzjBEl2hBNXeaPG+KLTZyreReyF8dlTFS5I6eqKiF+L1J1I6k1VTI4TgUuA5gIE2AX4PzgF/rWIuncXfsNInR94kRlmz5GjJkI01XmBUdWSy7zQTUJHcSwqdI7O1Oz0mwhmKo5NyUJgTpyF0kuLVElodAQrpVSdiE3TjgtYdKlxdNi0aVxHWEIG04ONqI3jJPAl2Ug/wWRNPb5jqprgpzCvjrf9dwHrX101QE0fp2uO7xOV53XBF8a5MC34AvwbQD+FL4FQGsfr9fxjwAUoA/+FIBIPv4z3qHr/b7ZrjoaHwDgFmVbtz8AXoJFeOHuCYnklTQ+Tpd8k30LWr8OeOM135JHsQeAmU75cQCuggfkQNF2Ba3PxVtxe8Y+KacfAWTP2CcB2u6QKi6+MkSFNkCPb/MPAZj/XLz9sGfsfv+b5b0rAKaLQtsVgNcPRVpY5El/I4XLKzqi315+RH+ujegNwzCqABvRG8Zosz8T7QfgoLinAebJcKhUeOVLN7ForYRXtmKbg9vm4CdC9Gb5Ef3ZtsOUYZwwL2eif4W7nwTgl7wok3WTZZh1SE4dbjL28vvYevsxT8ZefB8v3s5TAFwJF98HsPV2donjfh+seJkXPgyQKTMZux6AFcq2KWq/pzDdOgXeHHQy9uL7WHU7lJqMdYa9eLuf+L0WLr6PttsB/gOsBVI+gTt3c288BfCQk7H3rOHelV44IHMPGYjEoX/nGn68EuB5uFOUKTUZexWgvER3rgG8fijSwkCYjL2jspOx1tEbxijxciYCLm6jqxUb0WtsRH+SiXaU7+jPMh+9YRhGFWD56A1jNLDMw8apgI3ojVMU6+KNMYflujEMw6hyLLzSMEaQ/ZnicD4sXGpso63VC23BFLZg6qQTvVN+Mvb9gxi2YWW0/EGg5f6OZ2+dkyztXL242Tp645QjykX8S/F08/UACx6Bz/CabPEx/xEpfk16pfPlFDjT5yLYsoL5X/U7QPEH8E+ic6FoIqVBDluuZ/4jPHA9wBcfYcv1APO/Clfyk48C7IOZ8DFnZL9PW/N0HR99RJq9kIHrATL3KdtelwP4/67nf3pEhEvgndgTOHI9k+QGt1zPfNeH/FPyQXnDHuHR6wGue4Qt17uUyUy6zwe3bI13TVtgfryRef1srWWedFMum4v7EnTCzZF3LBQgUjupHJWkDjlY8BuArRcU4yBdy8DWWmb7GCjm/RBg66e9QvjK0EVaWDSy4G9kbmXDK3eW7+jPKGdY5+rFzVvuyq9ZRucDS5pfvSu/5hpdvGFlzaO/nzc3pXEKEd0bRfdGDEB/8S/rPvb9MMBkmBxO3V8O+qBPnQZJH1n8zlD0Q0HpBM1cvCm5Iv0SMx4MyEHW18gnbOiDPiaimu0nUyBTiNvWV6xSIFZX3y/99KjGs0BBGV/qyfgtqvrJQo8bZ/f5WomdTWvTe50OkAmXiOs7SVZOI+W60Gr+4eT8sW7ZNd4TJP3QH6uYTRVl440ES0ZiZ9YT5nj2jN22pf3m65YBzLn2xpbvProhVrrh0e/efN0yi7oxDMMYIxxHR9/Zsallvk8YOqd5IZs6OpOlHStrbDLWOFWI1kXOc9JVV9xfCbWPUled8wn7HaA8gy6Ysh2mbIepCtOxpT3pLYvTftuWu/LW0RvVT7QuAtjqO5HGvti+UT5XwTwa++iuA7fyXnfT3ssrp8BB76OfTLHzbcopnXnx/h0ld7Xmec8/88SArXA1rp/ZDzMR534/fACcnyc0O4+ZBWWPvi4A9Vq4JJkD4Qxl1WRnObFnkngyk5SpjWKtq9XREKtSl5I099HRQLPkQNheC9BXC3hhV9ZXGYBMg583zmc5ChPDQz7DX2uyar+5zwvPEGHzRTGTnP2JIi0sGpnzkgsGqCTH4Txqnt8yuELL/beb68YwDGOscByumznNC9u3bHPHnR2bWNg8J14K2MpYo+qJ7o38BFwW6qBO5jrlz3+OspAnwv9pBT8XqSV1UAu1cc14s4P8+VqOrLpirT8umuT++qE/ZVUe8nHbtJGUEspfDXFj8qWND4bpgxqXmCznayVmcKPUnG6hjv5wibi+k2RKzCXH/jL4q7pj3bJrvCY0mIVsrGImVZSJNxIs8Qfjj7nzWx78ygOdQOdjD7Xfcl0s5oZrrvvAbfc9bnH0RlXzeiYqiK95ssTb/aNss+dYdB/AS7fzwxAk/rzP54VtJWhbCZ5Eor3lwyubjiWOvnP1krlb7iq4QMsNK6Pl1tEb1Ynb9HXRw5BTjuxagK7raVSZGl9sBcteadkrx0L2yq7yHX3jCRlmrhvDMIwqxzp6owqxhGXGuOR4FkwNCwuvNKqNaF0Ui1nUUY8WXglYeGU1hVcOD+vojeohujcC32H5IX2IOUGiZQgZBjz+wyVRN5CMiomlR3H/1nl3fzESRjerQ18STWVLtB+F0lp/XIy6AZ+wIXGtrMxFatty8Q+0FsY/6DXKAB91Q4nOIBimD3x+tJyv1Z9NVknErhSjbijqe508QCbrTwtq+8A8xTCYWm+xn3TtVy07YU0Q6sZlhjZRFBMG8skbqTKq+uYMwzDGETaiN4zBie6N/PA2BzIwLEbBEw/xyxUr+h/rOajxtQpaOaeGxuHfPj/c9Zq5eLM5JUw0lVO1RFgIpf3+2Oc7C63VqlpBWCM3iBrXq/uKCbXcjaRFUkCaiutA0TB94IfmWV+rNlUlIYn6qK0LPwR8qX/mNf44CDOiGPmFCqKc93UHVPtRH0BtHfkgzMVMGlBXDEUxYaAmJakI1tEbRjn81Ku9lw2jDPbhMAzDGBvYiN4wSmKRlEb1YB29YaSxXt6oKqyjNwxN9O8RwGIlClORmwC6nHAx5LyExT68ckcdZ+0q1vOau9hRV0w1Q1BYLJ+SjXJKMQVCl5NsAngjyzm7RGeTsm2TairYtktmNXcpUw/7ZDU+k9e7APR5Y45cALuk2U1sjwDOXRy3baO/VD7cwmLYpW4HgG3QKJIueCMLiPGK8GTeUabuA6Bxsa8VT2zDL1OSHXXMlUsAZwESuOmEXeAS7XbDBDWFHPIFtQOH/LXaRdm17IRuV1rkloMB7UDQ31VsAdUIYomrVdkw+tHDOnrDMIyxwaiN6C2pmTH++EUmmpsShrU2bmDq9r7+L1AQyen4dadTfsKmjxcrOs3rYeFP+OrHAXJwvZSeLsPPd+QUmIxPf/Zd+M/S/sLH2bTM6+wSTdRg+nRl2/UUU0I6A1w7TwDQDz3wB3JfjjNgpzS7Cxb+DGDfx4q2vV/tAd4rN/t+2C4LxQILv8GmLxVv3+XsdMZrwpNxq2G7nPI3APZ9iRku0+fiWJVVU7nnYEzCC3AJPCunbuztRpjPArQto9VVcduBh1ScE6AboG0qnwTgIwdpmyrKrmXgEvgOfB5g81SABaLQNhXw+qFIC4s862+kMKWiSc2Olk9qNvGEDLOO3hhPbMlEwMK/Vv1oILhuJEMklr0Sy145rrJXjlpHb64bwzCMsYFNxhpGtC46fQUAC5UrIBBG9CFxGJbUzJKaWVIzsI7eGC9E6yK2eh9306YSvVLCdeP3dZpXynVTR+PWYj2vuZWuOnqCRHfT5V03hyl2vnuzNIULDeq6ce33SWnR1DfobQGoh6PAawDk4ByATItqdit7IoCZ85RtynXTG667ALYkXTc7KT6Bw7A3C4jxivBkjihTndmN83yt5rgPZFNDUtLVQLNy3Zybct0cHobrhp0AzeeWcN00a9fNqzGTYl4aKSrnumn2rpvkQ6gOzEdvjHX2ZiJgJuyBWU8AvH4156XUMgUGIoA9AKwB4J4CA5GXzJTZ2pq/YfefFCt+B4DPw6y/oe1PADKwUkpnQqYAsDvyp8BOmAjAarhbrjjrAXZ/kWDqTGlhj2oq2LYSvgvALWLq3bANfoA3oAecjZFMq07bz+7TfLN7YJar/78VbZtV8AfAq3ARXvh2lEzkMutf2P0Zf7wG7nkA8MZrnGEr3dcG1DvlfwHgWvg/AfhIrMqqS7nnl/FWXoJF8HM5vQSQ2eGfA7R9kVZdJewQGCZjL+UPAZj/S9ouVcovAbAIdvoh/95LAZpEoe1SwOuHIi0s8nN/I4VLKuqj7y3voz+xyQMb0RuGYYwNLLzSODX5aia6GoAdcBb8BIBL4EBKM4RXuiU5ewHohYJIdsi45tc+jMXjtvmeDk/ABQBMhSYp3SEBjmfKKZCB6QD0QLe0/yN8FKAzdYe0cJZqKtjWpIIDnandMAO2A5CFeu+wIeO9DuyEj0uzZ8G/AbBQ2fa2HAAb4BoRZuF98Wf17/AJOd4LLwNivMYZ1gRPAnAV7IVfAXAhbAHUHuiOmmLmec8hmKJ2SJ8CyLokJ7xUbjBXzIlJXu0VfpH8pHgdLhJl/CIqpkhYJvLrIvx4cL9pXo0XaWEgLM66tbKbg/eVH9HXWXilUaX8IBNdD5m1gHdP710B0PTjoSdjO1YANH+1xGQst4uLxIlXAMxbC5/lrTqAo9C8NhRLrTlyCpzhJ2M7PkdzCHm8AR4umjrIZGzHCprX8tcrAP5srZjaBlfzRAvAfpgNi10cez9cC/DiZC5eq+YbbgTga8q2TjmA76/gc2tFuMT7uANdK2iUG+xYITOlDyefqDdsLQ+vALhhLR0r/GRs41fhSwAdcadABzTHG2nuo6NOTcZmQRz9Trg16weyOk1xXnXfAzDvNYCOC71OaBnoqKNRFvE2/wCg47NewX1PZOJFWlg0Mudv5ILKdvT95Tv6WnPdGIZhVAEWdWOcUkTrIuCsFeyBWSqUxaJuLOqmmqNuzEdvnDrsy0QuRvsrWe56gL1fBGjKsTdL0+MAbyzjnNR2SOT9z37XAblYmrt7oUY6spy42/+KvXcU630bgD+Fpq/TdgdA5FfUSy3XK9XLKezJMgGA++HunLTfxt5WZapYuDdbbCrY9nkV7eNNzdGZ9VkHCtArNoQf803b2Xuub3ZvlqavA3Cbsq3XHwC/gg/ihbvqkwO6pkfZe50//g7c3QZ44zXBwjCN8R24+1EAloEz4MpYlS8v5e6fxltxKceektMrAIm6eQpgVSv3uCrOL5OKulm1lBsAmP9TVi0VZdcycDkcgGkA+5cCnCYKq5YCXj8UaWGRp/yNFK6qqOsmX951U2M+eqNa2JeJcGvx5wH8YgXnwSyVV2BTK1gKBCwFgm+ZKkuBUCjf0UfmozcMw6gCRs11Y/s2GGMF20XEMEYJG9EbY4JoXcQ8lXRlK8AOuIxYIhfvsLFcN5brxnLdHAvW0RsVxg/k01E0hnGKMTBIR39i30D2Y9kwDKPKsagbo5KEeLIQoYgLX7HwSguvPPXCK8mUH3mf2JJdc90YlcF7bP6y0nYYxphhEPfMCXrvzXVjGIZR5ZjrxqgA++QnapPOHSbeA1swZQumTs0FU4XyrpvIXDfG+CJaFzXdICfxGEQLr7TwylM5vHL0Lm4dvXHysEhKw6gI1tEbhmGMCWxEb4x7LMOBYQyOdfTGOMYll9eu6mJoeXzGspgYfR7AWVg+estHfwrlox+9jt4GWYZhGFWOdfTG6GIeG8MYJgPl/04Qi6M3RpH9mQgXzqjoqqPx63LyX4tC8PKuO3yVe+u482/o+hNcI111NP4EYMfHOSveJhRTILimvgVAazfUSON9sr7w23R9sVjPaX4BGh+gTeRfkNLGPu9n6Jogp7CvzkVXch+09hWN77pDmSoWdtUVmwq2fUGlXvCm9vF6HS6/QAGOwJ8CEInVM96h6/2+2a46Gh8A4BZlW7c/AF6CRXjh7gk+ZjHQ+Dhdy4q33yqPPUF4MnsAmOmUHwfgKnhADhRtV9D6XLwVlwLhSTn9CCApEJ4EaLtDqrgUCM7/loUGnCut7Qr+EID5z8XbDykQ9vv40feuAJguCm1XAF4/FGlhkSf9jRQur2QcfU/5OPqGE4ujt47eGBX8fvY/AJST2jEPQhz9V5UQkT/sqzzcyvWQWSuNzGPvCoCmHw8dR9+xAqD5qyXi6LldDHPiFQDz1sJneasO4Cg065VcrtYcdS9n+Dj6js/R3KaMf7ho6iBx9B0raF7LX68A+LO1YmobXM0TLQD7YTYsdolv+uFagBcnc/FaNd9wIwBfU7Z1ygF8fwWfWyvCJd7HHehaQaPcYMcKCXJ/OPlEvWFreXgFwA1r6Vjh4+gbvwpfAuiIT/N1QHO8keY+OupUHH0WxNHvhFuz/pvM9fM+jh6OwkSRz3sNoONCrxNaBjrqaIQuJ/kBQMdnvYIbCGfiRVpYNDLnb+SCE+tPT5DR6+htMtYwDGNMYFE3xngiujfC+QFcDEPiXVZP0ZWQVUJEXu/lA5CHbGikHj+qqSv1zhUPiWuqoIWhhTq5igquKCir8kESFLKpRoBaLyyIYbFms8VbINFUsK1ePtXhprJF4YDTcdeNfFMF9WSox20vGmnb6orXGtDCutj9Jm6wgHrspNRSphYfbA1IrsnAQErin6q84m5kn5ei4lWgoOYM85DTTR0uca1gQCEU1ZVSG7KI4ktfWayjN8YN0b1R8TPaC8iHO9CrPnA5JUTVykH4IR8a6ZWOoC/Vpm6tF4hrJszIiwSl2Qt5ampEEhQSFrrTft/LRKj28zFTkxWVpms/I6WRKnXCjNPpk2vllFVyj36z6D51ib7itTJa2Be738QNRsRfrIRaytTig80D1MS7kAzJ+QD/VOUVz4YvLbluJK6bKF49K8cD+GwGNcRcN96AGqJQq8+rESoG/b5iC6STQeaTN1JlmI/eGDEeyUTIR9H5ijcBMCGu1g2/I8fblBCRvyJV9sNieFEa6YZPAPCPcGbq6hn5DLumfgvAAhgQyQTpAs4Vw1Ca58FCeAzcLZwnpROkUziq7iUrXdUmmK+Mf0WZGu66WzUVbDtPXdcdzIeCHPdDQdng+vazYKc02w0fBsRj72ybKAfhEk64O2StET4GP1Oavw+I8Zpg4edyAN/P8lv4XwD4JVwAJAP0eU4SzgeaYC/MktO9gDw9J/xnaJHbDD+f+iAStXZ4HwDnQ7sou5Zdg2FC5F1AlF1FRD8UaWFgltzI71bUR3+kvI9+kvnoDcMwqgBz3RhjnQ2ZyPl43YjexWa48UnC99vrUjMC8KYSIvJIquyE01UjvTK4bijhTy6O6LUHoh4GRFIvQ+PGUutHImhUn7SgUC8X7VH3MuDX2BYNc8ZrU4OFvSlN135N6kKHoV8pBDeOG7qeDe9Ks70yHN4ety2s9zxbGTMl5anIxJ9AeOwlieDH2aJCRpp1tQ7FlQdSP+Aa4Yh6xd1vjry6LlLFDeFd3GoNKH+Pv9kJ6l/kPXYEMqoFUgaULErozEzdSJVhrhvjRLk34928zl3jOvp3APkRnXDa5tXv/TOUEJHPkipvw4DouI/9dgCaS82qhY7eFTknyRQ1PVgjHf1u9ese9TXwruyLt1F1BDXSteXUvRyQvL8HYbIyfpYyNdx1XjUVbJugnMXO1MkwFTYCMABnyigsAztE4TQVfbhHPcAQeh7mBbbDuSLs8hvtFQmlzhjXx80iSXB5bQfgXOgWH8gF8hIn5kr61KS4vuswoszLDSLCaRIcmdhJMLyg4Qt4H8VIStSjCBd1Dr23RcF9E3TFi7QwEB7dlyrqujlU3nUzxVw3hmEYVYC5bowxyl9lfPRHPT7j2IAapbpNpfvjVWqVh2GXEiLy96RKHgZUI8Fbsk/NcAbCINqpuRFcJ9SIpF9+v3dLm44wp9ctA+omNSrslzFmg7qXvCT8mgPblfHa1HDXtaqpYNvbasdtZ+p2eEumo38LkShk8Pt+HIWMNFsrPpC3lG09cuBqvSfCSfH7dQTJmfLjIK2TGCC/B2fKiP4tifdJjDPnw5ZUO6iBfyElnC5NafQvtjDMd0ub08qIcKc6Rn6KJYoycZ1gSUU3BR91rKM3jp9fZaJe+AwA/wqXANLhus/VuUAqbK9euieUl9a5np18tlT5FcxSjYTvkqZSrpvQ0bumngfgQyq6PZjRBe9XFd1i+CvgHYkzOeJ2tpJa7loTVSMF+U7aBh9Sxs9WpmrXPErubLsMfgLAImVqrwqdfweWAFAjwi1wgXI0vSmPItgWlpICr8h3xtFSPelu5ah5XuJnZpPEGXaZ5Bpwpi4EYDLsA9TXmONtWBCXuNclfO3pmQknfEOqaB99j/re+g0sB2Ar/Ea1n3byTAGUt/03gOiHIi0M9KdupCLYiN4YW0Q7IuDlswHmA/CUvJkimdhEPnIT43UPhe3o1PvPaYZarsph+IBq5BCcBcA7MCNlUr188p1+v9TKqxbcB2lWfOatX2rNkl4ApTBRemd9LwOqq8oq47WpE0s1FWw7JEPLQ2JAFiIZL/dBQX40ROJhb1TNHhLvfK26xBR1rYwSvpf6XtRPoD/+YmnCk5mgTA3fEM7+hLO7N/UDbhocUK+4/pHnhL+V04Kagu6HXvlWKKjnVlDtu8dyAGbLEKFXNY4M0vvjRYW4DmJJ4kYqgnX0hmEYVY519MYYItoQuRHRJWezQwZTYYmQ+6ke8h2SCmUD9svBtFLy09S6JOd5CI248eCUUm3m4vEebhgeKWEIdNkNU5Wm8w9MhIPyeahT4/EJ4kqeqMwI48F6ZYmkUExGAQaCpF79vAi/GCbAPjFsCvTIFSMxoAvOU424Qe4MZVukzK5RwsmpAWyP0qzH7ysS3z/cFyG/ioKpPXLsaiX6pgkpSQ4GRDnou3+dMKsCkAaUQkYF50wTYVa1nxNhTzyMJyjo06xqIW32kdHsZMcCFl5pHAMvZCJgB1wKwKPwHnwQgI0Ss6gnY90sYmKaK1JTqfVKiIrnc1X2QYOs6nS/613XMLHUxzL0C64ptwKzUblZCmpOb5Kq2CWaR5QLvklKCyqePdxLncyFHoQZyvgJytRw15FqKtjWJPMNp4upM6BbbrYf9snXZMgKsA+mSbNRPK5fZ39EXcIJ61Id/SHxWTtN9zTS30zOsCb1iPbKFSeLAYkXIu3sdlZllAIqchQ4KMbUqgnzgiyUq0jA4AAAIABJREFUdda6+YMDccvD10NiHj7cbHBe6SItDIQ3z60VDa/cVz68coaFVxqGYVQB5roxKk8mE7lf0G/LWM+tOz0IQCTDLjeCc8PVWVKkKUh11MyeG66EIaSr4ibxTpNGCmodVvqNGwZlrin342B2yAQJkQSf7IsvDnKj+NNhtwxX++F0KQ3Tg3l1LzlxJuRhtjK+SZmqB/IoubPtdPFuna5MfV38SJPhPbEh3NdutU64IEEvM5VteinpHqleA9lU4FNO3WC3DG+bSNItFh5VptaKtWHeWJN23WTjXjU33HY6Ia1Ng5zWlFoGFWaA3UxvCCGVzWWLT9s5svbFr9UQL9LCQF2ZqM2qwTp6Y1gUMtFr8iGcDQcA+dd5G2rVz/OMfK7cr3vtJ3HCGUoTpRlqTZKDRullJkGPdE9HxHGvSTR1RAmDGQW5it4a+rCy6rV4I65WohFgoty7W38fjA86PequdVPBtuBW7lGmTpJmXY/WI3LXB82JP5kZqRucrK41Twm7U73Y6UrziN97PPZMQhES6RhMvVCu6wzYH6+yS32FOKaphcTILIu7dyd8S75j6iAT777d9+5u+QKrh93qCylE3UyT59YtakhFRL9btUDqW+201I1UhOMd0W9YGS1/EGi5v+PZW+eU0jAfvTEEXxa/YclRzwTlcNcJVfRbNpFopWRR4i0enM5hFOwGtmG4l0wzO2j7icy0x/RxStRypzVqKX+mlPEDpS6X0Aw9WtrUkjYMlFJLFA0ydaHJplIXDE74GYGYnR9Ls5fHejuDcHdFffR7yvvoZ5Y1rHP14uYtd+XXLKPzgSXNr96VX3NNWsk2bjYMwxi/bNvSfvN1ywDmXHtjy3cf3ZBW6VxtHb0xGBszUY34eUs6MbWze6DMpvUDwyiilBzl4A4Jy/Lxium/ku0Mci1NptTvj3QjuXggYNr4kpdLaM5UC4ORny/BgMTPIx0XOJAqKtn+IMJjHf8W1AuRk9svySC/tEaPkRrOV5xhvqtjdHZsapk/1x3PaV7Ipo7OpMbqzz1kPnqjLN/LRIcl4ckB5cPNqS35dshxIozdKbuIvbRnI3h7dG4WVHhl8F/rRnSix/RnO/goQh5doFd5UQbUVXQoYdjKo1tu9rCavRxI9bmojs+5xRPNJvw52ocTbOuVycwJYmoPZMUFvw8K6nm655ArlfMyHRNJ6h51rGe6tF785ummwjMM2TG1cjrdEPFUycOkXFNaIbwfBldOhFeGV1YXaeFY43i8YR1b2v3i9NJ0PvDHD934PevoDcMwxi3N81sGKe1c/ccP3fi9Z+dYR2+U4D5JMV+v9gbRODdOHubKdoDBt6PHjyFuRKOHLX2lhD3xKqGRAfkFcKSUiyAMohPtD6hRcCjSo85ISdxwr0spBF9KmBNGxfP1yiXSpg4+GdsT/9WC7G7rQoBqROKuNSDCHvXrhLgOqZF7j7q7xIh+Qqlxd7mReAj1OSJ1KR+MWHvsI/oh3cdRfD+WNGGYnwgs0cr5UsKxxvGM6Oc0L2x/dJvbDaKzYxMLr1NhN52PPdTe3t4c3WZRN0aSL2eiEOAMfAyAAzADPvoSQM8iGjYB/NNCXoFPA/BTmAdIJ+ji6N22pYmfyRPUl0dijyEnP1uq7IQsnC2NTJA4+hmlPhJheaRr6mUA5kKN+rXuau2Pbznr9qS9GN6WlLy1siOrq+U+IZPj9+KE2+BCZbw2VW8Vi5I72z4MzwCwREy9ELIS5PcOHJHNV2ukhz0M71MeG2fqNGXbYbU37OtwvginpXpAnb3yZUkudjZJnGEfVo/oZXl0EyVcPRGQvju1gUl/vPd33xnOHifcI7nq3LdjeO+FF3QHXAVAJ+wQZVTG5oCL1wyBkm7DlrPiRVoYCJkyb6po1M3O8lE3ZwwadfPQjR3PfHHOIFE3NqI3ikQ9EZCZ4D+NrqN3b5ELXKjyFIAGyaxYgPNFrQeWAfLBdlP/n+gHyMcXxdf084RIPqGE4OWf6PdV/l84QzVS08/f1wIs7yWf2ksw9AuuqcO1AMvcGhxpwXUNT0ubjkNy0R/V+s7ugFLI1/peqU7dyx4VJr9MGa9NDXddo5oKtn2in5flusHUo7DeXQumwdXueynHtlqAberJ1PTzb7UAS5VtfbX+APitGNNXSwEa4jkQ/qG2eIOHa/3XzyfiOuEZfqKfXypTj8pdPwmkBshd8RlmpPsOvVQ6n8REldsSFUHbqxIWXXYU4ODEYmJU1A+1PXLRjwLwtCg42xrjRVoYKIyNkf5xRazOufXZ9Vui5prbXBz9NQCdq5fM3XJXodjlW0dvGIYxrrlmTSG/Rgvm3PpM/AeAdfSGJ3oxco6AgbkxuYu4mOBWqL4OwJnwAsAFV/AvatGjG0m5kbVbf/hcLaSctg21PuAENfhqqAUJRHmu1lc5CjtVIw21fiT4TH2JqJswondNvSnt14gkmLFb2nS8IRfdLR6nOqXQI36GSaqRvITlvAlPK+O1qeFyDaqpYNtztd4h81xt0dRDskj1ILwLz2YAamr9yPdt9WQaar2pTyvbjsiBa+E5Ec6A3viPqkjd4JuyzPi5uE54hs/VeveIM3WWXPddeVaa2lRi91w8NjdsWoAIwwbCkVoN26vcKW/B1okAXWrLGlQsUCQXfRavhlQMhCItDKS3uq0Io7cGzTp6AyB6IuL84t4TR1WRy0nZ796Fbjemgnfoup/kFwGwKZ6m2DmL9QxnoFu5krX/GlXLyQ/CWaqRbvlS6Y7nGXaEZIfdcgo0QDYVtTkl/t0T/MVTpBeIlEJIU6zvJSRZzCl3/OS4qYlbQ2nmVPs96qu0WxzfLnyzV4xxX29Z9WS6JQNB2I/JHafTLTTA0VT2Sp0EIhd/sTThxsNeiTl13cSkhaMvNR8wJT55EKbxEeFOOc2ri+bVmol8KSGSoucwzJTEQT2qHeLN9qgWSBk5eWwEXFpHb4wW0f8RAVwGdcXkrXq87OJA/KDMfTSn+/lB1324j25GfRpDoIv+6Ab06vmcEqJqOXkEA6qRfHzZToKSTeWUJKfCcvTnvJAyWyvkUo24Gw/J0HPqitrUhD2U0tRPI1fmxsOzKsS7vPQNZlSzkRJmUv1aXrmkEy+WJrFILR+/x2EuxSq3rCwIh5z9LMQfQqI6yvJEU4VSxyUvN3o97BjBOnrDMIwxweh934yFqWajYkSr7A1gGGOF40mBMDxsRH+KEnVFAP9ZvuuPwgvwGYCMWpeEJIB9n/tV/p8Ajs5g4m8BNsIscSvnJRDb+eh3AuCmddMLpl6S47OVEJHPlSqzZE0WstroKQCWyLyrJkzGuqacAedAjYrDC/no9XzzO3LRjTLfcEgpBL9HpO6lW+YM34WzlfHa1EEWTO2EudAu190pj+I0+HdAgt8vANTDPEM9mQE/He4nb8OqqPC9/YYYU0hlpwB+Ah+X451SKz4H74uc/JfKVDcJfJWEqyfCXDNlJMEw99YKHipgplSpkTV6yLOqF4Upcq2Zqv2QT2KXKLilQltEYaayMBRpYSBKSaoM6+gNwzDGBDYZa4wkL2SimbLzgptQPTyD2tnF9Sl6wupcAJ6HffD+GQA90HcewB7oldnaqRKy4saVbqjuJOkRfRjIv6WEqFquyiE4XzUyICFAr5bZeETP+20H4IJ4CgQ3kDw9HmP3plz0dPgJADlZVkqZFAjvyWh0O8xRxmtTBxnRb4fz5YfCW2LqHMhJWNEA7JefGhkJZDxPPZkBGaLuUrbpF65JRRA1xMOoUE/VGfN7YkmC7aJcUKZeKdd1L1ai5VflNQq4bQzCXuQ9ymAnfEXWVJcLr9wsYb3vg82ijFpDO1Nejl2ihlRE9HepFlCNIJYkbqQiWEdvjBjRD6NZnyl2lGH5a0YlDHEfoW7ISEffD69L93pA7YKNRDrm47t4u4WjLlwzvZVgCL0/oISodeqRtNkU78UmSrMhXC+Q2ErQHU+GAdWC6+inFeNIi/qTYIqqGxa+ltxKsEZtiTdZGb9fmTrIVoK6/UnK1D3xTa4nyX01qVrBS+NcEMFvRnwrwXekeg00pL5r6+JZhlz76S2WBsSMCcrUsHmhe90PpaqcFpdMULsVIv2ps3NKvEpWwlWBbtivvvbc6z453n7QDA/2iNip7T8tXqSFgSmpG6kyrKM/hdifiYADn+H81+BXIp0KMGsHfNh/8qILY73AzwCYBmfDrMcBZi2G7wF89L/wE+k9D8vKF+crd6NRl+smsWFpvfKizlJCRH6WVKmFjaqRevEaz1SrbALBR++acl3DO1ATT1kMdMRz3bj+4hC8ploIH/te6ZUmqnsJC6Ymqp1mt0gGlV5JRKxvDbWmf6Jq/5C0vBt6ZWn+AcjIt2CYjdgIc6XZetgkjyLYdlQNnN+TSxyF6akcZO+qYe9EWUyUSP+C+lrtVaa612gK/BpIef+nqe0YHUMumKqXKukR/bui4L7XX1PKqAVTAffds0NO66WWLtLCQFgwdRGVxEb0hmEYVY519MaJEr0QHbwU3G/Y1yQ1IjKsmgnr4b+Bygast1qe4AaYLhBkMlwKsBvOV+kDrwLiSc2WDJXU7ColRJKaLZHUXdthpmokJDW7ahhJzbpqffslk5otUUtF98pFf1Trx7kHlMLgSc0K7hJi/JLhJTXr+v/be/M4u6rj3vd7TvfpVmtoDSA1g0BCjQY0gM0skMAYzCAEHjC+L34yN04cx8ONEycOL7GJ+9PGQz6+ubkxyXUc573cXFt24oGEwcjCxsZgQIgZA6Il0WgWtFoSkrpbLbW6+7w/aled2meflmSBkHRUv48++qyuvVbt2mvvs/bav1WrqsC8fUnIgXn7SqZaULN9cApcng5q1uV6pmYfPdqKSkHN2tSYoYKa2QVu06Bm8zJBzbaphY84U4V7edcQQc1e04CahgMGNcunmyRxLGCYPntnwKzdAK8N5wxX+YBBzeRbcHT6kBcajpKgZocPMdBXP/LL8wC/ZOQFIETE2Rp5Fn3AN8A1sBKg+I7UYmyd1mpCV7JGJ1/7p1/MEreyKr6PMuBuBeCxAmSpm0KyYd2aiBDdyP5YIWnSDR1OSX0hGQUerq+wjbNE3RQA1qj+GpWYGVtVp2CNnnSrLkgWXIUSdeOUDOqi4hp42BnvTS1RN06V2faYxrp5rFAydS+cDMAb8Dr8WmPdmFOm9Ux9ISGvHna27dYCqhmjbtLv2lp3gWt0lHwsXcf3TNGZeqqet1P7yqNGlw0MA25JA71H8tCJcLOO+0Mtxq6HFcMBtsN6t95gKxO1elIZ4s0AWWbPpQ95oWFf5kKOCGJGHwgEAlWO2BkbOETE3tdA4FjB4dsZGxmmqhmv5PLAKwBMVP63BibCHK0jX6zD4UX1snhI/e2EfpFP6VEwHCQ75TNaOOMnfHlhSY/5DprHutA+ZWGk8o5FKTghzqXEAo3VpLd99qm12VR2RrPK/xJkcQvknYacnsXnRbL0THvU667Bff5bHvBady0WZqsJtjrj650l3quyzMJxsN1FlxRTt8KALoK0u1iYuOhmtY7Xlr6qcbb1uy90M6Yf6jMBy/a4HhinvkzZraFiWJmpm/WMZTHXBCepx3oZsrsKTDjyQK6No/S63jiQH6QlLbGGaH075IXeEjHsr45ohqnVQ2eYmvrmDAvqpjrxQC4PXP0cwDQZ4N8LswG64ZKfwnlaVYjSTZz+cwa+APDr9NAsFSfBazBjE8BppzLieYC7FjJF3wpPwWdVXwGWAPABIBuPHu7V8nucEJXfoE1+AKPhfapkGPxvAD5aKabugFvHQw34HXUkFw3yengYFrqG0j3Xw090aNsBH9OjpXj07lq26J/LYZEz/gZnaikevVNltn0M/haAP1VTF0EPPAjAaBgBnwSgoB6NbfABN+beDcBVzrYe50L+d/AnKhx0buyCf4Xf1fISdbW8gXKIYR+DrwNwKyzRIeNm3VxWFotiHZyflgy6UJo4Oh4V7tNukSmCvZX3uBt3Ux/AfXWpV5S98yx/oWzmsoUf6yt/yAsNA5kLOSIIjj4QCASqHDHQB34L5J/Jny7zcJmSXQ7AUppmg8xlpsByrS3uchvhcmqWAxQvSn3py2SwXxxCHgEYcUkyDT75HJbqXvxBeALQyZG4ZspJshumbLv5E05otizXJm9ArVNSr1tzH60UPL1sw9QqABrSG6YEe9zVW80TYY9uHSq4nKIVN0yZ180qGOmM96buZ8PUKjhBQ0csVwNGwl79POqEjRpKzGaa3a5nrA+fcLb5DVMj1ZiKG6bGuh5Ypa4+yymHGHaCzrLF1DP0vLu1rzyKmY8tmdHvx+tmvT5CsrPXZvGmpx1erkuE7VoZN6Mfq5Wf1LbWEK1vh7zQcJR43Rw+BEdfbViVy6+Hq1sBVrVAcPTB0QdHf4xw9G1Dc/QzgqMPCFbl8sCMVmYspK0FYMZXAaYJPfk4TZ8D2AnTVupsjdKM/vSLknH34YtSmefeBUAj1MK0HwDMuAS+C7D88wyAZJr/ue5JkdFZtqHLsJWd0a/V8slOiMona5NnYLxTUq+jzISDmNFLEOC5lWb0W9Lb/R8D4HLYoJ7jA/oVxIFm9G1wiTN+sjN1PzP6ZXA53K/nFVMvgb365tsO2zSSsM3oV8J09y6UOF8nOdv8jP4FfZdXnNG/rg3FmGmAGu+xTC0UJ/T5sAzOAaBRo1CUTYRfT1aCSjiYGf1sPVQ2o5dyO1wHwApod/rtnWdBfixWEtoQrW+HvNBwlMzow70yEAgEAoeIoG6qBKty+RmtAG0tBHVDUDdB3RyD1M2KoambmUHdHOcwxiaha4K6AYK6CeomqBuHoG4CgUCgyhEz+mMb+e/lGxYB0JzEIKQZxmtZ0oBKzKppiXAvMNZNquUjdg6Mh20AxfR8WT51k4SrMh06IZkEyvxUPvBrdDopM2vZuWNMi0ev29cz3AlxrUS+x7E0w6HXTffKog/iiAtRJZOjBqhViZ1rRHqGW9STjnD5m6zCcKWh/LXs0ylSUQ3zakVSdmm4ml5/nxow3G3Oyrk6FhEs73qmV7u9wdnW4NQOOOFeKAtM6WsW0zfLw3rGcjkV9YwN2qosMVO/dpdByBa743u1Gios6p8D7qr7nSqLAdDvKuNImJGww12CVSi6P/ucBjJGjqj2DFMxoz+GEXFsAoFqwuGLdRMz+mMS+e/lAZqhXSeMrsBZ9Ers4NWAEuGrkgp7gDccI7s7OcrYhKfMpbeDy4y/UaZ1ss9kesLRy+zMstlZlrgazdx2sjuDoV6P+kPycSDy8SqXZKd1WrNemfHGSvOvMo5e5ua9jqO3Jj06HcbVrIMeF+LcEiftTnP0osQ4+pwaJmrHO1PLLg1X0+uvUwN2u54vujo1bk27TtXWa0TlurRt9hlR4/ptTCYq73b3SZRz3V4G65k+Z6o8RQVtVUZt12bGlB6odXfcrx736Flq9dCAu1hTZWl7a11lay5fbLWuN6xCzv1Z5zRQycijgaM/fIiBPhAIBI4KHD7qJtwrjz2syuUbtNwL5lUZ7pXhXhnulRzL7pXPD+1eec6bMywG+mMJG3J5YNJiXRQFmlm3CBFagctYdzrApB8BcDYAK9l1I8AOOH0L/EA12GLstGQxtvWS0tgHfBDQAISzvwrAhcni13Mf5G74MAA/1kiKMjoLsSARY+y1JOjVtNc4Jz8hBEQ+Wps8ANP1rdMAvWpS3UEsxj4FwFy3GGtmdOrgLpAoKBdAt2bQqnEhGBv0pbXPKdmnfz4N853xo52p/mWMk4tt52ugyivU1PmwXcej1dCp7qo1LhPTJFXbqxzOaGdbwS26roCZKiSzGOvNe0rX7LNdKoadD88C8E54CmYB0KjzgDIOzTtuCka5BOi47OqocIXyiDLa2sKvvY3WamTNF5wPK27IHqOLsXIJ9nStBbS+HfJCw3C9kJuP6ED/3NAD/TvenGGxmhcIBAJVjuDojxnkv5fvFU9KH8EvFmMVsRgbi7HH+mJscPTHO4SXn6QsfHD0wdEHR5/Fsc7RPzM0dXNucPTVDc/Lr2uB4OiDow+OXjVTXRz94Rvog7oJBAKBowKRYeo4RS6X/w4AH1/EBPBlwRYntMJrmnz1IzeDzp7qpiet9sHaCZz2R4mGmn6AF2oZ7oiUbmeDTCQ2wQ6YfSsAu1g9DuAeAKY+CDDpihJ1M6izUZksl3ECPXoUWOmEqHyZNtkHC+EXqqRH43zdk0mCiqM1RFWXE/aoBumKC9Uwf7E9ablxyvX6I9nhrsXO1eX0z9coYGJqfSVVZluPOyqmDkINvAjoD16m/0bcz4XVbj1AQtWtcLadqAVgqp7iRJ16e1yqubTEGOv2MnTpuU5ypgprt1JbrU43kShvHmfAGv1oQL9FB9VI4BmYAgwd1KxN6aAp0KaV0VxXa2C6PkjbtJoZY39ucxp8HdSSsgs5IoiB/nhELpfPtXJLC0DTYmhPlQVNzdyyKBFaoekyxg5N3RQqUTdzHHVDmrqRj8nT4BQ0S/SFTN0BcOMHuRtWXwGwTn+3RchrvMMTnAbDCD2K+8yXsCciH61NCvATpW7yMEITUo+r5C5m1I2oGuWEI1SDtHoiTd2M1FZPKJlQ42Kz5NOBWSzs4oCeJe+MH+1MNQu9KrNthDs6So8OaFTF1dCl7IpRN8tgkqod4frKTtHv1Bp10w9jMtTNMldzVFqVxyhVLqFOT4ZR+j5o1FZl1M0YDatp6IQuR9x76kaEw7TJUNTNMH1OXnCVRTPQBdvT1I1VGOb+tEPD0nVQS+RCzst0QnUgBvqjERLhoGMRTc3OlyZdRsvhdYM7L+F1k24VXjfhdUMM9IFAIHCUINwrjyNYhIPvwC1pl8pwrwz3ynCvrGL3yieG9rq5MNwrqwaWK0oYko5FNC2mbREm9OUE7aXEUlZgIW3ngWaYQjNMdew3wxSaYepLF9FIUhO4BkhnmMJlmFp6dGeYOuegM0xd4jJMFfRNxkFkmDr3kDJMzdUMU9eoqee6DFMvv70ZpuQuTKYclqXLZ5iSJdD9ZJgqc688mAxTzXpoqAxT4l65Ik1bHmSGqeb0IS80WIap9x/RgX750AP9ReFeGQgEAlWA4OirHz5XlMw6JgDtLm9UuoyWGzKFyDBFZJiKDFORYcohgpodFYhcUYFA4PAhZvRHGLlcHii2pLY+yc6mXCvFlpQwNkzFhqnYMBUbpg4BMdAHAoHAUYGgbqoTwdgEAoG3AeFeeWSwXv2oJls0gDd49EaAS+9JfMHWzmfyr3l0fknoy9lWpebv4tFGgEvvBGAeAE+yYiHAWljQpd/8hj0wI1mza52YkDziP3c1AGNhEGYuB+i5iBE/A7jvap5SDuHXGmpBFmPlS1kWDMu2ZRZdfMExToiLXilN2uEEXXfNuYajMzpx7pWiStxHz4Qat99VCl1p5ke++qdDF2wCIJ8sNietBlRo15LT5c12mO6M96ZmPcdzzrZp8AIAc9TU6bBL1yeFKJul55VWG+E0t1y83XWg32og6NQ1+EEoZLxRfQ+sgiag0vq2GDYNXgbgLFilvpgjlNgpm4R2ZnbYCiFjvZFTq6y8SaJraPRNC7dp3qiblSZqh81aGecab+6VRvugDdH6dsgLDXZr/usRda98ZGj3ynnhXhkIBAJVgODoqwr5R/KTbDnyl1poZqIUdiXxvyfNh18ycX5J6MvZVqXmO5goW55kErAjKUuFkUAvPKEazL1yh87DJ6ZczUSTLLHKwuiILtgFcCo8RaK2AKcCOsOSR3aoePT22JXFoxe5xaN/GRqdkl6d4h1MPHqZGJ5yEPHoC2pqPXQAUKMn5UDx6NfDqc74g4xHX4ATXBeJAadCg85S5VvK/FblXD1wgotH36B9RaV49B16CRXj0Q935hW0ZrZLrWdGOFNl+j9UPPr6jDftAePRd2ghBwXnypnXaXitbtV+DWqdqgPGo691Z7FDXmgYfnS4Vx4+xED/diOXy6+dz+TW5M9iS1IIrxvC6ya8bsLr5vAgBvq3D+JJmWtl8rtL43tOR3wfcFh+DcUWIkxxhCmOMMXHT5ji8LoJBAKBwCEiZvRvEySLCFBsYR1B3SQI6iaom6BuDDGjP7YR/vKBQOCAKA79700iBqBAIBCocsSGqcOLft0Bcb8Gdr8f9sB7tcL9Wuh1+1NkFe4auB+ed8Ln0xXKWllht+aZmwS4KIOyNDoMtroo50IFvA5jNcLfKjeDqNX9L1thH5wPQB2sB2CzY2l+rYHjZa+QrHDKclnWvbJJy7ucEJV3aJNX4ELdxCS+iaP1WvbDOYoqaXWmk1jIxnFqGK7mqTBee7iorqJkvCTlz35dM1wFZzrjO5yp+3Gv3ASnauUmZ+puFQ6oSag9wHjY6twrZePSjrRtdq7z4SkV9mcWY336jk1wLqCn9rCe8YUZ2gPid/tGuskaZVQMo2Gn2+Envo9C3YjwGW0iQfwtiYolb1mjT3I+rd88Ji36v7mZmjFofTvkhYaxeiGLjuiGqV8OvWHq3bFh6qhFfy5fp1z8wpbEwWZhC+vAywVblGe8nBJHv7CF55zwuXSFslZWeE3zflgmEKBOW22EM1zmjSxHb7w50K+c7GuwA64XBxTl6P8NgHc/CPD6FUkyB9mfKWWhvH2CVmC7S+ux0glR+TJt0g5XKkcvWZmuAOAe/eV72LggqmRsmgaDKhmnXTEtzdHbKHaB9jBukDUPn63uWnLupXimM36ZM9WuertTZbadmhlGz3TcfRG6dPemcfRTYLXL+uQ5erHN0/Hb9BIOyNFvGpqjNwtlRaQRNsF1AKxUhr2M2l4DZ6UlB8PRS5OhOPo1ICmQR6X1D8XR2x4IGdPPSh/yQsMxztEv/UT++n8G5v7PVY9+piynCu13zJv62RjoDws6c3ngIbhcx/GmxXSoc+TkxqQsfyaFdm6WCpZhqoWmxdy8qCT05WwrKzRjsdwSAAAgAElEQVQtZNh5AJMrZZgagOZMhqk5LsNULp1hSoaeiZKBSGIqXMLUbwFc83mWwvIrAF5Nb5haCuiE2rakC8bqUVyGKZnciXyyNhkFP3RKxsISAE7J6MSFQBBVcgldUKMSa/JgOsPUKG31oM6XC24GajHivRLLMNXoEhAu1SRNIrHTeVVm21h3tFGPDsI5ALwM3TpjtW+vR2G6qh0LdwO6TDpWL9bUWoapLmjKxJr/pXNsbXTdXoZGVS5LxPOhUT1crVVZhqkT9WvPsBbyLoGCzzAlwgZtMtRA36DPyQpXGU2pkIdN6XtkFRrcn3aoIV0HtUQu5Bh0r2y/49LruW9g4DravzFv2iemDXzrWnd06Sem/vBDqweCow8EAoFjF6tXLPuD910H0LzwQ3O/fddSf/Cnd/3zx2/7THNQN4cB+e/lJywCaHAcC+1abod3p+QJfIqoyDAFRIYpIDJMRYap/aN91YtzZ/65lJunzeaHq9q51uib6741cB2E181bjvCkDAQCh4ZDca9ctSK7uJJFeN28ZdigK+a9IKEHcu4O3aI7m26BjS6U2S1a6IUZrQBtLcmk6TtwC0xywknpCmWtrLAerm4FWNUC8AoAE+E+AGpgohK4KDU5HF7U2eJDGnVAFvFk2j4Khuti6TNaOOMnfHlhSY9tmMrpjElmgmXuAnmd1uEYXum7vXpGaSIfBz4Gr02Zs9NPC9Ir/8uEd4uLfDCoE/a9+mUgsHgPxqE3uBXUQTWm1l3LoN7ZJtjqjK93lthV550qs227su3D1NStMKALEu0w0rH88hMtuo1a1oc1zrZ+51xhxvS7nWKGPa4HxsFrQGZfG27h15u6Wc9oVnmc5KIdeGSDNptwpHMBqohRel1vpP2FsrDAxdYQrW+HvNBbIob91RH1uvnZ0F43Vw9lWPsd827h/yRrsEs/kf/yzNWPZBZkl8b0MxAIBI4KHMqMvnna7GUrEpeh9lUvMnta2Sjffse8/PUxo39rsCqXn6H+M7QnTjUPweUqa1pMh8Usa6TjxpLcWrW1gMzQxetmEU2LaVtUEvpytlWp+ULazgOYUcnrZidMy3jd4LxuvpT2uhH3/0aohWkSCu0S+C7A8s+zFGSN/+fOm7NG/Ttlfmrzd0G9WyA42QlR+WRt8gyMd0rqdTo5ITNFxXndiCr5nj0HalRiZmxJe908pte0QX3vCs4BdG86pq4oMa+bNjjXGT/ZmWqnq3eqzLa5uoXiGjX1XNirnzgvwzaN92BeNythuvvokbBcJznbdruFAfO62Q1jMh9Arzuvm2V6FyZTDjFsrsYvkxgP4tTYqLGUy7xuvGaBfMfYIoH3urFQz816qKLXTTvcAMAKaHfPvn3c1KS9bsy1X1aUmtOHvNBgCU/ef0Rn9EuHntFfO6Rh7XdcOu2HH1r1yB83t39j3rSXb0t53bTfMW/qZ2ffNxAD/ZvFqlwemNHKOvWYDOomqJugbsoQ1M3B4JAGeir40bffMW/qitsGvzX1G/OmfXYZsWEqEAgEjhIc6oapa781OPAtL2j+zCPyXvjjRwb+GMLr5k0ifGwCgcBbhcMX1Cyom0PH13J53Ke0eebaZ3tuiDtUUe4ZBv+RVuuEtZkK+TSzIYXZujNebq0koFjhTpr9vs7BTHgJgCnK/9a6iyqzOa/C235F67sSoewqFOpG0knP0GoeRQ0Vi9uJLsrbtJU0ednl6RYDngHgkgzv7y2U/8233X4nFregXQ0TiLAO2uAmAJ5ykWpylfJZdyoTLUnyzHi73uJ+k4NLJByJtzPeufkXNQaDbJG1TpZuKUCz8/Rf6TrQHkJ7MJ7RCDZ56Mgk7LajYsw6IN0ndgho0G4/F3r1zp6vz1gZobAnTYuZVTYh8vdIhFOUXRzU7CtAjeuuM5U72ghnamXcA28nlXTqL2kF4cFeSR/ywjIjgT8/otTNkqGpmwVvzrAY6A8Ft+fywO8CMGkxwLpFCYEO0ExxEcBPQBnsJLi8FHg3xfklubVat0gVugxTQv2L0JezrUrNL2Pd0BmmdlTKMIXLMNWazjD1QcDyH8nq7oXJBpXnPsjd8GEAfgxXAUqsS1S1oXLGWlbPspyxIrecsQ/AdF0wkIhdYtLB5IyVkF5zDyJnrMS9uQC6dYyo0fBtHChn7NMw3xl/kDljn4Lz4UEArlBT58N2pZJXQ6cubtuGqc0wyQU12619RaWcsZZhqmLOWG/eU7qvLtulYtj5mqL4nfCUjphD5YzNLsYeMGfsCl0HHirD1FpdjH3BrXhzEDlj1wJa3w55ocFyxt58RAf6+4Ye6K+PoGaBQCBQBXjzFM1QiIH+t0b+x/nczWDswXotz9Qaq5Nv4R36UQxMnpmUJ8+E3Sm5tSopnAmwDibPTAnLKpS1Kh3dzt7TS7YlvoTrkgp9QDdMVw0yk1kL40vf4ebShwY1GyVzQDnv2Yl75T69TGBQp9syo5f4wxJ22OerAupddOKRTojKT9Amw6FLZ4jdUK+fHydndJJxr8yr8eZeaU12whjXUGrWuplg3v0wutPTz2692D1as9sZf4Iz1U5X71SZbaa/Vg3ohmLaDUnq1Oh19UGtqq1X5qfW2dbjJs72fdMDYzIjyC6Xoivvur0M1jP7nKm71CoplLlX1mWYup2uMmn3yl1azmlh0BFWOSdHhXmnf6ezJKcW4irk3Z+1TgMZI3dlLuSIIAb6owL5H+cBViT3I/kBn67lFVqvOQmi+4KPprtCyyvg3Sm5tSopXAHSdkVKWFahrFXp6LiUbYnP4KREWAeMdAGCjboZq7Fc0954FqWk16wdmVyY/DZk0MzrqC0Drnji92t1j1496g8Jp2GtRL4bTnNKenV46jsI6mZQjTfqxs41Oh3qZFCbj1bGadBVGJmmbkZq2dZLRjrjvalll4ar6fX3qwEjoc9xLIMu0osI61zP9CrtPsLZNsI1t/UVEZaFdml0ksH0zfKwnik4UxtVrRTKqJu+zGg1GrrcHffUjQg3apNBzUNJehHSOIuic2xFn4EuGK2zDesxb38xfWgwXQe15GiIdXP4EAN9IBAIHBU4fDP6WIw9WGzI5W07j+xm/YJuUCI2TMWGqdgwBcSGqTeHu4dejH1veN28DdiQy09qLTkqtrZAeN2E10143YTXzVuKwzfQB3VzACRboloOVC8QCATeHGIxNhAIBKocMdAfGeRyeeFbvgO3tSQ+keYwI5z7xxeVKgi2SFpXuKalxMwUtVxsYR1MdnJrZQotOXixJSUsq1DWygqvwb0AfORmsOTg05NW+2DtBE77o0RDNjk4aedF+WLcBDtg9q1AKTn4PQBMfRBg0hUljn5QaQdhRcrI3x49ivP9ESdOy1Jdr6Yu1OTg9dCjAR3vcd6BBuOvRVWXE/aoBumKC9PJwbu1lZf3aKFefyQ73LXYubqc/vka7lFMra+kymzrcUfF1EGo0eysfnOvJQefC6tVbU86ObjYdqIWgKl6igMmB+8aOjl4l57rJGeqLM+s1FZlObXb4JK05GCSg08Bho5e2aa8/xRo08oMnRzcKrS5P7c5Db4OasnRkBz88CE4+iFxu+PLatLLWeKYLqtSA5UqHBrKlpKywv1XKCtYTVuKRP+syPbl09vTy/RnzytNRqs39EDaN9l7QGcN9g5z/iy20mgLfTWqSq7L1njLSOeswf3OjLJ+y7kFUqtZq0MtznvSt/JKio77zjnj7XprMg1xt6bfuT/WOlPzaec/s8FWgGvdk+BdD7O33htTzEzo/LJtf3oRu6wa6ZP6kBjeP92wL+OQ7i8WPepdHgeVl/e3uNYtkNrCe53La+jNM6HFusA1LDvkhWVGAv/PEeXo/2Nojv4DwdG/5fhyLg8shPO2JJIlE1iwRWcFJ/DvE0CjhSx4AWDJHBZoZTp4eg7AZjhFZedt4ekJSYF6nh5dklurJXNUYRPA0xM4bwtLJpSEvpxtVWr+OEv+AGCBxCURF5+rWTEVJBjO47qmhj74m+DXDHwa4K/VNV5mgtdore1wwzaAnhMY0QZw5wy2q6PIk/D7gP5mHgfgCsCtRgoaNL456TVPXNxzafIjmKQx/WX98D8B+B0Xc9xgw5yoWgLAZVDjFjDlWl+ES11DcUa6Hh7Vwe4NWKBHe3WYMwd26RlZSFwOVzjj5ztT/RosTi62LYB/BOCTauoV0KsT7R1Qp4F3ajX00GaY5xaWHwb0IbTFc1tkvgduVOGwzO/8bnivlpfo59F8yiGGLYB/cKbKK/ZKXaEtGzFfdvGLBP06FqOXg/azCLdpQHx5H1gwH5vRvwa/A8Cv4DWXw8AGelsBng3oVxHqTXRy+pAXGvoyF3JEcPiomwi+GAgEAlWOmNGXI/+NfO5PQJyIOxNhbgJ0qnP1IP0TwL4Qt2lZK7MtcUBep57IAJ2Mm5AUGMu40SW5tSopzAOMmwCd5CaUhL6cbVVqvkPLMiffkZSHTwWhbrdrfiB0hrYbdlDTDTA4MkW2yLRrrMyXOwFGqG+dsA3iQ2k7Y2UGJw+W7Nv0ARWAYe7rfqQTovIR2qQAtU6JNazN7LbFzeiHuf9HakOcn2Kt2+5rNUc4s4e5CkbC+GvZozP6YTDSGe9Ntase5lSZbSP0Ho6oZKps9B+Rvtii6xmbpI90tvnCgCvXZ37nvgeGuW4vg/VMjTO1qOctqKke9RmJBE6wO+45PRHWOZapLBprXitYt9Q5/fZxUOs605tU5/4c4TRkzS4cHXPeWIx9O5B/KA9wZcJ7Dke/TuGMWfCsho45nTWzwGiHhwDOeFepMu0MfxdAv3Mf5lmGz0oKNDL8xpLcWp3xLlXYDDB8FjzLGbNKQl/Otio1X5isUyVUiHAfj9IwF2TgbnZBFCyV4OWJG3nxohSJbwPWOOB5AEYkxMGp83hBP4pzGmJXBtwud3HZVIIW+eQ5J0Tlz1IKy7PTKanXce2Jg0gluEH1Z1MJdrtet5rPQrdLJWgV9p9KcAM854z3pu4nleAGeFbfoM86U/eCvco7tXOs/3e4nrE+fM7Z5lMJDqoxFVMJjnIXuEF5jGcph/VMnTP1DD1vxVg3ZMIQDR4oleBmZbSyi7H9WqFNNW92VJit6IzSkz6XNkBW0RrSh7zQkF1aOCII6iYQCAQCh4iY0SfIfy/PO4FSjoYJlDagNkhZd8YOmpBkct1AarfqBK1gfpC0a7kd3p2SW6uSwvaSAV5YVqGsVenoeC1PdVc0LRHuBca6/Ny2M3Z8QkMV0/NlWalL3BXkvCck2UBkfir0VI1OJ2VmLV/KFrHAo9exBMOdENdK5HtctIPh0OumewfcGSuTowa3M9bONSI9wzXXjhFuJdAqDNd5pb+WfTpFKqphXq1Iyi4NV9Prtyhgw2GP9nzO1THGKe96ple7vcHZ1uDUDjjh3oyTkq9ZTN8sD+sZ261a1DM2aKuynbH9meBosjPW7vherYYKi/rngLvqfqeq6Lx0ik6/7YwdqfSktfX296cPFdN1UEuOhqBmQd0cXuTvytOhX+Pq4bwJJp6dlHcBZ+u+/rOTz8+Ef+hwFQSvskk1bVLZxLOT8sSzYU9Kbq1KCs8uGeCFZRXKWpWO7tGyuB9LiqbXEmE30OcCHffq1Z5eYr69o6RsKE/ijch5z4F/Ty4wp6sMA7qEIR/gFswAFzdGUHCb1Ouc0M41QpvUwS51AdoOBf0qH5vRiXMo9A58OxxBbE12pGPd2C9/h3NItYFgu1IE9U5JUR+Wohomxo9wptrpCk6V2eajVxb16B7V4P1i+12Y4n5VW9BFFh9I1TvvW7ihnnSsSkGXs6rour0M1jN7nKnm9S+tykaQfGa0klvgQ0Db/zv1LOb7aG199MqieyUUXZ0dqqpLhWXhd7y78IDTQMbInUfHUBhBzQ4XXsrlgUboVX62qInTii7i2KRW1rWU6NEIahZBzSKoWQQ1e8vx70P70f9fb86w4OgDgUCgynFcz+hX5fIz/gcATdCu3ikn0DoH4BPQtDip2baIGYtLXje3LwK4GfATaq1MOx0tAA+pSqBpMR2LkgKNdNxYklsr+zgQ3r9jEU2LaVtUEvpytlWp+ULazgOYIcEmZbP843R8DmAnTFsJq1SDed1clEywv3RRkhFCvqxlw1Qj1MI0iXl5SZJhavnnWQrXAvBz3XEj03Bx2JD5adbrZq2WT3ZCXGRBafIMjHdK6nU6OeEgvG5kQ/85lbxutujeZsFjek0bnNfNXD26f6+bNjjXGT/Zmbofr5tlMBfuB+AaNfVc2KufOC/DNo33YF43K2G6++iRXVQnOdu8180L+tFW0evGx5hcpndhMuUQw+a6jWzL1OumUbmjMmeVbPTKA3rdrNfHeagQCO0avXIFtLtn3z5uzAVWPshsM52sKDWnD3mhwbxu3n9EZ/T/NvSM/nciTPEh4D9zeeCml/jhLIDR0K2UX8GlGbpA61/zEvfPStj7EbAcgC++BHDnLETV/brXtBtkR6rfkn6BC4e71nnP2ym64SZVKGTCk3ABXOOE16QrlLWywi746EsA/zELlIUYoY+4OCmeqBqEXcnBFh0OOpXVX6X10agAQuxbTNer4G8c7+yfJP/AZqmb7KGyB9xYIFs2HHT8b77SQJ9zR9EL7IRihroZdPcF9y0/oGRFwcWHsVbWFvfneGWHyozfP3VzohI1QI2augVqYSIAa2CEi3WTVZuvxIEUMxny0MsvZoRWs2IwHDtE+uiJ7tEt6xBBllrJuwsvqy/Cer3SsuAcNnw3qP170luO/W31zLtdmu2I9oe8sMxI4LYjOtB/f+iB/sNB3QQCgUBgPzgeZ/RfzeX/AIBvw0cBOBW26Cx1O/wEgJbWUmjJL8NtlKJXfgnQl+THVdVtqt+iV/a1UPfbRK/8tir00Su/7IRfTlcoa2WFUvRKwKJXaquNUOO4i2z0yv+d7q4PAPAa7IAPyRqrRq/8NyhCy4MA378i+RyW6ZuU5SNmXFrhdkeMrHRCVL5Mm9wPf6HRK2UBU4Ln3OOT8SpsUiaqngBgAQyqZJx2xbR09EqpeSFcAP9LhRdqYZze6K3uWiw21i/gWmf8MmeqXfV2p8psu9CdVwrXQpeLXtmlO/IsaNoUWO0WSH30SrHNz83r1MIDRq98Av4YqBS90iwUf61GeAI+C8BK9d0tC/q4xIUJEhwweuU3tMlQ1M0SZQVHpfUPFb3SspxbrB5/yAsNFr3yvCM6o//e0DP6/zuCmv1W6Mzlz9KB+OYWJipH36Qc/Rkn8JM5AB0tJQL95kXkFjNZOfq84+hNVU4rN3mOXsfxpsWJsGkxkxuTMo6jb2rn5hZV2Fwy4OZFJaEvZ1tZoWkhw84DmFyJox+A5gxHP8dx9Lk0Ry+E1UThcO8E4BKmfgvgms+zFJZfAfAqnEqivwaWAkqvm/eIYKwexXH0QqGKfLI2GQU/dErG6q/0lIxOHFEmquQSuqBGJdbkwTRHP0pbPajf9QVtgvNt90qMo290vjpLleYWiZ3OqzLbxrqjjXp0EM4B4GXoVudUY5kehemqdizcDSghPlYv1tQaR98FTRmO/peOSW903V6GRlUur5/50KgOY9aqjKM/UVeyDGsh7zZ+eI5ehA3aZKiBvkGfkxWuMrqik4dN6XtkFRrcn3aoIV0HtUQu5LxMJ7ydCD/6twD57+WBCYsYidvlNAWA3yT7jCA2TMWGqdgwVd4zsWHqWMdxNNAHAoHA0YyY0b9Z5HL5DuU9nmzhHZoW6qPvBzj1/SmOPmniOPrvRIapyDDlriUyTEWGqcOBwzfQHxdeN0mC70AgEDgucbzM6AOBQOAoR8S6OUT4vK9+e4VlP7CAR3n9EqxVDmeD01OW5KylFeArLYlC0mGwzM8vN8Sdqyj3IVm8I1WtE9ZmKvhWVpitX/Rya3UJuXTSbECSHMxMwtEzRXdd+lSfZTbnVXjbr2h9VyIUjwXZjCNbxmZoNY+ifj7j0s6J8jZtJU1ehrzu3sq5MESXZHbbktk1YwukPmaWrMi1q2ECSyLaptn7nnLxynPpnUFS7lQ2owuGO+Ptev2OpKJTZbY16M6j8W6tuKhh08X9xjpZuqUAze4hXOk60B5CezCe0SyDeegoeRiUHxVjhJOcQTksbaF0+7nQq3f2fH3Gypz+9qTjCJlV9jv090iEUzQc06DbClfjuutM9TXaCGdqZdwDbyeVPYsvaQXx13olfcgLy4wE/vyIulf+69Dulb8bO2OHwtdyeeubW1wa0r96ibZZADNaeaEFIqhZBDWLoGYR1Aw40kHNDt9AX53UTULKtxyoXiAQCBw1CK+bQCAQqHIER/9bYKl+/uwmyRmF7AXtBGifx4PwsUcAlsxjgSRBHQadSfZXxnD7pwE+LK0AWDKdBSs1LN5YfjwdlOhYYKrMVbCT9nkAz+tGRzGgfbpaMob2JmeYtloyTxWOB2ifTvNKlkwvCX0526rU/CGWfAFgweOAhklcwJoLALphzl2ua4QE2Qn3wa0AX2pIcn3I5V6ZXDQ74RLxQWtOdsffeTG71UnuF/AXgAYClG2T1+8DGEjvnqzZx89UcqUTQiK/el/S5F/hVLhWldTs4/sFgA/vZSBDNNjOWFF1dwHgOqhRyUAh2Tf0a3iP20F0VwHgffv4eaGUQuS9WmGgoP6p7lo6dbPPL2GhM/5qZ6pddY1TZba9dx9/UwD43L7E1IWwW+m1TmiAj+nentUFgNWwQNXW7OP+AmhMCLGtr5AUgP+3wMdUOAjD0jum/q3A76jk7kLCBF5dtqtK+/C9+/jrAsBf7OPuQrKCdRP8Csj47bVluH6h8sqCwflljzdc0AIcx2Ls0zb4L7sBHhjONlfZll469QchHp+/1go+IoIdKguTgFoiF/LuI0rd/H9DUze/H9SNR/7R/JR5SflHzl+bJxMP3pHzGEvimz1lHjwAlIcplh4dScmtesp0eLIUpvjl6aAhEEqqzAe7nZHzAPpLKZvgSUZOV0saGXljSW6tEst/key2HTkdnmTK9JLQl7OtSs0XqpuwLAXIaPoQwy8A2SR5FvxGNViY4itVmE4OLo9IToh7Wb4YlmieeDFLdbG0CI8AOuDKS+LxAmTDFBdKDPgjTogy448Xkia7oN4pqS8kzPIj9fsNU1wAXU9+RMIUp83YoToFr+pJd7gwxVahFKbYKbEQCK/CI854b2opTLFTZbY9Xkj83B8vlEzdCzIB2A6vwaM5gBpt3ul6pr6QDFiPONt2awEoqjFJmOL0u7bBXeCrGmDg8XQd3zM1ztQz9LyW68oj55JJoX11wDDFtXqooh/9Onh5OMBOWOfGLFuuaNCTyuYAM0AWmWvTh7zQcJQkBz98qJ6BPv9oHuAXvKoD/Vlu40/TBUksj24Zg64EeBVmXAXAMJhcmtHn5yQ1mzSI8Ksw4wJdKRubuDrIz2DGla6CYHJy3tq0Ad1myZi0YdqqpHB8yQAvLKtQ1qp0dCmvSjx62V+kM3oZ0vcCL2dm9M2lGT3paZptQ98K08RX49Rk59JGmOImVtLxfkZ/8RAz+l0qmeeEkMgv1nnrK3CiU2Iz+nkHMaPvKCT6K87oL3az19f1pF2FJJbZTldh/zP6PjmFGn/xwc3oOwpcvI9H9Lxmqs3oB+BkuDQ9ox/veqZmHzu1FZVm9C+qMUPN6O0CO3RGf3FmRt+hFv7KmSpP0YIhZvSvZaJTHHBGn0s3saQCw3RGPwnO2g2waTiTXOXsjP5SwM3oZW/j6PQhLzQMZi7kiCA4+kAgEKhyxEB/AKzJ5YUt6Z7H0yochA+6SANbhCJv5d9a6JgH8DScfQ5ImOLp7J4HLhl0swuB8DRcN70UAkEcdWUKcN68UgXBluk0twKclg5T3OxCIDRnQyBMTyw/b14Sq0AM8MKyCmWtrPCalmdeDDBwMUAdnALARlj7vkphis8vpRCxHD3o3EdCICTf7TtYfTHAS1CEjzwI8OoVPKXdnoeLAPhpAbJhigulMMVPOSHqtfnTQtJkE3zEKdleSGIm/7B+v2GKC6AO6adImGLVIHPzd6pOgSzT1BdSYYrrtEIpTLG7FnPJfw5OcsZ7U0thip0qs61OqZufFhJTT4Iu9U0swmZYkgPIKXMyzfXM9kISAuF5Z9uJWgDOUGOSEAjpj6rrXQ88p2GKf5quY31YV0huu5j6WT2vfP2URQ7Ylpksv/kwxY/CL4cDjIZHK4UpvlB3FTwIOP5d6JoF6UNeaDhKQiAcPhzzi7FfyeWB8zR4xauwQMfQr7SUPhjDjz786MOP3hB+9EenH/23h16M/XgsxgYCgUAVIIKaVUZEKwsEAoED4himbjbk8sDpAORWlnJvr1Ou/F/hi1o5p1mZiov4GVyzGGDdIib9J6CJR8QRZTytEwC+SCnI8LoWJrWWosa0tgD8LgCTTJVWppniIoCfgBIbpaDHuVZ4N8X5Jbm1WrdIFbaDRjYWRkiEvpxtVWp+GetOB5gkFNLZAKxk140AO+D0LfAD1WCJR6YlPsatl5TYDOCDgAYhmS0pqy5M4uE+90Huhg8D8GMQDybhWyTBurClRjQJep0D3ElOiMpHa5MHYLrySJIVWkyqO4jEI8L+z3WJR8yMTu/26pK2d+vqS42m8ZNWQkPtc0r26Z9Pw3xn/GhnqqfXcHKx7Xxlja9QU+fDdmUYVkOnOgdb4pHNMMnltt6tfWW2FZwX4wp9WguugsGb95Ry5dkuFcPOh2cBeCc8pUFjGpXZK8vX8bq7p4JRGhFIIIxT4rGqpk4GlD+xJCfGL62FGwB4Ada6TFhGwozR+MxyCfZ0rQW0vh3yQoPlu7/5iFI33xqauvnEcRjrRoZ4GVhvb4Hg6IOjD44+OPpjn6P/x6EH+k++OcOC+ggEAoEqx7E3o1+ay0uwAYkHIDxNc6vu//wFS3RmHSEQIgRChECIEAjHUAiEbw49o//U8UPdSBCbBUZbB3UT1HOzlwAAACAASURBVI1DUDdB3Rzr1M3/Gnqg/3S4VwYCgUAVIKJXsjSXl6n6khaCugnqJqiboG6qj7r5h6Fn9P/teJjRJzEprwShaJ4EklCOiZ/clUkUSa4s5XeP6JURvTKiV0b0ymMoeuXxO6NflcujzLtM1YOjD44+OPrg6CviWOfo7xh6Rv+ZcK8MBAKBwH5wVM/oV+XyMkvlStrmMeMRgLZ5zFgM0LYIdArW9EiJumlziUc+rqqadDdpRwuPwE02C/4fcjhF3bTOAfiEtCI514zFJerm9kWg1E1pQq2VaaejBeAhVSkGdCxSSxrpuLEkt1b2cSDMTMcimhYnlylCX862KjVfSJvEo5eNrJcB8DgdnwPYCdNWwirVYNTNRckE+0sX0UhSE7gGgEaohWmyn/YS+C7A8s+zFK4F4OdKlMk0fAOg89Ny6kZ3J6KkATqdFPlkbfIMjHdK6nU6OSEzRcVTNwAsA+AcoW7SZmyhFLwTeEyvaYOjbiy+Zom6cUqMummDc53xk52pJerGqTLb5urSyTVq6rmwVz9xXoZtScz/5KKAlTDdffQIr3KSs22324D6gn60JdRNuq/8/tVlehcmUw4xbK6y3vNhmVI3jbAFyDAe2Z2xB0PdNOuhitRNu+6MXQHt7tm3j5savUcyGljsVdk/3pw+5IUGo27ef0Rn9N8Yekb/x9XK0W/K5X8Av9cC0A8/gFvmAfwA/nJRUgBaHgBYcxW180o1DU0PJIW+RXQAcNoDvHgVF7YkGj7yZwAToUOflT36LDY9QN+ipLmc1MIUD6oQtdCsAjrgtAcArrmK0c4AMaZvERug2cmtlSlscgZ4YVmFslZWeAMkssOHPw/GOcyjuQegbgSrpnOCBoqQo69BPQlrXEwPo/K7XQ/roem/AORh/edBV1wv6gZYOTJZ4JIftsSwlaXusoWvbg2rALQ5ISp/Spt0w6d0//0J0K3vkmc15LKHsSWiapcaP6iSE/S6LlTDBPYVvxD+BnBBC6RV3ikUw/aSrFc/oYaJ8U85U+2qu50qb5vBNHdosFy5EKM7xOxLYavrmQ8PAqzIpy7WNE/UwkmVuI4LddIibf8IcCGjywzDjd279M7u0PnDunSTxfriN5yiwRsEG90VifAr8A4A+qDODd820L+os5F3wItaGX0GNrvYEhu1GtrQ/tzoNPg6qCVlF3JEcHzFo8//Rx446SbQJ0weuEatsMnXnpocanQ1y45Kk4lOYg/uBHdUfoTDjF2dmjrRJm2VOrtTZfKJepZtMNoZcMbUpDAubViplSqUh/IMZ4AJy8q+lRUG00KZzTXA9hEAHTCK5LWHvtVOhn4l5Uk/FsJWF6Fe6de97rzAGyMB9uivUQb6Tqe8LBxKjR71h6SmtRJ5AbY6JTX6kVGX0Ymb1tU4Ya+b6e/WnulM18la5ZXsNmceZ3C/q7nbGe9NLbs0X5NKBTFPer4Pim7yK0sdO53aGlibTy7QlPQ6bbVOWMh8AO1M94C/WVn4k6K3oFdbZeeZZZJOGHR9K7+yende38QWP6QwmK4wmK5sBuzUfqtNV/Btyw5VNLKKcTQO9IFAIHAc4jjyuvlyLi9XexPMfpgXLwNYAzc8zL2XAdzwMO2XATQ/DNB6GcC5Lt3MDQ8nqr5yWWnmtUjnPnfCx5TFm/0wz14G6nWzFoBezVYzAEqr0Pww7ZeV6NF/B+ALDwOJhbPVKtFwJ+DmmGLAYi2sh4ed3FrNVoXiVLMYFullirA5XaGslRU2wHUPA6y4DJSUPEW9bnIwCS5WDTJb3AXrlMf8RaWZ0UgYCe8B4DkNYTj1UVovTSoMU8pVvG4815j1urEHuuiEOM9om7t5X7q8TqWHvTmvm4H0BMfcMPo1puOog/C6EUxwXjdF5x5+MF43FnLSvG7y6g++FRozforeIT2ftiTrddPvhFmvG+9NOE5dy3OUw7xuvKkdWrlG9Xuc7kghtGYxbZhBhI2ut7MNxQxptUu7ziqQHhzLvG7Mfn/IC70lcooj63Xzt0Nz9H9aNSEQ/jqXxy3a+MfuLHiZJCX3y244QMsz4WVXU/DFv6f1j0pKfP1cWihDgDho1ir/8xqptv55avl7gC855bgKuUN6M/vrzbqg4caRsnNlMVuJSOscYIWrcHL66tBqUmeiUkM1ldY8/Xlr4AvP0ap850xnp7xd5O1b9vAOusXYKU6IW4yVJq/qawkd2oTTPycTdxfX7aJKfrf1jgEw78wNapg/dR7WwPUA/MYNTHn3erBr2aYj8h5HK6/V9cx8Otm0H/HNtloXk9lcUQeVox+AsxxTL0vEtXCae2kJpzzF2ebvl21cEkqqLJuj39bUr7d7MuXo1/NKt8+AfvXWnaW3eD+vEH/VxguZDy4qPN15NufctdgNPUNXj7ZoYkLTLArtpLJu/4pWsMmfP+SFBuu6vzgmB/qln8hf/8/A3P+56tHPNFc6Gu6VgUAgcFSgOPS/odF+x6XXc9/AwODAqg/9cNonllY8erTM6G/XV1nLkmQ62nkr34RPAdAF34WPAPBdaOkHaK0FaLkX4NUbSh/g33Vq5SjAdXTWAoy/l9YbErXfhD8A4JR+OmuTWcNu3e7Uci9cl7RuraWlnw21AKf1J6cWmKoW/SjtrGX8vQA9NzDCGcBPk8LGWiZ6ubb6piocL6p+Ctcl5xKhL2dbWWEz3AXA7wE6+a2HU7cD9Ixjq26JQqdCm2FQJzXfd8kf8vBRANbDOrhJW4n3pJBRLdsBFo9L4j3IZFam4ZKj3M4l2Kqpw3EJVGQdWOTLtcnd8BmnZCucC8BD6Sm5wObLokrIsQ/CoEpO1AucQimDvNW8DM6DvwOgqC4l0koezS3uWnL6cXAvfMAZv9yZale91aky2y5z55XCB2C9fmnthN1wqZ5L7vap0OH2sokX8Spn2wQt4Ca5E6AzM1LMgRfc5f85oMZ7mIVi1cnwMHwagA69BW3pJveqK6ThTHjFfUDYJwsq/O/apNaxZ73OG/VezWnTkNZvk3Tji2R51uKDyI/shvQhLzTM0At5xxGd0f/N0DP6zw1p2NJP5O963+C3rgXa75g3dcVtUk4fPSoWY3fm8vbptHkBTQtA8+/I/5aJRnCfs/r5GwA2qGPA8NIRCnoU2Khff9NvKKlFw9s01tKj0Uu69If9/A0l5xY5qZziN7WprD0PugqCHj3LT2ChM2DiDUnhDZjj5NYKVTiiFmCiM8CEZWXfygrmniShbuSXsweuGgewErY4qtS+lAs6DOX1B1YLu/UHIHzOPQD0u+YFuH9cIhQyQT7AhQPtShspGNS9x2hiIHSMfkZb9ajlT+id7YFBja+yvdJOv3yGQAfWOZK9RwmZnZU8DrvgWa2wz1Xo0XP1u2vpdgTFOme8Xe+gu2qzp8eVuzKFddCjbvIT4FU9ZJ0p/lo9qlO8Id9wtvW6+7JSw0D1Os7K8HB6lcK6vSK69CHp0rbAMD11d6b+lvSfXS4fFs41CyfcojVHqal16TRY4n7WmNYv9sjqmvzw+yoZsKXSobI6uytdyLGB9lUvzp0pr2qap83mh6vauba5/OgRpm7yX8pH3tdAIBDg0KibVSuWHcTRo2JGHwgEAoFDca+cNnPuQRw9kgN9MZcfaAEY+3UGbk2Ep9yZcPRzW3hZd6L3w9Pq0ve0OvnJp+s5DwA0XlVyJjP6dZ8eBc66qrQz9q6rErUvKylctjP2FdV81lVJ86fgPZR2xvo9hKbqPSqxnbGT3c7Ys66i7oGk4HfG2ik61F9oroYwrHuAs64qJbluIlXOtrKCbQEXdxjzCxzTAzBtBGPdvs3szth1+sDJl7UwoeuhT5nNvFKi7bAPrukG+M7IxG1DOHr5dpQNqD4TN9DtUm+X7YwV+VPa5FmYrztjR0K3qj2YnbGCyW5n7Ei9xdPSO2MFE+Ac5/Y6QQsj0ztjxbC9KnwRJjvjn3Km2lV3O1VWnpApTIYOdcXZ6Vwt83qzJICahWC7OLMzttFtGCyos2wjdGVGkDInSPG1ze6MNQt3OlOl8g6NWlq2oXR5aS9gggPujP2FNhlqZ+xyve9j0vqH2hlru/lk1WFq+pAXGo7hnbHN02Yvu2u1dFH7qheZ/b7mCkePDG0SjE0gEAi8FZg6c+4/f/kb7UD7T3647OPvu7bS0aBuAoFA4KjAIe2Mbf7Mo/etyE+r+az40V8L3v0mOXoE3CvNk9L2UAxASytfbYHYMJVWRWyYig1TsWHquNkw9bWh3Sv/8tjaGXt7Lt+iiUGEoAdqvk6rcvQtytFvaOFf1BO8X7NwAN+BvwTga1Lfole6moYWCxLpOPrWqxK1/6K++WUc/Y8BuO0B+pRA/xr8pePov+ZOYar+UiXG0e/00SuVo+9Lc/R9jqP/F1VoHH3fVcm5RPi1dIWyVlYoRa8EHEcv0SvfGEHnfjn6Ox0RDPxXQKNXlnH00lEtytGfqyd6qzj6/4RbD5WjlyDANx8ER2/hgs/R6JWDGpyZA3H098BNvz1Hfz9c484rhZugQ0fPPnhDM4Pl1ezfiqN/46A5+vvhC0Aljt4slNfPVLhft54MxdEvdlE9BAcTvXKRXnVFjn4xSHqbMWn9B+ToLeiIP+SFBuPoz6jSgf7to24SUr7lQPUCgUDguMTxFaY4EAgEjkMc89ErSxEOWmk1xkaPBkdPcPRAcPTB0R/fHP1XhqZuvnD0c/S35/Itdybl1puw8maNnHLKEloXJD+nL7bCbIB9N/FV+DwAu+Af4L8B8A+OrAda1gFsm1RKPPIP7tQtyiCumZQwd4V1tE5K1H4V/gyATpd4ZC98WyxZxxplFmWFwMKM+GUAU3WLSiZCYR3Aq5OY4gw4Y11SGANjndxafVUVJolH1rFmUnIuEX4nXaGslRU2w78C8KeASzwi2e8k8YiRufJjOwn6NTvdEk2fJPEJhJfvgi1wtnaRdObfA/AZAH6kbvsyxp2k5yLz2ss5FnWDE6Lyjdrkl/ARTR8oo8BYreBjXdiFDDhVSwBYCAPuHS+/lVNc3hWruQCa4P/oJSzQo0XtIj9CWaSBn8N1zviNztTsy7jobFvgziuF62Cndshu2K7rBHk1oOjycuT0TbPN2ZZ3L5UdehPzbvQ0eI/+JZpxcyPlMAslRsK5sAQ+pBYKRV4WOeD7ujJkkKyKtmKxTYWo8NvaRN5/doEW6+b7ejvGpPVbhkiLyWwxf8wYtL4d8kKDdUjzER3ovzz0QH/b0czR58zui51Uy/YjP+Uy0BtcbCG3CeBnlP6XQ7bT1xoCbZMAXqkUnNqOAlt0uto8qaQW+CUAY9yEt1tVtU1KBcTYoI/7eFL4masgWKFn6UkbsHdSUlgFk5zcWplCCdC21xlgwrKyb2UFi2MuqT7lAanVYboPVrkRQcasZ2CUizwji2AycMiY2AN9buZlwciARwB4Xdf3ZDjbqOci85D1u7xa/emCyPe4REuPOvv71cJ61wOGXHoUlpu4FopuZir2bE7HfpFntB026NDT5S6wVlXtc9diYXNysN4Zv8eZmo3AXutsa9dJaLvqXw89Otz0QjEzde2DWhdl/jVnlcWLtw+dPh1SCy4Em2G7m+PndA/RHsqRUwtfd6ZKVKUm907yqIEn05JhsMe9mOUsck9NKE0k3pRdy16XblemIAOuMi67rL3VpG3Z+vCTlQ6VGTlcLyQT4/dtxeGjbmLXUiAQCFQ5DuOMPr82L9OGmtFsPjUR1lAqWwrm1SNLFF5uZ1LBRyKUN7+FsvPTyTV6KBu9stbNiV5X4mWNU2uFdjhRo1fa9/UancjYSSdowWZVXpVZdaKexSYaYkCPFkakDbNWZpKo6nFHB53+wfS5/IVIwTpBJgjSb3369bADChreFp0lNUKDo86N/SRNK8ucPefOW6PXlXMdktdAxzlt65FzbEa/E/r/+9X+vJun5xxHUaYTN6P333b97ob2u17yaxuDmVbFtIXeNjNsQAu+Qs41Kbs0r8qT0f5DxK6rzxmT17n/oFOb06mxv+R+dy5vTF1mLWewEqGUXe/JLmgV088DeqN9kzIJaQPkquvTmuv0kNUspFXJfHxYpTPK2oZpqFghe6isTq6S2dWEwzLQ5/fkASaxLg8w4ChCXzaGYboTrhudyion/293/6OjjAw647Rcn64D9Lvv036Xliyf+UQ6EYra1jIpj3O/1fVaRyr7OCFFp8Qk4/RivQHjtJBLy8v0FFWVGLDeCdenK5S1ssL2tFD+HK5dNBY63RtILrYPtqtz/TYdu4X1FpJkAzRoK59suuhewBZ4Hf1NykpA2U+oz519txPi1gasyVatUwd9+g7rVP0eNmJ6mmIrFF3kW3m09mSTyMMmdTm3P9FWMvTsddcy6N4ZW53xO5ypdglmT50rb1Jj7ERbHU0kyjfrdcmfo6HLpQqQR6jH2VavBeBU1VwP+9Khm4Fh6Rz3QvLsoDI2pQt5tVbue1lw48G0ZmAE9DiqTbgpuVMmlCYFF4B6ML2uu9PZafplZtYDBRVaP2TtLztUZuSooaM0v504fIkMg7oJBAKBKsdb73Vze2bhuKaShxnpaXWNOlRsdd9uA+pv4LMQ9JKCWD8yna9A4D+fa1RVnysPuElBjb7t+/Q70fuEmRLSXhxe6KtlU1fXOL+x4W76UFOpwoDqzB4dquwLw9yUHJ0D7nPLiXVu5c0+z4e5vcFGGgzXKViN23g5oH016PbT1qTpFE/+ZJF1PazY6gTYXsnZ9GDWrGzmW/Hs/rzmvJF36cUrNqx4FquZH/p6K7aSr5PdQ5g6LLM6Wubt6rPf/LbwptZlSEhvJ/qTQakVT9CRuRfFDAWUSws9UWa/9D7XtphpVafMw+50h+cyrcq4Jt+xdqhib5uGN+ncctTiLaZubs/lP65P8Og+ttYB/CN8Ufv19rpSebl+3J4Nf6Mf/p+GE/sAflTHCs2S0wBPa/lpuBCAJwD4pJZ7XU27tk9q+WkdsC6Hf1RVK9RFrBG61exu9Rr+pFP1BFzolgHE/UZ+dabqQq3cDZcD8DxJdj0x4DwtvOF2YD/tWq1QheL7cR48rZcpwifSFcpaWWGfumPLzgN7gYnr9NY0Wy2/29dgvH7evuTeEN3qXvkq1OvgUtRIyKuhCL8PwHe1K+SHZ44iVPK6GebKvmB+FLVqyfn6UhSXlRq96gN63cgixOxKXjfD0r9z6avpUKebMPbptTCE182AdsWLMMsZf5BeN6tgOqzU84qps6BHl0lEldkgnXkarHPrQ8OdVVmvmxH6Aq7odeOxSt2xsl43q9RCix+wSiMfFPRelHndrFaPdcMBvW6e16DBZV439jZarVs0umG1izBsT4vVtLZmDFrfDnmhYXjmQqoMb9lA/5RO5B/TEXBsXfJrHIQXdUz35U5tu8KNO6/A63WgA7f8LxNJ8332BOtzeqiQOdSvR0W+V+vn3azf8r3tzXgEP5fWtkuflf701LXXVTBVct52N9/ZpcJdzuyyVqZwMGNAr+siX/atrLAnLRTutUFHlkEXqgW9lrz61wM5XXbug3pdGJD/J2iFLtf8pyqUxd6i+5bKu/8Nfsq0/6l9HtpdedAtgZZ92OF+7X7O2On4dF/BmyQ1t6aXam29wab5/ooG3Vm2VDJ+/xumvH77hJXe8x+InVrZ8hfudgb4V6O9yezJ3KbPXi3sy8yyPf2dq+RBb4fEwl5nqrh1NqTXjT1eT/9ZqwvCOFOlH+rSTfqhy90s/0UrFu5J67d3Xs59g5K5Ha9XOlRmZN1Bf8MdowiOPhAIBKocbw1Hf3sun6WqDZ+Cf9TyJ13Z+Icu+BR8C0g7ftmLegzshNEA7HTvfLQ8VuetY5wzz59qMELTRprF89ReTj8kzSfHT/DL6NfPAfC3aYV+Prh/ErlihZyz0B/NO2E+U8G3soL1hnUOLu0U7rveY6zWGabfLg0Hypicgw/AnVo+3Vkis9EJKilrZR9zJzohKh+vTV6HgtYRtbJJp9l9ylS0Cu2B0U5iHdjhNk/iQj9uUVeWbucXVLbBtahNRutZRjvjxztT97MzVvxu7TE2Uxt0Z1AfnOh+IObucrJ73qSvTkzbZuda725Hl8sXL9gAp2l5h97usp2AuN+UfM+dDjv0Y26ifpSU3dzdlTYt+07AtRJhvU7V+/X7A2WcpNownbx3V1q6wG2tkqiu9ulv3I4/5IXePDnXn1UpR/9mB3pZel2oJNrX4eM6OkztTjKztS6g5etJ/dZbS+U3NDTxWBcAp2VJEqil81a+CZ8CoAu+qyGFvwst/QCttQAt9wK8ekNCO0hNgxwFuI7OWoDx99J6Q6L2mxpz9ZR+OmuTsW83/MjaXqdm19LSz4ZagNP6k1MLTFWLfjZ31jL+XoCeGxjhDEg4juvYWMtEL9dW31SF40XVT+G65Fwi9OVsKytshrsADaEsv5x6OHU7QM84troRVj7hNzti5PtuwSoPHwVgPawDDVqRDLgS8bVlO8DiccmChLAWsvNXViDsXIKtmsERtyFARg2RL9cmd8NnnJKtJJGQH0oHqxHYy1hUPQzAB2FQJSfqBU5xSyNW8zI4D/4OgCJcpkdP1NfDFnctxnrfCx9wxi93ptpVb3WqzLbL3Hml8AFYr8TITtgNl+q55G6fCh0uDsc8QGl0299vW6mNnJkAnZnheA684C7/zwE13sMsFKtOhofh0wB06C1oSze5Vxd1DGfCKy60jgXzQYX/XZsIyWM+FzZ83wsSzLshrd/C15yuryKbKJgxaH075IWGGXoh76jSgT6om0AgEKhyvKnF2PW5vHwQtema9Qh4Ak4G4NcjOX8BwDtgi07efdk2Z05s4R36ku9dQNcCgAkvcs5sxr8IMH4J59zKlK8DnHNr8nY+52SAFTcAzPx6EvdIagryehSY+RrjtyX1TyZRe85sTpGZQIHx2xj/cwBGU/i9pOZMjfF4zsnQyWkab/CdJ4POH02V8RHjtyXnfQlmOQMSU19jYmPKMGt1zmxVeELJALlMEfpytpUVxj/Amj8BOE1mKfcA8H7apgL0wHkPJZHjIJkuTu6C/6D/1qTfhO3ZDvs0FuAYOAXqOgB2N9G8EmD6dIrw8DiAtTrfl7mkBBG6vg+gN71jak4fP1PJfCeERH51X9JkAjwG16qSOX18pw7glj30ZnZM2fqnqNpWBzAZ8irprUs+bpbDe9zSm9S8oo+f1yXuUp1whVborUtu9Cx3LZv1IZ8FZzrjr3am2lXPcarMtiv6eFLPKwacCWfop1gfjIR3F5MOfbUOoA0WqNo5fdxfBzqvF9t66pIC8O06Pq7CqVCXXmr8bh0fUcm2uiRa0dWZ5Ujrma/VAXy4j211yXfpBzTGUdl9GJuh+5ZDvVZGPxblTonwTBdPsF75olrYox+jZ8KlvQAPN5Qqo58g9fBT/YK8AoDHtYJM+bvTh7zQ8EilLXjVhEOnbh7J5R/ULyDjISqyzy0v0TqrQrmszldnQTpMsYWUmw/PwjsBeNZ5UqKM4WUawGu+JiQC/nQtrZOTcq1z2DhDndhGOte9AfWPrNNUO73uNSgpbMxR/QtrAf52MrjsReYHMlCpE2qda4SvUOtaWfIgWfAQyrLBCRvSFcpaWeEKjdlkvYd+hgvOzATMGoB52odjlZfodZ4Psve9oF20W43/JNyhqiQg1KB6yzB0iChjbJoryU0o8cKmZCrMGiICl+92YdvfSAvlFqxLn1f4+i3QrqEZ16S3WGed97uUHxjuarYfdEiscbDdPQ9iqmwKXa113qmrTTklnXenCSvpijOdbb4HzJgc7Mp4o3pTx+nTkjXeDLP7YpXPdNGAPQ64QJXF2fCbA1WQ273qQJVlWcI2rovb6G/Sh7wwi8O3N/XI4lAGeuHlr4R599B6I8CN+kb9e/gjHaybXk9GmtYzadHgiq1Xl8q/uTopnK16gJZXkgWpFVfzI7gZgE74FbwLgF/BHwLwT1J/I8BDE5MfntQ0tKjj2OsTkx/kzI20TkzU/gjeD8AE2K78aa/OCFo28rqG0/0n+EPdkn6Knlpgqv5QJdth5kaAxyZyiTPgpI1JYQNc4OTW6keqUH5gJ23k9YnJuUT4T+kKZa2ssEbfVRZhGBitvsM7YYuLW1Cn/darUYIXu31kwMcAeAl26VaA3apTUs3J2smjuilBBnpZq1yvp/bY6Qas7U6IixUuTZbDVepaNxp26qD/vAt0bLAhRlQJBz0fii7yuzyZZUEspOYcOB3uU+EcLYyutLpbrz3wBMxzxq9xptpV73SqzLY5zo9eDJgHW9Xtbzfsg1l6XTLMTYbXnUuCDPFbnW1jXBgDcwwfA29k1kInuTiOL4D8+NZQDusZ68AXNKvfTv3sK4sl8IyuoxhOgtdBA1wldL+MpiK8R5+rnItv0+c4+qc1AnaN24+ChsJ+3eVXsFkX2hCtb4e80GARIxZU6UAfHH0gEAhUOd7u5OCBQCAQeJsRM/pAIBCocsRAHwgEAlWOGOgDgUCgyhEDfSAQCFQ5YqAPBAKBKkcM9IFAIFDliIE+EAgEqhwx0AcCgUCVIwb6QCAQqHLEQB8IBAJVjhjoA4FAoMoRA30gEAhUOWKgDwQCgSpHDPSBQCBQ5YiBPhAIBKocMdAHAoFAlSMG+kAgEKhyxEAfCAQCVY4Y6AOBQKDKEQN9IBAIVDlioA8EAoEqRwz0gUAgUOWIgT4QCASqHDHQBwKBQJUjBvpAIBCocsRAHwgEAlWOGOgDgUCgyhEDfSAQCFQ5YqAPBAKBKkcM9IFAIFDliIE+EAgEqhwx0AcCgUCVIwb6QCAQqHLEQB8IBAJVjhjoA4FAoMoRA30gEAhUOWKgDwQCgSpHDPSBQCBQ5YiBPhAIBKocMdAHAoFAlSMG+kAgEKhyxEAfCAQCVY4Y6AOBQKDKEQN9IBAIVDlioA8EAoEqRwz0gUAgUOWIgT4QCASqHDHQBwKBQJUj3oPXhwAAAFpJREFUBvpAIBCocsRAHwgEAlWOGOgDgUCgyhEDfSAQCFQ5YqAPBAKBKkcM9IFAIFDliIE+EAgEqhwx0AcCgUCVIwb6QCAQqHLEQB8IBAJVjhjoA4FAoMrx/wPsDlGq+9T1zQAAAABJRU5ErkJggg==" alt="plot of chunk unnamed-chunk-11"/> </p>
<p>The squared correlation coefficient relates to the linkage disequilibrium between two SNPs. A strong correlation, value close to 1, indicates two SNPs are seen together often, meaning they are inherited together. The blocks show regions of the chromosome that are inherited together. In this example, region3, there are three large blocks of SNPs that are inherited together. The further away two SNPs are, the more likely they will be separated during recombination leading to lower correlation between the two SNPs. </p>
</body>
</html>