# Gwas Result Result #27 — Analysis Report

Generated 2026-09-21 02:39:32 · 1 run(s) · CAT-G Data Analyzer

## 1. gwas result (#27)

2026-09-21 02:18:38 · gwas_result · status: completed

*GWAS tested 10 SNPs after QC (of 10 input); 0 passed significance.*

| Metric | Value |
|---|---|
| n_snps_input | 10 |
| n_snps_after_qc | 10 |
| maf_filtered | 0 |
| missingness_filtered | 0 |
| hwe_filtered | 0 |
| hwe_not_testable | 0 |
| ld_pruned | 0 |

![ld heatmap](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAnkAAAIpCAYAAADehy5WAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjksIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvJkbTWQAAAAlwSFlzAAAQ6wAAEOsBUJTofAAAXFtJREFUeJzt3Xl8VPW9//H3CUsSEmImETURMFiXUNQCt6YazRggoOCSsIq3FAahBS1pNWpQ/BVwA0ohtsVL0VZJpFB7WdS0ilAgMSAFbq0tXGlEFIFAXJoVRbIw5/cHNyNDFjJLZsvr+XicR5vvnOV9xiAfv9/v+R7DNE1TAAAACClh/g4AAAAA76PIAwAACEEUeQAAACGIIg8AACAEUeQBAACEIIo8AACAEESRBwAAEIIo8gAAAEIQRR4AAEAIosgDAAAIQRR5QAs++eQTGYah+fPnu32O+fPnyzAMffLJJ17LFWhausdQuO+W7iE/P1+GYai4uNjRVlxcLMMwlJ+f32abJ5KSkpSenu6VcwHoXCjy0Kk0/QX89NNP+zsKAAAdqqu/AwCB6LLLLtPXX3+trl35I+Kq//f//p8effRRhYeH+zuK29p7D1arVV9//bW6devWYVk++OADGYbRYecHELr4GwxogWEYioiI8HeMoNS1a9egL47bew9hYWEd8ntSX18vu92uiIiIoC6WAfgXw7VAC1qak3d228aNG3XDDTcoMjJSvXr10owZM/TVV1+d97ymaerxxx+XYRjKzs6W3W6X3W7XggULlJ6eroSEBHXv3l2XXnqppkyZoiNHjrR4jmeffVZXXnmlwsPD9a1vfUsLFy7U1q1bW5wLVl9fr8WLF+u6665TZGSkYmJilJGRoZKSknZ/HydOnNBPf/pTJSQkKDIyUoMHD9batWtb3Lel+WxVVVV65JFHdOWVVyoyMlIWi0XXXnutHnjggWbHFxUVaeTIkbJYLAoPD1f//v3185//XKdPn26271/+8henfw733nuv/v3vf8swDNlsNsd+bc2T82Re4fnm3/3mN79R//79FRERoaSkJM2fP1+NjY1O+9hsNhmGoYqKCv3oRz9yfMe7du2S1PKcvHPvr608TXMJt27dqgULFujyyy9XRESEvvOd72jjxo2SpP379+uOO+7QBRdcoNjYWNlsNn355Zdt3juAwBfc/7kN+MHGjRv13HPPacaMGbLZbNq6dateeOEFGYahFStWtHpcfX297r33Xq1Zs0a/+MUv9PDDD0uSTp06pZ///OcaM2aMbr/9dl1wwQXau3evXnrpJW3dulV79+5VXFyc4zy5ublasmSJUlJSdN9996murk4rV67Uq6++2uyajY2NGjVqlN5++23dc889mjlzpk6ePKnf//73Gjp0qF577TXdcccdbd5vY2OjRo4cqXfeeUejR4/WsGHDdOTIEd1777266qqr2vWdTZgwQUVFRfrRj36kgQMHqr6+Xh999JG2bNnitN9LL72k6dOna9CgQXr00UcVGxurd955R4899pjee+89vfLKK45933jjDWVmZqpXr16aPXu24uLitGHDBt12223tytTRnnvuOZWVlWnmzJmKi4vT66+/rieeeEIfffSRVq1a1Wz/jIwMxcfH69FHH5Xdbtcll1zi1TyPPfaY6urqdN9996lLly761a9+pczMTK1bt07Tpk3ThAkTdOedd+qvf/2rCgoKFB4erueff96rGQD4mAl0IkVFRaYk86mnnmpzv0OHDpmSzHnz5jVri4yMND/66COn/W+99VazW7du5pdffulomzdvninJPHTokFlVVWUOGTLEDA8PN//4xz86HWu3282vvvqqWYa//OUvpiRz8eLFjrYPPvjANAzDvOmmm8z6+npHe3V1tdmnTx9Tkrly5UpH+y9/+UtTkrlhwwanc9fX15uDBg0y+/Xr1+b3YJqm+eKLL5qSzJ/+9KdO7Tt37jQNw3DcY0v33ZRNkjlz5sw2r1NeXm5GRESYWVlZpt1ud/psyZIlpiSzuLjYNE3TPH36tJmUlGRGR0ebR44ccezX2Nhojho1ypRkTpkyxdHe9M/97O+mtbytta1cudKUZBYVFbV53qa2Hj16mJ988omj/fTp02ZWVlazc0yZMsWUZE6cOLHZfZumaV522WXmLbfc4tR27v21lacp93e+8x3z1KlTjvb33nvPlGQahtHsdzIzM9Ps1q2beeLEiWbXABA8GK4FXDR69GhdfvnlTm3Dhw9XQ0ODDh061Gz/I0eO6Oabb9Y///lPbdmyRRMmTHD63DAM9ejRQ5Jkt9tVXV2tf//73xo4cKAuuOAC7d6927Hva6+9JtM09eCDDzpN9r/gggt03333Nbv2qlWrlJSUpLS0NP373/92bDU1Nbrrrrt06NAhHThwoM37Xb9+vSRpzpw5Tu033nijhg0b1uaxkhQZGamIiAjt3r1bH3/8cav7rVu3TqdOndL06dNVUVHhlLept3HTpk2SpHfffVeffPKJJk+erD59+jjO0aVLl2Y5/WXSpEm67LLLHD+HhYXpsccek/TNd3q22bNnd+gDFj/+8Y+d5vcNHDhQMTExSkhIaPY7ecstt6ihoSGol8EBwHAt4LJzCzxJio+PlyRVVFQ0++zOO+/UqVOn9Pe//10DBgxo8ZyvvfaaFi9erHfffVf19fVOn1VWVjr+f1ORlJyc3Owc/fv3b9b2r3/9SydPnlSvXr1avZ/PPvuszWHXjz76SBdeeKEuuuiiZp8NGDCg2ZDrubp3765f//rXys7O1re+9S1dddVVSktL06hRo5SZmakuXbo4skpqc/j4s88+c2SSpG9/+9stZgoELWVrajt48GCzz9o79O2uln5vLRaLU5F8drvU8u8zgOBBkQe4qKkoaYlpms3avv/97+s3v/mNnnjiCa1evbrZchuvv/66Ro8ere9+97vKy8tT3759FRkZKUmaOHGi7Ha721ntdruuvvpqPffcc63uc80117h9/vb64Q9/qLvuuksbN25USUmJtmzZohdffFEpKSl6++23FRER4bjP3/3ud049YGdLTEx06/pt9ZCd+yCEvzT15nqirXtp7ffW1d9nAMGDIg/oYLm5uUpOTtYDDzygr7/+WmvXrnVadqOgoEARERF6++23nf6i/+qrr1RVVeV0rqbemNLS0mY9Vk09YWe76qqrdPToUaWnp7u9rMm3vvUtffDBB/r888+b9ea9//777T7PxRdfLJvNJpvNJtM0NWfOHC1atEivvPKKbDaboyfLYrEoIyPjvJmkM0+FnqulTE0PrpzdK9qkrSFkT7SUrantiiuu8OjccXFxPr0XAMGJOXmAD/zkJz/Rb3/7W7355pu64447dPLkScdnXbp0kWEYzXrsnnrqqWZtmZmZMgxDzz77rBoaGhztNTU1+s1vftPsupMnT1ZVVZWeeeaZFnM1DX+2ZcyYMZKkBQsWOLX/9a9/1datW897/MmTJ53uVzrTszZ48GBJ3wwJTpgwQREREZo/f36Ly3d8/fXXOnHihCRp8ODBuuyyy/Tyyy/r6NGjjn2alqM5V79+/dStW7dmQ8sffvhhi08le8Pvf/97HT582CnbwoULJX3znbrr6quv1l//+len7/XUqVNatmyZR+cFEFroyUOn9Pbbb7f62eOPP94h15w2bZp69OihyZMn69Zbb9Ubb7yhmJgYjR8/XuvWrdMtt9zi6OXatGmT9u/frwsvvNDpHFdffbUeeOABPfvss7r55pt19913q76+XitXrlRCQoKOHj3qNDT505/+VFu3btX8+fNVUlKiESNGKC4uTkePHtXOnTv18ccfn7f3Z8qUKXrxxRf1q1/9SkePHnUsofJf//VfGjRokP7+97+3efyBAwdktVqVlZWlAQMGqFevXvr444+1YsUK9ezZ01HwXHrppXr++ed177336uqrr9aUKVN0+eWXq7KyUqWlpdqwYYNee+01paenq0uXLvr1r3+t0aNHKyUlRTNnzpTFYtGGDRtaLBCjo6N177336vnnn9fdd9+toUOH6siRI1qxYoWuu+467dmzp73/GNutf//++t73vqf77rtPcXFxeu2117Rt2zZNnDhRQ4YM8ejcP/nJT3TPPfcoPT1dkydP1pdffqmXX35ZF1xwgZfSAwgFFHnolLZs2dLqAwOPPvpoh133nnvuUWRkpO6++24NGzZMmzZt0oQJE/Tll1/q2WefVW5urnr27Knhw4dr+/btuvnmm5udY+nSpUpMTNSKFSv02GOP6dJLL9UPf/hD9e/fX6NHj3bM55POvLnhT3/6k1544QUVFBTo6aefVmNjoy655BINHjxYixYtOm/mrl276q233tLjjz+utWvX6o033lD//v310ksv6f333z9vkdenTx9Nnz5dxcXF+vOf/6yTJ0/qkksu0Z133qnZs2erX79+jn0nT56sq6++WkuWLNFLL72kyspKxcXF6fLLL9fDDz+s6667zrHvXXfdpTfffFNz587VwoUL1bNnT91xxx36xS9+0eKDJkuXLpVhGFq3bp0KCwt1zTXXqKCgQH/72986pMibNWuWTp48qV//+tc6dOiQLr74Yv3sZz/Tz372M4/PPXHiRH366af69a9/rZycHPXt21czZ87U4MGD2/XEM4DOwTCZWQuEhF/84hfKzc3Vrl279L3vfc/fcfzKMAxNmTKl1TdRAEBnwJw8IMicO79NOjMnb9myZerVq5cGDRrkh1QAgEDDcC0QZNasWaMVK1bozjvvVGJioo4cOaKVK1fq2LFjeumll9S9e3d/RwQABACKPCDIDBw4UBdddJFWrFihiooKRUZGatCgQfrNb36jO++809/xAAABgjl5AAAAIYg5eQAAACGIIg8AACAEUeQBAACEoE754MWpU6e0b98+9erVy+33eQIAEIoaGxv1xRdf6Nprr3V6zzaCT6escPbt26eUlBR/xwAAIGDt2bNH119/vb9jwAOdsshreuXRdPVRzyD7CnJKO+a9qh0q8iJ/J3BPXZW/E7gn8sLz7xOITLu/E7jHCMJZL8H6XZ846u8E7onu7e8ELikv/7dSbri3xdcDIrgEV4XjJU1DtD3VVReom5/TuKb3pUH4F3iPi/2dwD2nuvg7gXuC9fsO1sKDIs93apu/7SUo9AzOP5NMZwp+QfhvJwAAAJwPRR4AAEAIosgDAAAIQRR5AAAAIYgiDwAAIARR5AEAAIQgijwAAIAQRJEHAAAQgijyAAAAQhBFHgAAQAiiyAMAAAhBFHkAAAAhiCIPAAAgBFHkAQAAhCCKPAAAgBBEkQcAABCCKPIAAABCkMtFnt1u14wZM5Samiqr1aqJEycqPz9fEREROnbsmCRpx44dstlskiSbzaZBgwbJarVqyJAhjn3OVVdXp8mTJystLU1ZWVmqqamRJFVXVyszM1NpaWmy2Wyqr6+XJC1btkxXXXWVrrjiCnfuGwAAIKS5XORt2rRJDQ0N2rlzp0pKSrR8+XJJUnJyshYuXNjiMcuWLVNJSYmGDx+uBQsWtLhPQUGBkpKStH37do0aNUp5eXmSpMWLFyszM1Pbt29XYmKiVq9eLUmaMGGC3n//fVfjAwAAdAouF3lRUVE6cOCA9u/fL9M0FRcXJ0kaO3asdu3a1WpPnSQNHjxYhw4davGzoqIijR49WpKUlZWloqKiNtsvvvhidevWrV2Za2trVVZW5tjKy8vbd7MAAABByuUiz2q1atq0aZo1a5b69eunpUuXSpIMw1Bubm6rPXWStHnzZg0YMKDFzyoqKmSxWCRJFotFFRUVkqTKykrFxsY2a3dFXl6e+vTp49hSUlJcPgcAAEAwcevBi6lTp2rbtm3at2+f1qxZ45gnN378eO3evVtlZWVO+2dnZ8tqtaqiokJz5sxp8Zzx8fGqqqqSdGYeXnx8vCQpLi5O1dXVzdpdkZOTo6NHjzq2PXv2uHwOAACAYNLV1QOOHz+u6OhoxcTEKDo6WtHR0TJNU9KZ3rzZs2frmWee0cCBAx3HLFu2TDfffHOb501PT1dhYaEGDRqkwsJCpaenO7VPmTLFqd0VMTExiomJcfk4AACAYOVyT15ZWZlGjhyptLQ0paamKiMjQ+Hh4Y7Px40bJ7vd7nIQm82mgwcPymq1qrCwUDk5OZKk3NxcrV+/XlarVYcPH9akSZMkSa+88ooyMjJ0/PhxZWRkaPv27S5fEwAAIFQZZlM3XCdSVlamPn366EH10wVq38MbgWLeicX+juC6Hgn+TuCeU//2dwL39LjY3wncY7r+H4cBwQjC5UaD9buubfnBvYDXM8nfCVxSVvaZ+lx2l44eParevXv7Ow484PJwracqKys1ZswYp7ahQ4dq7ty5vo4CAAAQsnxe5MXFxam4uNjXlwUAAOhUgnCcAQAAAOdDkQcAABCCKPIAAABCEEUeAABACKLIAwAACEEUeQAAACGIIg8AACAEUeQBAACEIIo8AACAEESRBwAAEIIo8gAAQKd28OBBzZw5UwMHDlTXrl11zTXXtOs40zS1aNEi9e3bV5GRkbrxxhu1a9euDk7bfhR5AACgU3v//ff1xhtv6IorrtC3v/3tdh/385//XPPmzdODDz6oP//5z0pISNCIESP08ccfd2Da9qPIAwAAndqdd96po0ePat26dRo8eHC7jjl16pQWLlyohx56SA8++KCGDRumV155RXFxcVqyZEkHJ24fijwAANCphYW5Xg7t3LlTtbW1mjBhgqOte/fuGjNmjN58801vxnMbRR4AAICLSktLJUnJyclO7f3799eRI0f09ddf+yOWk67+DuBPOaWPq/elF/o7hkue6Jnr7wgum1f3vL8juKfHxf5O4J7T9f5O4J6GWn8ncE9EcP07RJJkBOl/3/fs6+8EbjEr/unvCC4xq6r8ct3q6mp9+eWXHp/Hbrc365mLiYlRTEyMx+c+W1VVlcLDwxUREeHUbrFYZJqmqqqqFBkZ6dVruqpTF3kAAMD/qqurlWiJ19eye3yuyMjIZr1o8+bN0/z58z0+d7ChyAMAAH715Zdf6mvZNUN9Fa0u7p9Hp/X810e0Z88eJSQkONq93Ysnnemxq6ur06lTp5x686qqqmQYhiwWi9ev6SqKPAAAEBBi1EUXqJvbxzcN0iYkJKh3797eCdWKprl4H3zwgb7zne842ktLSx3r5vlbkE7MAAAA8J/U1FTFxMRo7dq1jraGhgZt2LBBo0aN8mOyb9CTBwAAAkKYPOt9cvfYkydPOpY9OXz4sGpra7Vu3TpJ0i233KJevXpp2LBhOnz4sA4ePChJioiI0GOPPab58+erV69euvbaa7V8+XJVVFTo4Ycf9uAuvIciDwAABAR/FXmff/65xo8f79TW9HNRUZHS09N1+vRpNTY2Ou0ze/ZsmaapJUuW6IsvvtDAgQO1adMmXX755W4m8S6KPAAA0KklJSXJNM029ykuLm7WZhiGHnvsMT322GMdlMwzFHkAACAgGPKsJ8/wVpAQQZEHAAACgiHPCjWKPGc8XQsAABCC6MkDAAABwV8PXoQqijwAABAQKPK8iyIPAAAEBIo87+L7AAAACEH05AEAgIDAEireRZEHAAACAsO13sX3AQAAEILoyQMAAAGBnjzvosgDAAABgSLPu/g+AAAAQpDLRZ7dbteMGTOUmpoqq9WqiRMnKj8/XxERETp27JgkaceOHbLZbJIkm82mQYMGyWq1asiQIY59zlVXV6fJkycrLS1NWVlZqqmpkSRVV1crMzNTaWlpstlsqq+vlyT96Ec/0g033KAbbrhBixYtcufeAQBAAGl6utbdjadrnblc5G3atEkNDQ3auXOnSkpKtHz5cklScnKyFi5c2OIxy5YtU0lJiYYPH64FCxa0uE9BQYGSkpK0fft2jRo1Snl5eZKkxYsXKzMzU9u3b1diYqJWr14tSXr44Ye1a9cu7dy5U6+//ro++uijVjPX1taqrKzMsZWXl7t62wAAoIN5UuB5OtQbilz+PqKionTgwAHt379fpmkqLi5OkjR27Fjt2rWr1Z46SRo8eLAOHTrU4mdFRUUaPXq0JCkrK0tFRUVttl911VVnbiAsTF27dlWXLl1avW5eXp769Onj2FJSUly8awAA0NEo8rzL5e/DarVq2rRpmjVrlvr166elS5dKkgzDUG5ubqs9dZK0efNmDRgwoMXPKioqZLFYJEkWi0UVFRWSpMrKSsXGxjZrb/L73/9effr0UVJSUqvXzcnJ0dGjRx3bnj172nu7AAAAQcmtonfq1Knatm2b9u3bpzVr1jjmyY0fP167d+9WWVmZ0/7Z2dmyWq2qqKjQnDlzWjxnfHy8qqqqJJ2ZhxcfHy9JiouLU3V1dbN2Sdq4caNefvllvfTSS23mjYmJUe/evR1bQkKCO7cNAAA6ED153uXy93H8+HHV1tZKkqKjoxUdHS3TNCWd6c2bPXt2swchmubkFRQUOHrrzpWenq7CwkJJUmFhodLT09tsLykp0dNPP621a9cqIiLC1dsAAAABxvDChm+4XOSVlZVp5MiRSktLU2pqqjIyMhQeHu74fNy4cbLb7S4HsdlsOnjwoKxWqwoLC5WTkyNJys3N1fr162W1WnX48GFNmjRJkjR58mSdOHFCmZmZSk9PZwgWAADgLC4vhpySkqJ33nmn1c8Nw9DevXsdP+fn57frvOHh4Vq1alWzdovF4ujJO9snn3zSrvMCAIDg0LSEiifH4xs+f+NFZWWlxowZ49Q2dOhQzZ0719dRAABAAOGNF97l8yIvLi5OxcXFvr4sAABAp8K7awEAQECgJ8+7KPIAAEBAoMjzLoo8AAAQECjyvIvvAwAAIATRkwcAAAICS6h4F0UeAAAICAzXehffBwAAQAiiJw8AAAQEevK8iyIPAAAEBIo87+L7AAAACEH05AEAgIBAT553UeQBAICAwTIo3kORBwAAAgI9ed7F9wEAABCC6MkDAAABgZ487+rcRV7kRVKPi/2dwiXz6p73dwSXPRE+w98R3DKvboW/I7inrtLfCdxiHtji7whuMQZO9ncE13Xp7u8E7jHt/k7gFqPXYH9HcIlR95nfrk2R5118HwAAACGoc/fkAQCAgGFIMjx4vNYwvRYlJFDkAQCAgGAYpsI8qNQMmRKFngPDtQAAACGInjwAABAQDMPD4VqJnryzUOQBAICAYMizN17wtgxnFHkAACAgnOnJ82ROHs7GnDwAAIAQRE8eAAAICF6ZkwcHijwAABAQKPK8i+FaAACAEERPHgAACAhhHi6GHMb6KU4o8gAAQEBgCRXvYrgWAAAgBNGTBwAAAoOHD17AGUUeAAAICDxd610UeQAAICAYMj1744XJgxdnY04eAABACKInDwAABIQw48zmyfH4BkUeAAAICB7PyaPIc+LycK3dbteMGTOUmpoqq9WqiRMnKj8/XxERETp27JgkaceOHbLZbJIkm82mQYMGyWq1asiQIY59zlVXV6fJkycrLS1NWVlZqqmpkSRVV1crMzNTaWlpstlsqq+vlyTNmjVLt9xyi1JSUpSbm+vOvQMAAIQsl4u8TZs2qaGhQTt37lRJSYmWL18uSUpOTtbChQtbPGbZsmUqKSnR8OHDtWDBghb3KSgoUFJSkrZv365Ro0YpLy9PkrR48WJlZmZq+/btSkxM1OrVqyVJeXl5evvtt7Vnzx7t3r1b77//fquZa2trVVZW5tjKy8tdvW0AANDBDJkeb+4oLS3V8OHDFRUVpUsuuUS5ubmOTqW2VFRUaObMmerbt6+ioqJ0zTXXaMWKFW5l6AguF3lRUVE6cOCA9u/fL9M0FRcXJ0kaO3asdu3a1WpPnSQNHjxYhw4davGzoqIijR49WpKUlZWloqKiNtu7d+8uSaqvr1ePHj2UmJjY6nXz8vLUp08fx5aSkuLiXQMAgI7WNFzryeaqqqoqDR06VPX19dqwYYMWLFigF154QTk5Oec9dvz48SosLNSTTz6pP/3pT7rtttt033336be//a0bd+99Lhd5VqtV06ZN06xZs9SvXz8tXbpUkmQYhnJzc1vtqZOkzZs3a8CAAS1+VlFRIYvFIkmyWCyqqKiQJFVWVio2NrZZuyTNmDFDl19+uRITE3XBBRe0et2cnBwdPXrUse3Zs8elewYAAKFpxYoVqq2t1auvvqpbb71V9957rxYvXqwVK1bo+PHjrR736aefqqioSAsWLJDNZtPQoUO1ZMkSWa1WvfLKKz68g9a5tYTK1KlTtW3bNu3bt09r1qxxdGmOHz9eu3fvVllZmdP+2dnZslqtqqio0Jw5c1o8Z3x8vKqqqiSdmYcXHx8vSYqLi1N1dXWzdkl6/vnndejQIX3++ed66623Ws0bExOj3r17O7aEhAR3bhsAAHQgf/Tkbdy4URkZGY6RSUmaMGGC7Ha7Nm/e3OpxDQ0NktSsk+mCCy6QGSDr9blc5B0/fly1tbWSpOjoaEVHRztuxjAMzZ49W4sWLXI6pmlOXkFBgaO37lzp6ekqLCyUJBUWFio9Pb3N9lOnTkmSunXrpujoaPXo0cPVWwEAAAHEMEyFebC5s5ByaWmpkpOTndpiY2OVkJCg0tLSVo/r06ePRowYoQULFmj//v06ceKE/vu//1ubN2/Wj3/8Y5dzdASXl1ApKyvTgw8+qLCwMDU2NmrUqFEKDw93fD5u3Dg99dRTLgex2WyaPn26rFarLBaLCgoKJEm5ubmaMmWKXnzxRSUlJelnP/uZpDNzAL/66ivV19fLarU6ij8AABCcvLWEyrkPWMbExCgmJqbFY6qqqhzTws5msVhUWVnZ5vU2bNigu+++2zEVrUuXLlq2bJnGjh3revgO4HKRl5KSonfeeafVzw3D0N69ex0/5+fnt+u84eHhWrVqVbN2i8Xi6Mk72xtvvNGu8wIAgM7l3Acs582bp/nz53v1GqZpaurUqfrwww+1Zs0aJSQk6C9/+YseeOABWSwWTZw40avXc4fPF0OurKzUmDFjnNqGDh2quXPn+joKAAAIIMb/bZ4cL0l79uxxmn/fWi+edKYzqWlt3rNVVVU5zdM71xtvvKG1a9dq7969uvbaayWdmWL2+eef66GHHuqcRV5cXJyKi4t9fVkAABDgDDfn1Z19vCQlJCSod+/e7TomOTm52dy7mpoalZeXN5urd7b9+/erS5cuuuaaa5zaBw0apN/97nc6efKk358XcOvpWgAAgFAwcuRIbdmyxbGShyStXbtWYWFhGjFiRKvHXXbZZTp9+rTTFDVJevfdd3XRRRf5vcCTKPIAAECA8McSKjNnzlTPnj2VlZWlzZs3a+XKlXrkkUc0c+ZMpxctDBs2TFdccYXj51GjRqlv374aN26cfv/732vr1q2aPXu28vPzlZ2d7Y2vw2M+H64FAABoSZikMA8m5bnTc2WxWLR161ZlZ2crKytLPXv21PTp0/XMM8847Xf69Gk1NjY6fu7Zs6e2bt2qxx9/XLNnz1Z1dbX69eunvLw8zZo1y/2b8CKKPAAA0Kn1799fW7ZsaXOflp4nuOKKK/THP/6xg1J5jiIPAAAEBG89eIEzKPIAAEDA8GQJFTijyAMAAAHBW2+8wBk8XQsAABCC6MkDAAABgTl53kWRBwAAAkKY4eESKgzXOmG4FgAAIATRkwcAAAICD154F0UeAAAICBR53sVwLQAAQAjq3D15dVXSqS7+TuGaHhf7O4HL5tWt8HcEtzwRPtPfEdwyr+55f0dwi/Ef0/0dofMw7f5O4J6GE/5O4J4u4f5OEDQMmTLkwdO1Hhwbijp3kQcAAAIGw7XeRZEHAAACg2HI8GQdFDtV3tmYkwcAABCC6MkDAAABwQg7s3lyPL5BkQcAAAICc/K8i5oXAAAgBNGTBwAAAgMvr/UqijwAABAQmJPnXXwdAAAAIYiePAAAEBAMw5DhwdMTnhwbiijyAABAQDAMD4drqfGcUOQBAIDAwBoqXsWcPAAAgBBETx4AAAgIPF3rXRR5AAAgIBhhhgwP1rrz5NhQRM0LAAAQgujJAwAAAYHnLryLIg8AAAQGD5dQEUWeE4ZrAQAAQhA9eQAAIDCEGWc2T46HA0UeAAAICMzJ8y6Xh2vtdrtmzJih1NRUWa1WTZw4Ufn5+YqIiNCxY8ckSTt27JDNZpMk2Ww2DRo0SFarVUOGDHHsc666ujpNnjxZaWlpysrKUk1NjSSpurpamZmZSktLk81mU319vdNxkydPVkZGhqu3AQAAAkzTEiqebPiGy0Xepk2b1NDQoJ07d6qkpETLly+XJCUnJ2vhwoUtHrNs2TKVlJRo+PDhWrBgQYv7FBQUKCkpSdu3b9eoUaOUl5cnSVq8eLEyMzO1fft2JSYmavXq1Y5j/v73v6u2tva8mWtra1VWVubYysvLXb1tAACAoOJykRcVFaUDBw5o//79Mk1TcXFxkqSxY8dq165drfbUSdLgwYN16NChFj8rKirS6NGjJUlZWVkqKipqs12SnnzyST3++OPnzZyXl6c+ffo4tpSUlPbdLAAA8JmmN154suEbLn8dVqtV06ZN06xZs9SvXz8tXbpUkmQYhnJzc1vtqZOkzZs3a8CAAS1+VlFRIYvFIkmyWCyqqKiQJFVWVio2NrZZ+5///GcNHjxYvXr1Om/mnJwcHT161LHt2bOn3fcLAAB8wzAMjzd8w62ad+rUqdq2bZv27dunNWvWOObJjR8/Xrt371ZZWZnT/tnZ2bJaraqoqNCcOXNaPGd8fLyqqqoknZmHFx8fL0mKi4tTdXW1U/vp06eVl5enhx56qF15Y2Ji1Lt3b8eWkJDgzm0DAAAEDZeLvOPHjzvmwUVHRys6OlqmaUo6U4HPnj1bixYtcjqmaU5eQUGBo7fuXOnp6SosLJQkFRYWKj09vdX2Tz/9VFVVVRo7dqymTJmif/zjH3ryySddvRUAABBIDJ2pTNzd6Mhz4nKRV1ZWppEjRyotLU2pqanKyMhQeHi44/Nx48bJbre7HMRms+ngwYOyWq0qLCxUTk6OJCk3N1fr16+X1WrV4cOHNWnSJF166aV677339NZbb6mgoEADBw7U3LlzXb4mAAAIHE1LqHiy4Rsur5OXkpKid955p9XPDcPQ3r17HT/n5+e367zh4eFatWpVs3aLxeLoyWtJUlKStmzZ0q5rAAAAdBY+Xwy5srJSY8aMcWobOnQoPXEAAHRyZ56Qdb87jqdrnfm8yIuLi1NxcbGvLwsAAAKcp8ugUOQ547VmAAAgMPBeM6+i5gUAAAhB9OQBAICAYBgeDtfSkeeEIg8AAASGMMOjBy/kybEhiOFaAACAEERPHgAACAg8d+FdFHkAACAgGB4O13o01BuCGK4FAAAIQfTkAQCAwGD83+bJ8XCgyAMAAAGBN154F0UeAAAICIbh4Zw8nrxwQs0LAAAQgujJAwAAAYElVLyLIg8AAAQEllDxrs5d5EVeKPW42N8pXHO63t8JXFdX6e8EbplX97y/I7jlifAZ/o7glnnmGn9HcI9p93eCziM8zt8J3BNs/94+3eDvBPCSzl3kAQCAwBEmz54W4EkDJxR5AAAgMIQZZzZPjocDNS8AAEAIosgDAACBwdA3Q7bubG525JWWlmr48OGKiorSJZdcotzcXNXXt28u5bFjxzRlyhT16tVLkZGR6t+/v1avXu1eEC9juBYAAAQGPwzXVlVVaejQobryyiu1YcMGHTt2TDk5OTp58qSee+65No8tLy/XjTfeqKuvvlovvPCCYmJi9P7776uurs7dO/AqijwAABAY/PDgxYoVK1RbW6tXX31VcXFnnuBubGzU/fffrzlz5igxMbHVY3Nzc9WnTx+99dZb6tKliyRp2LBhbkXvCAzXAgCATmvjxo3KyMhwFHiSNGHCBNntdm3evLnV42pra/Xf//3fuv/++x0FXqChyAMAAIGhabjWk81FpaWlSk5OdmqLjY1VQkKCSktLWz3u73//u+rr69WtWzfdcsst6tatmy655BLNnj1bDQ2BsdYgw7UAACAwhMnDOXln/qe8vNypOSYmRjExMS0eUlVVpdjY2GbtFotFlZWtL+b/6aefSpKmT5+uH/7wh5o/f7727NmjuXPnKiwsTAsXLnTvHryIIg8AAISUlJQUp5/nzZun+fPne/UadvuZt91kZGRo6dKlkqQhQ4boxIkTWrJkiebOnavIyEivXtNVFHkAACAwNC2h4snxkvbs2aOEhARHc2u9eNKZHruamppm7VVVVU7z9Fo6TpKGDh3q1D5s2DA988wzOnjwoK699lpX0nsdRR4AAAgMXlpCJSEhQb17927XIcnJyc3m3tXU1Ki8vLzZXL2zffvb327zvKdOnWrX9TsSD14AAIBOa+TIkdqyZYuqq6sdbWvXrlVYWJhGjBjR6nGXXXaZrr32Wm3ZssWp/S9/+YsiIyPPWwT6AkUeAAAIDJ687cLNNfZmzpypnj17KisrS5s3b9bKlSv1yCOPaObMmU5r5A0bNkxXXHGF07HPPPOMCgsL9cADD+gvf/mLFixYoCVLlignJ0dRUVGuh/EyhmsBAEBg8MMbLywWi7Zu3ars7GxlZWWpZ8+emj59up555hmn/U6fPq3GxkantjvvvFN/+MMf9NRTT+k3v/mNEhIS9MQTT+jRRx91/x68iCIPAAAEBkNuv3/Wcbwb+vfv32zY9VzFxcUttt999926++673btwB2O4FgAAIATRkwcAAAKD4eFwreFJN2DoocgDAACBwQ9z8kIZw7UAAAAhyOUiz263a8aMGUpNTZXVatXEiROVn5+viIgIHTt2TJK0Y8cO2Ww2SZLNZtOgQYNktVo1ZMgQxz7nqqur0+TJk5WWlqasrCzH6tPV1dXKzMxUWlqabDab6uvrnc6bnp6uMWPGuHPvAAAgkPhhCZVQ5vLXsWnTJjU0NGjnzp0qKSnR8uXLJZ1ZMbq1l/EuW7ZMJSUlGj58uBYsWNDiPgUFBUpKStL27ds1atQo5eXlSZIWL16szMxMbd++XYmJiVq9erXTeYuLi7VhwwZXbwMAAASapuFaTzY4uFzkRUVF6cCBA9q/f79M03S8123s2LHatWtXqz11kjR48GAdOnSoxc+Kioo0evRoSVJWVpaKiorabJeknJwcpaWlac2aNW1mrq2tVVlZmWMrLy9v/w0DAAAEIZcfvLBarZo2bZpmzZqljz/+WNnZ2YqPj5dhGMrNzdWCBQt0zz33tHjs5s2bNWDAgBY/q6iocLzs12KxqKKiQpJUWVmp2NjYZu1LlizRhRdeqOrqag0bNkzXX3+9rrzyyhbPnZeXpyeeeMLVWwUAAD5kGJLhwZArD9c6c+urnDp1qrZt26Z9+/ZpzZo1jnly48eP1+7du1VWVua0f3Z2tqxWqyoqKjRnzpwWzxkfH6+qqipJZ+bhxcfHS5Li4uIc75M7u/3CCy+UJMXGxurWW2/VP/7xj1bz5uTk6OjRo45tz5497tw2AADoSAzXepXLRd7x48dVW1srSYqOjlZ0dLRM05QkGYah2bNna9GiRU7HNM3JKygocPTWnSs9PV2FhYWSpMLCQqWnp7fZ3lT4NTQ0aMeOHbrqqqtazRwTE6PevXs7toSEBFdvGwAAdDQevPAql7+OsrIyjRw5UmlpaUpNTVVGRobCw8Mdn48bN052u93lIDabTQcPHpTValVhYaFycnIkSbm5uVq/fr2sVqsOHz6sSZMmSTrzGpGbbrpJN910k7KysvSd73zH5WsCAACEKpfn5KWkpOidd95p9XPDMLR3717Hz/n5+e06b3h4uFatWtWs3WKxOHryzrZp06Z2nRcAAAQJFkP2Kp+/8aKysrLZunZDhw7V3LlzfR0FAAAEEoo8r/J5kRcXF6fi4mJfXxYAAKBT4d21AAAgMBjy7OEJOvKcUOQBAIDAECYPh2u9liQk8HUAAACEIHryAABAYPB0rTu6rpxQ5AEAgMDA07VeRZEHAAACAw9eeBUdmwAAACGInjwAABAYGK71Koo8AAAQGHjwwqv4OgAAAEIQPXkAACAwMFzrVRR5AAAgMFDkeRXDtQAAACGIIg8AAASGpnXy3N2CrCPvq6++0n/913/pn//8p6Pt+eef99r5KfIAAEBgaBqu9WQLIlOnTtW6det0xx13aNWqVZKkP/7xj147f+eek2faz2zBpKHW3wlcZh7Y4u8IbjH+Y7q/I7hlnrnG3xHc8oTxn/6O4JZ59S/4O4Lrwrr7O4F7Tn3h7wTu6ZHg7wSuCevix2urUy2h8tFHH+ndd9/Vp59+qttuu00xMTFePX/nLvIAAAD8JDY2VpJ0ySWXaPPmzcrIyFBtrfc6c4Ks5gUAACHLMDzfgki/fv30xRdfaMOGDfryyy+1fv16devWzWvnp8gDAACBwfDCFkR+97vfqVevXvrP//xPHTlyRFdeeaU+/PBDr52fIg8AAMCPkpOTdeTIEa+flyIPAAAEhk42XNtk4cKFevrpp/W3v/3Nq+flwQsAABA4grNO80hubq4qKir0ve99T/Hx8br44otlnFWwGobhtJZee1HkAQAA+NF//Md/6Lvf/a7Xz0uRBwAAAoOnQ65BOlybn5/fIeelyAMAAIGhky2G3NEo8gAAQGAw5GFPnteShARqXgAAgBBETx4AAAgMni5oTE+eE4o8AAAQGDrpgxcdheFaAACAEERPHgAACAwM13oVRR4AAAgMDNd6FcO1AAAAIYiePAAAEBgMedb9REeeE4o8AAAQGBiu9SqKPAAAEBh48MKrXO4UtdvtmjFjhlJTU2W1WjVx4kTl5+crIiJCx44dkyTt2LFDNptNkmSz2TRo0CBZrVYNGTLEsc+56urqNHnyZKWlpSkrK0s1NTWSpOrqamVmZiotLU02m0319fWSpK+//lozZ85URkaG0tPT9cEHH7hz/wAAACHJ5SJv06ZNamho0M6dO1VSUqLly5dLkpKTk7Vw4cIWj1m2bJlKSko0fPhwLViwoMV9CgoKlJSUpO3bt2vUqFHKy8uTJC1evFiZmZnavn27EhMTtXr1aknSk08+qbvuuktbtmxRcXGxrr766lYz19bWqqyszLGVl5e7etsAAKCjNQ3XerLBweUiLyoqSgcOHND+/ftlmqbi4uIkSWPHjtWuXbta7amTpMGDB+vQoUMtflZUVKTRo0dLkrKyslRUVNRm++bNm1VUVKT09HTl5OSosbGx1evm5eWpT58+ji0lJcXV2wYAAB2MGs+7XC7yrFarpk2bplmzZqlfv35aunSpJMkwDOXm5rbaUyedKcwGDBjQ4mcVFRWyWCySJIvFooqKCklSZWWlYmNjm7X/7//+r1JTU1VcXKyGhga9/PLLrV43JydHR48edWx79uxx9bYBAACCilsPKk+dOlXbtm3Tvn37tGbNGsc8ufHjx2v37t0qKytz2j87O1tWq1UVFRWaM2dOi+eMj49XVVWVpDPz8OLj4yVJcXFxqq6ubrF95MiRkqTbb79d//jHP1rNGxMTo969ezu2hIQEd24bAAB0JLryvMrlIu/48eOqra2VJEVHRys6OlqmaUo605s3e/ZsLVq0yOmYpjl5BQUFjt66c6Wnp6uwsFCSVFhYqPT09Dbbhw0bpr/97W+SpN27d+uqq65y9VYAAEAgMbywwcHlIq+srEwjR45UWlqaUlNTlZGRofDwcMfn48aNk91udzmIzWbTwYMHZbVaVVhYqJycHElSbm6u1q9fL6vVqsOHD2vSpEmSpEWLFunpp59Wenq63n//fU2fPt3lawIAAIQql9fJS0lJ0TvvvNPq54ZhaO/evY6f8/Pz23Xe8PBwrVq1qlm7xWJx9OSdrXfv3nrrrbfadW4AABAEDEMKYzFkb/H5YsiVlZUaM2aMU9vQoUM1d+5cX0cBAACBhMWQvcrnRV5cXJyKi4t9fVkAABDoeK2ZV3nyGmAAAICgV1paquHDhysqKkqXXHKJcnNzHSuHtNcvf/lLGYahO+64o4NSuo531wIAgMDgh+HaqqoqDR06VFdeeaU2bNigY8eOKScnRydPntRzzz3XrnN8+umneuKJJ3TRRRe5HqADUeQBAIDA4Ifh2hUrVqi2tlavvvqq4y1ejY2Nuv/++zVnzhwlJiae9xy5ubm66667dPjwYZev35EYrgUAAJ3Wxo0blZGR4SjwJGnChAmy2+3avHnzeY/fsWOHXnvttWZrBAcCijwAABAY/LAYcmlpqZKTk53aYmNjlZCQoNLS0jaPPX36tGbNmqXHH388IN+mxXAtAAAIDGHybJ28/+u6Ki8vd2qOiYlRTExMi4dUVVUpNja2WbvFYlFlZWWbl1u+fLm++uorPfjgg27F7Wj05AEAgJCSkpKiPn36OLa8vDyvX+Pzzz/X3LlzlZeXp+7du3v9/N5ATx4AAAgMXnq6ds+ePU7Dp6314klneuxqamqatVdVVTnN0zvX3Llzdd111yktLU3V1dWSzjyw0djYqOrqakVHR6trV/+WWRR5AAAgMHjp6dqEhAT17t27XYckJyc3m3tXU1Oj8vLyZnP1zlZaWqqSkhJZLJZmn1ksFm3cuFG33XabC+G9jyIPAAAEDh+/tGLkyJFasGCBqqurHXPz1q5dq7CwMI0YMaLV4375y186evCaPPDAA4qMjNTChQt13XXXdWDq9qHIAwAAndbMmTO1bNkyZWVlac6cOTp27JgeeeQRzZw502mNvGHDhunw4cM6ePCgJGngwIHNzhUbG6vo6Gilp6f7KH3bKPIAAEBg8MNiyBaLRVu3blV2draysrLUs2dPTZ8+Xc8884zTfqdPn1ZjY6P72fyAIg8AAAQGP7zWTJL69++vLVu2tLlPcXHxec/Tnn18iSVUAAAAQhA9eQAAIDCEGR4uhuzjpzYCXOcu8oywM1swibjQ3wlcZgyc7O8InYtp93cCt8yrf8HfEdzyRPcf+TuCy+Y1vOjvCO6JutTfCdxiHjn/+08DiVnefM04n/HDnLxQFmQVDgAAANqjc/fkAQCAwEFPnldR5AEAgMBAkedVFHkAACAweDpXPtjm2Xcwvg0AAIAQRE8eAAAIDIY8WwaF0VonFHkAACAwMCfPqxiuBQAACEH05AEAgMDAgxdeRZEHAAACA8O1XkXJCwAAEILoyQMAAIEhzPDs6VpPjg1BFHkAACAwMCfPqyjyAABAYGBOnldR8gIAAIQgevIAAECA8LAnj1deOKHIAwAAgcEwPJyTR5F3NoZrAQAAQhA9eQAAIDCwhIpXUeQBAIDAwNO1XuXycK3dbteMGTOUmpoqq9WqiRMnKj8/XxERETp27JgkaceOHbLZbJIkm82mQYMGyWq1asiQIY59zlVXV6fJkycrLS1NWVlZqqmpkSRVV1crMzNTaWlpstlsqq+vlySNGTNG6enpSk9PV2xsrP70pz+5c/8AAAAhyeUib9OmTWpoaNDOnTtVUlKi5cuXS5KSk5O1cOHCFo9ZtmyZSkpKNHz4cC1YsKDFfQoKCpSUlKTt27dr1KhRysvLkyQtXrxYmZmZ2r59uxITE7V69WpJ0oYNG1RcXKy33npLsbGxGjFiRKuZa2trVVZW5tjKy8tdvW0AANDRmhZD9mSDg8vfRlRUlA4cOKD9+/fLNE3FxcVJksaOHatdu3a12lMnSYMHD9ahQ4da/KyoqEijR4+WJGVlZamoqKjN9iavvfaabr31VoWHh7d63by8PPXp08expaSktP+GAQCAbzQN13qywcHlIs9qtWratGmaNWuW+vXrp6VLl0qSDMNQbm5uqz11krR582YNGDCgxc8qKipksVgkSRaLRRUVFZKkyspKxcbGNmtvsmrVKv3gBz9oM3NOTo6OHj3q2Pbs2dOuewUAAD7U9OCFJxsc3OrXnDp1qrZt26Z9+/ZpzZo1jnly48eP1+7du1VWVua0f3Z2tqxWqyoqKjRnzpwWzxkfH6+qqipJZ+bhxcfHS5Li4uJUXV3drF2SPvvsM3344Ye66aab2swbExOj3r17O7aEhAR3bhsAACBouFzkHT9+XLW1tZKk6OhoRUdHyzRNSWd682bPnq1FixY5HdM0J6+goMDRW3eu9PR0FRYWSpIKCwuVnp7eZrsk/eEPf9DEiRNl0D0LAEAIMDycj0c9cDaXi7yysjKNHDlSaWlpSk1NVUZGhtN8uHHjxslut7scxGaz6eDBg7JarSosLFROTo4kKTc3V+vXr5fVatXhw4c1adIkxzGrVq1y+hkAAAQx5uR5lcvr5KWkpOidd95p9XPDMLR3717Hz/n5+e06b3h4uFatWtWs3WKxOHryzvXuu++269wAAACdjc8XQ66srNSYMWOc2oYOHaq5c+f6OgoAAAgkLIbsVT4v8uLi4lRcXOzrywIAgEBHkedVrBoIAAAQgnh3LQAACAxhYWc2T46HA0UeAAAIDAzXehVFHgAACAwUeV5FvyYAAEAIoicPAAAEBuP/3njhyfFwoMgDAACBwTCkMIZrvYXhWgAAgBBETx4AAAgMPHjhVRR5AAAgMBhhHs7JY4DybHwbAAAAIYiePAAAEBgYrvUqijwAABAYKPK8iiIPAAAEBt5d61V8GwAAACGoc/fkmfYzWzAJxieHunT3dwL3BNvvRrALC87fk3kNL/o7gsue6DbN3xHcMu/0Kn9HcIvRd4S/I7jECPtM0s/9mcCP1w4tnbvIAwAAgYM5eV4VhN1CAAAAOB968gAAQIDwcDFk+q6cUOQBAIAAYcizOXkM156NkhcAACAE0ZMHAAACAw9eeBVFHgAACAyGh3PygnGZsQ5EkQcAAAIEc/K8iZIXAAAgBNGTBwAAAoMhD+fkeS1JSKDIAwAAASJMng0yMkB5Nr4NAADQqZWWlmr48OGKiorSJZdcotzcXNXX17d5THl5uXJzczVw4ED17NlTvXv31n/+53/q8OHDPkp9fvTkAQCAwOCHJVSqqqo0dOhQXXnlldqwYYOOHTumnJwcnTx5Us8991yrx7377rvasGGD7r33Xt1www3697//raeeekopKSn63//9X/Xq1cv9+/ASijwAABAY/FDkrVixQrW1tXr11VcVFxcnSWpsbNT999+vOXPmKDExscXjbr75ZpWWlqpr129KqdTUVPXt21cvv/yyHnroIffuwYsYrgUAAJ3Wxo0blZGR4SjwJGnChAmy2+3avHlzq8fFxsY6FXiS1Lt3b/Xq1UvHjx/vsLyuoMgDAAABwvDC5prS0lIlJyc7tcXGxiohIUGlpaUunevAgQP6/PPP1b9/f5dzdASGawEAQGDw0hsvysvLnZpjYmIUExPT4iFVVVWKjY1t1m6xWFRZWdnuS5umqZ/85CdKTEzUPffc0/7MHYiePAAAEBia5uR5sklKSUlRnz59HFteXl6HR58/f762bt2ql19+WVFRUR1+vfagJw8AAISUPXv2KCEhwfFza7140pkeu5qammbtVVVVTvP02vLb3/5WTz75pF588UUNGzbM9cAdxOWePLvdrhkzZig1NVVWq1UTJ05Ufn6+IiIidOzYMUnSjh07ZLPZJEk2m02DBg2S1WrVkCFDHPucq66uTpMnT1ZaWpqysrIcX3h1dbUyMzOVlpYmm83mWLdm48aNuv7665WWlqZ77rlHDQ0N7tw/AAAIKJ7Px0tISFDv3r0dW1tFXnJycrO5dzU1NSovL282V68lr776qu677z49+eSTuvfee1250Q7ncpG3adMmNTQ0aOfOnSopKdHy5cslnfmSFi5c2OIxy5YtU0lJiYYPH64FCxa0uE9BQYGSkpK0fft2jRo1ytG1unjxYmVmZmr79u1KTEzU6tWrJUk/+9nPtG7dOm3fvl1RUVF66623XL0VAAAQSJrm5HmyuWjkyJHasmWLqqurHW1r165VWFiYRowY0eaxxcXFuueee/TDH/5QP/vZz1y+dkdz+duIiorSgQMHtH//fpmm6ejKHDt2rHbt2tVqT50kDR48WIcOHWrxs6KiIo0ePVqSlJWVpaKiojbbBwwYoOrqapmmqerq6jYXHaytrVVZWZljO3dCJgAA6Jxmzpypnj17KisrS5s3b9bKlSv1yCOPaObMmU5r5A0bNkxXXHGF4+d//etfysrK0pVXXqkf/OAH2rVrl2P76KOP/HErzbg8J89qtWratGmaNWuWPv74Y2VnZys+Pl6GYSg3N1cLFixo9amSzZs3a8CAAS1+VlFRIYvFIunM+HhFRYUkqbKy0vHUy9ntkyZN0m233aaYmBhde+21uuGGG1rNnJeXpyeeeMLVWwUAAD5kyJDhwWLIhhtLqFgsFm3dulXZ2dnKyspSz549NX36dD3zzDNO+50+fVqNjY2On3fv3q2amhrV1NTopptuctp3ypQpys/Pd+sevMmtp2unTp2qbdu2ad++fVqzZo1jntz48eO1e/dulZWVOe2fnZ0tq9WqiooKzZkzp8VzxsfHq6qqStKZeXjx8fGSpLi4OEcX6tntM2fO1K5du/TBBx/owgsv1B//+MdW8+bk5Ojo0aOObc+ePe7cNgAA6FC+XydPkvr3768tW7bo5MmT+uyzz/SLX/xC3bt3d9qnuLhYn3zyieNnm80m0zRb3AKhwJPc6Mk7fvy4oqOjFRMTo+joaEVHR8s0TUmSYRiaPXu2nnnmGQ0cONBxzLJly3TzzTe3ed709HQVFhZq0KBBKiwsVHp6ulP7lClTnNq7du3q6Pnr1auX01j6udpaHwcAACAUudyTV1ZWppEjRyotLU2pqanKyMhQeHi44/Nx48bJbre7HMRms+ngwYOyWq0qLCxUTk6OJCk3N1fr16+X1WrV4cOHNWnSJEnS008/rWHDhslqteof//iHvv/977t8TQAAEED88OBFKHO5Jy8lJUXvvPNOq58bhqG9e/c6fm5vl2V4eLhWrVrVrN1isaiwsLBZ+/jx4zV+/Ph2nRsAAAQD94dcvzkeTXy+GHJlZaXGjBnj1DZ06FDNnTvX11EAAEAgOeutFW4fDwefF3lxcXEqLi729WUBAAA6FV5rBgAAAoNheDavjp48JxR5AAAgQDAnz5t4DAUAACAE0ZMHAAACAw9eeBVFHgAACAyernXHOnlO+DYAAABCED15AAAgQPDghTdR5AEAgMDAnDyvosgDAAABwpBnM8ko8s7GnDwAAIAQRE8eAAAIDAzXehVFHgAACAwUeV7FcC0AAEAIoicPAAAEiDB51v9E39XZKPIAAECA8HC4lqdrnXTuIu/EUan2pL9TuKZnX38ncJ1p93cC9zSc8HcC94TH+TuBe0594e8E7om61N8JXDbv9Cp/R3DLE11+4O8Ibpl3arm/I7imrtbfCeAlnbvIAwAAAYQ3XngTRR4AAAgMRtiZzZPj4UCRBwAAAoMhD5dQ8VqSkEDJCwAAEILoyQMAAAGCOXneRJEHAAACA3PyvIpvAwAAIATRkwcAAAIEw7XeRJEHAAACg+HhGy88eltG6GG4FgAAIATRkwcAAAJEmDzrf6Lv6mwUeQAAIDAwXOtVFHkAACAwGIaHS6hQ5J2Nfk0AAIAQRE8eAAAIECyh4k0UeQAAIDAwJ8+rGK4FAAAIQfTkAQCAAMESKt5EkQcAAAIDw7VeRckLAAAQglwu8ux2u2bMmKHU1FRZrVZNnDhR+fn5ioiI0LFjxyRJO3bskM1mkyTZbDYNGjRIVqtVQ4YMcexzrrq6Ok2ePFlpaWnKyspSTU2NJKm6ulqZmZlKS0uTzWZTfX29JOmvf/2rI0NeXp479w4AAAJKmBc2NHH529i0aZMaGhq0c+dOlZSUaPny5ZKk5ORkLVy4sMVjli1bppKSEg0fPlwLFixocZ+CggIlJSVp+/btGjVqlKNwW7x4sTIzM7V9+3YlJiZq9erVkqSf/OQneuWVV/T2229r06ZN+uijj1rNXFtbq7KyMsdWXl7u6m0DAIAOZ3wzZOvOxhIqTlwu8qKionTgwAHt379fpmkqLi5OkjR27Fjt2rWr1Z46SRo8eLAOHTrU4mdFRUUaPXq0JCkrK0tFRUVttldXV6tv374yDEODBg3S22+/3ep18/Ly1KdPH8eWkpLi6m0DAICOZsizIo8az4nLRZ7VatW0adM0a9Ys9evXT0uXLpUkGYah3NzcVnvqJGnz5s0aMGBAi59VVFTIYrFIkiwWiyoqKiRJlZWVio2NbdZ+4YUX6p///Kfq6+tVVFSkysrKVq+bk5Ojo0ePOrY9e/a4etsAAABBxa2na6dOnaqpU6fqxIkTSk9P14wZMyRJ48eP1+LFi5WWlua0f3Z2tnr27Kl+/frpl7/8ZYvnjI+PV1VVlZKSklRdXa34+HhJUlxcnKqrq2WxWJzaX3jhBT388MMyDENXXnmlEhMTW80bExOjmJgYd24VAAD4DEuoeJPL38bx48dVW1srSYqOjlZ0dLRM05R0pjdv9uzZWrRokdMxTXPyCgoKHL1150pPT1dhYaEkqbCwUOnp6W22X3vttdq0aZMKCwtVXV2tW2+91dVbAQAAgcSjoVoPl18JQS4XeWVlZRo5cqTS0tKUmpqqjIwMhYeHOz4fN26c7Ha7y0FsNpsOHjwoq9WqwsJC5eTkSJJyc3O1fv16Wa1WHT58WJMmTZJ0Zp7dkCFDdNttt2nWrFmOHj4AAAC4MVybkpKid955p9XPDcPQ3r17HT/n5+e367zh4eFatWpVs3aLxeLoyTtbTk6OoxAEAAChwNMnZOnJO5vP33hRWVmpMWPGOLUNHTpUc+fO9XUUAAAQSIywM5snx8PB50VeXFyciouLfX1ZAACAToV31wIAgADBcK03UeQBAIDAYBgeDtdS5J2NIg8AAAQIevK8iRmKAAAAIYiePAAAEBg8XdCY4VonFHkAACAwsISKV/FtAAAAhCB68gAAQIDgwQtvosgDAACBgTl5XsVwLQAAQAiiyAMAAAEizAub60pLSzV8+HBFRUXpkksuUW5ururr6897nGmaWrRokfr27avIyEjdeOON2rVrl1sZOgJFHgAACAyGvhmydWtz/ZJVVVUaOnSo6uvrtWHDBi1YsEAvvPCCcnJyznvsz3/+c82bN08PPvig/vznPyshIUEjRozQxx9/7HqQDsCcPAAAECDc74375njXrFixQrW1tXr11VcVFxcnSWpsbNT999+vOXPmKDExscXjTp06pYULF+qhhx7Sgw8+KElKS0vTVVddpSVLlmj58uXu34aX0JMHAAA6rY0bNyojI8NR4EnShAkTZLfbtXnz5laP27lzp2prazVhwgRHW/fu3TVmzBi9+eabHZq5vejJAwAAAaH80wqPnpAt/7TizP+Wlzu1x8TEKCYmpsVjSktLde+99zq1xcbGKiEhQaWlpa1eq+mz5ORkp/b+/fvryJEj+vrrrxUZGenyPXhT5y7yontLPS/2dwqXmBX/9HcElxm9Bvs7gnu6hPs7gXtOn3+ycEDqkeDvBG4xj7T+X/qByug7wt8R3DLvlP+Hv9zxRMT9/o7gkho1+Pya0dHRiouLU0rKFI/PFRkZqZSUFKe2efPmaf78+S3uX1VVpdjY2GbtFotFlZWVrV6nqqpK4eHhioiIaHacaZqqqqqiyAMAAJ1bbGysPvroI3355Zcen8tutysszHk2Wmu9eKGOIg8AAPhdbGxsiz1qHc1isaimpqZZe1VVldM8vZaOq6ur06lTp5x686qqqmQYhiwWS4fkdQUPXgAAgE4rOTm52dy7mpoalZeXN5tvd+5xkvTBBx84tZeWljrWzfM3ijwAANBpjRw5Ulu2bFF1dbWjbe3atQoLC9OIEa3PX01NTVVMTIzWrl3raGtoaNCGDRs0atSojozcbhR5AACg05o5c6Z69uyprKwsbd68WStXrtQjjzyimTNnOq2RN2zYMF1xxRWOnyMiIvTYY49pyZIl+tWvfqVt27bpnnvuUUVFhR5++GF/3EozzMkDAACdlsVi0datW5Wdna2srCz17NlT06dP1zPPPOO03+nTp9XY2OjUNnv2bJmmqSVLluiLL77QwIEDtWnTJl1++eW+vIVWUeQBAIBOrX///tqyZUub+xQXFzdrMwxDjz32mB577LEOSuYZhmsBAABCEEUeAABACKLIAwAACEEUeQAAACGIIg8AACAEUeQBAACEIIo8AACAEESRBwAAEIIo8gAAAEIQRR4AAEAIosgDAAAIQRR5AAAAIYgiDwAAIAS5XOTZ7XbNmDFDqampslqtmjhxovLz8xUREaFjx45Jknbs2CGbzSZJstlsGjRokKxWq4YMGeLY51x1dXWaPHmy0tLSlJWVpZqaGknSunXr1L9/f0VERDjtv3fvXt18881KTU3V0qVLXb0NAACAkOZykbdp0yY1NDRo586dKikp0fLlyyVJycnJWrhwYYvHLFu2TCUlJRo+fLgWLFjQ4j4FBQVKSkrS9u3bNWrUKOXl5UmS0tPT9d5776l3795O+//4xz/WypUrtWPHDr3++uv66KOPWs1cW1ursrIyx1ZeXu7qbQMAAAQVl4u8qKgoHThwQPv375dpmoqLi5MkjR07Vrt27Wq1p06SBg8erEOHDrX4WVFRkUaPHi1JysrKUlFRkSTpwgsvbNaLV1dXp9raWl155ZUKCwvTHXfcobfffrvV6+bl5alPnz6OLSUlxaV7BgAACDYuF3lWq1XTpk3TrFmz1K9fP8dQqWEYys3NbbWnTpI2b96sAQMGtPhZRUWFLBaLJMlisaiioqLV85y9b3v2z8nJ0dGjRx3bnj172rxHAACAYOfWgxdTp07Vtm3btG/fPq1Zs0b19fWSpPHjx2v37t0qKytz2j87O1tWq1UVFRWaM2dOi+eMj49XVVWVJKm6ulrx8fGtXj8uLk7V1dWOn8+3f0xMjHr37u3YEhIS2nurAAAAQcnlIu/48eOqra2VJEVHRys6OlqmaUo605s3e/ZsLVq0yOmYpjl5BQUFTj1wZ0tPT1dhYaEkqbCwUOnp6a1miIiIUM+ePfXxxx/LNE29+eabslqtrt4KAABAyHK5yCsrK9PIkSOVlpam1NRUZWRkKDw83PH5uHHjZLfbXQ5is9l08OBBWa1WFRYWKicnR5JUXFysjIwMHT9+XBkZGVq7dq0k6de//rUmT56sm266SaNGjdIVV1zh8jUBAABClWE2dcN1ImVlZerTp4+OHi5U794X+zuOS8yKf/o7gsuMXoP9HaFzOV3v7wTuCevq7wRuMY9u8XcElxl9R/g7gnvqqv2dwC1PRNzv7wguqVGDntUhHT16tNnKFgguPv+3amVlpcaMGePUNnToUM2dO9fXUQAAAEKWz4u8uLg4FRcX+/qyAAAAnQqvNQMAAAhBFHkAAAAhiCIPAAAgBFHkAQAAhCCKPAAAgBBEkQcAABCCKPIAAABCEEUeAABACKLIAwAACEEUeQAAACGIIg8AACAEUeQBAACEIIo8AACAEESRBwAAEIK6+juAPzQ2NkqSysv/7eckrjOrqvwdwWVG3Wf+jtC5nG7wdwL3hHXxdwK3mOU1/o7gMiMsSP9M1tX6O4FbahRcfyZP6MzfkU1/VyJ4dcoi74svvpAkpdxwr5+TAIA//NzfARAEvvjiCyUlJfk7BjxgmKZp+juEr506dUr79u1Tr1691LWrd+vc8vJypaSkaM+ePUpISPDquTtSMOYOxswSuX0pGDNL5Pa1YMzdkZkbGxv1xRdf6Nprr1VERIRXzw3f6pQ9eREREbr++us79BoJCQnq3bt3h16jIwRj7mDMLJHbl4Ixs0RuXwvG3B2VmR680MCDFwAAACGIIg8AACAEUeR5WUxMjObNm6eYmBh/R3FJMOYOxswSuX0pGDNL5Pa1YMwdjJnhe53ywQsAAIBQR08eAABACKLIAwAACEEUeQAAACGIIg8AACAEUeQBAACEIIo8AACAEESR14Heeecdf0c4r0OHDunIkSNObX/729/8lKZ9/vWvf6mmpkaS9NZbb6mkpMTPidzz9NNP+zuCSz766COtX79e+/bt83eUNn344YeSJLvdrvz8fD3wwANavny56uvr/ZysbW+88YYaGhr8HcNlp0+f1htvvKE9e/ZIktauXavf/va3OnnypJ+Tnd+hQ4eUn5+vvLw8vfLKK6qqqvJ3pPNqyrx06VLl5+fr448/9nckBDDWyetAI0aM0ObNm/0do1WLFy/Wa6+9pq5du6pv375auXKlunXrpqFDh2rbtm3+jteiRx99VLt27VJUVJQSExP12WefKTo6WklJSVqwYIG/47VqxIgRTj+bpql//OMfGjRoUED/jowbN07r1q3TH/7wBy1dulS33nqrduzYoTFjxuinP/2pv+O1qOn396GHHlJDQ4Nuv/127d69W/v379crr7zi73itslgsslgsuvXWW/WDH/xAqamp/o7ULt///vfVo0cPVVdXyzAMxcfHq1evXvrnP/+p119/3d/xWvXss89q06ZNGjRokIqKinTVVVepurpas2bNavbnNVDk5eXp9ddf1x133CGLxaLq6mq9+eabuv322/XQQw/5Ox4CkQmPXXnllc22K664woyKivJ3tDbdeOONjv+/atUq89ZbbzVra2vN9PR0P6Zq2/e+9z3TNE2zrq7O7N27t9nY2GiapmnedNNN/ox1Xj/4wQ/MrKws829/+5v5ySefmIcOHTJvueUW85NPPvF3tDYNGTLENE3TtFqtZnV1tWmapllfX+/45xCIhg8fbpqm2ez3uOleAlV6erp56tQp85VXXjFvv/128+qrrzbnzZtnfvjhh/6O1iar1Wqapmna7XYzOTnZ0X7LLbf4KVH7nJ2vrq7OHDlypFlfX28OGzbMf6HO46abbjLtdrtTm91uD/h//8F/uvq7yAwFhmHo/fffV9euzl/n8OHD/ZSofcyzOnEnTZqk+Ph4jRw5Ul9++aUfU7UtIiJCktS9e3fdfvvt6tKli+PnQPbyyy/rvffe05NPPqn/+I//0EMPPaQePXrosssu83e0NtXU1Gjnzp2qq6vTBRdcIEnq1q2bunXr5udkrfvWt76lt956SykpKSosLNRtt92m//mf/1F4eLi/o7XJMAyFh4fr7rvv1t13363PP/9cf/jDH3TPPffof/7nf/wdr1Vff/21ysrKVFlZqZqaGpWVlSk2Njbgh2sbGhp07NgxXXrppfrggw9kt9vVrVs32e12f0drlWEYOnTokC6//HJH2+HDh2UYhh9TIZBR5HlBdna2qqqq1KtXL6f26dOn+ylR+1x33XU6ePCgrrjiCknSyJEjFRUVpSlTpvg5WeuuvfZaNTY2qmvXrlqxYoUkqa6uTj169PBzsvMbNGiQXn31Vb3xxhu64447VF1d7e9I53Xdddfpt7/9rfr376/KykrFxcWppqZGUVFR/o7WqmeffVZPPvmktm7dql/96leKj4/XsGHD9Lvf/c7f0dpksVicfr7ooov005/+NGCHxZs88cQTGj16tK655hq9/PLLSk9PV7du3QJ+zunSpUs1evRoVVdXq1evXnrhhRckSXfeeaefk7Vu+fLlstlsOnHihCwWi6qqqhQTE6Nly5b5OxoCFHPyAD+w2+369NNPlZiY6O8oAIJMfX29KioqFBcXF/A91PAvijwvaGxs1OrVq1VUVKSKigrFx8crPT1dkyZNajaEG0iCMXcwZpa+yb1t2zZVVlYGXe5g+r6DMbNEbl8L1twtWb9+vcaOHevvGAhAFHlecO+99yohIUFZWVmOJ54KCwt15MgR5efn+zteq4IxdzBmlsjtS8GYWSK3rwVr7pYE+koO8CM/PvQRMpqeLjtXoD9dFoy5gzGzaZLbl4Ixs2mS29eCMXdkZGRQruQA/2ExZC+wWCzKz893LNBbU1Ojl19+2fE0YqAKxtzBmFkity8FY2aJ3L4WjLmTkpK0b98+HThwwLF9+OGHuvHGG/0dDQGK4VovqKmp0eLFi1VUVOQ03+qRRx5RbGysv+O1KhhzB2Nmidy+FIyZJXL7WjDmXr9+vYYMGaK4uDin9rffflu33HKLn1IhkFHkAQAAhKDgeoQoQJ04cUJLly5VcXGx01NaOTk5iomJ8Xe8VgVj7mDMLJHbl4Ixs0RuXwvG3E2Zz+19DOTM8C/m5HnBlClTdOmll2rt2rV67733tG7dOvXt2zegFxWWgjN3MGaWyO1LwZhZIrevBWPupszr1q0LmszwM/8+9xEa0tLSWmxv7emtQBGMuYMxs2mS25eCMbNpktvXgjF3MGaGfzFc6wX9+vXT/PnzNXr0aMerZl5//fWAfy9pMOYOxswSuX0pGDNL5Pa1YMwdjJnhXzx44QX19fUqKChwrJx+4YUX6pZbbpHNZlP37t39Ha9VwZg7GDNL5PalYMwskdvXgjF3MGaGn/m7K7EzeOqpp/wdwS3BmDsYM5smuX0pGDObJrl9LRhzB2NmdCwevPCBbdu2+TuCW4IxdzBmlsjtS8GYWSK3rwVj7mDMjI5FkecDZpCOiAdj7mDMLJHbl4Ixs0RuXwvG3MGYGR2LIs8HDMPwdwS3BGPuYMwskduXgjGzRG5fC8bcwZgZHYsizweC9b+ugjF3MGaWyO1LwZhZIrevBWPuYMyMjsXTtV5QXl6uZ555Rp988okGDBigOXPmOL3kuqqqShaLxY8JWxaMuYMxs0RuXwrGzBK5fS0YcwdjZvgXPXleMGnSJF166aXKzs5WY2OjHnjgAafPA/UPXTDmDsbMErl9KRgzS+T2tWDMHYyZ4We+fJQ3VJ272vgtt9zinyAuCsbcwZjZNMntS8GY2TTJ7WvBmDsYM8O/eOOFF5w4cUI7d+5s8efU1FR/xTqvYMwdjJklcvtSMGaWyO1rwZg7GDPDv5iT5wVTp05tsd0wDL300ks+TtN+wZg7GDNL5PalYMwskdvXgjF3MGaGf1HkAQAAhCCGa71g6tSpLa5PZBiGXnzxRT8kap9gzB2MmSVy+1IwZpbI7WvBmDsYM8O/KPK8YPr06U4/l5WVadGiRQG/ZlEw5g7GzBK5fSkYM0vk9rVgzB2MmeFn/njaI1QdP37c/PGPf2x+97vfNdeuXevvOO0WjLmDMbNpktuXgjGzaZLb14IxdzBmhn9Q5HnB8ePHzezsbPP6668Pqj9wwZg7GDObJrl9KRgzmya5fS0YcwdjZvgXD154QWRkpPr3768xY8YoLMx5fek5c+b4KdX5BWPuYMwskduXgjGzRG5fC8bcwZgZ/sWcPC9YsWKF088NDQ0yTVPh4eF+StQ+wZg7GDNL5PalYMwskdvXgjF3MGaGf1HkecGmTZu0fPlyxcbGKi8vT2vWrFF0dLSGDx/u72htCsbcwZhZIrcvBWNmidy+Foy5gzEz/Mw/o8Sh5exXy/Tr18/88ssvTbvdbqampvovVDsEY+5gzGya5PalYMxsmuT2tWDMHYyZ4V/05HlBY2OjJOlf//qX+vXrp6ioKElS166B/fUGY+5gzCyR25eCMbNEbl8LxtzBmBn+xW+GF1xzzTWaPn26Dhw4oGnTpkk6807B06dP+zlZ24IxdzBmlsjtS8GYWSK3rwVj7mDMDP/i6VovME1TmzZtUlRUlNLS0iRJ5eXlKisr0/XXX+/ndK0LxtzBmFkity8FY2aJ3L4WjLmDMTP8iyIPAAAgBIWdfxcAAAAEG4o8AACAEESRBwAAEIIo8gAAAEIQRR4AAEAIosgDAAAIQRR5AAAAIYgiDwAAIARR5AEAAISg/w+yqDkq08QK9AAAAABJRU5ErkJggg==)

*ld heatmap*

![missingness plot](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAA2QAAAF+CAYAAADtOtkcAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjksIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvJkbTWQAAAAlwSFlzAAAQ6wAAEOsBUJTofAAAVupJREFUeJzt3Xd8jff///HnETJkCJUQYseqTcVobaWUFq1d4lO0VtWoVb4EtWq2tOrTGv10fKyido1GUbVaqsMoQkNTq5EYsZLr90c+zq+nSZDj5FzOyeN+u1239LyvcZ7XRfPyOtc4FsMwDAEAAAAAnC6b2QEAAAAAIKuiIQMAAAAAk9CQAQAAAIBJaMgAAAAAwCQ0ZAAAAABgEhoyAAAAADAJDRkAAAAAmISGDAAAAABMQkMGAAAAACahIQMAAAAAk9CQAQCyhMjISFksFodvd9u2bbJYLNq2bZvDtw0AcH80ZACQCS5cuKA333xTFStWlL+/v7y9vVW8eHF17dpVUVFRZsezWr9+vSIjI82OAQBAlmUxDMMwOwQAuJP9+/fr2WefVXx8vNq3b6/w8HB5e3srOjpaq1ev1k8//aT169erWbNmZkdVv3799N577ykrlII7d+7ozp078vb2duh2k5OTdevWLXl6eipbNj7nBABkTHazAwCAO7l8+bJatWqlbNmy6eDBgypTpozN/LfeeksrVqyQn5+fSQmzruzZsyt7dseXvWzZsjm8yQMAZB18lAcADvTBBx/o7NmzmjVrVqpm7K42bdqoTp06NmOxsbHq0aOH8ufPLy8vL5UtW1Zz5861WebuvUr//e9/NXHiRIWGhsrb21uNGjXS8ePHU73PF198oSeeeEI+Pj567LHH1LFjR/3+++/W+d26ddN7770nSbJYLNbp1KlTevLJJ1WxYsU081etWlU1atSQJJ06dUoWi0WTJ0/W7NmzVaxYMfn4+Kh27drat29fqnUfZD/TU7RoUT3zzDPatm2bdb/Kly+vr7/+WpK0atUqVaxYUd7e3qpataq+//57m/XTuofshx9+UPPmzRUUFCRvb28VKVJEXbp0UWJionWZpUuXqnr16goICJC/v7/Kli2r8ePHW+endQ9Zt27d5O3trbNnz6pVq1by8/NTUFCQ3njjDSUlJdlkuHTpkrp06aKAgAAFBgYqIiJCBw8elMVi0aJFi+zapmEYmj17tipUqCBvb28FBwere/fuunjxosP3HwDwcDhDBgAOtGbNGvn4+KhNmzYPvM758+dVs2ZNJSUlqU+fPgoODtbWrVvVp08fXbp0SaNGjbJZ/u2335aHh4feeOMNxcfH6+2331bnzp21Z88e6zKffvqpunTpomrVqmnSpEm6cOGC3n33Xe3cuVMHDhxQ3rx59eqrr+qPP/7Q5s2b9cknn1jXDQoKUkREhF599VUdOnTIpjE7fPiwDhw4oDlz5thk+vzzzxUXF6c+ffooOTlZ7733nho1aqQffvhBYWFhdu1nWk6ePKmOHTvqlVdeUZcuXTR16lQ999xz+ve//63hw4erT58+ypYtmyZOnKi2bdvqt99+k4eHR5rbunDhgp5++mnlzZtXQ4cOVe7cuRUTE6PVq1fr2rVr8vHx0ZYtW9ShQwc1bNhQkyZNkoeHh44ePaqdO3feN2tycrKeeeYZhYeHa9q0adqyZYumT5+uEiVKqHfv3tZlWrZsqT179qh3794qW7asVq9erYiICLu3KUm9e/fW/PnzFRERoX79+ikmJkazZ8/W3r17tW/fPnl7e2f6/gMAHpABAHCY3LlzG5UrV041npCQYFy4cME6XblyxTqvZ8+eRr58+Yzz58/brNOjRw/Dx8fHiIuLMwzDMKKiogxJRpkyZYybN29al3vnnXcMScZPP/1kGIZh3Lp1y8iXL59RtmxZ4/r169bl7q4/ePBg61jfvn2NtEpBXFyc4e3tbQwZMsRmfMSIEUaOHDmMCxcuGIZhGNHR0YYkw8vLy4iOjrYud/ToUSN79uxG586dM7yf6SlSpIghydi+fbt17KuvvrK+/8mTJ63j8+bNMyQZmzdvto6NGTPGZl9XrVplSDL27duX7nsOGDDACAgIMO7cuZPuMnePa1RUlHUsIiLCkGSMHTvWZtkqVaoY1apVs77+4osvDEnGtGnTrGNJSUlGw4YNDUnGwoULM7zNb7/91pBkfPzxxzbL7dixw5BkzJs3z6H7DwB4OFyyCAAOlJCQkOb9YT179lRQUJB16tevn6SUS8uWL1+uZ599VhaLRRcvXrROTZo0UWJios2ZL0nq2rWrPD09ra/vXv548uRJSSkPFTl37px69+4tHx8f63L169dXtWrVtG7duvvuR2BgoJ577jl9/vnnSk5Otmb9/PPP1axZM+XNm9dm+ZYtW6po0aLW16VKlVLTpk2t72XPfqalVKlSNpd73r10sn79+ipWrFiq8bvHJC25cuWSJK1du1a3b99Od5lr165p06ZN982Wlp49e9q8rlOnjk2mjRs3ysPDQ6+++qp1LFu2bOrbt6/d21y6dKn8/Pz0zDPP2BznMmXKKF++fNanfDpj/wEA90dDBgAO5O/vrytXrqQaHzVqlDZv3qzNmzfbPADiwoULiouL04IFC2watqCgILVr105SyqV+f1e4cGGb17lz55YkxcXFSZJOnz4tSSpdunSqHGXLltWpU6ceaF8iIiJ09uxZ6z/gd+zYodOnT6tLly6pli1ZsmSqsVKlSuny5cu6fPmyXfuZln/u+92molChQmmO3z0maalXr55efPFFjR07Vo899phatmypDz/8UNeuXbMu06dPH5UuXVrNmzdXwYIFFRERodWrVz/QUylz5MihkJAQm7HcuXPbZDp9+rTy5cuXqom/e5mnPds8duyYrl69qnz58qU61ufOnbMe58zefwDAg+EeMgBwoLJly+rgwYO6ffu2cuTIYR0vX768ypcvL0k29zTdPfvUsWNHvfzyy2lus1y5cjav07snytH/SG7atKny5cunTz/9VI0aNdKnn36qwMBAtWzZMsPbsmc/05LevttzTCwWi5YtW6a9e/dq7dq12rx5s1555RVNmjRJu3fvVnBwsIKDg3XgwAFt2bJFGzZs0MaNG/Wf//xHLVq00OrVq+/5RdOZ8Qj8B9lmcnKyHnvsMS1evDjN+Xcb+MzefwDAg6EhAwAHatmypb777jstX75cHTt2vO/yQUFB8vf31507d9S4cWOHZChSpIgk6ejRo2rSpInNvCNHjthcWnivf1B7eHioc+fO+vDDDzVz5kwtX75cbdu2lZeXV6plf/vtt1Rjx44dU2BgoAIDA5WUlOTw/XSU8PBwhYeHa9y4cdqwYYOaN2+uDz/8UCNHjpQkeXp6qnnz5mrevLkMw9CIESM0ZcoU7dq1S08++eRDvXeRIkW0detWXb161eYsWVpPzXxQJUqU0ObNm1WzZs0H+noFM/cfAMAliwDgUL169VJISIgGDRqkI0eO3Hd5Dw8Pvfjii1q1apV+/PHHVPMvXLiQ4QxPPPGE8uXLp3nz5unGjRvW8R07dmj//v1q0aKFdczX11dS+pf2RURE6MqVK3r11VcVFxenrl27prncmjVrbC6FPHbsmL766is1b9480/bzYcXFxaU6g1a1alVJKd8nJ6U8kv7vLBaLqlSpYrPMw2jatKmSkpI0b94869jdp1Taq3379kpOTta4ceNSzUtKSrL+WT8K+w8A4AwZADhU7ty5tWrVKj377LOqXLmyOnTooPDwcHl6eiomJkYrVqzQtWvXbO6Fmjx5srZt26ZatWqpZ8+eKleunOLi4nTw4EGtXLnSpql6EDly5NDUqVPVtWtX1alTRy+99JL1sfcFCxbUsGHDrMs+8cQTkqR+/fqpWbNmyp49u1q2bGlt1CpWrKhKlSpp6dKlKlasWLpnRO4+bKNv375KTk7WnDlz5O3trTFjxmTafj6sjz/+WO+9955at26tEiVKKDExUQsXLrQ2j5LUo0cPXbx4UY0aNVKhQoV09uxZzZkzRyEhIapbt+5DZ2jVqpXCw8M1dOhQRUdHWx97/9dff0m69xnM9NStW1d9+/bV1KlTdejQITVt2lReXl46fvy4li9frnHjxqlbt26PxP4DAGjIAMDhwsPD9csvv2jGjBlau3atli1bpqSkJBUoUEBPPvmk3nnnHTVs2NC6fHBwsPbs2aPx48dr1apVmjt3rvLkyaOyZctq+vTpdmXo0qWLcubMqUmTJmnYsGHKmTOnmjdvrilTptg8IbFNmzYaMGCA/vvf/+q///2vDMNQdHS0tSGTUs6SDRo0SC+99FK6DUKnTp3k6+urGTNmKDY2VlWqVNGsWbNUqlSpTN3Ph1GvXj3t27dPS5cu1Z9//qmAgABVqVJFs2fPtj6l8aWXXtJHH32kDz74QHFxccqXL5+aN2+uMWPGyN/f/6EzeHh4aN26dXr99df1n//8R9myZVPr1q01evRoPfXUUzYPgMmIOXPmqGrVqvrggw80cuRIZc+eXYULF1a7du2sf/cehf0HAEgWg0clAQDu4b333lO/fv109OhRmwZLkk6dOqVixYpp0qRJGj58uEkJ3c+qVavUunVr7dy5k/u0AMDNcQ8ZAOCePvroI9WqVStVMwbHSExMtHmdlJSk2bNnKyAgwHpPFwDAfXHJIgAglWvXrmn16tX65ptvdPDgQS1fvtzsSG7rtddeU2JiomrVqqWbN29qxYoV2rVrlyZOnGjzxd4AAPdEQwYASOXChQvq1KmTAgMDNXToUL3wwgtmR3JbDRs21PTp07V27VrduHFDYWFhmj17tvr162d2NACAE3APGQAAAACYhHvIAAAAAMAkNGQAAAAAYBIaMgAAAAAwCQ0ZAAAAAJiEhgwAAAAATEJDBgAAAAAmoSEDAAAAAJPQkAEAAACASWjIAAAAAMAkNGQAAAAAYBIaMgAAAAAwCQ0ZAAAAAJiEhgwAAAAATEJDBgAAAAAmoSEDAAAAAJPQkAEAAACASWjIAAAAAMAkNGQAAAAAYBIaMgAAAAAwCQ0ZAAAAAJiEhgwAAAAATEJDBgAAAAAmoSEDAAAAAJPQkAEAAACASWjIgEdM/fr19cYbbzz0drp166YXX3zRAYkAAFlRZGSknnjiCae/76JFi5Q3b96H3s62bdtksVh09epVB6QCMg8NGdzauXPn1KNHDxUqVEje3t4qWLCgWrZsqSNHjliXsVgs8vHx0dmzZ23W/WdjVL9+fVksFlksFnl7e+vxxx/X+++/7/DMK1as0JgxYx56O++8844++ugjByQCAGSWB6lTWU379u3166+/PvR2ateurdjYWPn6+jogFZB5spsdAMhMbdq0kWEY+uyzz1SkSBGdPXtWGzdu1OXLl1MtGxkZqQ8//PCe2+vdu7dGjx6tGzduaPXq1erXr59y586tjh07Oixznjx5HLKdXLlyOWQ7AIDMk5E6lVX4+PjIx8fnobfj6emp/PnzOyARkLk4Qwa3FRcXp127dmnq1KmqW7euihQpotq1a2vcuHGqWbOmzbJ9+/bVokWLdOzYsXtuM2fOnMqfP7+KFi2q/v37q1GjRlq9enWay546dUoWi0XLly9X7dq15ePjo1q1aun3339XVFSUKlSoIH9/f7Vv317Xrl2zrvfPM3Pvv/++wsLC5OXlpfz586tbt27WecuXL1f58uXl7e2tvHnzqmnTpkpOTpaU+pLF+vXra+DAgRo4cKACAwMVGhqq9957zybzt99+q4oVK8rb21s1a9bUsmXLZLFYdOrUKUn//zKSdevWqVSpUgoICNCLL76o+Ph46zaSk5M1ceJEFS1aVDlz5lTVqlW1bt06mz+XTp06KW/evPLx8VGZMmW0bNkySdKtW7fUt29f5c+fX97e3ipevLhmz559zz8TAHBVD1qn3njjDZUsWVI+Pj4qUaKEJkyYYP1dL/3/3/fjxo1TUFCQ8uTJo2nTpunmzZvq1auXAgICVKxYMa1du9a6zt3L+TZs2KAKFSrI29tbderU0W+//XbPzB9++KFKly5tvVJk4cKF91y+fv36GjBggPr166dcuXIpJCREn332meLi4tS2bVv5+fmpXLly2rNnj3Wdf16y+OOPP6p+/fry9/dXQECAwsPDdfjwYUnS6dOn1bJlSwUGBsrX11cVK1bUN998Y7OPdy9ZfJAaduXKFXXq1Em+vr4qWLCgPvjgAz3xxBOKjIy0LmOxWDR//ny1aNFCOXPmVPny5bVz506b/d6xY4eefPJJ+fj4qEiRIho6dKhu3rxpnW9vbYd7oiGD2/L395evr69WrVqlW7du3XPZ+vXrq1GjRho5cmSG3sPHx+e+2x49erTGjh2rffv26caNG+rQoYMmTJighQsX6quvvlJUVJRmzZqV5rr79+9X//79NX78eB07dkzr1q1TjRo1JEmxsbHq2LGjunfvriNHjujrr79Ws2bN7pllwYIFypcvn/bt26eBAwfqtddesxa1hIQEPffcc6pWrZoOHDigkSNHavjw4am2ceXKFc2ZM0dLly7Vxo0b9d1332ny5MnW+ZMmTdJnn32mDz/8UD///LN69eqlNm3a6MCBA5Kk//u//9Phw4e1ceNG/frrr5o1a5b1bN67776rNWvWaPny5Tp69KgWLFigAgUK3HOfAMBVPWidypUrl/7zn//o8OHDmjp1qmbMmJHqkvRNmzbp0qVL2rlzp8aPH68hQ4aoVatWKl++vL7//nu1bNlSXbt2TXU/1bBhw/TOO+9o79698vb2Vps2bdL9x/9nn32msWPHasqUKTp8+LAiIyP1+uuv68svv7znfi5YsEChoaHav3+//vWvf6l79+7q3LmzWrdurQMHDujxxx9X165dZRhGmut37txZhQoV0v79+7V//34NGDBA2bKl/BO2b9++unnzpnbs2KFDhw5p3Lhx8vb2TjfL/WrYoEGDtGfPHq1bt04bN27UmjVrdPTo0VTbGTdunLp166Yff/xRFStWVKdOnXT79m1J0okTJ9S8eXN16NBBP/30kz799FOtX7/eWlMdXdvhBgzAjS1ZssTIlSuX4ePjY9StW9cYM2aMcezYMZtlJBlr1qwxfvjhByNbtmzGvn37DMMwjHr16hmDBw+2Lvf310lJSca6desMT09PY86cOWm+d3R0tCHJWLRokXXsww8/NCQZBw4csI717dvXaNSoUZrv88UXXxgBAQFGQkJCqu1///33hiTj1KlTab5/RESE8cILL9hst0GDBjbLhISEGHPnzjUMwzDmzp1rBAcHGzdv3rTOnzlzpiHJiI6ONgzDMBYuXGhIMk6fPm1dZsSIEUaNGjUMwzCMGzduGDlz5jT27t1r8z7t27c3evfubRiGYbRs2dL417/+lWbm1157zWjYsKGRnJyc5nwAcDcPUqf+acyYMTZ1IyIiwihevLiRlJRkHStYsKDx/PPPW19fuHDBkGTs2LHDMAzDiIqKMiQZK1assC7zxx9/GDly5DC++uor6/tUq1bNOr9EiRLG0qVLbbIMGzbMaNasWbpZ69WrZ9SvX9/6+saNG4aHh4fx+uuvW8f27dtnSDJiYmIMw0ipNY899ph1vr+/v00t/bsKFSoYkZGRac67u49XrlyxbvdeNSwhIcHIkSOHsXLlSuv8c+fOGTly5DDGjBljHZNkjB071vr66NGjhiTj8OHDhmEYRvfu3Y0+ffrYZNmwYYORM2dOIzk5+aFqO9wTZ8jg1tq1a6c//vhDy5cvV926dbVq1SpVqFBBX331Vaplq1SporZt22rEiBHpbu/dd9+Vn5+fvL299cILL6hv377q3bv3PTNUrFjR+t/58uWTJJUrV85m7Pz582mu27hxY4WGhqp48eKKiIjQ4sWLrZc8VKpUSfXr11eFChXUvn17LViwQAkJCffMUqFCBZvXISEh1vc+evSoKlWqJE9PT+v88PDwVNsICAhQ4cKF09zG8ePHdf36dTVo0EB+fn7WacWKFTpx4oQk6ZVXXtF///tfValSRcOHD9e+ffus2+ratat++OEHlSlTRgMGDNDXX399z/0BAFf3IHVqyZIlql27tvLlyyc/Pz9NmjRJv//+u812ypcvbz1rJKXUlr/Xmrx588rDwyNVvbl7ZkZK+X1epEgR65UTf3ft2jWdOHFCERERNr/fZ86caf39np6/10EvLy/lypUrVR2UlG4tfP3119W9e3c9/fTTmjp1qk6fPm2d169fP7311lt66qmnNHbs2DSz/929atjJkyd1+/Ztm9oXHBysokWLptrO3+tpSEiITf4ff/xRH330kc1xatOmja5fv64///zT4bUdro+GDG4vZ86cat68ucaPH68DBw6ofv36mjBhQprLvvXWW9q2bVu6jUDXrl118OBBRUdH69q1a5oxY4ZNAUxLjhw5rP9tsVjSHEvv8pCAgAAdPHhQn332mfLmzathw4YpPDxciYmJ8vDw0NatW7Vu3TqFhYVp6tSpevzxx3Xu3LkHynK/97ZnG3cvhdm4caMOHjxonX799VctWLBAktSiRQudOnVKr732mk6dOqWnnnrK+ufxxBNPKDo6WpGRkUpISNDzzz+vnj17ZigfALiae9Wp7777Tp07d9Zzzz2n9evX68CBAxo4cGCqSxzT+t38zzFJdt+LdPf3+8KFC21+v//yyy/atGnTPde9X7a7tTG9bOPHj9fPP/+sZ555RuvXr1eZMmW0detWSSkf8h0/flwdO3bU/v37ValSJf3nP//JUBZ7jsm98l+9elV9+/a1OU6HDh3Sb7/9pqCgIIfXdrg+GjJkKRaLRaVKlbJ5iMbfhYWFqXv37mneOyVJgYGBCgsLU8GCBe/biDlKjhw51KRJE02fPl179+7VoUOHrGeVsmXLpjp16mjChAn68ccfdfXq1TTP/j2I0qVL68cff7Qp8n8/e/UgHn/8cXl6eiomJkZhYWE2U8GCBa3L5cuXTy+//LIWL16scePGaf78+dZ5gYGB6tixoxYsWKCPPvpIixYt4mZmAFnGP+vUrl27VKJECQ0fPlzVqlVTyZIlrQ9acoS/P0zjzz//1OnTp1W2bNlUy+XLl08hISE6efJkqt/vRYoUcVie9JQpU0aDBw9WVFSU6tevr88++8w6r0iRIurbt6/WrFmj7t27a9GiRXa9R/HixZUjRw6b2nfhwoUMH+8qVarol19+SXWcwsLClD17ygPOnVXb4Rp47D3c1qVLl9SuXTv16NFDFSpUkI+Pj7755hstWLBAQ4cOTXe90aNHKywsTIZhmPKFmH+3du1aRUdHq27dusqVK5eWLl2qHDlyqHjx4tqzZ4+2bt2qJk2aKCgoSNu3b9fVq1dVunRpu96rU6dOGjlypHr37q033nhDJ06csD7h8O6nf/fj7++vgQMH6vXXX9edO3dUu3ZtXb58Wdu3b1doaKjatm2rMWPGqFq1aipXrpyuXbumr776SmXKlJEkzZw5UwUKFFDlypVlGIZWrVqlkiVLOq35BQBnepA6VbJkSUVHR2vp0qWqVq2aVq5cqY0bNyowMNAhGSIjI5U7d27lyZNHQ4YMUcmSJdW4ceM0lx09erTeeOMN+fv7q0mTJrpx44Z2796tpKSk+16+b6/ExEQNHTpUbdu2VZEiRXTq1CkdOnRIDRo0kCQNGDBAzZs3V8mSJXXhwgV9++23euqpp+x6L39/f0VERGjw4MEKDAxUnjx5NGLECHl5eT1wHZSkoUOHqlatWnr99df18ssvy8fHRz///LO+++47TZ061am1Ha6Bhgxuy8/PT9WrV9fUqVN14sQJJScnq0iRItaCkp4CBQqof//+mjJlihPTpi0wMFDLly/XmDFjdPPmTZUpU0YrVqxQaGiorly5ou3bt2vmzJm6evWqihUrpg8//NDmfoCMCAgI0OrVq9W7d29VrlxZVapU0ejRoxUREXHPJ1b906RJkxQcHKy33npL0dHRyp07t6pXr67Ro0dLSvlUcNiwYTp16pR8fX3VoEEDvfPOO5IkX19fTZo0ScePH1f27NlVs2ZNffHFF3btDwA86h6kTj333HMaMGCA+vTpo1u3bun555/XsGHDNG/ePIdkmDhxovr27avjx4+revXqWrFiRbofgvXq1Uu+vr6aNm2aBg8eLH9/f1WqVOmeH3I+LA8PD128eFEvvfSSzp07p+DgYHXq1EmDBg2SJN25c0e9e/fW2bNnFRgYqBYtWjxU/Z4xY4ZeeeUVNWvWTLlz59b//d//KTo6OkN1sHLlyoqKitKoUaNUu3ZteXh4KCwsTBEREZKcW9vhGiyGkc4zRgFkefPmzdOQIUN0+fJlzlIBgBvZtm2bGjRooCtXrsjPz8/sOI+sy5cvq0CBAvrkk0/0wgsvmB0HboozZACsFi1apJIlS6pAgQL6/vvvNXr0aHXt2pVmDACQJXz//ff67bffVL16df31118aPXq0AgIC9Mwzz5gdDW6MhgyAVWxsrMaMGaM///xToaGh+te//qWxY8eaHQsAAKcwDENvv/22jh07Ji8vL4WHh2v79u3y9fU1OxrcGJcsAgAAAIBJuA4JAAAAAExCQwYAAAAAJqEhAwAAAACT0JABAAAAgEmy5FMWb9y4oZ9++klBQUHKnj1LHgIAeGTcuXNHFy5cUIUKFTL05auwRW0DgEdHRmpblvyN/dNPPyk8PNzsGACAv9m7d6+qV69udgyXRW0DgEfPg9S2LNmQBQUFSUo5QCEhISanAVzYRx+l/OzRw9wccGmxsbEKDw+3/m6GfahtgINQ2+AAGaltWbIhu3spR0hIiEJDQ01OA7iwLVtSfkZGmhoD7oHL7B4OtQ1wEGobHOhBahsP9QAAAAAAk9CQAQAAAIBJaMgAAAAAwCSPTEN2/Phx9erVS5UrV1b27NlVvnz5NJebP3++SpUqJW9vb1WqVElr1651clIAAO5v2bJlev755xUaGipfX19VrlxZCxYskGEYNstR1wAga3tkGrJffvlF69atU1hYmB5//PE0l1m8eLF69uyp9u3ba8OGDapVq5Zat26t3bt3OzktAEnSmjUpE4BUZsyYoZw5c2r69Olas2aNmjVrpp49e2rcuHHWZahrwCOI2gYnsxj//KjOJMnJycqWLaU/7Natm/bv36+ff/7ZZpnSpUurWrVq+vzzz61jtWvXVmBgoNavX//A73XmzBkVKlRIMTExPIkKAEzmrr+TL168qLx589qMvfLKK1qyZIni4uKULVs2h9U1yX2PIwC4ooz8Tn5kzpDdbcbSc/LkSR07dkzt2rWzGe/QoYO2bt2qmzdvZmY8AGn55JOUCUAq/2zGJKlKlSpKSEjQtWvXqGvAo4raBid7ZBqy+zly5IgkqUyZMjbjZcuW1a1btxQdHW1GLCBrmzcvZQLwQHbu3KmCBQvK39+fugY8qqhtcDKX+RbOuLg4SVJgYKDNeO7cuSVJf/31V7rrJiQkKCEhwfo6NjbW8QEBALiHnTt3avHixZo+fbqkh6trErUNANyFyzRkD2PGjBkaO3as2TGATNFyxJKHWn/NpPYOSgIgPWfOnFH79u3VoEED9e/f3yHbpLbBnZ3t8HC1qeDih6uNgDO5zCWLdz8xjI+Ptxm/+wljnjx50l130KBBiomJsU579+7NvKAAAPzN5cuX1axZMz322GP64osvrPdMP0xdk6htAOAuXOYM2d1r7I8cOaLSpUtbx48cOSJPT08VL1483XUDAgIUEBCQ6RkBAPi7xMREtWjRQvHx8fruu++UK1cu67yHqWsStQ0A3IXLnCErXry4SpUqpWXLltmML1myRI0aNZKnp6dJyYAsbNGilAlAKnfu3FG7du10+PBhbdy4UQULFrSZT10DHlHUNjjZI3OG7Pr169bvXDl9+rQSEhK0fPlySVK9evUUFBSkyMhIde7cWSVKlFCDBg20ZMkS7dmzR9u3bzczOpB1hYWZnQB4ZPXp00dr167V9OnTlZCQYPNlz1WqVJGXlxd1DXgUUdvgZI9MQ3b+/Hm1bdvWZuzu66ioKNWvX18dO3bU9evXNXnyZE2ePFmlS5fWypUrVatWLTMiA7j7xbXNm5ubA3gEbdq0SZI0ePDgVPOio6NVtGhR6hrwKKK2wckemYasaNGiMgzjvst1795d3bt3d0IiAPc1cWLKT4oWkMqpU6ceaDnqGvCIobbByVzmHjIAAAAAcDc0ZAAAAABgEhoyAAAAADAJDRkAAAAAmOSReagHABc0fbrZCQAAcCxqG5yMhgyA/WrUMDsBAACORW2Dk3HJIgD77dmTMgEA4C6obXAyzpABsN/dL7zdudPcHAAAOAq1DU7GGTIAAAAAMAkNGQAAAACYhIYMAAAAAExCQwYAAAAAJuGhHgDs9+abZicAAMCxqG1wMhoyAPZr3tzsBAAAOBa1DU7GJYsA7Hf8eMoEAIC7oLbByThDBsB+3bql/OS7WgAA7oLaBifjDBkAAAAAmISGDAAAAABMQkMGAAAAACahIQMAAAAAk/BQDwD2e/VVsxMAAOBY1DY4GQ0ZAPt16WJ2AgAAHIvaBifjkkUA9ouLS5kAAHAX1DY4GWfIANivZcuUn3xXCwDAXVDb4GScIQMAAAAAk9CQAQAAAIBJaMgAAAAAwCQ0ZAAAAABgEh7qAcB+L75odgIAAByL2gYnoyEDYL8BA8xOAACAY1Hb4GRcsggAAAAAJqEhA2C/evVSJgAA3AW1DU7GJYsA7JeUZHYCAAAci9oGJ+MMGQAAAACYhIYMAAAAAEzicg3Z6tWrVaNGDfn7+yskJETt2rXTyZMnzY4FAAAAABnmUg3Ztm3b1Lp1az3++ONauXKlZs2apR9//FFNmjRRYmKi2fGArKdBg5QJAAB3QW2Dk7nUQz0WL16sIkWKaMGCBbJYLJKk4OBgNWzYUPv371edOnVMTghkMePHm50AAADHorbByVzqDNnt27fl7+9vbcYkKVeuXJIkwzDMigUAAAAAdnGphqxbt2769ddf9f777ys+Pl4nT57Um2++qSpVqujJJ59Md72EhASdOXPGOsXGxjoxNeDGWrdOmQAAcBfUNjiZSzVkderU0cqVKzV8+HAFBgaqRIkSOnfunDZs2CAPD49015sxY4YKFSpkncLDw52YGnBjFy6kTAAAuAtqG5zMpRqyXbt2qUuXLurZs6e+/vprLVu2TMnJyXr22Wfv+VCPQYMGKSYmxjrt3bvXiakBAAAAIG0u9VCP/v37q2HDhpo+fbp1rGbNmipcuLA++eQTvfLKK2muFxAQoICAAGfFBAAAAIAH4lJnyH799VdVrlzZZiw0NFR58+bViRMnzAkFAAAAAHZyqTNkRYoU0Q8//GAzdvr0aV28eFFFixY1JxSQlVWoYHYCAAAci9oGJ3OphqxXr14aMGCAXn/9dbVs2VKXLl3SW2+9peDgYLVr187seEDWM3eu2QkAAHAsahuczKUasv79+8vLy0tz587V/Pnz5e/vr1q1amnZsmV67LHHzI4HAAAAABniUveQWSwW9erVSz/++KOuXr2q2NhYrVixQmXKlDE7GpA19eiRMgFI0/Hjx9WrVy9VrlxZ2bNnV/ny5VMtU79+fVksllTTkSNHTEgMgNoGZ3OpM2QAHjH8gxG4p19++UXr1q1TjRo1lJycrOTk5DSXe/LJJzVt2jSbMe6NBkxCbYOT0ZABAJBJWrZsqeeff16S1K1bN+3fvz/N5QIDA1WzZk1nRgMAPCJc6pJFAABcSbZslFkAwL1RKQAAMNk333wjX19feXt7q169etq+fft910lISNCZM2esU2xsrBOSAgAcjUsWAdivYEGzEwAur169euratatKliypP/74Q9OmTVPjxo31zTffqFatWumuN2PGDI0dO9aJSYEsgtoGJ6MhA2C/JUvMTgC4vH82VS1atFC5cuU0fvx4rV+/Pt31Bg0apB5/exJcbGyswsPDMy0nkGVQ2+BkNGQAADxCfH199eyzz2r58uX3XC4gIEABAQFOSgUAyCzcQwbAfsOGpUwAALgLahucjDNkAOz37bdmJwDczrVr17R27VpVr17d7ChA1kRtg5PRkAEAkEmuX79uvQ/s9OnTSkhIsF6KWK9ePR05ckRTp05V69atVbRoUf3xxx+aPn26/vzzTy1btszM6AAAJ6EhAwAgk5w/f15t27a1Gbv7OioqSqGhobp165befPNNXbp0Sb6+vqpdu7Y++OADHtABAFkEDRkAAJmkaNGiMgzjnsts3LjRSWkAAI8iGjIA9vP2NjsBAACORW2Dk9GQAbDfli1mJwAAwLGobXAyHnsPAAAAACahIQNgv7ffTpkAAHAX1DY4GQ0ZAPutXp0yAQDgLqhtcDIaMgAAAAAwCQ0ZAAAAAJiEhgwAAAAATEJDBgAAAAAm4XvIANjvq6/MTgAAgGNR2+BkNGQA7Ofra3YCAAAci9oGJ+OSRQD2W7AgZQIAwF1Q2+BkNGQA7EfRAgC4G2obnIyGDAAAAABMQkMGAAAAACahIQMAAAAAk9CQAQDwP6dPn9ZPP/1kfX3z5k1NmDBBL730khYtWmReMACA2+Kx9wDs9/nnZicAHKpnz56qXLmy3n77bUnSsGHD9P7776tixYpaunSprl+/rj59+picEkCmorbByThDBsB+hQunTICbOHjwoOrUqSNJunPnjj7++GNNmTJF+/fvV2RkpObOnWtyQgCZjtoGJ6MhA2C/L79MmQA3ceXKFeXKlUuStGfPHiUkJKhDhw6SpKeeekonT540Mx4AZ6C2wcloyADYb+rUlAlwE6Ghodq9e7ckacWKFXr88ccVEhIiSYqLi1POnDnNjAfAGahtcDLuIQMA4H+6d++uUaNGadmyZTpw4IBmzpxpnbd7926VLVvWxHQAAHdEQwYAwP8MHz5cBQoU0L59+9SnTx9169bNOi8uLk49evQwLxwAwC25ZEP28ccfa9asWTp8+LD8/PxUvXp1rVixQj4+PmZHAwC4uK5du6pr166pxj/44AMT0gAA3J3LNWQTJkzQlClT9Oabb6pWrVq6ePGitm7dqqSkJLOjAQDcxMaNG7Vv3z7FxMRo1KhRKly4sLZv366wsDAVKFDA7HgAADfiUg3Z0aNHFRkZqdWrV6tZs2bW8RdeeMHEVEAWNnu22QkAh7pw4YJatWql3bt3q1ChQoqJiVGvXr1UuHBhLViwQL6+vnrvvffMjgkgM1Hb4GQu1ZAtXLhQxYoVs2nGAJioShWzEwAONWDAAF24cEE///yzSpYsKU9PT+u8xo0b66233jIxHQCnoLbByVzqsfe7d+9WhQoV9NZbbyk4OFienp568skntWfPnnuul5CQoDNnzlin2NhYJyUG3NzOnSkT4CbWrVunCRMmqGzZsrJYLDbzChUqpDNnzpiUDIDTUNvgZC51huzPP//U999/r59++knvv/++cubMqYkTJ6pJkyb67bffFBwcnOZ6M2bM0NixY52cFsgChg9P+Unhgpu4c+eOfH1905wXFxdnc8YMgJuitsHJXOoMWXJysq5evarly5frxRdfVPPmzbV69WoZhqE5c+aku96gQYMUExNjnfbu3evE1AAAV1GjRg0tWLAgzXmLFy/Wk08+6eREAAB359AzZKdOndLx48dVtWpV5cmTx5GbliTlzp1bjz32mCpWrGgdy5Mnj6pUqaJffvkl3fUCAgIUEBDg8DwAAPfy1ltvqUGDBqpbt65efPFFWSwWrVq1SpMmTdK6deu0k0/MAQAOZvcZssGDB2vAgAHW1ytXrlTp0qXVpEkTlSxZUt9//70j8tkoV65cuvNu3Ljh8PcDAGQttWrVUlRUlCwWiwYPHizDMDRhwgTFxsZq69atqlq1qtkRAQBuxu6GbOXKlXriiSesr9988001b95chw4dUnh4uEaNGuWQgH/XokULXbp0SQcPHrSOXbp0ST/88IOqVavm8PcDAGQ9tWrV0jfffGN9INSVK1e0c+dO1apVy+xoAAA3ZPcli7GxsSpcuLAk6cSJEzp69Kg+/fRTlS9fXq+99poiIiIcFvKuVq1aqXr16nrxxRc1YcIE+fj4aNKkSfLy8lKfPn0c/n4A7mPMGLMTAJnGx8dHPj4+ZscA4GzUNjiZ3Q1Zrly5dP78eUnS5s2blSdPHutZKi8vLyUmJjom4d9ky5ZN69ev18CBA/Xqq6/q1q1bqlOnjrZv3678+fM7/P0A3MfTT5udAHho/fv3f+BlLRaL3nnnnUxMA8B01DY4md0NWd26dTV69GidO3dO06ZNU6tWrazzjh49aj175mh58+bVJ598kinbBpBBhw+n/Cxb1twcwENYs2bNAy9LQwZkAdQ2OJndDdnMmTPVpUsXDR8+XFWrVtWECROs8z755BPVqVPHIQEBPMJ69kz5yZPn4MKio6PNjgDgUUJtg5PZ3ZAVLFhQX3/9dZrzvvrqK3l7e9sdCgAAAACyggw/ZXHz5s167rnnVKFCBT3zzDOaP39+qmUCAgLk6enpkIAAADjTxYsXNXLkSDVu3FjlypVT48aNNWrUKF24cMHsaAAAN5Shhmz16tVq2rSpduzYIV9fX/3000965ZVXFBkZmUnxAABwnj179qhkyZKaM2eOcuXKpXr16ilXrlyaPXu2wsLCtGfPHrMjAgDcTIYuWZw8ebIaNGigL7/8Un5+fkpOTtZrr72madOmafTo0cqWze6vNQMAwHR9+/ZVuXLltH79egUEBFjH4+Pj1axZM/Xr10/79u0zMSEAwN1kqCE7fPiwPv30U/n5+UlKeQz9qFGjNHfuXJ0+fVrFihXLlJAAHlH9+pmdAHCoX375RcuWLbNpxqSUr3oZPny42rdvb1IyAE5DbYOTZaghi4+PV1BQkM1YcHCwJCkuLo6GDMhqOnQwOwHgUGFhYbp8+XKa8+Lj41W8eHHnBgLgfNQ2OFmGn7KYkJCgv/76y/r6zp07aY5LUp48eR4yHoBH2v++HF7/+2AGcHVTp05V3759VahQIdWrV886vm3bNkVGRmrOnDkmpgPgFNQ2OFmGG7KmTZumOd6oUaNUY0lJSRlPBMB1tGmT8pPvaoGbGDJkiOLj49WwYUPlypVLQUFBunDhguLj45U7d24NGzZMw4YNk5TyJdE//vijyYkBOBy1DU6WoYZs4cKFmZUDAADTVatWTRaLxewYAIAsJEMNWURERGblAADAdIsWLTI7AgAgi+E59QAAZJLjx4+rV69eqly5srJnz67y5cunudz8+fNVqlQpeXt7q1KlSlq7dq2TkwIAzJKhM2Qvv/zyAy9rsVg0f/78DAcCAMBMu3fv1rJlyxQTE6MbN27YzLNYLPryyy8feFu//PKL1q1bpxo1aig5OVnJycmpllm8eLF69uypkSNHqmHDhlqyZIlat26tHTt2qGbNmg+9PwCAR1uGGrIDBw7cd5lbt27p8OHDNGRAVtCxo9kJAId65513NHDgQAUHB6tEiRLy9PR8qO21bNlSzz//vCSpW7du2r9/f6plxowZow4dOmj8+PGSpAYNGujQoUMaN26c1q9f/1DvD8AO1DY4mcMashs3bujf//633n77bXl4eKgjf5kB99e3r9kJAIeaNm2a+vXrp1mzZilbtoe/qv9+2zh58qSOHTumKVOm2Ix36NBBQ4YM0c2bN+Xl5fXQOQBkALUNTpbhx97/U2JioubOnatp06bpr7/+0ksvvaQRI0aoRIkSjsgH4FH2v+8hVPaH/lUCPBKuXbum559/3iHN2IM4cuSIJKlMmTI242XLltWtW7cUHR2dat5dCQkJSkhIsL6OjY3NvKBAVkJtg5PZXXGuXbumKVOmqGjRonrzzTf1/PPP69ixY/roo49oxoCson79lAlwE+3bt9eGDRuc9n5xcXGSpMDAQJvx3LlzS5L++uuvdNedMWOGChUqZJ3Cw8MzLSeQpVDb4GQZbv2vXLmid999V7NmzdK1a9fUs2dPDR06VAULFsyMfAAAOM2sWbPUo0cPderUSY0bN07VKElSm7tfGmuyQYMGqUePHtbXsbGxNGUA4IIy1JBFRkbq3Xff1Z07d9SrVy+98cYbCg4OzqxsAAA41ZEjR/Ttt9/q1KlTWrx4car5FotFSUlJDnu/u2fC4uPjlT9/fuv43TNnefLkSXfdgIAABQQEOCwLAMAcGWrIxo0bJ4vFovDwcB05csTmk7l/yuijgQEAMNvLL78sPz8/rVmzRqVKlXropyzez937w44cOaLSpUtbx48cOSJPT08VL148U98fAGC+DDVkdevWlcVikZRy6SIAAO7k8OHDWrFihZ555hmnvF/x4sVVqlQpLVu2zPp4fElasmSJGjVqlOkNIQDAfBlqyLZt25ZJMQC4pCZNzE4AOFTlypV17tw5h23v+vXr1u8SO336tBISErR8+XJJUr169RQUFKTIyEh17txZJUqUUIMGDbRkyRLt2bNH27dvd1gOABlAbYOT8TxPAPYbPdrsBIBDvf/+++rWrZtCQkLUsGFDZX/Ix16fP39ebdu2tRm7+zoqKkr169dXx44ddf36dU2ePFmTJ09W6dKltXLlStWqVeuh3huAnahtcLIMVZqLFy/qjz/+UMWKFW3GDx06pHHjxunw4cPKnz+/BgwYoJYtWzo0KAAAma1OnTq6ffu2mjVrpmzZssnHx8dmvsViUXx8/ANvr2jRojIM477Lde/eXd27d89wXgCA68tQQzZixAh9//33+uGHH6xjp0+fVp06dXT9+nVVqlRJP//8s1q3bq2vv/5adevWdXhgAI+QFi1Sfq5da24OwEEGDx5svVcaQBZFbYOTZagh+/bbb1N9gjdz5kxdvXpVGzZsUJMmTZSYmKinn35aU6ZMoSED3N3ly2YnABwqMjLS7AgAzEZtg5Nly8jCZ8+eVfny5W3G1qxZo8qVK6vJ/26A9PHxUb9+/XTo0CHHpQQAAAAAN5ShM2QWi8XmUo5z584pOjpaAwYMsFkuNDRUFy9edEhAAACc6fjx41q0aJGOHTumGzdupJq/evVqE1IBANxVhhqy0qVLa8uWLdazYWvXrpXFYrG+vis2NlZBQUGOSwkAgBPs27dP9erVU5EiRXTs2DFVrFhR8fHxOnXqlEJDQxUWFmZ2RACAm8lQQ9a/f3917dpVcXFxyp8/v+bOnauwsDA1btzYZrmvvvpKFSpUcGhQAI+gJ54wOwHgUEOHDlW7du00f/585ciRQ/Pnz1fVqlW1a9cudezYUcOGDTM7IoDMRm2Dk2WoIevcubPOnj2r2bNnKy4uTtWqVdP7779v8z0t58+f15o1azR27FiHhwXwiJk1y+wEgEP9+OOPGj58uLJlS7nF+u4li7Vr11ZkZKSGDx+upk2bmhkRQGajtsHJMvyNl0OHDtXQoUPTnR8cHKxz5849VCgAAMxgsVjk6ekpi8Wi4OBgnT59WrVr15aUcn/0sWPHTE4IAHA3GXrKIgDYiIhImQA38fjjj+vEiROSpFq1amn69On6+eefdfToUU2ePFklSpQwOSGATEdtg5Nl+AwZAFj97x+ugLt45ZVXdPr0aUnSxIkT1aRJE1WqVEmS5Ovrq+XLl5sZD4AzUNvgZDRkAAD8T5cuXaz/XbZsWR0+fFi7du3SjRs3VLNmTQUHB5uYDgDgjly2Ibt69arKlCmjs2fPat++fXqCJ+IAABzMz88v1Ve7AADgSC57D9n48eN1584ds2MAAFzcxYsXdejQoVTjhw4d0osvvqhy5cqpcePGWrNmjQnpAADuziUbsiNHjui9997j0fqA2YoWTZkAFzZixAh169bNZuz06dOqU6eOvvzyS/n4+Oinn35S69attX37dnNCAnAeahuczCUvWXzttdfUq1cvlS5d2uwoQNb26admJwAe2rfffqvu3bvbjM2cOVNXr17Vhg0b1KRJEyUmJurpp5/WlClTVLduXZOSAnAKahuczOXOkC1fvlw//fSTRo8e/cDrJCQk6MyZM9YpNjY2ExMCAFzJ2bNnVb58eZuxNWvWqHLlytb7x3x8fNSvX780L20EAOBhuFRDdv36dQ0aNEgTJ05UQEDAA683Y8YMFSpUyDqFh4dnYkogCxk0KGUCXJjFYpHFYrG+PnfunKKjo1WvXj2b5UJDQ3Xx4kVnxwPgbNQ2OJlLNWRvvfWW8uXLp3/9618ZWm/QoEGKiYmxTnv37s2khEAWs3dvygS4sNKlS2vLli3W12vXrpXFYkn1dMXY2FgFBQU5Ox4AZ6O2wclc5h6y06dPa/r06Vq5cqXi4+MlpTz6/u7Pq1evys/PL811AwICMnRGDQCQdfTv319du3ZVXFyc8ufPr7lz5yosLEyNGze2We6rr75ShQoVTEoJAHBXLtOQRUdH69atW3r22WdTzWvQoIFq1Kih3bt3m5AMAODKOnfurLNnz2r27NmKi4tTtWrV9P777yt79v9fIs+fP681a9bwdF8AgMO5TENWuXJlRUVF2YwdPHhQAwcO1AcffKDq1aublAwA4OqGDh2qoUOHpjs/ODhY586dc2IiAEBW4TINWWBgoOrXr5/mvGrVqqlq1arODQRA8vc3OwEAAI5FbYOTuUxDBuARtGGD2QkAAHAsahuczKUbsvr168swDLNjAAAAAIBdXOqx9wAeMRMmpEwAALgLahucjIYMgP02bODSDgCAe6G2wcloyAAAAADAJDRkAAAAAGASGjIAAAAAMAkNGQAAAACYxKUfew/AZF9/bXYCAAAci9oGJ6MhA2A/T0+zEwAA4FjUNjgZlywCsN8HH6RMAAC4C2obnIyGDID9Pv00ZQIAwF1Q2+BkNGQAAAAAYBIaMgAAAAAwCQ0ZAAAAAJiEhgwAAAAATMJj7wHYb9kysxMAAOBY1DY4GQ0ZAPuFhJidAAAAx6K2wcm4ZBGA/ZYt45NEAIB7obbByThDBsB+77yT8rNtW3NzAADgKNQ2OBlnyAAAAADAJDRkAAAAAGASGjIAAAAAMAkNGQAAAACYhIYMgP0++CBlAmC3RYsWyWKxpJqGDx9udjQga6K2wcl4yiIA+5Uvb3YCwG1s3LhRuXLlsr4uWLCgiWmALIzaBiejIQNgv6+/TvnZsKG5OQA3UK1aNeXNm9fsGACobXAyGjIA9hs9OuUnRQsA4C6obXAy7iEDAOARUK5cOXl4eKh48eKaNGmSkpKSzI4EAHACzpABAGCikJAQjR07VjVq1JDFYtHq1as1atQonT17VnPmzEl3vYSEBCUkJFhfx8bGOiMuAMDBaMgAADBR06ZN1bRpU+vrJk2ayMfHRzNnztTIkSMVEhKS5nozZszQ2LFjnRUTAJBJuGQRAIBHTLt27ZSUlKSDBw+mu8ygQYMUExNjnfbu3eu8gAAAh+EMGQD7TZhgdgIgywoICFBAQIDZMQD3Q22Dk9GQAbBfvXpmJwDc0uLFi+Xh4aEqVaqYHQXIeqhtcDIaMgD2O3Qo5WfFiubmAFxY06ZN1bBhQ1WoUEGStHr1av373//W66+/rvz585ucDsiCqG1wMhoyAPbr0yfl586d5uYAXFiZMmU0f/58nTlzRsnJySpVqpRmzZql1157zexoQNZEbYOTudRDPZYtW6bnn39eoaGh8vX1VeXKlbVgwQIZhmF2NAAA7PLOO+/o2LFjun79um7cuKFDhw6pf//+slgsZkcDADiBS50hmzFjhooWLarp06crKChImzdvVs+ePRUTE6MxY8aYHQ8AAAAAMsSlGrI1a9Yob9681tcNGzbUpUuXNGPGDP3f//2fsmVzqRN+AAAAALI4l+pg/t6M3VWlShUlJCTo2rVrJiQCAAAAAPu51BmytOzcuVMFCxaUv79/usskJCQoISHB+jo2NtYZ0QD3N2iQ2QkAAHAsahuczKUbsp07d2rx4sWaPn36PZebMWOGxo4d66RUQBbSpo3ZCQAAcCxqG5zMpS5Z/LszZ86offv2atCggfr373/PZQcNGqSYmBjrtHfvXielBNzcmTMpEwAA7oLaBidzyTNkly9fVrNmzfTYY4/piy++uO/DPAICAhQQEOCkdEAW0qFDyk++qwUA4C6obXAyl2vIEhMT1aJFC8XHx+u7775Trly5zI4EAAAAAHZxqYbszp07ateunQ4fPqwdO3aoYMGCZkcCAAAAALu5VEPWp08frV27VtOnT1dCQoJ2795tnVelShV5eXmZmA4AAAAAMsalGrJNmzZJkgYPHpxqXnR0tIoWLerkRAAAAABgP5dqyE6dOmV2BAB/17Wr2QkAAHAsahuczKUaMgCPmFdeMTsBAACORW2Dk7ns95ABeAQkJqZMAAC4C2obnIwzZADs9/TTKT/5rhYAgLugtsHJOEMGAAAAACahIQMAAAAAk9CQAQAAAIBJaMgAAAAAwCQ81AOA/Vq0MDsBAACORW2Dk9GQAbDf8OFmJwAAwLGobXAyLlkEAAAAAJPQkAGwX5MmKRMAAO6C2gYn45JFAPa7ft3sBAAAOBa1DU7GGTIAAAAAMAkNGQAAAACYhIYMAAAAAEzCPWQA7FerltkJAABwLGobnIyGDID9pk41OwEAAI5FbYOTcckiAAAAAJiEhgyA/Tp2TJkAAHAX1DY4GZcsArBfTIzZCQAAcCxqG5yMM2QAAAAAYBIaMgAAAAAwCQ0ZAAAAAJiEe8gA2K9UKbMTAADgWNQ2OBkNGQD7LVhgdgIAAByL2gYn45JFAAAAADAJDRkA+732WsoEAIC7oLbBybhkEYD9DhwwOwEAAI5FbYOTcYYMAAAAAExCQwYAAAAAJqEhAwAAAACTcA8ZAPs99pjZCQAAcCxqG5yMhgyA/b780uwEAAA4FrUNTsYliwAAAABgEhoyAPaLjEyZAABwF9Q2OJnLNWRHjhzR008/LV9fX+XPn19Dhw7VrVu3zI4FZE1btqRMAB4KtQ14hFDb4GQudQ9ZXFycGjZsqJIlS2rFihU6e/asBg0apOvXr2vOnDlmxwMAIMOobQCQtblUQ/bBBx8oISFBK1euVJ48eSRJd+7cUZ8+ffTmm2+qQIECJicEACBjqG0AkLW51CWLGzZsUOPGja0FS5LatWun5ORkbdq0ycRkAADYh9oGAFmbS50hO3LkiF5++WWbscDAQIWEhOjIkSPprpeQkKCEhATr65iYGElSbGxs5gQFnCgx4dJDrX/mzBn7V7558+5GHioDsra7v4vv3LljchJzUNuA1GKvX3+o9Q1qG0yWkdrmUg1ZXFycAgMDU43nzp1bf/31V7rrzZgxQ2PHjk01Hh4e7sh4gEsq9H5fB2yk0MNvA1nehQsXVLRoUbNjOB21DcgEa9Y+/DaobXCAB6ltLtWQ2WvQoEHq0aOH9fWNGzcUExOjYsWKKXt21zwEsbGxCg8P1969exUSEmJ2HFNwDDgGWX3/Jfc4Bnfu3NGFCxdUoUIFs6O4FGqbe+IYcAyy+v5L7nEMMlLbXOo3du7cuRUfH59qPC4uzuba+38KCAhQQECAzVhYWJjD85khJCREoaGhZscwFceAY5DV919y/WOQFc+M3UVtS83V/z47AseAY5DV919y/WPwoLXNpR7qUaZMmVTX08fHxys2NlZlypQxKRUAAPajtgFA1uZSDVmzZs20ZcsWXb582Tq2bNkyZcuWTU2aNDEvGAAAdqK2AUDW5lINWa9eveTv769WrVpp06ZNWrhwoYYMGaJevXplue9pCQgI0JgxY1JdrpKVcAw4Bll9/yWOgTugtv1//H3mGEgcg6y+/1LWOwYWwzAMs0NkxOHDh/Xaa69p165d8vf3V9euXTVhwgR5enqaHQ0AALtQ2wAg63K5hgwAAAAA3IVLXbIIAAAAAO6EhgwAAAAATEJDBgAAAAAmoSEDAAAAAJPQkLmQNWvWqFKlSvL29lapUqW0cOHCDG+jVatWslgsmjZtWiYkzHz2HIN9+/bp5ZdfVlhYmHLmzKmSJUtqxIgRunbtmhMS2+fIkSN6+umn5evrq/z582vo0KG6devWfdczDEOTJ09W4cKF5ePjo1q1amn37t1OSOx49hyD2NhYDR06VJUrV5a/v79CQ0PVqVMnnT592kmpHcvevwd/N2vWLFksFrVo0SKTUgIPh9pGbbsfahu17Z/crrYZcAk7duwwPDw8jFdffdX4+uuvjVGjRhkWi8VYtmzZA29j/fr1Rr58+QxJxtSpUzMxbeaw9xgMHjzYeOqpp4x58+YZUVFRxpw5c4w8efIYDRo0cFLyjPnrr7+MkJAQo27dusbGjRuN+fPnG7ly5TL69u1733UnTZpkeHp6GjNmzDC2bNlitG7d2vD39zdOnDjhhOSOY+8xWLNmjVGiRAljwoQJxtatW40lS5YY5cuXN4KDg43z5887Kb1jPMzfg7tiY2ONwMBAIzg42Hj22WczMS1gH2obtY3aRm2jthkGDZmLaNKkiVG7dm2bsY4dOxply5Z9oPVv3LhhhIWFGQsWLHDZomXvMUjrl9Vnn31mSDL279/v0IyOMHHiRMPX19e4dOmSdWzevHmGh4eHcfbs2XTXS0xMNAICAowRI0ZYx27evGkUKVLE6N27d6ZmdjR7j0FcXJxx+/Ztm7GYmBjDYrEY06ZNy7S8mcHeY/B3Xbp0Mbp27WrUq1fPbYoW3Au1jdpGbaO2UdsMg0sWXcDNmzcVFRWltm3b2ox36NBBhw8f1qlTp+67jWnTpil37tzq1q1b5oTMZA9zDIKCglKNValSRZL0xx9/ODSnI2zYsEGNGzdWnjx5rGPt2rVTcnKyNm3alO56u3btUkJCgtq1a2cd8/T0VJs2bbR+/fpMzexo9h6DwMBAZc+e3WYsNDRUQUFBj+Sf9b3Yewzu2rlzp1atWqXJkydnZkzAbtQ2ahu1jdomUdsk7iFzCSdOnNDt27dVpkwZm/GyZctKSrkW915+//13TZo0Se+++64sFkum5cxMD3sM/mnnzp2SlGp7j4IjR46kyhUYGKiQkJB77ufdeWkdo99//12JiYmOD5tJ7D0GaTl27JjOnz9v/bviKh7mGCQlJalfv34aOXKkQkJCMjMmYDdqG7WN2kZtk6htEg2ZS4iLi5OU8hf273Lnzi1J+uuvv+65/sCBA9WmTRvVrFkzU/I5w8Meg7+7ePGiIiMj9fzzz6tkyZIOy+gocXFxqfZTStnXe+1nXFycvLy85O3tnWo9wzCsx9AV2HsM/skwDPXv318FChRQx44dHZgw8z3MMXj//fd17do1DRw4MJPSAQ+P2kZtk6htErVNorZlv/8iyAzx8fGKjY2973LFixd/qPfZtGmTNm3apKNHjz7UdjKDs47B392+fVsdOnSQJM2dO9dh28WjKTIyUlu3btXGjRvl6+trdhynOH/+vEaPHq3//Oc/8vT0NDsOshhqG7UNmY/a5n61jYbMJMuWLVPPnj3vu9zhw4etn5TFx8fbzLv7qdDfr8P9p/79+6t///7KmTOnLl++bB2/ceOGLl++nOanFM7irGNwl2EYevnll7V3717t2LHjkT3dnTt37lT7KaXs6732M3fu3Lp586Zu3Lhh80liXFycLBaL9Ri6AnuPwd99+OGHGjdunObPn69GjRo5OmKms/cYjB49WhUrVlSdOnWs/8/fuXNHd+7c0eXLl+Xn55fqXgTAUaht1Lb0UNuobRK1LV1mPlEED+bGjRtGjhw5jFmzZtmMr1692pBkREdHp7uupHtOiYmJmZzeMR7mGNw1aNAgw9PT09iyZUsmpXSMOnXqGK1atbIZu3z5smGxWIyFCxemu97WrVsNScbBgwdtxgcNGmQUKVIkE5JmHnuPwV0rVqwwPDw8jHHjxmVSwsxn7zGoV6/ePf+f37BhQyYnBx4MtY3aRm2jthkGtc0weOy9y2jSpInx1FNP2Yx17tz5vo/FjYqKSjVJMnr16mVERUUZSUlJmRnboew9BoaR8h0mFovF+PzzzzMrnsNMnDjR8PPzM+Li4qxjH3744QM/GnjkyJHWsVu3bhlFixZ1yUcD23MMDCPl77yXl5fRq1evTE6Zuew9BgcOHEj1/3ylSpWMmjVrGlFRUTaPGgbMRm2jtlHbqG3UNhoyl3H3iyN79+5tREVFGaNHjzYsFouxdOlSm+U8PDyMl19++Z7bkot+V4u9x+Du97K89NJLxnfffWczPYpfqHj3SxPr1atnfPXVV8aCBQuMwMDAVF+a2LBhQ6NEiRI2Y5MmTTK8vLyMWbNmGVu3bjVeeOEFl/7yzIweg19//dXIlSuXUb58eePbb7+1+bM+fvy4s3fjoTzM34N/cqfvaoF7obZR26ht1DZqGw2ZS/nyyy+NChUqGJ6enkZYWJgxf/78VMtIMiIiIu65HVctWoZh3zGIiIhI9xT3g1wiYIZff/3VaNSokeHj42MEBwcbb7zxhnHz5k2bZerVq5fqco3k5GRj4sSJRmhoqOHl5WXUqFHD2LVrlxOTO449x2DhwoXp/lnf7/+LR5G9fw/+yZ2KFtwPtY3a9nfUNmqbYWS92mYxDMOw/w40AAAAAIC9+B4yAAAAADAJDRkAAAAAmISGDAAAAABMQkMGAAAAACahIQMAAAAAk9CQAQAAAIBJaMgAAAAAwCQ0ZAAAAABgEhoyuLzIyEhZLBYVLFhQycnJqeY/+eSTslgs6tatm806fn5+Dsuwbds2WSwW7d+/32HbzEpmzZql9evXmx0DAB4Z1DbXR23Dg6Ihg1vIkSOHLl68qO3bt9uMnz59Wt99912qAtWjRw9FRUU57P2rVq2q7777TmXLlnXYNrMSihYApEZtc23UNjyo7GYHABzB09NTjRs31n//+1/Vr1/fOr548WKVK1dOHh4eNsuHhoYqNDTUYe8fEBCgmjVrOmx7rs4wDN26dUteXl5mRwEAl0Vte7RQ25BZOEMGt9GxY0ctX75ct2/fto59/vnn6tSpU6pl/3lZx+3btzVkyBAVLlxYXl5eCgkJUcuWLRUfH/9A89O6rMNisejtt99WZGSk8uXLp7x58+pf//qXrl27ZpNl586dqlKliry9vVWxYkVt3rxZlStXtrkMpVu3bipfvry2bdumKlWqyNfXV+Hh4fr+++9ttmUYhqZNm6ZSpUrJy8tLxYsX18yZM22WOXPmjNq1a6d8+fLJ29tbxYoV08CBAx94flru5lu/fr0qVaokLy8vrVmzRteuXVO/fv1UunRp5cyZU0WLFlWvXr2sx02SihYtqtOnT+u9996TxWKRxWLRokWLrPMXLVqkihUrytvbWwULFtTIkSOVlJR0zzwA4C6obdQ2uD/OkMFttGzZUt27d9emTZv07LPP6tdff9WhQ4e0atUqLVmy5J7rTpo0SR988IGmTJmicuXK6eLFi9q0aZNu3rz5QPPTM2fOHNWpU0cff/yxjh07piFDhihfvnyaPHmyJCk2NlbPPPOMqlatqqVLlyo+Pl69e/dWfHy8KleubLOtP//8U/3799fw4cOVK1cujRgxQq1bt9aJEyeUI0cOSdLrr7+ujz76SCNHjlSNGjW0a9cuDRs2TD4+PurVq5ckqWvXrvrjjz/07rvvKl++fPr9999tiu395qfnjz/+UP/+/TVq1CgVLlxYhQsX1vXr15WUlKQJEyYoKChIMTExmjBhglq1amW9rGblypVq3ry5nnrqKQ0ePFiSVKJECUnSjBkzNHToUA0cOFDTp0/X4cOHrUXr7jEEAHdGbaO2IQswABc3ZswYw9fX1zAMw+jUqZPx0ksvGYZhGKNGjTJq1aplGIZhVKpUyYiIiEhzHcMwjGeffdZo06ZNuu9xv/lRUVGGJGPfvn3WMUlGeHi4zXIRERFGiRIlrK+HDBli5MqVy0hISLCO7dixw5BkkzciIsKwWCzGzz//nOo9d+zYYRiGYRw/ftywWCzGvHnzbN5z2LBhRv78+Y2kpCTDMAzD19fXePfdd9Pdl/vNT0tERIQhydi9e/c9l7t9+7axc+dOQ5Jx9OhR63iRIkWMvn372iybkJBg+Pn5GSNGjLAZnzt3ruHj42NcvHgxQxkBwJVQ26htyDq4ZBFupWPHjvryyy+VmJioxYsXq2PHjg+0XtWqVbV+/XpFRkZq3759qZ5odb/56Xn66adtXj/++OM6c+aM9fW+ffvUoEED+fv7W8eeeuop5cmTJ9W2ChQooHLlytlsS5J1e1u2bJEkvfDCC7pz5451aty4sf7880/FxMRY92XatGmaO3eujh8/nuaxuNf89Dz22GOqUaNGqvFPPvlEVapUkZ+fn3LkyKGnnnpKknTs2LF7bm/Xrl26evWq2rZtm2p/EhMT9fPPPz9wNgBwZdQ2ahvcGw0Z3ErTpk2VI0cOjR49WtHR0WrXrt0DrTdy5EgNGzZMH3/8scLDw5U/f36NHTtWhmE80Pz0BAYG2rz29PS0uRQkNjZWQUFBqdYLDg5+oG1J0o0bNyRJFy9elGEYyps3r3LkyGGd7hbOu0VryZIlatSokUaOHKmSJUuqTJkyWrFihXW795ufnnz58qUaW7lypbp27arw8HAtXbpUu3fv1sqVK21yp+fixYuSUoro3/enZMmSNvsDAO6O2kZtg3ujIYNbyZEjh1544QXNmDFDDRs2TPMXaVq8vLwUGRmp6Oho/fbbb+rRo4ciIyP16aefPtB8e4WEhOjChQupxs+fP5/hbeXJk0cWi0Xffvut9u3bl2qqVKmS9T0XLFigixcvau/evSpdurTat2+vkydPPtD89FgsllRjy5YtU+XKlTVv3jw1b95cNWrUUO7cuR94fyRpxYoVae5Ps2bNMnJ4AMBlUduobXBvPNQDbqdHjx46f/68evbsadf6YWFhmjhxoubNm6fDhw9neH5GVK9eXfPmzdOVK1esl3bs2LFDf/31V4a31ahRI0nSpUuX1LJly/suny1bNlWvXl1vvfWWVq9erePHj6t48eIPPP9BJCYmWj/tvOuzzz5LtZynp2eqTxVr1aqlnDlz6syZM2rdunWG3hcA3A21jdoG90VDBrcTHh6uVatWZWidVq1aqVq1atbH7q5Zs0ZxcXFq2LDhA82318CBA/X+++/r2Wef1ZAhQ3T58mWNHTtWefPmVbZsGTuBXapUKfXt21ddunTRkCFDVKNGDd2+fVvHjh1TVFSUVq1apfj4eDVt2lRdunRR6dKldevWLc2ePVuBgYGqWrXqfedn1NNPP62+fftq/PjxqlWrltavX6+tW7emWq5s2bL6+uuvtXnzZuXOnVvFihXTY489pnHjxmno0KE6c+aM6tevLw8PD508eVJffvmlvvjiC+XMmTPDmQDAFVHbqG1wXzRkgKQnn3xSS5cu1fTp03Xnzh2VLl1an332mRo3bvxA8+0VEhKiDRs2qH///nrxxRdVokQJvfPOO+rXr59y5cqV4e29++67Kl26tObNm6dx48bJz89PpUuXVtu2bSVJ3t7eqlChgmbPnq3ff/9dPj4+euKJJ7Rp0yblzZtXN2/evOf8jHr11Vd18uRJzZ49W1OnTlXTpk31+eefp/qi0YkTJ6p379564YUXdOXKFS1cuFDdunXT4MGDVbBgQc2YMUOzZ89Wjhw5VKJECbVo0SLVp5MAAFvUNmobXIPFuN+dmwCc6rffflOZMmW0YMECRUREmB0HAICHRm0D0kdDBphsxIgRqlixogoUKKCTJ09q4sSJSkxM1JEjR+Tn52d2PAAAMozaBjw4LlkETHbr1i0NGzZM586dk4+Pj+rXr6+pU6dSsAAALovaBjw4zpABAAAAgEn4HjIAAAAAMAkNGQAAAACYhIYMAAAAAExCQwYAAAAAJqEhAwAAAACT0JABAAAAgEloyAAAAADAJDRkAAAAAGASGjIAAAAAMMn/A89zj0nXtKQIAAAAAElFTkSuQmCC)

*missingness plot*

![genotype pca](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAhoAAAGtCAYAAABQnt12AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjksIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvJkbTWQAAAAlwSFlzAAAQ6wAAEOsBUJTofAAAR3ZJREFUeJzt3X1ck/X+P/DXgLFxL3cmqIjoUTQhycAbUsC8Ce81u1ErsSNlVpZ6LI/fEj2VdbRI8ybzpkzK7GiWRy3LWzIsJUuPKchRIVBIBRRwY2Nsn98fHvZzMoEhF9vY6/l48FCu67Nr72tj22vX53N9LpkQQoCIiIhIAk7WLoCIiIhaLgYNIiIikgyDBhEREUmGQYOIiIgkw6BBREREkmHQICIiIskwaBAREZFkGDSIiIhIMgwaREREJBkGDSIiIpIMgwYRNdhTTz2Fu+++GwaDwdql3LGkpCTIZDLJti+TyZCUlCTZ9ql56fV6dO3aFc8884y1S7E7DBoEANBqtVizZg2GDBmCu+66C66urvD29kbPnj3x/PPP48iRI9Yu8bYWLFiAr7/+2tplWCwvLw8ymczkR6lUomvXrpgzZw6uXr1a6zYGgwGff/45RowYgaCgILi6usLHxwcxMTGYP38+CgsLzd5XcXExFAoFZDIZPvroo0bVe+zYMWzYsAGvv/46nJz41pGXl4cFCxbg+PHj1i7ljn399ddYsGCBtcswunbtGhYsWICDBw9auxQjZ2dnLFy4EOvXr8d//vMfa5djXwQ5vNzcXBERESEAiPvvv1+8/vrrYv369WLFihXi2WefFe3atRMAxP79+61dqlkAxOTJk61dhsVyc3MFABEfHy/S0tJEWlqaWLp0qRg0aJAAICIjI4VWqzW2Ly0tFfHx8QKA6Nmzp0hJSRHr168X77//vnjiiSeEp6enaN26tdn7euedd4STk5MICwsTffv2bVS9I0eOFGFhYcJgMDTq9rZm8uTJ4k7eAg8cOCAAiI8//tjs+srKSlFVVdXo7TenO30smlrNayMlJcXapZjQ6/WiXbt24qGHHrJ2KXbFxVoBh2yDRqPBiBEjcObMGfzrX//Cww8/XKvN+++/j08++QRubm5WqLDl69SpEx5//HHj7y+++CJGjBiBXbt2Yfv27cbn5LHHHsPBgwfxz3/+Ey+//HKt7ZSWlt72W+m6deswcOBAjB49Gi+88AJOnTqFu+++u8E15uXlYdeuXUhJSZG0u6ElUSqV1i5BMnq9HlqtFu7u7tYupVEaW7+TkxOeeOIJLFmyBIWFhQgODpaowhbG2kmHrGv58uUCgJg7d67FtzUYDGLNmjUiOjpauLu7C3d3d9G3b1/x1Vdf1WqL/x11OHLkiEhISBAeHh7Cx8dHPProo+LSpUu12l+9elXMnDlThIaGCldXV9G6dWvx2GOPiZycHGObmm+U5n6EEGL48OFCqVSK0tLSWtv/888/hVwuN/lmUlPj/v37Rb9+/YS7u7vw9/cXkydPNlujVqsV//znP0VERIRQKpXCy8tLPPDAAyI9Pb1Bj1/Nt7a//vWvtdbVPC+LFi0SQgjxzTffCACN+ib1ww8/CADis88+E6WlpUKhUIiXXnrJom0sWbJEABC//PJLrXVxcXGiQ4cOIi8vT4wbN060atVKuLu7i0GDBoljx47Vaq/X68WyZctEZGSk8XFLSEgQ33//fa22HTp0EHFxceL48eNi8ODBwtPTU3h7e4uxY8eKs2fPmrT9+OOPBQBx4MCBWtsx943d3LKsrCwxffp0cffddwtvb2+hVCpFjx49xJIlS0R1dbWxXUpKitm/u7i4OGMb3OZI26effipiYmKMr5nevXuLzz///LaPa1FRkXj88ceFn5+fUCqVon///iIzM7NWe3MMBoN4//33Rc+ePYW3t7fw8PAQYWFhYsKECaKwsFAIceMxNrcvNUdqavb11KlTYs6cOSIkJEQ4Ozsb199uP293xKeqqkqkpqaKe++9V7i7uwtPT08REREh5s+fL4T4/8/jrT8dOnQQQtR9tMPc30B99VvyPiaEED/++KMAIN5///2GPAUkeETD4W3duhUAkJycbPFtp0yZgo0bN2L06NGYNGkSAGDbtm0YO3YsPvjgA0ybNs2k/YkTJ5CYmIgnn3wSjz76KI4dO4Z169bh2rVr2L17t7FdRUUFYmNjcfr0aUyYMAH3338/zp07h1WrVmH37t3IyMhA9+7d0a1bN6SlpeGJJ55A//798fTTT5vc3zPPPINdu3YhLS0NM2bMMFm3YcMG6HS6Wrf57bffsHXrVkyZMgWPP/44jh49ik8++QRHjhxBZmYmPD09AQDV1dUYNmwY0tPTMWHCBEybNg1qtRqffvopBg4ciK+//hojRoyw+DGtkZOTAwAIDAwEAPzrX/8CgFqPaUOsXbsWPj4+GDt2LNzc3DB69GikpaXh7bffhkKhaNA2Dhw4ADc3N9xzzz1m16tUKsTFxSEqKgpvvPEGCgoKsGrVKgwYMAAZGRkmt0tKSkJaWhpiY2OxaNEiXL9+HevWrcPQoUOxceNGk6M7AHDhwgUkJCRg1KhRWLx4MbKysrB69WocPnwYx44dQ9u2bS1+TG7n4MGDOHDgAEaMGIGOHTtCo9Hgm2++wZw5c3D+/HmsWrUKADBu3DjodDosWrQITz/9NPr37w8AuOuuu+rc/vz58/H6668jIiICKSkpEELg008/xYQJE3D+/HnMmzfPpL1KpUL//v3Rq1cvvP7667h06RLee+89JCYm4vz58/Dy8qrz/hYtWoRXX30Vw4YNw9SpU+Hq6or8/Hzs3r0bhYWFCAoKwtKlS5GamopDhw4hLS3NeNt+/fqZbGvSpElwcXHBc889B09PT3Tt2rXBj2sNnU6HxMRE7Nu3D3FxcZg/fz68vb2RlZWFLVu2YOHChRgwYADee+89zJw5E2PHjsW4ceMAwPjaa6zb1W/p+9h9990HV1dXHDhwAC+88MId1eQwrJ10yLr8/f2Ft7d3reUGg0FcuXLF5KeiosK4/uuvvxYARGpqaq3bjhw5Unh7e4vy8nLjMgBCJpOJjIwMk7bPPPOMACDOnDljXPbaa68JAOLNN980aXvw4EEBQDzwwAMmy3Gbb1TV1dWiffv2okePHrX2rXPnzqJjx44m4w3wv29OW7ZsMWmfmppa6xvU0qVLBQCxbds2k7ZVVVUiKipKdOzYsVY9t6r5Zvb4448bH+MzZ86IJUuWCBcXF9GqVStx+fJlIYQQvXr1EgBESUlJvdu92dWrV4Wbm5t45plnjMt2794tAIhNmzY1eDsdOnQQ4eHhZtfFxcUJAOK5554zWf7LL78IJycnk2/5+/btEwBEYmKiyRGCy5cvi9atW4tWrVqZ/J3VfNtesmSJyba3bdtW63lviiMa169fN7uPEydOFM7OzqKoqMi4rL4xGrfWl5OTI5ycnMQ999wjVCqVyX326NFDODs7i9zcXOPymse15qhWjc8//1wAEB9++KHZ+71ZVFSU6NatW73t6hqjUXNE4P777zc75uR2rz9zj0/NkbEZM2bUGuuj1+uN/6/rqEVjj2iYq9/S97EanTp1Ep07d661nMzj0HEHV1ZWBm9v71rLL126hMDAQJOfl156ybg+LS0Nbm5uePTRR1FcXGzyM2bMGJSXl+Onn34y2Wbfvn1rfUsaPHgwgP//DR4AvvzyS3h7e2PWrFkmbePi4pCQkID9+/ebPSPjVs7OzkhOTsbvv/+On3/+2bj84MGDOHv2LKZOnVprvEGXLl0wfvx4k2XPPfccWrVqhS+//NJk/0NDQ9G/f3+TfS8rK8OoUaOQm5trsk91+fTTT42Pcc0ZJ5GRkdizZ4/xiEZZWRkAmH2u6tt2ZWUlpkyZYlw2ePBgtG/fHuvWrWvwdq5cuQJ/f/8629z6bbxXr14YOnQo0tPTUVxcDADGx/C1116Ds7OzsW1gYCCee+45XLt2Dfv27TPZjpeXV61vjmPHjkW3bt3w1VdfNempth4eHsb/a7ValJaWori4GA8++CD0ej1++eWXRm/766+/hsFgwCuvvGIyNsDDwwNz5syBXq/H9u3bTW7j5OSEmTNnmiwz95q5nVatWuHixYtIT09vdN01Zs+eDblcfkfb+PTTT+Hh4YFFixbVeu1JfSaTufob8z4GAP7+/rh8+bKk9bYk7DpxcD4+PigvL6+13M/PD3v27AFwI3Tcejg7KysLlZWVdR62vnTpksnvYWFhtdrUfHiVlJQYl50/fx5333232cF0EREROHDgAHJzc+Hr61vHnt0wdepU/OMf/8DatWvRp08fADe6ElxcXPDUU0/Vat+9e/day1xdXdGpUyecOnXKuCwrKwtqtdoYBMy5dOkSunTpUm+NiYmJmDVrFmQyGRQKBUJDQ9GuXTuTNj4+PgCA8vJy+Pn51bvNGmvXrkX79u3h7++Ps2fPGpc/+OCDWLduHc6dO4dOnTo1aFtCiNuua9WqldmBcd27d8e3336Lc+fOISAgAOfPnwdw43m8Vc2yc+fOmSzv1KmT2S6e7t27IysrC1euXKm3y6Kh1Go1Xn/9dWzevBl5eXm11peWljZ6243Z9+Dg4FqvA3Ovmdt56623MHbsWMTHx+Ouu+5C//798cADD2DChAnGv6mGasjfcn1ycnIQHh5uEuiai7n6G/M+Btx4LXBQdMMxaDi4Hj16ID09HefPnzcJAq6urhg0aBAAmH3DNRgM8PHxMY7xMOfWsxpu/gZ7q7o+xO5EUFAQRo0ahS+++AJLly5FVVUVtm3bhlGjRqFNmzaN3q7BYEDXrl2xYsWK27bp0aNHg7YVHBxsfKxvJyIiAseOHcOvv/5ab9saR48eNZ7v/5e//MVsm/Xr12PRokX1biswMLBBH2zWVNcbf3V1dYO2MWnSJGzfvh1Tp07FgAEDEBAQABcXFxw7dgxz585t9onK7vQ107t3b5w9exZ79+7FgQMHkJ6ejq1bt2L+/Pn44YcfEB4e3uBaLD1Do6GPuaUa+zybq78x72PAjZDXunXreiqlGgwaDm78+PFIT0/H2rVr8dZbbzX4dl26dEF2djaioqLqPaRuqU6dOuHs2bPQarW1vsn+/vvvkMlk6NixY4O3N23aNGzbtg2bNm1CZWUltFptrUGgNU6fPl1rWVVVFc6dO4fOnTsbl3Xp0gUFBQWIj4+Hi4v0L6OHH34YGzZswJo1axocNNauXQuZTHbbU5PffvttfPzxx/jHP/5R7z706NED+/btg06nM3v4/Nq1a2ZP96t5PGuOmtT8e+rUKfTu3duk7e+//27Spsa5c+fM/i2cPn0a3t7exqNKNUd6zB11qDmaUJeysjJs374djz/+ONasWWOy7r///W+t9pZ+o715328Nobfb96bg7u6OUaNGYdSoUQCA3bt3IzExEW+//TY2bNgAwPJ9uZmfn1+DH/MuXbogJycHKpWqzqMaddVzp8/zrfVY+j6m0Whw4cKFOxrs7Wg4RsPBTZ06Fd27d8c777yDLVu2mG1j7pvTk08+CQB4+eWXza43d7ixocaNG4eysjIsX77cZPmhQ4ewf/9+DBw40KTbxNPTs85D2oMGDUKnTp2wdu1arFu3DqGhoRgyZIjZtjk5ObW+3axcuRLXrl0zjn4Hbuz/1atX8eabb5rdzp3svznDhg3D4MGDsWXLFqSmppptc/XqVbz44osAgOvXr2Pz5s2IjY3FE088gfHjx9f6SUpKwp9//omdO3fWe//x8fHQaDT49ddfb9vm1iMjx44dw3fffWc8MgDA+BguWrTI5OhAcXExVq5ciVatWuGBBx4w2U5FRUWtv4WvvvoKWVlZGDNmjLFvv+bQ+N69e03aHjp0yGSMzu3UbOfWv+eKigqzj3nNWRAN7U6pqfWdd96BRqMxLler1ViyZAmcnZ0xevToBm2roa5cuVJrWa9evQCYdr1Yui8369q1K3766Seo1WrjMo1GU+s5A4DHH38cKpUKr732Wq11N/891FWPl5cXgoKCsH//fpPnqqSkxOJZbxvzPnbs2DFUVVUhISHBovtyZDyi4eCUSiV27dqFkSNH4pFHHsH999+PoUOHIjg4GCqVCufOnTOeWhkaGmq83UMPPYTk5GSsXbsWJ06cwJgxY9CmTRsUFhbi2LFj+Oabb6DT6RpV05w5c/Dll19izpw5OHHiBPr162c8vdXHxwfvv/++Sfs+ffpg7969+Oc//4mQkBDIZDI89thjxvUymQxPP/00XnnlFQDAG2+8cdtvTBEREUhKSsIPP/yAbt26ITMzExs2bECXLl3wt7/9zdjuxRdfxL59+7BgwQL88MMPGDJkCPz8/FBQUIDDhw/j/PnzFn+7qs8XX3yBcePGYfbs2fj0008xZswYtG/fHiqVCr/++iu+/PJLuLm5YdmyZfj8889x/fp1sxOw1XjooYfw4osvYt26dRgzZkyd9z1+/Hi88sor2LlzZ60jEQAQEBCAnTt34uLFixg8eDAKCgqwcuVKKJVKLF261Nhu4MCBeOKJJ5CWloaEhASMHTvWeHrr5cuXsXHjxlqnMXbq1AmLFi0yHgXJysrCBx98gMDAQLzxxhvGdl27dsXQoUOxevVq6PV69OrVC1lZWdiwYQMiIyNx4sSJOvfRy8sLDz74ID777DMoFAr07t0bRUVFWL9+vdkxIN27d4eXlxdWrVoFd3d3tGrVCq1bt8bAgQPNbr9z5874v//7P7z++uvo06cPJk2aZDy99eTJk3jzzTdNXmNNoVu3bujduzdiYmLQrl07lJaW4pNPPgEATJ482diuT58+WLFiBaZPn47hw4dDLpejd+/eDTpyOGPGDEyYMAHx8fF48skncf36dWzcuNHsGJAXX3wRu3btwnvvvYfffvsNiYmJ8Pb2Rk5ODr7//nvjkR1/f3907twZmzdvRqdOnXDXXXfBw8MDI0eONN7n3//+dwwdOhRjx47FlStXsHbtWnTs2NGikN+Y97GdO3fCxcXF5IsH1cNap7uQbdFoNGL16tVi0KBBIjAwULi4uAhPT08RGRkppk+fLo4cOWL2dps2bRLx8fHCx8dHuLq6ivbt24vExETxwQcfmLSDhZP6lJaWipdeekl06NBByOVyERAQIB577DGT02Br5OTkiMGDBwsvLy+TCbtudvnyZeHq6ipcXFyMExXdqqbGmgm73NzchK+vr3jiiSfEn3/+Wat9dXW1WLVqlejdu7fw9PQUSqVShIaGinHjxokvvvjC7H3crK4Ju26nurpafPbZZyIxMVG0bt1auLi4CC8vLxEdHS0WLFhgPP0yOjpayGQyceHChTq3179/f+Hs7CwKCgrqve+RI0eKkJCQWqcl1kwslZubK8aNGyd8fHyEu7u7GDhwoNmJpfR6vVi6dKmIiIgQCoVCeHp6ioSEBPHdd9/Vamtuwi4vLy8xevRo8d///rdW+0uXLonHHnvMWMOAAQPE4cOHG3x6a0lJiXjmmWdE27ZthUKhEF27dhWLFy8We/fuNft3umvXLhEVFSUUCkWDJ+xKS0sTMTExws3NTbi5uYnevXubPdW45nE153bbvtVbb70l4uLiROvWrYVcLhdt2rQRDz74YK3J0fR6vZg9e7Zo27atcHJyMjth182n3t7qvffeEx07dhRyuVx06tRJLFmyxHgq862PmbmJ7iIjI8WCBQtM2h05csQ4cR5umrBLiBuvg3nz5ong4GDh6uoq7r77bvHxxx/XeXprXfU39H2MU5A3jkwIiUbhEdmQq1evIjg4GImJidi2bZvZNjKZDJMnTzb2W5OpY8eOITo6Gl988YXJkZL4+Hjk5eWZHTR8p0JDQxEaGmpTF9cix7V582Y8/vjj+PXXXxEZGWntcuwGx2iQQ1i7di00Gk2jZtakG3r16oWkpCSkpKS0iMvEE1lCr9cjJSUFf/3rXxkyLMQxGtSibd68GRcvXsTrr7+O6Ojo2w4CpYZp7CXmieyds7Mzzpw5Y+0y7BKDBrVoEyZMgEKhQN++fS2aCZOIiJoGx2gQERGRZDhGg4iIiCTDoEFERESSYdAgIiIiyTj8YFCNRoOTJ08iMDCwWa5ZQUREZO+qq6tx5coVREREmL3S9s0c/pP15MmTiImJsXYZREREdufo0aOIjo6us43DB42aKz8ePXoUQUFBVq6GiIjI9hUVFSEmJsb4GVoXhw8aNd0lQUFBaNeunZWrISIish8NGXLAwaBEREQkGQYNIiIikgyDBhEREUmGQYOIiIgkw6BBREREkmHQICIiIskwaBAREZFkGDSIqMXTGwxQa3XQGwzWLoXI4Tj8hF1E1HKpNFVIP56PzOxCqDQ6eCjliA4PRlzPEHgoXa1dHpFDYNAgohZJpanC2h2/ISu/GE4yGRRyF1SotdiecQY5BSVIHhnFsEHUDNh1QkQtUvrxfGTlF8PX0w1B/l7w877xr6+nG7LzS5B+PN/aJRI5BAYNImpx9AYDMrML4SSTwV0pN1nnrpRDJgMysws5ZoOoGTBoEFGLo9XpodLooJCb7x1WyF2g1uig1embuTIix8OgQUQtjkLuDA+lHFpdtdn1Wl013JVyKOTOzVwZkeNh0CCiFsfZyQnR4cEwCAG1RmeyTq3RQQggOjwYzk58CySSGs86IaIWKa5nCHIKSpCdX4IylQYKuQu0umoIAYSH+COuZ4i1SyRyCAwaRNQieShdkTwyyjiPhlqjg6+XkvNoEDUzBg0iarE8lK4Y1qczhsaEQavTQyF3ZncJUTOz21fcli1bMHr0aLRr1w4eHh7o2bMnPvroIwghrF0aEdkYZycnuCvkDBlEVmC3RzRSU1MRGhqKd999F4GBgdizZw+Sk5NRUFCAlJQUa5dHREREsOOgsWPHDgQEBBh/HzhwIEpKSpCamorXXnsNTvzmQkREZHV2+2l8c8ioERUVhfLycqhUKitURERERLey26Bhzo8//oi2bdvCy8vL2qUQERER7Ljr5FY//vgjNm/ejHfffbfOduXl5SgvLzf+XlRUJHVpREREDqtFBI0LFy7g0UcfRUJCAmbMmFFn29TUVCxcuLCZKiMiInJsMmHn54Neu3YN/fv3h0wmw6FDh+Dj41Nne3NHNGJiYlBQUIB27dpJXS4REZHdu3DhAtq3b9+gz067PqJRWVmJESNGoKysDD/99FO9IQMAvL294e3t3QzVERERkd0GjerqajzyyCPIysrCoUOH0LZtW2uXRERERLew26Axffp07Ny5E++++y7Ky8vx888/G9dFRUVBoVBYsToiIiIC7DhofP/99wCA2bNn11qXm5uL0NDQZq6IiIiIbmW3QSMvL8/aJRAREVE9WtSEXURERGRbGDSIiIhIMgwaREREJBkGDSIiIpIMgwYRERFJhkGDiIiIJMOgQURERJJh0CAiIiLJMGgQERGRZBg0iIiISDIMGkRERCQZBg0iIiKSDIMGERERSYZBg4iIiCTDoEFERESSYdAgIiIiyTBoEBERkWQYNIiIiEgyDBpEREQkGQYNIiIikgyDBhEREUmGQYOIiIgkw6BBREREkmHQICIiIskwaBAREZFkGDSIiIhIMgwaREREJBkGDSIiIpIMgwYRERFJhkGDiIiIJMOgQURERJJh0CAiImogvcEAtVYHvcFg7VLshou1CyAiIrJ1Kk0V0o/nIzO7ECqNDh5KOaLDgxHXMwQeSldrl2fTGDSIiIjqoNJUYe2O35CVXwwnmQwKuQsq1FpszziDnIISJI+MYtioA7tOiIiI6pB+PB9Z+cXw9XRDkL8X/Lxv/Ovr6Ybs/BKkH8+3dok2jUGDiIjoNvQGAzKzC+Ekk8FdKTdZ566UQyYDMrMLOWajDgwaREREt6HV6aHS6KCQmx9poJC7QK3RQavTN3Nl9oNBg4iI6DYUcmd4KOXQ6qrNrtfqquGulEMhd27myuwHgwYREdFtODs5ITo8GAYhoNboTNapNToIAUSHB8PZiR+nt8OzToiIiOoQ1zMEOQUlyM4vQZlKA4XcBVpdNYQAwkP8EdczxNol2jQGDSIiojp4KF2RPDLKOI+GWqODr5eS82g0EIMGERFRPTyUrhjWpzOGxoRBq9NDIXdmd0kDMWgQERE1kLOTE9wVDBiW4KNFREREkmHQIGrheBEoIrImdp0QtVC8CBQR2QIGDaIWiBeBIiJbwa4TohaIF4EiIlvBoEHUwvAiUERkSxg0iFoYXgSKiGwJgwZRC8OLQBGRLWHQIGpheBEoIrIlPOuEqAXiRaCIyFbY9Veas2fPYtq0aejZsydcXFzQo0cPa5dEZBNqLgI1KrYL2vh5QiF3Rhs/T4yK7cJTW4moWdn1EY1Tp05h165d6N27NwwGAwwcRU9kxItAEZEtsOt3nZEjR6KgoABbt27Fvffea+1yiGzSjYtAyRkyiMgq7Pqdx4lvnERERDaNn9REREQkGbseo9EY5eXlKC8vN/5eVFRkxWqIiIhaNocLGqmpqVi4cKG1yyAiInIIDtd1MmvWLBQUFBh/jh49au2SiIiIWiyHO6Lh7e0Nb29va5dBRETkEBzuiAYRERE1H7s+oqFWq/HNN98AAP744w+Ul5dj69atAIC4uDgEBgZaszwiIruhNxg4sRtJwq6DxuXLl/Hwww+bLKv5/cCBA4iPj7dCVURE9kOlqUL68XxkZhdCpdHBQylHdHgw4nqGcKp6ahJ2HTRCQ0MhhLB2GUREdkmlqcLaHb8hK78YTjIZFHIXVKi12J5xBjkFJbwuDjUJHh8jInJQ6cfzkZVfDF9PNwT5e8HP+8a/vp5uyM4vQfrxfGuXaDf0BgPUWh30vOZWLXZ9RIOIiBpHbzAgM7sQTjIZ3JVyk3XuSjnKVBpkZhdiaEwYx2zUgV1P9WPQICJyQFqdHiqNDgq5+Y8BhdwFao0OWp0e7goGDXPY9dQw/OshInJACrkzPJRyaHXVZtdrddVwV8qhkDs3c2X2g11PDcOgQQT2r5LjcXZyQnR4MAxCQK3RmaxTa3QQAogOD2a3yW3U1/UkkwGZ2YV8TwG7TsjBsX+VHFlczxDkFJQgO78EZSoNFHIXaHXVEAIID/FHXM8Qa5dos9j11HAMGuSw2L9Kjs5D6YrkkVHGsK3W6ODrpWTYboCarqcKtdbseq2uGr5eSnY9gUGDHNjN/as3H/pUa3TG/tVhfTpbsUIi6XkoXTGsT2cMjQnjzKAWqOl62p5xBmqNrtZ7CLue/j8+AuSQ2L9KZMrZyQnuCjk/GC0Q1zME3UICcO26BkUlFSgtr0RRSQWuXdew6+kmPKJBDon9q0R0p9j11DAMGuSQ2L9KRE2BXU/146NBDomn9hFRU2LX0+3xiAY5LJ7aR0QkPQYNcljsXyUikh6DBjk09q8SEUmLQYMINf2rDBhERE2N76xEREQkGQYNIiIikgyDBhEREUmGQYOIiIgkw6BBREREkmHQICIiasH0BgPUWp3VLhLJ01uJiIhaIJWmyjghoUqjg4dSbpUJCe84aKhUKly6dAmVlZXw9/dHmzZtmqIuIiIiaiSVpgprd/yGrPxiOMlkUMhdUKHWYnvGGeQUlCB5ZFSzhY1GBY2TJ09iw4YN2LNnD06fPg0hhHGdj48P+vXrh4cffhgPP/ww3N3dm6xYIiIiql/68Xxk5RfD19MN7kq5cblao0N2fgnSj+djWJ/OzVKLRWM0fvrpJ8THx+Oee+5BRkYGBg0ahPXr1+Pf//43vvvuO3zxxReYN28ePD09MXv2bLRt2xZvvvkmVCqVVPUTERHRTfQGAzKzC+Ekk5mEDABwV8ohkwGZ2YXNNmbDoiMaI0aMwIwZM7Bx40aEhNR9Zcvq6mp89913SE1NhcFgwGuvvXZHhRIREVH9tDo9VBodFHLzH/EKuQvUGh20On2zXHrBoqDxxx9/wNPTs2EbdnHB8OHDMXz4cB7RICIiaiYKuTM8lHJUqLVm12t11fD1UkIhd26WeiyKMg0NGbfy8PBo1O2IiMi+WfvUSkfk7OSE6PBgGISAWqMzWafW6CAEEB0e3GxXqm6S01uFEFi3bh327NkDIQQGDRqE5ORkOPFy20REDslWTq10VHE9Q5BTUILs/BKUqTRQyF2g1VVDCCA8xB9xPese/tCUmiRozJkzB9u2bcNDDz0ElUqFuXPnIisrC0uXLm2KzRMRkR2xpVMrHZWH0hXJI6OMYU+t0cHXS2n782gUFhYiODi41vLPPvsMv/32m3EOjYSEBEyfPp1Bg4jIAdnSqZWOzEPpimF9OmNoTBi0Oj0Ucudm6y65mUX3GBERgbfffhs6nWmfj4eHB/Ly8oy///HHHxyXQUTkgGzt1Eq6MWbDXSG3SsgALAwaP//8M3788Ud0794dO3fuNC6fN28e4uPjERMTgx49emDu3Ll49dVXm7xYIiKybZacWkmOwaKuk7/85S/YuXMndu3ahVmzZmHVqlVYtmwZnnrqKURHR+PgwYMAgPj4eEREREhRLxER2TBbO7WSrK9Rg0GHDx+OIUOGIDU1FX379sVTTz2F+fPnM1wQETm4mlMrt2ecgVqjqzVGo7lPrSTra/QzLZfL8corr+DkyZMoKipC165dsXHjxqasjYiI7FBczxB0CwnAtesaFJVUoLS8EkUlFbh2XdPsp1aS9VkUNIqLi/Hkk0+iTZs28PX1xdChQ1FaWoq0tDRs2bIFy5cvR9++fXHs2DGp6iUiIhtXc2rlqNguaOPnCYXcGW38PDEqtgtPbXVAFnWdJCUloaCgAO+//z7c3d3x4Ycf4sEHH0ReXh769euHo0ePYt26dRg+fDhGjBiBdevWSVU3ERHZMFs5tZKsz6Jn/dChQ3jnnXfwyCOPYMSIEdi4cSMuXryI8+fPAwBkMhmSk5ORnZ3d6OnKiYio5bD2qZVkfRbPo5GWlobS0lKo1Wp8+OGH8Pb2rnUl11atWnGyLiIiIrKs6+Tjjz9GUlISAgICIJPJEBYWhi1btkChUEhVHxEREdkxi+fRyMjIgEqlQlVVFXx9faWqi4iIiFqARnWaeXh41AoZQgjk5ORAo9E0SWFERERk/5psdE55eTm6deuGzMzMptokERER2TmLuk5mzJhx23VVVVUQQmDp0qXYsmULZDIZli1bdscFEhERkf2yKGisWLECPj4+aNWqVa11BoMBMpkMhw8fhlKpZNAgIiIiy7pOZsyYAYPBgKSkJGRlZSE3N9f4c+LECQgh8MUXXyA3N9c4twYRERE5LouCxtKlS/Hjjz/iwIEDCA8Px9atW43rZDJZkxdnr/QGA9RaHfQGg7VLISIisiqLr94aERGBgwcP4rPPPsOLL76IFStWYPny5bUm7XJEKk0V0o/nIzO7ECqNDh5KOaLDgxHXM4Rz+xMRkUNq9FknkyZNwpkzZ9CrVy/07t0bM2bMcOijGipNFdbu+A3bM87gz9LrqNLp8WfpdWzPOIO1O36DSlNl7RKJiIia3R2d3urp6Yl3330XmZmZuHTpEkJCQqBUKpuqNruSfjwfWfnF8PV0Q5C/F/y8b/zr6+mG7PwSpB/Pt3aJREREzc7irhNz7r77buzevbspNmWX9AYDMrML4SSTwV0pN1nnrpSjTKVBZnYhhsaE8cJCRETkUOz6Uy87OxuDBw+Gh4cH2rRpg5dffhlVVc3fRaHV6aHS6KCQm89tCrkL1BodtDp9M1dGRERkXU1yRONWgwcPhsFgwL59+6TYPADg6tWrGDhwIP7yl79g27ZtuHjxImbNmgW1Wo0VK1ZIdr/mKOTO8FDKUaHWml2v1VXD10sJhdy5WesiIiKyNkmCRm5uLgwSn9q5evVqlJeX46uvvoKfnx8AoLq6GtOnT8e8efMQHBws6f3fzNnJCdHhwdiecQZqjc6k+0St0UEIIDo8mN0mRETkcCT55Dt79qzkE3Z9++23GDRokDFkAMAjjzwCg8GA77//XtL7NieuZwi6hQTg2nUNikoqUFpeiaKSCly7rkF4iD/ievL0XyIicjx2+xU7Ozsb4eHhJstatWqFoKAgZGdnN3s9HkpXJI+MwqjYLmjj5wmF3Blt/DwxKrYLkkdGcR4NIiJySBZ1nRw9ehT33XcfnGygC+Dq1atmr7ni6+uL0tLS296uvLwc5eXlxt+LioqarCYPpSuG9emMoTFh0Or0UMid2V1CROTA9AaDw38eWBQ0+vTpg9atW+ORRx7BpEmT0Lt3b6nqkkxqaioWLlwo6X04OznBXeGYf1BERMSZom9m8adheHg4Pv74Y/Tr1w+dOnXCa6+9htOnT0tRW518fX1RVlZWa/nVq1dNxm3catasWSgoKDD+HD16VMoyiYjIwXCmaFMWB43Fixfj0qVLSEtLQ7du3bB48WJEREQgKioK77zzDi5cuCBFnbWEh4fXGotRVlaGoqKiWmM3bubt7Y127doZf4KCgqQulYiIHAhnijbVqOP77u7umDhxInbu3ImioiKsWLECXl5eeOWVVxAaGoq4uDisWbOmqWs1kZiYiL179+LatWvGZVu2bIGTkxOGDBki6X0TETkCR7sSdVPsb30zRctkQGZ2ocM8pkATzKPh5+eHZ599Fs8++ywKCgqwadMmbNq0Cc8++yyefvrppqjRrGnTpmH58uUYM2YM5s2bh4sXL2LOnDmYNm1as86hQUTU0jja+IKm3F9LZop2lLF8TbqX7du3xyuvvIITJ07g999/b8pN1+Lr64t9+/bBxcUFY8aMwdy5czF16lSkpqZKer9ERC2Zo40vaOr9rZkpWqurNrteq6uGu1LuUDNFW3REo0OHDlAoFA1q261bt0YVZIlu3bph7969kt8PEZGjuHl8wa2zHNeMLxjWp/Md3YctnfLZ1PvLmaJrsyho5ObmSlUHERFZmdRXora1Lhmp9jeuZwhyCkqQnV+CMpUGCrkLtLpqCAGHnClakmudELVEtvQtjEgKUo4vqOmiyMovhpNMBoXcBRVqLbZnnEFOQYlVZlCWan9rZoquCVVqjQ6+XsoWPc6lLpIEjaysLGzZsgXz58+XYvNEzcrWvoURSUXKK1E3R5eMpaTcX84U/f9JstenT5+WfPZNoubgaAPjyLHVjC8wCAG1Rmey7k7GF9jqKZ9S7e+t9+GukDtsyAAsPKJR1zVEblZRUdGoYohsjS1+CyOSkhTjC2z5lE+Op5CeRUEjICAAMpms3nZCiAa1I7JlUg+MI7JFUowvkLKL4k5xPIX0LAoa7u7uGDBgAB599NE622VmZuKDDz64o8KIrM2Wv4URSampxxfY+imfHE8hLYuCRq9evWAwGDB58uQ623l6ejJokN2z5W9hRM2hKa9EbQ9dFLzytjQsekR79+6NzMzMett5eHggJMT6fzREd6I5BooROYqaLopRsV3Qxs8TCrkz2vh5YlRsF6uc2krNRyaEEA1trNPpoFar4ePjI2VNzerChQto3749CgoK0K5dO2uXQzam5qyT7PwSyGSo9S2Mb5BEluOcNPbPks9Oi7pO5HJ5iwoZRPXhQDGipscuCsfCmUGJ6sGBYkREjWfxGI2vv/4ahgZOqlJQUIC//e1vvKIqtQiceIeIyHIWHdF48sknMX36dDz99NMYPXo0YmNjERkZicDAQCgUCly7dg25ubk4duwYvv32W/z8888YNWoUnn32WanqJyIiIhtm0WBQAKisrMTmzZuxceNGHD58GNXV1SbrhRAICgrC+PHjMXXqVERERDRpwU2Ng0GJiIgsI9lgUABwc3PDlClTMGXKFGg0Ghw/fhxFRUXQaDTw8/ND165dERoa2tjaiYiIqAW5o8GgSqUSffr0aapaiIiIqIXhqDYiIiKSDIMGERERSYZBg4iIiCTDoEFERESSYdAgIiIiyTBoEBERkWQsDhppaWm45557EBgYiAEDBmDHjh212hw5cgTOzs5NUiARERHZL4uCxr///W9MnjwZbdq0wV//+lcYDAaMGTMGycnJ0Ov1UtVIREREdsqiCbveeustPP3001i9erVx2aZNm/Dss88iPz8fX375JTw9PZu8SCIiIrJPFh3ROH36NB5++GGTZRMnTsShQ4dw6tQpxMfH4/Lly01aIBEREdkvi4KGu7s7rl+/Xmt5ZGQkMjIycP36dcTGxuLs2bNNViARERHZL4uCRmRkJL755huz6zp06ICMjAy0atUKTz31VJMUR0RERPbNoqAxbtw47N69G6WlpWbX+/v74+DBg4iLi4OFV58nIiKiFkgmHDwRXLhwAe3bt0dBQQHatWtn7XKIiIhsniWfnRbPo1FRUQGNRnPb9RqNBhUVFZZuloiIiFogi4LG/v374efnhyNHjty2zZEjR+Dv748ffvjhjosjIiIi+2ZR0Fi5ciUeffRRxMXF3bZNXFwcJkyYgPfff/+OiyMiIiL7ZlHQyMjIwEMPPVRvu7Fjx+LHH39sdFFERETUMlgUNK5evYrAwMB62wUEBODq1auNLoqIiIhaBouCRkBAAM6dO1dvu/PnzyMgIKDRRREREdkivcEAtVYHvcFg7VLshkXXOomLi8PKlSsxadIkuLiYv2l1dTVWrlyJhISEJimQiIjI2lSaKqQfz0dmdiFUGh08lHJEhwcjrmcIPJSu1i7Ppll0ROPvf/87Tp48ieHDh+P06dO11mdlZWHEiBH4z3/+g7lz5zZZkURERNai0lRh7Y7fsD3jDP4svY4qnR5/ll7H9owzWLvjN6g0VdYu0aZZdEQjIiICmzdvxuTJkxEREYHg4GCEhIRAJpMhPz8fFy9ehJeXF7744gv06NFDqpqJiIiaTfrxfGTlF8PX0w3uSrlxuVqjQ3Z+CdKP52NYn85WrNC2WTxh1+jRo3HmzBksWLAA4eHhuHr1KkpLS9G1a1f84x//wJkzZzBq1CgpaiUiImpWeoMBmdmFcJLJTEIGALgr5ZDJgMzsQo7ZqINFRzRq3HXXXXjttdeauhYiIiKbotXpodLooJCb/7hUyF2g1uig1enhrrD4u7tDsDhonD59GqtXr0Zubi7atm2L8ePHY9CgQVLURkREZFUKuTM8lHJUqLVm12t11fD1UkIhd27myuyHRfHrxx9/xL333ouVK1ciMzMT69evx9ChQ7F69Wqp6iMiIrIaZycnRIcHwyAE1BqdyTq1RgchgOjwYDg78WjG7Vj0yKSkpCA8PBx5eXn4888/UVJSgjFjxuDVV1+Vqj4iIiKriusZgm4hAbh2XYOikgqUlleiqKQC165rEB7ij7ieIdYu0aZZ1HVy8uRJrF69Gu3btwcAeHt7491330VYWBgKCgqMy4mIyDHoDQZodXoo5M4t9lu9h9IVySOjjPNoqDU6+HopOY9GA1kUNIqLi2tdd74mXBQXFzNoEBE5CEebwMpD6YphfTpjaExYiw9WTc3iwaAymUyKOoiIyE7UTGCVlV8MJ5kMCrkLKtRabM84g5yCEiSPjGqRYQO4MWaDZ5dYxuKgkZCQACczKa5///4my2UyGcrKyu6sOiIisjmcwIosYVHQSElJkaoOIiKyA/VNYFWm0iAzuxBDY8LYtUAAGDSIiMgCnMCKLMW/AiJyeLz0d8PVTGCl1VWbXa/VVcNdKecEVmTUqCnIiYhaAkc7c6Ip1ExgtT3jDNQaXa0xGpzAim7FoEFEDsmRz5y4U3E9Q5BTUILs/BKUqTRQyF2g1VVDCHACK6rFbiPnnj17MHHiRHTq1AkymQzPP/+8tUsiIjty85kTQf5e8PO+8a+vp5vxzAkyr2YCq1GxXdDGzxMKuTPa+HliVGwXBjSqxW6PaOzevRsnTpxAXFwcSktLrV0OEdkRnjlx5ziBFTWU3f5VLFmyBKdOncJHH30EHx8fa5dDRHbEkjMnqG43JrCSM2TQbdntX4a5ScOIiBqCZ04QNR9+WhORw+Glv4maj92O0Wis8vJylJeXG38vKiqyYjVEZC08c4KoedhM0CgrK2vQh35YWBhcXRs/ojk1NRULFy5s9O2JqGXgpb+JmofNBI0tW7YgOTm53nZZWVkIDw9v9P3MmjULU6dONf5eVFSEmJiYRm+PiOwXz5wgkp7NBI2pU6eaBACpeHt7w9vbW/L7ISL7wUt/E0mHrywiIiKSjM0c0bDUH3/8gczMTACAWq3GuXPnsHXrVgDA+PHjrVkaERER/Y/dBo0DBw5gypQpxt93796N3bt3AwCEENYqi4haGL3BwPEbRHfAboNGUlISkpKSrF0GEbVQvLIrUdOw26BBRCQVXtmVqOnwOCAR0S14ZVeipsOgQUR0k/qu7CqTAZnZhdAbDFaqkMi+MGgQEd2EV3YlaloMGkREN+GVXYmaFoMGEdFNeGVXoqbFs06IiG7BK7sSNR0GDSKiW/DKrkRNh0GDiMgMXtmVqGkwaBAR1YFXdiW6M3z1EBERkWQYNIiIiEgyDBpEZFP0BgPUWh1n3iRqIThGg4hsAq+WStQyMWgQkdXxaqlELRe7TojI6ni1VKKWi0GDiKyKV0slatkYNIjIqni1VKKWjUGDiKyKV0slatkYNIjIqni1VKKWjWedEJHV8WqpRC0XgwYRWR2vlkrUcjFoEJFN4NVSiVomBg0isim8WipRy8JXMxEREUmGQYOIiIgkw6BBREREkmHQICIiIskwaBAREZFkGDSIiIhIMgwaREREJBkGDSIiIpIMgwYRERFJhkGDiIiIJMOgQURERJJh0CAiIiLJMGgQERGRZBg0iIiISDIMGkRERCQZBg0iIiKSDIMGERERSYZBg4iIiCTDoEFERESSYdAgIiIiyTBoEBERkWQYNIiIiEgyDBpEREQkGQYNIiIikgyDBhEREUmGQYOIiIgkw6BBREREkmHQICIiIskwaBAREZFk7DJo6PV6LF68GAMGDEBAQAD8/PyQkJCAQ4cOWbs0IiIiuoldBo3Kykq89dZb6NWrFz755BNs2rQJvr6+SEhIwP79+61dHhEREf2Pi7ULaAw3NzecP38evr6+xmWDBw9Gjx498N5772HgwIFWrI6IiIhq2OURDWdnZ5OQUbMsMjIShYWFVqqKiIiIbmWXQcOc6upq/Pzzz+jWrZu1SyEiIqL/scuuE3MWL16MixcvYubMmXW2Ky8vR3l5ufH3oqIiqUsjIiJyWDYTNMrKyhr0oR8WFgZXV1eTZXv27EFKSgrmz5+PXr161Xn71NRULFy48I5qJSIiooaRCSGEtYsAgHXr1iE5ObnedllZWQgPDzf+/uuvvyI+Ph5jx47FJ598Uu/tzR3RiImJQUFBAdq1a9e44omIiBzIhQsX0L59+wZ9dtpM0GiMs2fPIjY2FlFRUdixYwfkcrnF27DkwSIiIiLLPjvtdjBoUVERhgwZgpCQEGzdurVRIYOIiIikZTNjNCxRWVmJxMREFBcXY9myZfj999+N6xQKBaKioqxYHREREdWwy6Bx6dIlnDhxAgAwatQok3UdOnRAXl6eFaoiIiKiW9ll0AgNDYUdDy0hIiJyGHY7RoOIiIhsH4MGERERSYZBg4iIiCTDoEFERESSYdAgIiIiyTBoEBERkWQYNIiIiEgyDBpEREQkGQYNIiIikgyDBhEREUmGQYOIiIgkw6BBREREkmHQICIiIskwaBAREZFkGDSIiIhIMgwaREREJBkGDSIiIpIMgwYRERFJhkGDiMjO6A0GqLU66A0Ga5dCVC8XaxdAREQNo9JUIf14PjKzC6HS6OChlCM6PBhxPUPgoXS1dnlEZjFoEBHZAZWmCmt3/Ias/GI4yWRQyF1QodZie8YZ5BSUIHlkFMMG2SR2nRAR2YH04/nIyi+Gr6cbgvy94Od9419fTzdk55cg/Xi+tUskMotBg4jIxukNBmRmF8JJJoO7Um6yzl0ph0wGZGYXcswG2SQGDSIiG6fV6aHS6KCQm+/tVshdoNbooNXpm7kyovoxaBAR2TiF3BkeSjm0umqz67W6argr5VDInZu5MqL6MWgQEdk4ZycnRIcHwyAE1BqdyTq1RgchgOjwYDg78S2dbA/POiEisgNxPUOQU1CC7PwSlKk0UMhdoNVVQwggPMQfcT1DrF0ikVkMGkREdsBD6YrkkVHGeTTUGh18vZScR4NsHoMGEZGd8FC6YlifzhgaEwatTg+F3JndJWTzGDSIiOyMs5MT3BUMGGQf+JdKREREkmHQICIiIskwaBAREZFkGDSIiIhIMgwaREREJBkGDSIiIpIMgwYRERFJxuHn0aiuvnGRoqKiIitXQkREZB9qPjNrPkPr4vBB48qVKwCAmJgYK1dCRERkX65cuYLQ0NA628iEEKJ5yrFNGo0GJ0+eRGBgIFxcHD53SaqoqAgxMTE4evQogoKCrF0O/Q+fF9vE58U28Xm5obq6GleuXEFERASUSmWdbR3+k1WpVCI6OtraZTiUoKAgtGvXztpl0C34vNgmPi+2ic8L6j2SUYODQYmIiEgyDBpEREQkGQYNajbe3t5ISUmBt7e3tUuhm/B5sU18XmwTnxfLOfxgUCIiIpIOj2gQERGRZBg0iIiISDIMGkRERCQZBg0iIiKSDIMGNTu9Xo/FixdjwIABCAgIgJ+fHxISEnDo0CFrl+bw9uzZg4kTJ6JTp06QyWR4/vnnrV2Sw8nOzsbgwYPh4eGBNm3a4OWXX0ZVVZW1y3J4Z8+exbRp09CzZ0+4uLigR48e1i7JbjBoULOrrKzEW2+9hV69euGTTz7Bpk2b4Ovri4SEBOzfv9/a5Tm03bt348SJE4iLi0OrVq2sXY7DuXr1KgYOHIiqqips27YNixYtwpo1azBr1ixrl+bwTp06hV27dqFz587o3r27tcuxKzy9lZqdXq9HeXk5fH19TZb16NEDnTt3xo4dO6xYnWMzGAxwcrrx/SM0NBQjRozAihUrrFyV43jrrbfw5ptvIj8/H35+fgCANWvWYPr06cjPz0dwcLCVK3RcN782kpKS8Msvv+D333+3clX2gUc0qNk5OzubhIyaZZGRkSgsLLRSVQTA+EZK1vHtt99i0KBBxpABAI888ggMBgO+//57K1ZGfG00Hh85sgnV1dX4+eef0a1bN2uXQmQ12dnZCA8PN1nWqlUrBAUFITs720pVEd0ZBg2yCYsXL8bFixcxc+ZMa5dCZDVXr141OzbG19cXpaWlzV8QURNw+MvEU9MoKytDUVFRve3CwsLg6upqsmzPnj1ISUnB/Pnz0atXL6lKdEh38rwQETUFBg1qElu2bEFycnK97bKyskwODf/666946KGHMHHiRMyfP1/KEh1SY58Xsg5fX1+UlZXVWn716lWTcRtE9oRdJ9Qkpk6dCiFEvT83f5idPXsWiYmJ6NevH9atW2fF6luuxjwvZD3h4eG1xmLUHJXic0T2ikGDrKKoqAhDhgxBSEgItm7dCrlcbu2SiKwuMTERe/fuxbVr14zLtmzZAicnJwwZMsR6hRHdAXadULOrrKxEYmIiiouLsWzZMpNz0RUKBaKioqxYnWP7448/kJmZCQBQq9U4d+4ctm7dCgAYP368NUtzCNOmTcPy5csxZswYzJs3DxcvXsScOXMwbdo0zqFhZWq1Gt988w2AG6+T8vJy42sjLi4OgYGB1izPpnHCLmp2eXl56Nixo9l1HTp0QF5eXvMWREYbNmzAlClTzK7jW0XzyMrKwgsvvIDDhw/Dy8sLTz75JN58800O1rWyut63Dhw4gPj4+OYtyI4waBAREZFkOEaDiIiIJMOgQURERJJh0CAiIiLJMGgQERGRZBg0iIiISDIMGkRERCQZBg0iIiKSDIMGERERSYZBg4iIiCTDoEHkABYsWACZTGb8CQwMxMCBA3Ho0KFabU+ePImJEyciODgYrq6uuOuuuzBu3Djs27fP2OaXX37BlClT0K1bNzg5OWHEiBEW1bNy5UpER0cbf6+qqsLLL7+MAQMGwMPDAzKZDMXFxbVuFxoaarIfN//8/PPPdd7nH3/8gQkTJiAoKAheXl6Ijo7Gl19+adImJycHsbGx8Pb2xvDhw3H58mWT9f/973/h5+eHCxcumCyvqKiAn58fMjIyLHociBwBgwaRg3Bzc8NPP/2En376CR988AFKSkrwwAMPmFzUbvv27YiOjkZOTg7efPNN7N27F6tWrYKbmxuGDBmCsrIyAEBGRgYOHTqEe++9FyEhIRbVoVar8cYbb2Du3Lkmy9auXQulUon+/fvf9rZfffWVcR9qfgYMGIDAwEDcd999t72dVqvFgw8+iOPHj2PZsmXYtm0bunXrhocffhjfffedsV1SUhJCQ0OxZcsWFBQUYNasWSbbeemllzB79my0a9fOZLmXlxdeeOEFzJs3z6LHgsghCCJq8VJSUoSHh4fJsj/++EPIZDLx3HPPCSGEKCoqEt7e3uKBBx4QWq221jb2798vVCqVEEIIvV5vXB4XFyeGDx/e4Fo++ugj4e/vL3Q6nclyg8EghBDi448/FgDElStX6t3W9evXhYeHh5g+fXqd7X766ScBQBw4cMC4TK/Xi9DQUPHUU08JIYSoqKgQAMTly5eFEEJs3rxZBAYGGtvv3LlThIWFCY1GY/Y+8vLyBABx/PjxeusmciQ8okHkoEJCQhAYGIjc3FwAwNq1a1FeXo733nvP7JVCExIS4O7uDgBwcmr8W8cnn3yC0aNHw8XFxWS5TCazeFvbt2+HSqXCpEmT6myn0+kAAD4+PsZlTk5O8PLyMl6VtqqqCsCNIz8A4O7uDq1Wa1w3c+ZMpKamQqFQmL2PDh06ICYmBhs2bLB4P4haMgYNIgdVXl6OkpISBAcHAwDS09MRHByMiIgIye6zsrIShw8fRmxsbJNsb9OmTQgNDUW/fv3qbNe3b1/cfffd+L//+z/k5ubi2rVrWL58OXJycpCcnAwA8PPzQ1hYGJYvX47S0lKsWbPGOI7kvffeQ1hYGEaPHl3n/fTr1w979uxpkn0jaikYNIgcSHV1Naqrq5GXl4cpU6ZAr9dj/PjxAICLFy9aPN7CUsePH4dOp0NkZOQdb6ukpATff/89JkyYUG9bFxcX7N+/H6WlpQgLC4Ovry9eeeUVbN68GX379jW2++CDD/D222/D398fv/76K1JTU1FUVITFixdj6dKl9d7PPffcg9OnT6OiouJOdo2oRWHQIHIQKpUKcrkccrkcHTt2xIEDB7BixQoMHTrU2KYx3ReWKCoqAgAEBgbe8bb+9a9/QafTYeLEifW2raysxPjx4yGEwFdffYV9+/Zh8uTJmDhxItLT043thgwZgj///BPZ2dnIy8tDZGQkXn75ZSQlJSE8PBwfffQROnXqhLvuuguzZ8+GXq83uZ+AgAAIIXDp0qU73j+ilsKl/iZE1BK4ubnhhx9+gEwmQ0BAANq3b28y1qJt27bIzs6WtAaNRgMAtx3nYIlNmzYhMjISPXr0qLft+vXrcfToUVy4cAEBAQEAgIEDB+Ls2bP4+9//jsOHDxvburm5oWvXrgCAw4cPY+/evThz5gxOnjyJZ599Fvv370fHjh0RGxuLLl264JlnnjHetma/Kisr73j/iFoKHtEgchBOTk6477770KtXL3To0KHWgM74+HhcvHgRp06dkqwGPz8/AMC1a9fuaDv5+fnIyMho0NEMADh9+jTatm1rDBk1oqKicO7cObO3MRgMmDFjBhYtWgRvb28cOHAAERERiI2NRXBwMB566KFa4zFq9svf39/ynSJqoRg0iAgAMHXqVHh7e2PmzJnGszRudvDgQajV6ju6j5ojBTVnujTW559/DgANGp8B3Dgj5MKFC7hy5YrJ8mPHjiE0NNTsbdavXw9nZ2ckJSUZl928/yqVynjGSo28vDz4+PigTZs2DaqLyBGw64SIAABt2rTBxo0b8cgjjyA2NhbPPfccwsLCUFxcjK+//hqfffYZSkpKAABXrlwxjm24cuUKrl+/jq1btwIAhg0bZjwN9lYdO3ZEUFAQjh07hsTERJN13377LVQqFX755RcAwI4dO+Dl5YXu3buje/fuJm03bdqE2NjY2w5e7dy5Mzp06GCczXTixIlYtGgRhg0bhrlz58LLywtbtmzB/v37kZaWVuv2165dw6uvvoodO3YYx63Ex8fjpZdewpIlSxAaGorPP/8cixYtMrndL7/8gn79+t3R6b9ELY51p/EgouZgbsKu2zlx4oSYMGGCaNOmjXBxcRGBgYFi7NixYv/+/cY2Bw4cEADM/uTm5ta5/RdeeEH069ev1vIOHTqY3V5KSopJu1OnTgkAYtWqVbe9jw4dOoi4uDiTZceOHROJiYmidevWwsvLS9x7770iLS3N7O1nzJghkpKSai3/8MMPRfv27YWfn594/vnnTSYdq6qqEn5+fmL9+vV17D2R45EJccuxPyIiCf3nP/9BVFQUzp8/jw4dOli7nCaza9cuTJw4ERcvXoSnp6e1yyGyGTy+R0TNKjIyEqNGjcKyZcusXUqTevfddzF79myGDKJbMGgQUbNbvHixcUbSluD69euIi4vDzJkzrV0Kkc1h1wkRERFJhkc0iIiISDIMGkRERCQZBg0iIiKSDIMGERERSYZBg4iIiCTDoEFERESSYdAgIiIiyTBoEBERkWQYNIiIiEgy/w93mPYWbgoVfwAAAABJRU5ErkJggg==)

*genotype pca*

### Method & Provenance

- **Method:** GWAS — association analysis
- **Engine:** Python · scikit-learn/statsmodels 1.8.0
- **Input datasets:** gallery — example_gwas.csv
- **Started:** 2026-09-21 02:18:38.023727
- **Permanent link:** https://cat-g-omxpert.wishartlab.com/results/by-token/ce9e69bd9b5841ff84c3245d1bf700a6

| Parameter | Value |
|---|---|
| target_column | Trait_value |
| encoding_model | additive |
| significance_threshold | 0.005 |
| n_snps_input | 10 |
| n_snps_after_qc | 10 |
| n_samples | 20 |
| n_significant | 0 |
| lambda_gc | 1.620086 |

**Command history**

```python
# app.services.gwas_qc_service — SNP QC + allele encoding
```

```python
from sklearn.linear_model import LogisticRegression
```

```python
for snp in snp_matrix.columns:  # per-SNP association test
```

```python
    model = LogisticRegression(max_iter=1000).fit(X_snp, y)
```

```python
    p = association_pvalue(model)
```

```python
fdr = statsmodels.stats.multitest.multipletests(p_values, method='fdr_bh')
```