{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "8e2619df",
   "metadata": {},
   "source": [
    "# Sequence Search with Malva\n",
    "\n",
    "This notebook demonstrates how to query arbitrary nucleotide sequences against\n",
    "a Malva index and visualize the results at single-cell resolution.\n",
    "\n",
    "We will:\n",
    "\n",
    "1. Load the clustered PBMC data from the previous step\n",
    "2. Query a CD3E transcript isoform sequence against the index\n",
    "3. Project query results onto the UMAP and compare with gene-level quantification\n",
    "\n",
    "**Prerequisites:** Run `1_run_malva.sh` and `2_simple_analysis.ipynb` first."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8ed107f2",
   "metadata": {},
   "source": [
    "## Load clustered data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "id": "4c9a59f6",
   "metadata": {},
   "outputs": [],
   "source": [
    "import scanpy as sc\n",
    "import matplotlib.pyplot as plt\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import dnaio\n",
    "\n",
    "from malva.indexes import MalvaIndex"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "beaa7c79",
   "metadata": {},
   "outputs": [],
   "source": "# Load the clustered h5ad from the previous notebook\nadata = sc.read_h5ad(\"quant/pbmc_1k_v3/pseudoquant_clustered.h5ad\")"
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "16b57a00",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAV0AAAFxCAYAAADQ5nv6AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAAOLpJREFUeJzt3XecFfW9//HXKXu2sb3B0ot0KTZUBAGxAkZSjBRjRE1i4s8Yo9GbeNUYYyP3agyK10iJPRbEgo1QJIKgsqL03nfZvudsPXV+f2xYXNll2zlz9uy+n4/HPh67M9+Z+Rxd3gzf+X6/YzEMw0BERExhDXcBIiKdiUJXRMRECl0RERMpdEVETKTQFRExkUJXRMRECl0RERMpdEVETKTQFRExkUJXwm7x4sVYLBYOHDjQouPuv/9+LBZLs9paLBbuv//+lhcnEmQKXREREyl0JWLdc889VFdXh7sMkRaxh7sAkday2+3Y7foVlsiiO11plz744APGjRtHfHw8CQkJTJkyha1bt9Zr01Cfrtvt5je/+Q0ZGRkkJCRw5ZVXcuTIkQavcfToUebMmUNWVhbR0dEMGzaMhQsX1muzevVqLBYLr732Gn/+85/p0aMHMTExXHTRRezZsye4H1o6Bd0mSLvzwgsvcN1113HppZfy6KOPUlVVxfz587ngggv46quv6NOnT6PH3njjjbz44ovMnDmT888/n5UrVzJlypST2uXn53PuuedisVi45ZZbyMjI4IMPPuCGG27A5XJx22231Wv/yCOPYLVaueOOO3A6nTz22GPMmjWLDRs2BPnTS4dniITZokWLDMDYv3+/UV5ebiQnJxs33XRTvTbHjh0zkpKS6m2/7777jG//Cm/atMkAjF/+8pf1jp05c6YBGPfdd1/dthtuuMHo1q2bUVRUVK/tNddcYyQlJRlVVVWGYRjGqlWrDMAYMmSI4Xa769r99a9/NQBj8+bNbf780rmoe0HaleXLl1NWVsaMGTMoKiqq+7LZbIwZM4ZVq1Y1euz7778PwK233lpv+3fvWg3D4M0332TatGkYhlHvOpdeeilOp5OcnJx6x1x//fU4HI66n8eNGwfAvn372vJxpRNS94K0K7t37wZg0qRJDe5PTExs9NiDBw9itVrp379/ve2DBg2q93NhYSFlZWU8++yzPPvssw2eq6CgoN7PvXr1qvdzSkoKAKWlpY3WI9IQha60K4FAAKjt1+3atetJ+4MxWuH4NWbPns11113XYJsRI0bU+9lmszXYztDbrqSFFLrSrhy/S83MzGTy5MktOrZ3794EAgH27t1b7+52586d9dodH9ng9/tbfA2RtlKfrrQrl156KYmJiTz00EN4vd6T9hcWFjZ67OWXXw7Ak08+WW/7E088Ue9nm83GD37wA9588022bNnSomuItJXudKVdSUxMZP78+Vx77bWcccYZXHPNNWRkZHDo0CGWLVvG2LFjmTdvXoPHjho1ihkzZvD000/jdDo5//zzWbFiRYPjaR955BFWrVrFmDFjuOmmmxg6dCglJSXk5OTwr3/9i5KSklB/VOmkFLrS7sycOZPs7GweeeQR5s6di9vtpnv37owbN47rr7/+lMcuXLiQjIwMXnrpJZYuXcqkSZNYtmwZPXv2rNcuKyuLzz//nAceeIAlS5bw9NNPk5aWxrBhw3j00UdD+fGkk7MYehIgImIa9emKiJhIoSsiYiKFroiIiRS6IiImUuiKiJgoLEPGAoEAubm5JCQkNPsdVyIi7ZlhGJSXl5OdnY3V2vj9bFhCNzc396RxkyIiHcHhw4fp0aNHo/vDEroJCQlAbXGnWjVKRCRSuFwuevbsWZdvjQlL6B7vUkhMTFToikiH0lSXqR6kiYiYSKErImIiha6IiIkUuiIiJlLoioiYSKErImIiha6IiIkUutL+BHxQVQw+d7grEQk6ha6E1Zo1a7jiiivIyMjAYrFgsVh45o7vw7q/wOo/whfP4q2u4I9//CP9+vXD4XDQo0cPfvOb31BRURHu8kVaTKErYZWTk8Py5ctJTU09sdHv+c83Bjj3M2faudx///0cPHiQfv36UVBQwBNPPMHUqVMJBAJhqVuktRS6ElbXXnstLpeLjz76qMH9ObuP8eKKrQD89eaL2LFmCW+++SYAn3zyCUuXLjWrVJGgUOhKWKWlpREbG9vo/g8+31f3/Q/GDYLNLzFlyhRiYmIA+PDDD0Neo0gwKXSlXTtc6Kr7PjM5DowA1rIDpKenA3Do0KFwlSbSKgpdiTzuMgzDCHcVIq2i0JV2rWfGiaU/C8qqAAikDqK4uBiAXr16haUukdZS6Eq7dtnZ/eq+f/PfOyHzdJZ9vJKampra/ZddFq7SRFrFYoTh32kul4ukpCScTqcWMe/klixZwu9+9zt8Ph8HDx4EICM9jcQ4B2OG9uKlOy9h5oNv8sqqbVitVgYOHMjevXvxer2MGzeO1atXn/J9VCJmaW6u6bdVwsrlcrF37966wAUoLCpm76E8jlbHwKQH+MdHm7j33nvp1asXe/fuJSMjg1tvvZVly5YpcCXi6E5XRCQIdKcrItIOKXRFREyk0BURMZFCV0TERApdERETKXRFREyk0BURMZFCV0TERApdERETKXRFREyk0BURMZFCV0TanYDPwO/umAvV28NdgIiIx2XgLgN7nMHmv8PhVQYBL6QOg/hMKD8CMSnQb5qV7hdYwl1umyh0RSRsfNUGGx8PcHglGP6T95dsrf067tjnAYbfYGHI7Mj9R3rkVi4iEW/9nwIcWt5w4DZm2z8MakoCoSsqxHSnKyKmKvzG4PAqA9chg8Kclh8f8MGKXxlMnm8QnRx5XQ0KXRExzaan/Ox+o+3nqToG3/xfgLPvsrX9ZCZT94KIhJy3yuDzR4ITuMcd+BAOLo+8bgbd6YpISPmqDVbeEsC1P/jn/uIRg6R+Bsn9I6ebQXe6IhJSWxaEJnABjADsXxZZ43kVuiISUvuWhfb8NSWhPX+wKXRFJGTK9hj4a0J7jS7doSDHwFMRGXe86tMVkZBpyfjb1rBGwY6XDXa8bGCLgUHXWBh2Xfu+l2zf1YlIREs+DeKyWnfs/269hqtX27l6tZ0nts5ssE3Ae+J7fw1sW2xw8OP2PaJBoSsiIWOxWjjz9pbHzKq8xawvbN34sp3/NHDub79dDQpdEQmprLPB0oI5DMeq97Joz20MTDyXtOgeLb6ecx98PCfAR9f7Kd3Z/sJXoSsiIWWxWMge27y2/oCPv227FgtWbh3yAlZaP+PMdQBW/TqA62D76m5Q6IpIyI34mZWY9KbbvX7wAXaXf86NA+eRGdu3zdf1u+HjOUa7mrmm0BWRkOvS3cKlC6wM/Wnjbfa6vmTpwUcZlzWLcVkNPzhrDSMAXz5m4Ha2j64GDRkTEVM4Ei3Ed7UADYff4cqtBPCzvvBNPi9cCoAnUAXAhqIlXLsmif87/xBx9qQWXzvgg6NrDNLPLWT/po1ERUfT/6xziY6La+3HaTWFroiYpqa46TbewMmzKfyGD7/hwzBaf7e6/eMd7HzpbgyjtqshKiaWKbfeQf8zx7T6nK2h7gURMU3q0MYXppnQ7Tpem+Cr95UR3RuA8zOu5rUJPuKjklt1XQODfUcW1AUugLemmvf++hjVFeWtOmdrKXRFxDSZoyykjzD/us74F/E4dp603ed2s3v9WlNrUfeCiJhq/GNWPrwuQFV+022fOm9vm69X3OVZKuMaX3WnqtyJu6qS6Lj4Nl+rOSxGWzpJWsnlcpGUlITT6SQxMdHsy4tImJUfNvjk9gDVRaG9To19BwWpd52yjcVqxQgE6D54GBOvu4msfgNada3m5pq6F0TEdAk9LVz2gpUz74S4rqG5hse2H1f86022MwK1/bxHd2zln/ffjauoMDQF/YfudEUkrAzDIHdt7at3DB/YYuHwihYcT4CA1YklEIeVaADKY5dR2uXvYGl5vFmsNs6aNp1zvvdDYuK7NPu45uaaQldE2oWA38/mlR+z5V9rKNkZIMP5e2ycOvQqY1ZRFv8Kfls+lkAM8TUXkVh5FbnpN4PF16Z6Mnr3ZcaDfyHKEd2s9s3NNT1IE5GwCAT85O/dg2EYdB1wGsv++hi7NvxnJEE0FKTcQ7fSJxo81uqACjZQnHhiv2GtoSJuGb64PUDbAheg8OB+tn+6mhGTLm3zub5NoSsipju4eRMfPv0EFSWNP0nzRu3nWPLviHOPJb5mEjYjAYAB0y30uKSSl/74ODSwpEINJw8Na61NHy1T6IpIZKsoKWbp3D/hc7ubbOtx7MTj2Ell6uucMfIBhlwxgLRhFt7+y1/xB6pPeexHW3axfNvuBvc9+sPLsVmbHkdQknsEIxDA0oy2zaXQFRFTbf1kRbMC99u83nIOlDyFddvZbH16Ba6igmYfGx/tIC2+/hoLFpr3yna/x4PX48YRE9uiek9FoSsipiovacYCDA3I37eH/H17WnzckG6ZXHPOyFZdM6ZLl6AGLih0RcRk2acN4uuPQ/xe9m/ZfCSPTYdziY2KokdKEpcNH0j3lOatVDbiosuCXo8mR4iIqQaeN46kzBDNiPgOq8VCQkw0qXGxlNe42Z5XwN9WrONoqbPJYzP79GfczJ8GvSaN0xUR01WXu/jHnbdQWVoSsmsUllcQ73AQF+0AYOexQv6+5nMAzunbk6vPbnzlnfiUVG6atxCbvfmdAZoGLCLtVmxCIr945nn6jj47ZNfISOhSF7gAg7pmEOeIAqCs6tQjH7qkpLYocFtCoSsiYfP9u+/j4p/9P2xRUS0+tqlhXCu376W08kS47jpWSJXHC0BK/Kkfjg06b1yL62muTvUgze1xUVNTQlxsJlFRcfj9XvwBN46o5s+vFpHgGnHRpQwZN4EDmzayc/1adn3277pFaBpjtdq45OZb+fCpxxtt89neg3yweQdJcbE47DYKXRUAOOw2xp3W+Esve48YzejLr2zdh2mGDt2nW1i8md3738VVcQiMAAGjdmqgxWLFYrETCHgAiIvtxuD+36dr5pkhq0VEmsdbU8Pa115k65oVeKqqCfh9gAWL1YIRCJDZtz/jZlxHn5Fn8MmLC/ny3SUNnmf93kN8fSSPfGc5VR4vCbHR9E1LYfLQ08hMPPlGa9iEyQwZO4Fep4/EYmneON5v6/QL3uQVfMlXW55ucF/+sXLmzHqj0WPvuPPnPPjgY0Q79JBPJNz8Ph+1L7O04Pd6cMTWn+hwaMvXfPbmqxzZtrlV57fabFx47Y2ccfm0NtXZ6Re82b771Ub3RUXZGDQko962ygoPRw7XDiOp8X7DqrV30KfnxQwe8KOQ1ikip/btB1oNPdzqNXwkvYaPZP9XX/LBU/9LdbmrRee/4tY7GXTuBW2us7k6ZOh6fVXUuBsfipKaFsf/zptab9v8Jz/jyGEnXRIcTJjcn4DhY9+hD4iPy6Jn9vhQlywibdR39FlMmvMLlv31sRYdF/C1fUWyluiQoxds1mgsluZ/NJezhuUf1U4vvGLaYGJjTzxJPXR0dbDLE5EQGThmLCnZPVp0TFqPXiGqpmEdMnStVhvJSc1/z9Gyd3bgrvERFWVj2vQh9fbVuMuCXJ2IhIrVZuPqex9i0HnjmrUyWJ9RZ5LZp58JlZ3QIbsXAM4Y9ktWr78bv7/mlO28Hj/vvb0DgImT+5GaWr+T3jBOPXRFRNqXLimpTL3tLnxeL1Vlpez54jPyD+wjd+d2KkqK8Hk82KOjGTZ+EuNnzzG9vg4butHRiYwf82dWrbuD2iefDVuxfA9lpdVYLDD9R8NP2u/xunB7yol2JISwWhEJNntUFIkZmZxxxffqbfdUV2GLcoRsxllTOmT3wnGxMSl073p+o/sNw+Ct17cCcPaYnvTqndxgu/yinFCUJyJh4IiNC1vgQgcPXYChp80gNXlQg/s2fHa4bpjYD3588l3ucW63i6KSbXi9VSGpUUQ6jw7bvXBcVFQc555xF6XOvVRU5mLBwqHcTyhz7WXJa1sAGDg4neEjGl9qbvf+twCwYKNPr0sY0g7G7tb4A3xWVoHVYiHaAuvKKukWE8XUjGQS7LZwlycijeiwM9KaUl1TO4535943yM1f36JjkxP7c86o32K3x4SitHq+cFbyP/uPsT6/EM8Lz+Jbuxpn/jF8VhvWrG7ETJ5C3I9/UjdtMc5q4fkR/bggRX3QImbq9NOAm6vGXcr6nEepqm7+O5cA0lKGEBebiddbSVrKELp3Ox+7LTqotW1yVfG9r3bjDhg4H72Xmo/eBcDWpz9GZQWBwnwAEm75HXHfn1F3XLLdxjdjh+EI4sv0ROTUtJ5uM8VEpzDunAcYMWQOPbIvxGZtXnAWl27ncO4nHCv8kq27XmD9xofx+oLb5/u3Q/m4A7V/J3o3bwLAcc75pC98g/Tnl4KjtlZ/fl6948p8flaXlAe1FhEJjk4fugA2m4Me3S5gxODrcLRyaJir4hAHD68Ial2ri08EZ9TpowHwfL6Oojk/pOgnV4HHTdTpo4m7+tqTjq3ya3yxSHvU4R+ktVRW+igOHPlXs9s7jlQSvdNJ3LYy3Idz2Jv2FEk/+AExQ4fgfPsdfKUlJEyaRMqMGfWWizMMg8pPP6X662+I6ppFwmWXY+sSX7d/e14+lX4//OeYxN/8AYwANR+/h//A3tpGUVHY+w3EmlD/nzJWYHyq+nRF2qNO36f7XW6Pi/U5j1BZdQyO/6dpYG1Na5WPzOd2Ebuj6RfcAdjS0uj79lKi0tMJVFZy6Gc/p3rjxhP7U1Lo+fe/EzPwNI489DBXjLiAwpS0uv2VLy+kYsE8ooaNJPmB/yFQVkrJbTdiuMqI+/5MEm65s67tTd3T+dPAls0/F5G26fRLO7ZWtCORsWf9NwdeeIKKhf/EXu6jcmQK5RO74ekRj73UTZd/5xO7w0nMocpmn9dfXMyRn/+cvm++SdH8+fUCF8BfWsqhG24gKiuL5XFJFF54YhaNUVNNxaL5YBhEj7sIa3Iq1uRUHMNH4l73Ce6cDSQANuC/+nXjlt5ZQfqvISLBptBtgP9QHt7HXiD6P3e6iRuKSNxQ1Obz1mzdxrEHH6TsjTcb3B9wOnE7nWybPrPedsNdA/7a5ed8u7fXbvO48R3YB4AlJpYYC3w9djhJUfpfKtKe6U/od1Rv2Urevfee6FoIstIXX2qyTbqztN7P1qQUokacgfebHGr+9T7e7ZsxqqoIlBYDkHbFlTw/or8CVyQCaPTCt5S89BIHfvQj3Nu2hbWOi9f/m7jq+sPPkv/0OHHX/BRHz95ElRaTSICBZ5zFXfOf5eDcP+nBmUiE6PQP0jwHD1Lx708xjAAFj80Frzes9Rz3zYBBPHzdLzmWnglASsDHraf14mc9M7C14qV5IhJampHWDPlz51KycFHIuhKOsyQkYJS3fLJCwGZj37iJdLnqe4y75CLsVoWtSHul0QtNKF+1ipIFC025llFejjU9nUBR8x/GJU6dSvZjjzJMU3lFOpROG7rOd95pss3zJSW85XKS6/VSYxik2myMjInl5rQ0BsW0bLGblgRu8owZdLvv3hadX0QiQ6e9jQpUNb1OwhfVVZT4/PSIiqJnVBSFPh8fV5Tz08OHqAqEZpqtY8AAsu76XUjOLSLh16HvdGt27cL1/vsYHg/2zEwsFguOPn2IHzeO6H79qfxkzSmP/0u3bKK/9c/7J4sKeaa4GGcgwH6Ph2EtvNs9idVa++XzYYmPJ+l736Pr3XdhcTjadl4Rabc6ZOgGqqvJu+9+XI10IdjSUvGXlDa479uirVb+VV7OcyXFVAQCHPB4AEi12ejjiGri6KY5+ven7+uvEaiowJaWVm9tBhHpmDpc6Ja+9hr5jzyKcYruA39xSbPPV+T38U3NiTcK94iK4qnuPYi3tv3tDDGDB2ONicHa1jtmEYkYHaZPt3rrVvZOu5Jj9953ysBtqWuSU9g6cBD/6tefyxMSOOL18tvco1QG/G0+d/L06UGoUEQiSYcIXb/LxaE5N+DZvTsk57dYLGRHRfGz1NpVv/Z4PCxztW2R8KiePYk//7xglCciEaRDhG7ePf9NwNm8JRabq8zv5x2nE8+3Jk6sqTyxqlh1G0YvWGJi6PbAH9tUn4hEpojv063Zto3yjz8O+nkrA37uPpbH/fnH6BnloDzg55ivdqWveKuVyQmtW+sg+eqrSZtzPY4+fYJYrYhEiogPXefbJ0YoOBK9BLwWfNVt/1gJVhtXJCSwuaaGw14PPsOgq93OuOwofj/ewWldnJTti8N1MK7Z5+wy+SLd4Yp0chEfuv7KCqz2ANnnlZLQ3Y0RgLL9cZTujsNXZSO+q5v4rjXkb0oi4Gn+iINEm42/ZHev+9liNeh1USFxaT7AADzEZ3mISfZS8HXSKc9lsQeIG9id7F9c1boPKSIdRuSG7qEN8NULZGRsJ3VyETHJtf/0t1ghpX8VKf3rj2CIz/KSuyGJqoJooOXjYTNGuP4TuPWlDqqkZFcXfNWNB3pyvyq6jtgIb/wAdv0YrnqmdlKEiHQ6kRm6G56FD2rfCRYFRCU3fUhUvJ/ek0o4lpNA6a4EEntV0aV7DTXFDrAY2GMDFG/vgt9dG56x6W6yRruwx/qpKbPj6NLwEDGLFeIy3bgOxtJgmFuM+n8BfPNPOO0SOP2HLfvMItIhRF7oVpfB8tYvBtP1jHISetQQl+HFYoGk3icmPsSmeTi4IgNrVIBeE0qw2mtHLkTFeU55Tr/bStZoF9FJXnI3pNTd9docfrLOdBKd9J075JUPwpArwa7pviKdTeSF7v414Ktu0yniMxteqDwuw0tsuptuY8rqArcpvhorKQMqSejhBmDAtHyqixwYAYjN8NDgxLXS/fD3SXD9+xDTvt6GLCKhFTkdixUF8MHddd0KodL1rDKiE5o/28weE6gLXDje3eAhvmsjgXtc/mZY97c2VCoikSgy7nSry2DhpVCyL+SXiklu+/TeZtv2Nkz6AxTvhS+eg8KdkDEIzr4R0vqbV4eImCYyQverF0IWuGsO+njkUw9f5PopqqrtUpg/JYZfnGVCf6vFAofWwwvfB+9/ZrvtXQEb/wHXLoFe54a+BhExVWR0LxxaH7JT5+T5Wb7PR2psGJZVHHIlfPT7E4F7nLcSPvqD+fWISMhFRujGZ4Ts1NeOiMJ1dwIfzW7+zLKg6H4mjL4Wjm5seP/RL6Gq+UtQikhkiIzQPfO6kJ06Lc5KbJTJd7kX/hfM+Qji08DaSA+P1Q72aHPrEpGQi4zQTexBa2aRtT9WmPgHmHg32KLAEQ+DpzTcdPCU2v0i0qFExoO0qmJq1zuIcDf+C3qcWX/b5XNrRy/kbzmxLXNY7XYR6XAiI3TT+kNsGlQXh7uStrE18F61hCz4+b9hz3Io2gXpA2HAxVqbQaSDiozQPfQZeIP3Cp6wsNig6+kN77NaYeCltV/SYrvzy3nu3/tZv7+YKo+fxBg7EwdlcsGAdJw1XmxWCxU1PrKTY7lgQDpWq4WSSg+lVR56p8axv6iS3QUV9E6No7TKi8fv59x+acQ5IuOPh0QWi2EYpv+73eVykZSUhNPpJDGxiWmw+/8Nz38PjNBMWliy3cvvltfgC8BBZ+1/iow4C4nRMKaHjZe+H6RRDYnd4fZtwTmXALD5iJMFa/fz3te5+ALN+zXumx5Pn/Q41uwqwh8wsFstDR4bY7dy60WncfOE/npLszRLc3Ot/f9VvmZuyAIXwOU22Fta/w9dYZVBYRX0SAzi30djbwveuYS5H+3gqVV7W3zc/qJK9hedGBfdWFjX+AI89tFOnv/sIH+6ahgXD+3a6lpFvq393+k+3AvcwX3/meni0uD27RoCFiSbjziZNu9TU6+ZmRDNPVOHcOXI7g3u9/oDrN5ZSGG5m7P6pDAwq3Wvc5LI1XHudJN6QEGEh25VMax+GCbfH+5KOoT3t+SZfs2Ccje3vrKJf6w7wC2TBjBxUFbdvl355Vy/6AuOlp1Y/e57o7L5nx+NxG6z4qz2snJHPv4ATBiUQXoX/eXbmbX/0D3nJnjvtnBX0XY5zyt0g2TzkfD9JbzxYBnXL/qSrokxLP3VWLomxfCrl3LqBS7A25tyARh/Wga/f2szbl/t26NtVgt3XTaIn43XgkadVfsP3bOur13Wcd2T4KkIdzWtV10a7goilrPKy+sbD7Mt18Xh0iq+OBD+/5bHXDVc9dSnPDXrTHYXNPx7+fam3LrwPc4fMHjo/R0M7ZaIw27jzY1HKHd7GTsgnR+c0YOYqOa/x08iU/sPXYAJd8F5v4TCXZD7FXx4FwROfl9Zu9bngnBXEJF2Hitnxt/XU1J56rd3hMMxl5tPdxe26thrF37Ot5+mvL/5GK99cZiXbzqX+OjI+GMprRM5I/CjE2pnc51zY+1kgphTv4G3XbHHwqT/DncVEedP723lir+uaZeBe9zHW/NbdVxDj6+/PuLkpQ0H21iRtHeRE7rfljW0/YZYTDL85G0YNKV2Ou/IGbXTf3ueE+7KIsarGw4x5L8/ZMGnB/C389nfW/NcQT3f0q9ym24kES1y/x1zzk21Sx9+Ng/crtpVuQZPqX0bQ7hYrPCDBdBvQu2XNFsgYLB001EefG8rJVUR1nUURAeLK5tuJBEtckMXavt6z7+l9q0SCd0gPh1evgZ2fRDa6/abVLtmQvZoKNoDBVsg7TQY84vau3Bpthqvn02Hy/jbit2s3Rvha2sEQaXHz77CCvpldAl3KRIikR26ULv84bfXNLjySXj5x5Cb0/pzWh0Qnwrlx+pvT+lT261x+g9bf26p88rnh3jkgx04qxt+O3Nn9diHO3h61plYrZp+3BFFfuh+V5dM+Nmq2lf8fP0q7HwfKpr5sKNLFlz6UG3XQHw6HP689iuha+3KX7ER9PCunduwr5jfL9ncERbsDLoPt+azcO1+bhzXL9ylSAi0/2nAbRXwg/MwlOfDxkWQuwmiYiA+CwhAZSH4fdB/Apx3S23ASkgdLavm6mc+O2lCQahVbvsE1+dL8BYfxmJ3ENNrBMkTricqpZupdTRH9+RY1t49KdxlSAt0nGnAbWW11XYLpPSBXmPCXU2n9+7Xufz61a9o5qJgQVP+9ceUfPgkAPakLPzV5VTtWkfNkW1kX/83bF1SzC2oCXlOc/9CEvO0eMhYdXU1n376Kdu2nbxMYU1NDc8//3xQCpOOwzAM3t50lB8+s47/94r5gWv4vZR9shiAuIHn0/0XC8i+cT4WRyyBqjKc618zt6BmCBhQ6e68ozg6shaF7q5duxgyZAjjx4/n9NNP58ILLyQv78TiI06nk+uvvz7oRUpku/ftrfz61U18Gabpu+683QSqa8fTxg06HwB7QhrR2YMBqN7XyBuZw8ymB2kdUotC96677mL48OEUFBSwc+dOEhISGDt2LIcOHQpVfRLh9hSU88L68M6y8rtOTNW1xiXXfW+Lr/3e52rdVN5QGtUzWeswdFAtCt1169bx8MMPk56ezoABA3j33Xe59NJLGTduHPv27QtVjRLBPt1dFO4SGmf+M+RmSYix83+zzwh3GRIiLQrd6upq7PYTz94sFgvz589n2rRpXHjhhezatSvoBUpkO1gS/nfb2RIz6r4PVJ1YFtL/n+/t39ofTj2SY7jrssF89d8Xk5UUG+5yJERaNHph8ODBfPnllwwZMqTe9nnz5gFw5ZVXBq8yiXhVHh9vbDwS7jKI7nYa1thEAtUuqnauJX7ohfjKi3Hn7gAgtt+ZYa4QRvRI4u1fjdX72DqBFoXu9OnTeeWVV7j22mtP2jdv3jwCgQDPPPNM0IqTyPbJzkLKa8L/BN5iiyJ5/E8o+WgeVbvWcfSZG/BXl2N4qrHGJpI4xvwZhpkJ0QzrnojNYmHWmN5MHJxpeg0SHh1/coSEzdubjvLrVzeFu4w6FVtX4fr8rbrJEbG9R5J84XVEpTb83rNgio+2MXFQJnarhUuGdeXSYV01OqGDCdnkiAMHDrB8+XI8Hg8XXnghw4cPb1Oh0nGNPy2DaLu17lU14dZl2ES6DJtoyrVS4x116wCf0zeVe6cOZXh3TSOXFobuqlWrmDp1KtXVtbNl7HY7CxcuZPbs2SEpTiJbSryDP0wZwn3vbG2vAwVCIjMhmg2/v4hyt49AwCA5zhHukqQdaVH3wgUXXEB6ejrz588nJiaGe+65h7feeovc3JYtvKzuhc5l8xEnr315mPX7iht9n1hHYbXAZ/91EVmJMeEuRUzW3FxrUegmJyezbt06hg6tXTO2qqqKxMRE8vPzSUtLC3px0vH8+P8+Y8P+knCXETIvzjmHCwa2jyFoYq7m5lqLxum6XC7S09Prfo6LiyM2NhanM3yvxJbIctnwjruKW9fEGAWuNKnFD9I++ugjkpJOPBAIBAKsWLGCLVu21G3TeF1pzJTTu/HQ+9vxtveXn51CrMNGtcdfb5vdauHpWZpFJk1rUfeC1dr0jbHFYsHv95+yjboXOp+8ihJe2L6KTQX7OVbm4eCenvirUxtpbWDtcgxrbAmGPxq/qzv42s8MrXP7pXLVqO68sP4g+a4aRvVM4ZZJAxjVMzncpUkYhWTIWCDQPob+SGTZV3aMme//hTL3iZcu2nsexlqZji9/JMa3A9Xiw9FjA9a4EyuS2dN24c09k0BllpllN+q68/pw+enduOacXuEuRSJQUF/BHggEeO+994J5SukAnvzq3XqBC2CxgK1LEVjrvx/NlrqvXuACWKwBorp+DZj/l35mQjTHZ+Ymx0Vx37ShXH56+3vThESOoLw5Ys+ePSxcuJDFixdTWFiI16sXDcoJa45sbXSftUs+/pIT/xSzJeQ12M5i92CNKyZQFfwHVRYgxm7FZrOQ3iWagAFdou18/4zu/PT8PhSUuymp9DAgs4uWW5Q2a3XoVldX8/rrr/Pcc8+xdu1axo0bx7333sv06dODWZ9EsL1lefz2k4VU+dyNtrGn7iNQkYXhaU7ffmjudE/vkcQ7t1zQ6P7s5Fiyk9tPn7JEthaH7hdffMFzzz3Hq6++Sv/+/Zk1axbr1q3j6aefrhu/K+IL+PnZ8qc4WlF8ynYWmxdH77X4nT0JVKVheGMgurzBtrbk/QSq08E49d3mmcNLydlla2aQQ/+MLs1qJxIMLQrdESNG4HK5mDlzJuvWrWPYsGEA3H333SEpTiKLL+Dn3rUv8eGBHKp97ma/Xt1i9WNPOQApB07ZzpZQhDX+YwLlXfEWDQRf/MmNrF62uD8jum8Ab/5w/GV9mrz+FPXRiola9CBt586djB8/nokTJ+quVuoxDIPvv/MQS/Z8RlULArelLFY/tqSjOLpvBMvJy0ba03dhsdZ2Q0RnbSMmuvb5Qs/UWCYOOrk/+MqR2Vw0RMsqinladKe7b98+Fi9ezM0330x1dTUzZsxg1qxZWni5E3O6K3lx22q+KdzPrtKWrcHRFtYYF47en+Ir6Y9Rk4zFXoMt+QC2hPy6NgEC3H9NOpO6n0V6l2isVguf7y/hvW9y8foDXDK0KxMGZej3V0zV6vV0V65cycKFC1myZAk1NTXccccd3HjjjQwcOLDJYzU5IvJVeGt4aMNrvLX7s5Dd1QbDwxf8hOmnnRfuMqQTCMnaC982adIkXnzxRfLy8pg3bx4rV65k8ODBjBgxorWnlAixIW8XE//5e5a088C1AIVVrnCXIVJPUN8csWnTJhYuXMiTTz55yna6041cHr+Xia/9geKahkcYNFfJh19T8fUBPMecBCprsCXFETcom7Qrz8KRGdzfiZtHXs6vz9B6IBJaIb/TbcioUaOaDFyJbGtzd7Q5cAFKV2ymelcetjgH9uR4fMUVuNbt4tDDS/FXe4JQ6Qn/982HbC48ENRzirRWix6kTZo0qck2FouFFStWtLogad+qvY1PdGiJ5PFDSDx/IFFpCQAUvLKW0uWb8TurqNp2lIQz+wblOgABw+Cm5fP41w//RBeHJjlIeLUodFevXk3v3r2ZMmUKUVFRoapJ2rFzswfhsNnx+Nv2lt+0afVfex47sBulyzcDYIlq+B9gVouFQCt7w8rclbyz93NmDrmwVceLBEuLQvfRRx9l0aJFvP7668yaNYs5c+boxZSdTGpMAr8efSVzv1wStHMagQBln2wHICojkbghPertj7HYWTL9DyRHx/POns+Z++US/EbLpwTvd+Y33UgkxFrUp3vnnXeybds2li5dSnl5OWPHjuWcc87hmWeeweXSU+LO4obTL+bqgY2vVdASAbeXo/M+omrLYWxJcXS/9TKs31lUJiEmjr6JWaTGJPDT4RfxsxGXtupa/ZM180zCr1UP0s477zz+/ve/k5eXx69+9SsWLlxIdna2grcTufOs6dgsbXsO63NWceiRd6jcdJCorCR6/ddVRHc/eWHzwmoX6/N21v2c5Ghg+m8TusWnMK3/2W2qVyQY2vSnJicnh08++YTt27czfPhw9fN2IgnRcdw88vJWH+8+WsLBB5fgPlhI7MBu9P7D9FMOFdtWfKjue5enqkXXSnbE8ftzfkR8lN7QK+HX4lXGcnNzWbx4MYsXL8blcjF79mw2bNigtRg6oVtGTyXJEc9Dn7+O0cJpEkfnfYSvuPZ17IEaD0eeeL9uX9L4ISSPH1KvfUbsiffyJUXHtehaZZ4qfrP6OZ6YeBOTe49q0bEiwdai0L3iiitYtWoVl1xyCXPnzmXKlCnY7UFZB10i1LXDJuIJeJn75VstOs7wnXiPnvtQ/eUf44f3PKn9/2xcys7So6w5spUan6fFIxl8RoC5Xyzhol4jtdaChFWLX0zZrVs3MjMzT/mLm5OTc8rzaEZax/PittUs3LKc3MoS7FYbvsCpX07aEIfVjifQtqFoTVn+wwfomaDXpEvwheTFlPfee6/uEqRBs4dOYPbQCXj8PkrdFTyy4Q2WH9yEz2g6fBMdcdx/3gz+9tW77HcVhLTO361ZzOjMfsweMpHsLo29jVgkdIK69kJz6U63c/D4fQSMANU+D2/tXs/6vB1sLNhLpbemrk1GbCKLLruNAcndePSLN1i0xZzZjMnR8bx0xW81jEyCprm51qLQTUlJafBONykpiYEDB3LHHXdw8cUXB6046XgqvTW8t+8L9pTm0Scpk+/1H1M3Nbe42sWUtx446c3BoXJZnzN4YuJNplxLOr6QdC888cQTDW4vKytj48aNTJ06lTfeeINp06a1qFjpPOKjYvjxoHEN7kuLTeTNK/+LP657lTVHt4S8ln8f3Rbya4h8V4tC97rrrjvl/lGjRvHwww8rdKXVundJ49lLfsVL21bz589fa/VaC83RReN2JQyCurTj1KlT2bFjRzBPKZ3UrKETWHvNY/RNzArZNa4acG7Izi3SmKCGrtvtxuFwBPOU0omlxHTh7av+wG1nXEm/pK5kxCQSbQvOrMfx3Ye1aUadSGsFdWbDggULGDVqVDBPKZ2cwxbFL0Zezi/+E5DVPg/v7N3AxmN7qPF7sVmsxNij2Oc8xjeFB5ucGTem60BuP+sqRmYEb71ekZZoUejefvvtDW53Op3k5OSwa9cu1qxZE5TCRBoSa3fw40HjGnwYV1Dl5Eh5EY/nvM0Xx3bX25cWk8AvR13BrCETTKpUpGEtCt2vvvqqwe2JiYlcfPHFLFmyhL59dQch4ZEZl0RmXBIvXH47O0qOcLS8mEGpPeiRkBbu0kTqaHKEiEgQhOXFlCIicmoKXREREyl0RURMpNAVETGRQldExEQKXREREyl0RURMpNAVETGRQldExEQKXREREyl0RURMpNAVETGRQldExEQKXREREyl0RURMpNAVETGRQldExEQKXREREyl0RURMpNAVETGRQldExEQKXREREyl0RURMpNAVETGRQldExEQKXREREyl0RURMpNAVETGRQldExEQKXREREyl0RURMpNAVETGRQldExEQKXRERE9nDXYBIqHiPVeI5WoEtOZrofklYLJZwlySi0JWOx/AFKH5lBzVbi+u2RXWLJ+2nw7AnRYexMhF1L0gHVPrm7nqBC+DNq6T09V1hqkjkBIWudCi+omqqvipocJ97Txnuwy6TKxKpT6ErHYrzowOn3F/41NcUzN9EwOM3pyCR71CfrnQIht+g5OXtVH+nW6EhnoPl5P15AzGDU4k/uysxA5JDX6DIf+hOVzqEyg15zQrc4wy3n+qvCyl6bjPlnx4NYWUi9elOVzqEii/yTtpWXFXGE2sXs3zPWgoqiukSHc/QzP48etnv6J2cXdfOuWwf8WdmYY3VHwcJPf2WScQy/AF8xTWUvrkbX15VvX0lVWVc+fzPOeTMw2GLom9qTzAMNh7dSn5FUb3QxYCSN3aSfu0wkz+BdEYKXYk4hjeA88P9VH6Zj+Fu+IHY3H8/xyFnHgPT+/Lyj/+HrC7pAHj8XgzDOKl9zdYSCp77hqRJvYnulxTS+qVzU+hKxCl5cxfVmwob3W8YBu/uWAVAdkImM//5Ww478+iT3J1fnjuLq4ZObvA4zx4nhXu+IeVHA4k/MysktYvoQZpEFF9JDdVfNx64UNuX66wpB2D1/g24aipIik5ge+Fe/t+7D7Bsx+pTHl/2zl4MXyBYJYvUo9CViOIrrIKTewfqtwmc6HI4La03a3/xKmt/8SqnpfUGYHHOklMeb7j9VH1z6mAXaS2FrkQUa6KjyTZpcck4bFEADMkcgMMWhcMWxZDMAQAccR5r8hw1u8vaVKdIYxS6ElF8BdVNtomy2RnTYyQAOwr24vX78Pp97CjYC0DflB5NX8iv7gUJDT1Ik8jSzNUZ7xx/Ixte/ppdxQc4/5kfA3CsohCbxcYt581u8vjovhrBIKGhO12JKDEDU8De9K/t6OyhvHrNE5zXazROdzluv4dxvc/irdlPcX7vM055rD0rjjiNXpAQsRgNDVoMMZfLRVJSEk6nk8TERLMvLxGucmM+pW/savKBWmvEnp5OyvQBWOOign9y6dCam2vqXpCIE39mFo5eCVTlFFCzoxjvd2ajNcWWGk3mr0bjza2g4tOj+IqqsWfFkzC+O9F91K0goaU7XYlohi/Asblf4ne6m9Xemuwg8+ZReoOEBJ3udKVTsNitZP16NK4Vh6j6pggsED0wBavVQs3uUgJuP1gt2BMdxJ6eTpdzs7WwjYSVfvsk4lnjokie1p/kaf3DXYpIkzR6QUTERApdERETKXRFREyk0BURMZFCV0TERApdERETKXRFREyk0BURMZFCV0TERApdERETKXRFREyk0BURMZFCV0TERApdERETKXRFREyk0BURMZFCV0TERApdERETKXRFREyk0BURMZFCV0TERApdERETKXRFREyk0BURMZFCV0TERApdERETKXRFREyk0BURMZFCV0TERApdERETKXRFREyk0BURMZFCV0TERApdERETKXRFREyk0BURMZFCV0TERApdERETKXRFREyk0BURMZFCV0TERApdERETKXRFREyk0BURMZFCV0TERApdERETKXRFREyk0BURMZFCV0TERApdERETKXRNVlFRwX333cfgwYOJjY2la7dsfnL9zzh6rCjcpYmICezhLqCzmTZtGqtXr8Zms9Gz7yDycw/xwuK/s27D57y8dCVnn5aCxWIJd5kiEiK60zXR1q1bWb16NQBzbvsT//vCSuYu/hiAvdu/5pVX/snuvKowVigioabQNUmB082/NhXW/Wyx1v6nt1pO/C/45os1HCyoNr02ETGPuhdMUOPx89mOMrr3GUivfoM5tG8Hz/3vH/h46QsU5B6qa1dSdIwajz+MlYpIqOlO1wQHCqoJGGCz2bjn8ZcZf+kPSExKJf/oQYaMHEPX7n0AsNnspCc6wlusiISU7nRNUFnjq/s+PTOb2+5/qu5nj7uGOVNHANC9d3+G9Oxien0iYh7d6ZogKT6q7vu9O7+hurICAL/fzz/mPUBVhQuAW3/2E5K/1VZEOh7d6ZqgV0YsO49W4vEZrHz3FVa8+wpde/ShrKQAV1kJALfddhsXTTg/zJWKSKjpTtcEDruVC4enkpHoYMDQ0WR170V+7kE87mrOPPNMnnvuOR5//PFwlykiJrAYhmGYfVGXy0VSUhJOp5PExESzLx9WAcPAApoAIdLBNDfX1L1gMqvCVqRTU/eCiIiJFLoiIiZS6IqImEihKyJiIoWuiIiJFLoiIiZS6IqImCgs43SPz8dwuVzhuLyISNAdz7Om5puFJXTLy8sB6NmzZzguLyISMuXl5SQlJTW6PyzTgAOBALm5uSQkJGg6rIh0CIZhUF5eTnZ2NlZr4z23YQldEZHOSg/SRERMpNAVETGRQldExEQKXREREyl0RURMpNCVsJswYQK33XbbSdsXL15McnIyAPfffz8Wi4XLLrvspHZz587FYrEwYcKEk/YdOXIEh8PB8OHDG7y2xWKp+0pKSmLs2LGsXLmybv+aNWuYNm0a2dnZWCwWli5d2pqPKFJHoSsRo1u3bqxatYojR47U275w4UJ69erV4DGLFy/m6quvxuVysWHDhgbbLFq0iLy8PNauXUt6ejpTp05l3759AFRWVjJy5Eieeuqp4H4Y6bQUuhIxMjMzueSSS/jHP/5Rt23dunUUFRUxZcqUk9obhsGiRYu49tprmTlzJgsWLGjwvMnJyXTt2pXhw4czf/58qqurWb58OQCXX345Dz74INOnTw/Nh5JOR6ErEWXOnDksXry47ueFCxcya9YsHA7HSW1XrVpFVVUVkydPZvbs2bz66qtUVlae8vyxsbEAeDyeoNYtcpxCVyLK1KlTcblcrFmzhsrKSl577TXmzJnTYNsFCxZwzTXXYLPZGD58OP369eP1119v9NxVVVXcc8892Gw2LrzwwlB9BOnk9DZgiShRUVHMnj2bRYsWsW/fPgYOHMiIESNOaldWVsaSJUv49NNP67bNnj2bBQsW8NOf/rRe2xkzZmCz2aiuriYjI4MFCxY0eE6RYFDoStglJibidDpP2l5WVtbgak1z5sxhzJgxbNmypdG73JdffpmamhrGjBlTt80wDAKBALt27WLgwIF12x9//HEmT55MUlISGRkZQfhEIo1T94KE3aBBg8jJyTlpe05OTr1wPG7YsGEMGzaMLVu2MHPmzAbPuWDBAn7729+yadOmuq+vv/6acePGsXDhwnptu3btyoABAxS4Ygrd6UrY3XzzzcybN49bb72VG2+8kejoaJYtW8Yrr7zCu+++2+AxK1euxOv11o3j/bZNmzaRk5PDSy+9xODBg+vtmzFjBg888AAPPvggdnvTv/4VFRXs2bOn7uf9+/ezadMmUlNTGx2mJnIqutOVsOvXrx9r1qxhx44dTJ48mTFjxvDaa6/x+uuvNzgZAiA+Pr7BwIXau9yhQ4eeFLgA06dPp6CggPfff79ZtX355ZeMHj2a0aNHA3D77bczevRo7r333uZ9OJHv0Hq6IiIm0p2uiIiJFLoiIiZS6IqImEihKyJiIoWuiIiJFLoiIiZS6IqImEihKyJiIoWuiIiJFLoiIiZS6IqImOj/A+AMse5ih1iBAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 400x400 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Quick look at the clusters\n",
    "plt.rcParams[\"figure.figsize\"] = (4, 4)\n",
    "sc.pl.umap(adata, color=[\"leiden\"], wspace=0.4, ncols=3, legend_loc = \"on data\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3c031c1b",
   "metadata": {},
   "source": [
    "## Load the Malva index"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "14b03642",
   "metadata": {},
   "outputs": [],
   "source": "# Open the index we built in step 1\nmindex = MalvaIndex(\"indices/pbmc_1k_v3\")"
  },
  {
   "cell_type": "markdown",
   "id": "6bb82780",
   "metadata": {},
   "source": [
    "## Query a custom sequence\n",
    "\n",
    "We query `sequence.fa`, which contains one isoform of the CD3E gene.\n",
    "\n",
    "This demonstrates how you can search for any sequence of interest:\n",
    "transcript isoforms, viral sequences, circular RNAs, etc."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "49191a9d",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Load the query sequence from FASTA file\n",
    "query = []\n",
    "with dnaio.open(\"sequence.fa\") as fasta:\n",
    "    for record in fasta:\n",
    "        query.append([record.sequence])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "77b0d422",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Search the index for cells containing k-mers from the query sequence\n",
    "mindex.open()\n",
    "results = mindex.where(\n",
    "    sequence = query,\n",
    "    sliding_size = 64,        # window size for k-mer matching\n",
    "    pct_threshold = 0.65,     # minimum fraction of k-mers that must match\n",
    "    count_at_most = 100000,   # upper count threshold\n",
    "    count_at_least = 0,       # lower count threshold\n",
    ")\n",
    "mindex.close()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "91b009a5",
   "metadata": {},
   "source": [
    "## Project results onto cells"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "32bcdfb1",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Convert results to a dataframe and merge with cell annotations\n",
    "df_results = pd.DataFrame({\n",
    "    \"cell\": results[0][0],\n",
    "    \"expression\": results[0][1]\n",
    "})\n",
    "df_results.index = df_results[\"cell\"].astype(str)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d227deed",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Add query expression to the AnnData object\n",
    "adata.obs[\"query_expression\"] = np.nan_to_num(\n",
    "    adata.obs.merge(df_results, left_index=True, right_index=True, how=\"left\")[\"expression\"].values,\n",
    "    nan=0\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "027167d6",
   "metadata": {},
   "source": [
    "## Compare sequence query with gene-level quantification\n",
    "\n",
    "The query sequence is a CD3E isoform. We can compare our sequence-level\n",
    "search results with the gene-level CD3E expression from `malva quant`.\n",
    "\n",
    "The sequence query results correlate with CD3E expression,\n",
    "as both capture T cell populations. Sequence-level queries\n",
    "can detect isoform-specific expression patterns not visible at the gene level."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "id": "2dc1e022",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABW8AAAFxCAYAAAAMOVp9AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAA5FdJREFUeJzs3XWcVNX7wPHPnZ7tYAOWpbtEOgQDFBURVAxUUH62GNjxVQkDW7FQDLAwsFEREAnp7u4Ftjum7/39sbIwbLA9u8Pzfr3u98vce885zww4Z+5zzz1H0TRNQwghhBBCCCGEEEIIIUSdovN1AEIIIYQQQgghhBBCCCGKk+StEEIIIYQQQgghhBBC1EGSvBVCCCGEEEIIIYQQQog6SJK3QgghhBBCCCGEEEIIUQdJ8lYIIYQQQgghhBBCCCHqIEneCiGEEEIIIYQQQgghRB0kyVshhBBCCCGEEEIIIYSogyR5K4QQQgghhBBCCCGEEHWQJG+FEEIIIYQQQgghhBCiDpLkraiTZs6ciaIoHDp0qELlJk6ciKIo5TpXURQmTpxY8eCEEEIIISqgsr9rhBBCCCGEkOStEEIIIYQQQgghxFlo//793HXXXbRo0QKLxUJISAj9+/dn6tSp2Gw2AJo1a4aiKCiKgk6nIywsjM6dO3PnnXeyevXqEut96KGH6NatGxEREQQEBNC+fXsmTpxIXl6e13knbnCWtq1atarGPwMh6jqDrwMQojo988wzPPnkk74OQwghhBCiyOjRo7nhhhswm82+DkUIIYQo8scff3DttddiNpsZM2YMnTp1wul0smzZMh577DG2b9/O9OnTAejatSuPPPIIALm5uezcuZPZs2fz8ccf89BDD/Hmm2961b127VoGDBjA2LFjsVgsbNy4kZdffpm///6bpUuXotN5jyWcPHkyzZs3LxZjq1ataujdC1F/SPJW+BWDwYDBIP+shRBCiMqy2+2YTKZiF1X+oqCggICAgFptU6/Xo9fra7VNIYQQoiwHDx7khhtuoGnTpvzzzz80bNiw6Ni4cePYt28ff/zxR9G+uLg4br75Zq86XnnlFW688UbeeustWrduzT333FN0bNmyZcXabNmyJY8++ihr1qyhT58+Xscuu+wyevToUV1vTwi/4p+/yoVfmjt3LgMGDCAwMJDg4GCGDh3K9u3bvc4pac5bh8PBQw89RFRUFMHBwVx55ZUcPXq0xDaOHTvG//3f/xETE4PZbKZjx4589tlnXucsXrwYRVH4/vvvefHFF2ncuDEWi4VBgwaxb9++6n3TQgghzhrLli2jZ8+eWCwWWrZsyUcffeTVrx06dAhFUZg5c2axsiXN416RPu3bb7/lmWeeIS4ujoCAADZt2oSiKLz11lvF2lqxYgWKovDNN9+U+72pqsrbb79Nx44dsVgsxMTEcNddd5GZmVl0zoQJE9DpdCxcuNCr7J133onJZGLz5s1eMX/33Xc8/fTTxMbGEhgYyJVXXklCQoJX2QsuuIBOnTqxfv16Bg4cSEBAAE8//TRQ+PtgwoQJtGrVCrPZTHx8PI8//jgOh8OrjgULFnDeeecRFhZGUFAQbdu2LarjhHfffZeOHTsSEBBAeHg4PXr0YNasWUXHS5vz9oMPPqBjx46YzWYaNWrEuHHjyMrKKvE97NixgwsvvJCAgADi4uJ49dVXy/35CyGEEKd79dVXycvL49NPP/VK3J7QqlUrHnzwwTLrsFqtfPnll0RERPDiiy+iaVqZ5zdr1gygWF8nhCibDFEU9cKXX37JLbfcwpAhQ3jllVcoKChg2rRpnHfeeWzcuLGoEyjJ7bffzldffcWNN95Iv379+Oeffxg6dGix85KTk+nTpw+KonDfffcRFRXF3Llzue2228jJyWH8+PFe57/88svodDoeffRRsrOzefXVV7nppptKnfNHCCGEKM3WrVu55JJLiIqKYuLEibjdbiZMmEBMTEyl6qton/b8889jMpl49NFHcTgctGvXjv79+/P111/z0EMPeZ379ddfExwczPDhw8sdz1133cXMmTMZO3YsDzzwAAcPHuS9995j48aNLF++HKPRyDPPPMOcOXO47bbb2Lp1K8HBwcybN4+PP/6Y559/nnPOOcerzhdffBFFUXjiiSdISUnh7bffZvDgwWzatAmr1Vp0Xnp6Opdddhk33HADN998MzExMaiqypVXXsmyZcu48847ad++PVu3buWtt95iz549/PLLLwBs376dK664gi5dujB58mTMZjP79u1j+fLlRfV//PHHPPDAA4wcOZIHH3wQu93Oli1bWL16NTfeeGOpn8nEiROZNGkSgwcP5p577mH37t1MmzaNtWvXFn0mJ2RmZnLppZdy9dVXc9111/HDDz/wxBNP0LlzZy677LJy/z0IIYQQJ8yZM4cWLVrQr1+/KtUTFBTEVVddxaeffsqOHTvo2LFj0TG3201WVhZOp5Nt27bxzDPPEBwcTK9evYrVk52dTVpamtc+RVGIjIysUnxC+AVNiDpoxowZGqAdPHhQy83N1cLCwrQ77rjD65ykpCQtNDTUa/+ECRO0U/9Zb9q0SQO0e++916vsjTfeqAHahAkTivbddtttWsOGDbW0tDSvc2+44QYtNDRUKygo0DRN0xYtWqQBWvv27TWHw1F03tSpUzVA27p1a5XfvxBCiLPLiBEjNIvFoh0+fLho344dOzS9Xl/Urx08eFADtBkzZhQrX9U+rUWLFkX7Tvjoo480QNu5c2fRPqfTqTVo0EC75ZZbyv3e/v33Xw3Qvv76a6/9f/31V7H9W7du1Uwmk3b77bdrmZmZWlxcnNajRw/N5XIVnXMi5ri4OC0nJ6do//fff68B2tSpU4v2nX/++Rqgffjhh15tf/nll5pOp9P+/fdfr/0ffvihBmjLly/XNE3T3nrrLQ3QUlNTS31/w4cP1zp27FjmZ3Dq7xpN07SUlBTNZDJpl1xyiebxeIrOe++99zRA++yzz4q9hy+++KJon8Ph0GJjY7VrrrmmzHaFEEKIkmRnZ2uANnz48HKd37RpU23o0KGlHj/RX/76669e+1euXKkBRVvbtm21RYsWeZ1zoo8saTObzRV9a0L4JZk2QdR5CxYsICsri1GjRpGWlla06fV6evfuzaJFi0ot++effwLwwAMPeO0/fcSRpmn8+OOPDBs2DE3TvNoZMmQI2dnZbNiwwavM2LFjMZlMRa8HDBgAwIEDB6rydoUQQpxlPB4P8+bNY8SIETRp0qRof/v27RkyZEiF66tMn3bLLbd4jVYFuO6667BYLHz99ddF++bNm0daWlqxOe/KMnv2bEJDQ7n44ou9YunevTtBQUFe/XinTp2YNGkSn3zyCUOGDCEtLY3PP/+8xPnsx4wZQ3BwcNHrkSNH0rBhw6K+/wSz2czYsWOLxdS+fXvatWvnFdNFF10EUBRTWFgYAL/++iuqqpb4/sLCwjh69Chr164t92fy999/43Q6GT9+vNfcwnfccQchISFecwxC4aimUz9zk8lEr1695DeHEEKISsnJyQHw6kerIigoCChcyOxUHTp0YMGCBfzyyy88/vjjBAYGkpeXV2Id77//PgsWLPDa5s6dWy3xCVHfybQJos7bu3cvQNEF1elCQkJKLXv48GF0Oh0tW7b02t+2bVuv16mpqWRlZTF9+vSi1TRPl5KS4vX61AtsgPDwcACv+fuEEEKIM0lNTcVms9G6detix9q2bVssGVme+irap5W0unNYWBjDhg1j1qxZPP/880DhlAlxcXGl9skl2bt3L9nZ2URHR5crlscee4xvv/2WNWvW8NJLL9GhQ4cSy53+eSmKQqtWrYrNKxsXF+d1s/VETDt37iQqKqrMmK6//no++eQTbr/9dp588kkGDRrE1VdfzciRI4uSrk888QR///03vXr1olWrVlxyySXceOON9O/fv+QPhMLfJ1D894jJZKJFixZFx09o3LhxsTn9w8PD2bJlS6ltCCGEEKU5cQ19erK1sk4kZE9PBoeEhDB48GAAhg8fzqxZsxg+fDgbNmwoNh1Sr169ZMEyIUohyVtR550Y6fLll18SGxtb7HhJo3Eq28bNN9/MLbfcUuI5Xbp08Xpd2qrR2hkmaRdCCCEq4/Tk3Qkej8frdWX6tNNH3Z4wZswYZs+ezYoVK+jcuTO//fYb9957r9do0TNRVZXo6GivEbynOj2BeuDAgaIbt1u3bi13O6Up6b2pqkrnzp158803SywTHx9fVHbp0qUsWrSIP/74g7/++ovvvvuOiy66iPnz56PX62nfvj27d+/m999/56+//uLHH3/kgw8+4LnnnmPSpElVjh/kN4cQQojqFRISQqNGjdi2bVu11HeinlatWpV53tVXX83o0aP59ttviyVvhRClk+StqPNOjJqNjo4uumtXXk2bNkVVVfbv3+81umX37t1e50VFRREcHIzH46lwG0IIIURVREVFYbVaixKWpzq1vzrxhMfpKzSfPkqzOvu0Sy+9lKioKL7++mt69+5NQUEBo0ePrlAdLVu25O+//6Z///6lJolPUFWVW2+9lZCQEMaPH89LL73EyJEjufrqq4ude/rnpWka+/btK5aYLi2mzZs3M2jQoFKT4ifodDoGDRrEoEGDePPNN3nppZf43//+x6JFi4o+38DAQK6//nquv/56nE4nV199NS+++CJPPfUUFoulWJ1NmzYFCv9+W7RoUbTf6XRy8OBB+S0ihBCixl1xxRVMnz6dlStX0rdv30rXk5eXx88//0x8fDzt27cv81yHw4GqqmRnZ1e6PSHORjLnrajzhgwZQkhICC+99BIul6vY8dTU1FLLnliB+Z133vHa//bbb3u91uv1XHPNNfz4448l3n0sqw0hhBCiKvR6PUOGDOGXX37hyJEjRft37tzJvHnzil6HhITQoEEDli5d6lX+gw8+KFZfdfVpBoOBUaNG8f333zNz5kw6d+5cruToqa677jo8Hk/R1AunOrEK9QlvvvkmK1asYPr06Tz//PP069ePe+65p9jq0wBffPGF1+OeP/zwA4mJiUV9/5liOnbsGB9//HGxYzabjfz8fAAyMjKKHe/atStQeAEKkJ6e7nXcZDLRoUMHNE0r8XcLwODBgzGZTLzzzjteo2c//fRTsrOzGTp06BnfgxBCCFEVJ+agvf3220lOTi52fP/+/UydOrXMOmw2G6NHjyYjI4P//e9/RTdEs7KySuwDP/nkEwCZHkGICpKRt6LOCwkJYdq0aYwePZpu3bpxww03EBUVxZEjR/jjjz/o378/7733Xollu3btyqhRo/jggw/Izs6mX79+LFy4kH379hU79+WXX2bRokX07t2bO+64gw4dOpCRkcGGDRv4+++/S7yAE0IIIarDpEmT+OuvvxgwYAD33nsvbrebd999l44dO3rNa3r77bfz8ssvc/vtt9OjRw+WLl3Knj17itVXnX3amDFjeOedd1i0aBGvvPJKhd/b+eefz1133cWUKVPYtGkTl1xyCUajkb179zJ79mymTp3KyJEj2blzJ88++yy33norw4YNA2DmzJl07dqVe++9l++//96r3oiICM477zzGjh1LcnIyb7/9Nq1ateKOO+44Y0yjR4/m+++/5+6772bRokX0798fj8fDrl27+P7775k3bx49evRg8uTJLF26lKFDh9K0aVNSUlL44IMPaNy4Meeddx4Al1xyCbGxsfTv35+YmBh27tzJe++9x9ChQ0tdCCYqKoqnnnqKSZMmcemll3LllVeye/duPvjgA3r27FmhBeGEEEKIymjZsiWzZs3i+uuvp3379owZM4ZOnTrhdDpZsWIFs2fP5tZbby06/9ixY3z11VdA4WjbHTt2MHv2bJKSknjkkUe46667is5dvHgxDzzwACNHjqR169Y4nU7+/fdffvrpJ3r06FFiPzd37lx27dpVbH+/fv28nlIR4qykCVEHzZgxQwO0gwcPFu1btGiRNmTIEC00NFSzWCxay5YttVtvvVVbt25d0TkTJkzQTv9nbbPZtAceeECLjIzUAgMDtWHDhmkJCQkaoE2YMMHr3OTkZG3cuHFafHy8ZjQatdjYWG3QoEHa9OnTveIAtNmzZ3uVPXjwoAZoM2bMqLbPQQghxNljyZIlWvfu3TWTyaS1aNFC+/DDD4v1awUFBdptt92mhYaGasHBwdp1112npaSkVHufdrqOHTtqOp1OO3r0aKXf3/Tp07Xu3btrVqtVCw4O1jp37qw9/vjj2vHjxzW326317NlTa9y4sZaVleVVburUqRqgfffdd14xf/PNN9pTTz2lRUdHa1arVRs6dKh2+PBhr7Lnn3++1rFjxxLjcTqd2iuvvKJ17NhRM5vNWnh4uNa9e3dt0qRJWnZ2tqZpmrZw4UJt+PDhWqNGjTSTyaQ1atRIGzVqlLZnz56iej766CNt4MCBWmRkpGY2m7WWLVtqjz32WFEdmlby7xpN07T33ntPa9eunWY0GrWYmBjtnnvu0TIzM8v1Hm655RatadOmZX7mQgghxJns2bNHu+OOO7RmzZppJpNJCw4O1vr376+9++67mt1u1zRN05o2baoBGqApiqKFhIRoHTt21O644w5t9erVxerct2+fNmbMGK1Fixaa1WrVLBaL1rFjR23ChAlaXl6e17kn+sjSNrm+FkLTFE2TlQ6EEEIIIeqiiRMnMmnSJJ8vTHXuuecSERHBwoULfRoHFI7mufDCC5k9ezYjR470dThCCCGEEELUKJnzVgghhBBClGrdunVs2rSJMWPG+DoUIYQQQgghzjoy560QQgghhChm27ZtrF+/njfeeIOGDRty/fXXex33eDxnXPwsKCiIoKCgmgxTCCGEEEIIvyYjb4UQQgghRDE//PADY8eOxeVy8c0332CxWLyOJyQk0LBhwzK3119/3UfRCyGEEEII4R9kzlshhBBCCFFhdrudZcuWlXlOixYtZIVoIcrh2LFjPPHEE8ydO5eCggJatWrFjBkz6NGjBwCapjFhwgQ+/vhjsrKy6N+/P9OmTaN169Y+jlwIIYQQNU2St0IIIYQQQgjhI5mZmZx77rlceOGF3HPPPURFRbF3715atmxJy5YtAXjllVeYMmUKn3/+Oc2bN+fZZ59l69at7Nixo9ioeCGEEEL4F58kb1VV5fjx4wQHB6MoSm03L4QQPqNpGrm5uTRq1AidTmauEXWf9NlCiLNVSX223W7H6XSWq+zp35lmsxmz2Vzs3CeffJLly5fz77//llpXo0aNeOSRR3j00UcByM7OJiYmhpkzZ3LDDTdU9K0JPyV9thDibFWVPttkMtX5G6E+Sd4ePXqU+Pj42m5WCCHqjISEBBo3buzrMIQ4I+mzhRBnuxN9tt1up3nTIJJSPGcsExQURF5ente+CRMmMHHixGLndujQgSFDhnD06FGWLFlCXFwc9957L3fccQcABw4coGXLlmzcuJGuXbsWlTv//PPp2rUrU6dOrdL7E/5D+mwhxNmuMn12bGwsBw8erNMJXIMvGg0ODgYKP9SQkBBfhCCEED6Rk5NDfHx80fegEHWd9NlCiLPV6X220+kkKcXD4fXNCAku/emZnFyVpt0PFfveLGnULRQmZ6dNm8bDDz/M008/zdq1a3nggQcwmUzccsstJCUlARATE+NVLiYmpuiYECB9thDi7FXVPtvpdEry9nQnHuEICQmRTkUIcVaSR9lEfSF9thDibHd6nx0UrBAUXHo/rlKx701VVenRowcvvfQSAOeeey7btm3jww8/5JZbbqlC5OJsI322EOJsV9k+u66TCReFEEIIIYQQopw8mnrGrSIaNmxIhw4dvPa1b9+eI0eOAIWPcwIkJyd7nZOcnFx0TAghhBDFVXef7SuSvBVCCCGEEEKIclLRzrhVRP/+/dm9e7fXvj179tC0aVMAmjdvTmxsLAsXLiw6npOTw+rVq+nbt2/V35AQQgjhp6q7z/YVn0ybIIQQQgghhBD1kYpKWeN0yj5a3EMPPUS/fv146aWXuO6661izZg3Tp09n+vTpQOEjoOPHj+eFF16gdevWNG/enGeffZZGjRoxYsSIyr8RIYQQws9Vd5/tK5K8FUIIIYQQQohy8mgaHq30kTplHStJz549+fnnn3nqqaeYPHkyzZs35+233+amm24qOufxxx8nPz+fO++8k6ysLM477zz++uuvOr24ihBCCOFr1d1n+4okb4UQQgghhBDCh6644gquuOKKUo8risLkyZOZPHlyLUYlhBBCiLpAkrdCCCGEEEIIUU5nmiOvvsyfJ4QQQvg7f+mzJXkrhBBCCCGEEOWkouHxgwtBIYQQwt/5S5+t83UAQpSLqwAKMkCrH5NJCyGEEGer/JwCju9Pwu1y+zoUIWqEv6xcLYQQtnw7x/cn4bQ7fR2KEDXCX/psSd6KOmfp0qVcfvnlREVFoSgKiqLw4aNXw4rXYPEk2Dcfl8vFpEmTaNGiBSaTicaNG/PQQw+Rl5fn6/CFEEKIs5K9wMEbt0/j2pjbGdN+PMNb3M/Ux770dVhCVLsTi5+UtQkhRF3mcXv48JHPuS72dsa0f5BhbR5g8t3T0eT7S/gZf+mzZdoEUeds2LCBBQsW0KJFC9LS0rwPepxwaBH/d/fDfLVwOzqdjtatW3PgwAHefvttNm7cyD///INOJ/clhBBCiNr01p0f8s+sZWhNYyEuBodex5w1Cczp9j+6xgTy0s+PYrKYfB2mEFWm/reVdVwIIeqyj5/4ih/f+h13ixhc7ePBqGd+Sg4LLn6eDh4Pb/32JAHBVl+HKUSV+UufLRkuUeeMHj2anJwc5v1Y8midDXuT+GrhdgCm3juYXVs38uOPPwKwZMkSfvnll9oKVQghhBBA2vEMFn+3Ai06Apo0BP0pPzGDA9mUnM8VQTdz/ECS74IUQgghBLZ8O39MX4AnKhRXl2Zg1Bcd00IC2GEycmXYGLYt3+m7IIUQXiR5K+qcyMhIrFYr2DNLPD53zYGiP19zXhtYO42hQ4disVgA+Ouvv2olTiGEEEIUSj6UiupRoWFkyScEB6IGWHigz/9qNzAhaoDnv8VPytqEEKKuykrOxp7vwN08usTjWqAFNTqMJy95oZYjE6L6+UufLclbUXcFNChxd0JqTtGfo8MCoCAVnU5HgwaF5x85cqRWwhNCCCFEocZtGoICmIyln2Q0kp2Wg73AUWtxCVETPNqZNyGEqKsi4yLQG3RgLr3P1ixGHDYnR/ccq8XIhKh+/tJnS/JW1F0BpYzeKYVMri6EEEL4RmiDEMJjQiEnv+QTVBXyCgAoyLXVYmRCVD+1HJsQQtRVJrORxm0aocvILfUcXXrhsdSjGbUVlhA1wl/6bEneinonPiqk6M8pWQVgsKCqKunp6QA0adLEV6EJIYQQZ62bnhkJCcng9hQ/eDwVxeXGYDQQERNW67EJUZ1UFDxlbCqKr0MUQogyjX1hFIb9yWB3FjumT0hDl2cHoPOA9rUdmhDVyl/6bEneinrn0p4tiv7847+7od1w/vjjD+z2wg7m0ksv9VVoQgghxFlr6B2D6XNhB9i8B1IyCi8Ic/NhzxGUg8cBuGXydT6OUoiqU7Uzb0IIUZf1G96TS67vi3npDvSHUlAKHCjZ+Ri3Hsa4fj8AV9x9CQajwceRClE1/tJny3+Jos756aefePzxx3G73UX7nvt6La//soXe7Rrx9eOXMurCDnyzaAcPTlvI+wtHsX9/YQczYMAARowY4aPIhRBCiLOXwWjg+d+eZN28Tayfv5mdq/awd+tB3E4PodEh3PnaGC4efb6vwxRCCCHOeoqi8PiM+7h07EWs/G0dezceYOeqvbgcLgLDAxk94VqufmCor8MUQvxHkreizsnJySlKxp6Qmp5Bajo0bt4GBr3A5wMn0PqFF/jiiy/Yv38/UVFRjBw5khdeeAGdTgaUCyGEEL6gKAo9Lz2Xnpee6+tQhKgxJx61LOu4EELUB10GdqDLwA6+DkOIGuMvfbYkb0Wdc+utt3LrrbeWeY7RaGTSpElMmjSpdoISQgghhBAC/7kQFEIIIfydv/TZkrwVQgghhBBCiHJSNQVVK/1ir6xjQgghhKg9/tJnS/JWCCGEEEIIIcrJX0bxCCGEEP7OX/psmRxUCCGEEEIIIYQQQggh6iAZeSuEEEIIIYQQ5eRBh6eMMTCeWoxFCCGEEKXzlz5bkrdCCCGEEEIIUU7aGebP0+rJ/HlCCCGEv/OXPluSt0IIIYQQQghRTv4yf54QQgjh7/ylz5bkrRBCCCGEEEKUk0fT4dHKeARTq8VghBBCCFEqf+mzJXkrhBBCCCGEEOWkoqCWMX+eSj25EhRCCCH8nL/02aW/AyGEEEIIIYQQQgghhBA+IyNvhRBCCCGEEKKc/GX+PCGEEMLf+UufLclbIYQQQgghhCinM8+fVz8ewRRCCCH8nb/02TJtghCiXlLdGh5H/fiiFUIIIc5mHrcHe4HD12FUm8L588rehBBCiPpIVVVs+XZfh1Ft/KXPlpG3Qoh6wZ6h4coHY5DGlo8g4R8N1QWRnSA4HrL3gyEAml6i0OxSBUWpH1/CQgghhL/JSc8lOy2H8JhQZjzzLQu+WIItz07r7i1oeU4zDm07gt6o54Lr+3PFXRdjMNavSxIVHR4/WPxECCGEyM/OJyMpi4iG4Xw75Wf+mL6A3Mx8mnWKp12vVhzZdRzVo9J/RC9G3H8ZlgCzr0OuEH/ps+vXLyUhxFnHnqGx7jWVxNVQ0vdq+rbC7YTUTRrp26HHo5K8FUIIIWpTfnY+b98znSXfr0RTi3fae9cfYO/6A0Wvty/fzfoFm3n+1ydrM8wq85dHMIUQQpy9nHYnH4yfyV8z/sHj8hQ7fmhbAoe2JRS93rV6LyvnrOP1fyZgNBlrM9Qq8Zc+W6ZNEELUaUseUUlcRYmJ29Ic/EMjY5daYzEJIYQQorjJ177B4m9XlJi4Lc2qOetZ8evaGoxKCCGEEKd7++7p/DF9QYmJ29LsWLGbuZ8srMGo6r5p06bRpUsXQkJCCAkJoW/fvsydO7fMMrNnz6Zdu3ZYLBY6d+7Mn3/+WeF2ZeStEKJOUV0aR5doJG+ArAMaOYcqV8+ShzUuek8jtIWMwBVCCCFqgsfjYdWc9az+YwOJB5PZ9M+2MxcqwQs3vMkr85+j84D21RxhzVDRofrBI5hCCCHOLuvmb2bZT6tJO5bO6j82VKqO9x+cQYO4SPoN71nN0dWM6u6zGzduzMsvv0zr1q3RNI3PP/+c4cOHs3HjRjp27Fjs/BUrVjBq1CimTJnCFVdcwaxZsxgxYgQbNmygU6dO5W5XkrdCiDrDbddYPF4lc3c11GWDNVNULv5YX/XKhBBCCOHF4/bwzLCXWTdvU5XrcjncvHTj23x9eBo6Xd1/MNCjKXi00m8Ol3VMCCGE8IWXx7zDwq/+rXI9qkfl5dHv8F3ix1gDLdUQWc0qb5+dk5Pjtd9sNmM2F5/fd9iwYV6vX3zxRaZNm8aqVatKTN5OnTqVSy+9lMceewyA559/ngULFvDee+/x4Ycflvt91P1fR0KIs0LOYY0Fd1RP4vaErH2w9nUPWj2Zx0YIIYSoD5IOpfBg//9VS+L2hLRjGTw3/BU87vI/vukrnv8WPylrE0IIIeqCjKRMnrrsxWpJ3J5gy7Pz+ODJ2PLt1VZnTSlvnx0fH09oaGjRNmXKlDPX7fHw7bffkp+fT9++fUs8Z+XKlQwePNhr35AhQ1i5cmWF3oeMvBVC+FzaVo3FD6to7uqv+9AfENFGo+WVMgpGCCGEqKqD244w/rxnKcgpqPa6V/+xga+e/4FbJl1f7XVXJ1XToZax+IkqN42FEELUASkJadzX+ykyk7Kqve5dq/fy0cOfM/6ju6q97upU3j47ISGBkJCQov0ljbo9YevWrfTt2xe73U5QUBA///wzHTp0KPHcpKQkYmJivPbFxMSQlJRUkbcht4WFEL63cnLNJG5PODBHLqKEEEKI6jDl5qk1krg94Y/pC2qs7uoiI2+FEELUB2/cNq1GErcn/P3VUhw2R43VXx3K22efWIDsxFZW8rZt27Zs2rSJ1atXc88993DLLbewY8eOGn0fMvJWCOFTKZtU7Gk124Y9o2brF0IIIc4GyYdTObjlSI22kZWSg6ZpKIo8MSOEEEJUVn52PhsXbqnRNhw2J/nZBZitpSc6/ZHJZKJVq1YAdO/enbVr1zJ16lQ++uijYufGxsaSnJzstS85OZnY2NgKtSm3hYUQPnV8Wc2Pio1oBykbNHKOyAhcIYQQorKW/rCqxtto26sVmxdvZ//mQzXeVmWpnFwApaRN9XWAQgghznqrft9ATc/i06hlLAm7jrN73f6abagKaqPPVlUVh6PkEch9+/Zl4cKFXvsWLFhQ6hy5pZGRt0II31IUoAZ7FT0kroHjKwq/lqO6Qp9ndVgiZESPEEIIURGqp2bTkopO4eDWwzw2aBIArbu34JlvH6JRy4qNTqlpKjrUMsbAlHVMCCGEqA013WcDZKZk8ehFEwGIbxfH07MepFXX5jXebkVUd5/91FNPcdlll9GkSRNyc3OZNWsWixcvZt68eQCMGTOGuLi4ogXPHnzwQc4//3zeeOMNhg4dyrfffsu6deuYPn16hdqVXxZCCJ9q1L/ySdQ3t9/AdYsNXLfYwNvbbyz5JA9e8+mmboJVz8uYGCGEEKKi+l7Zo0br11QNR4Gz6PXe9Qd49sqX0erYAmAeTXfGTQghhPClnpd1xWDU11wDCthy7UUvE3Yd43+Xv4TT7iyjUO2r7j47JSWFMWPG0LZtWwYNGsTatWuZN28eF198MQBHjhwhMTGx6Px+/foxa9Yspk+fzjnnnMMPP/zAL7/8QqdOnSrUroy8FUL4VHRXhZiekLy2YuUWJc5kVeoPlWozdRMc+UcltoeCKURG4AohhBDl0aRdHJffPog/P1l45pOryZGdx5jz4XwGXN2b8JiwWmu3LCoKKqX/fijrmBBCCFEbwqJCGT3hWmY8823NNFDCfdWMpCy+f+03hoy9kKjGkTXTbgVVd5/96aeflnl88eLFxfZde+21XHvttRVq53RyW1gI4XOtRlTsqyjJtp8Z+8bTJqQPkebGlWpz9fMac65V2fS+iuapWyN6hBBCiLrq8jsvrvU23x33CaPi7+b1//ugTozokZG3Qggh6oOLx1xQ621+PuE7bm52D5NGvk5+dn6tt386f+mz60eUQgi/Fn0uGALLd65HdfPujtEo6Hig/ZfoqPyjIKoT9v6gsepFFVUSuEIIIcQZtTq3GVHxtT+axuP2MG/mIp4Z9nKdSOBWp4kTJ6IoitfWrl27ouN2u51x48YRGRlJUFAQ11xzTbGVq4UQQojTRTWOpHX3FrXerqpqLPtpNY8Omowtz1br7fsjSd4KIXzOYFU4934FpRzfSLMPT2Zv7hpub/Me0dbqmQz96CJYcIeKI0sSuEIIIURZ9Ho9D067E6PZ6JP2Ny7cyth2D5J40HfJSw+6M24V1bFjRxITE4u2ZcuWFR176KGHmDNnDrNnz2bJkiUcP36cq6++ujrfkhBCCD9137u3ERBi9Unb+zYc4Obm49i/+ZBP2oea6bN9oX5EKYTwe82G6Lj4Ex2xvUs/Z3/OOn45/AoDYm5iQEwpC5RVUs5B2PKhLGQmhBBCnEnvy7vxybY3GTx6oE/aTzmSxtS7K7ZKc3VSNeWMW0UZDAZiY2OLtgYNGgCQnZ3Np59+yptvvslFF11E9+7dmTFjBitWrGDVqlXV/daEEEL4mQ592vDp9re4ZvxQn7Sfk57Ly6Pf8UnbUDN9ti9I8lYIUWeENlewNij9yzMhfzsqHlal/sjopaGMXhpKmuMIAKvTfmL00lAK3NmVbv/wf+uvJO3fy6Z5f7Bv7SpUj6fS9QkhhBD+qlHLWGKaRvms/fULtpCXlc/hHQnMmTaPJbNX1tp0CuoZRvCo/11i5eTkeG0Oh6PUOvfu3UujRo1o0aIFN910E0eOFP6+Wb9+PS6Xi8GDBxed265dO5o0acLKlStr9o0KIYTwCw3iIolr3dBn7R/alkDC7mMkHkhmzofzWfj1v7U2nUJ5++y6zuDrAIQQ4lT29DNPXeBS7cX2eTQ3Hs2NplV+6gPNrfHzq89zYP3qon0hUTFc/dREIuPiK12vEEII4Y/Sj2f6tP137v2YRd8uL3odHhPK5F+foF2v1jXarqrpUMtY4OTEsfh4798OEyZMYOLEicXO7927NzNnzqRt27YkJiYyadIkBgwYwLZt20hKSsJkMhEWFuZVJiYmhqSkpCq/FyGEEGcHX/fZXz3/A4u/XY6qFl6vB4YG8Mx3D9PjknNqtN3y9tl1nSRvhRB1SkQHhcRVJSdgL2h4Cxc0vMVr37iVLUl1HKZf1HWM7zirSm279UkcOSVxC5CTmsyf77zO6FemVqluIYQQwt+079OGvz77x2ftn5q4BchMzmbytW/w5f730Rsqv6DpmXhQ8FD6k0InjiUkJBASElK032w2l3j+ZZddVvTnLl260Lt3b5o2bcr333+P1eqbeQqFEEL4l3a9a/bG5pn8M2uZ1+v87AKev+4Nvkn4iIDgmuvryttn13X1I8UshDhrtBymYAqt/XY95JMSOqnEYymH9pN65FDtBiSEEELUcRfdeB4NW0SX+3xNUdBMBmpyedDUhHQ2LdpWgy2UX0hIiNdWWvL2dGFhYbRp04Z9+/YRGxuL0+kkKyvL65zk5GRiY2NrIGohhBD+qNfl59K2Z8tyn68poFoN1OSUsAU5Npb/sqbmGvAjkrwVQtQp5jCFQe8r6EzlO//9vvv5/gJ3lUbdqthIbHAXbkNiqefkZ2bgtBVUug0hhBDC31gCzLyz8iWCwgPLPE9TFJyd4im4sjsFw3tgu6wrrhqcLzcrNYf87Pwaq//EI5hlbVWRl5fH/v37adiwId27d8doNLJw4cKi47t37+bIkSP07du3qm9FCCHEWUKn0/HawglEN2lQ5nkakNcrjpS7epB8by9S7uhOXreamy83Jz2X3My8Gqu/pvvs2lI/ohRCnFWC4nSc96IOvaV22ssMmoGqyy31uKIo/DhlAu+OvZ4fXnyW9GMJtROYEEIIUceFRYXy8rxnCY4IKvUcZ9emuNrHgalwxjYtyIKzV0vcjSOrPyAFXr3lPUaE38pDA59l99p91d6Eh5OPYZa8Vcyjjz7KkiVLOHToECtWrOCqq65Cr9czatQoQkNDue2223j44YdZtGgR69evZ+zYsfTt25c+ffpU+3sTQgjhv6xBVl6Z/yxRZfS/eb0bkzugKWpg4WgqNdhM7oXNyT+3ZhK4nzzxNVdHjuXeHo+zYeHWaq+/uvtsX5HkrRCiTorpoTD0Gx3n3A3msJppQ8NNvnkpDtOWss/TNPhvO7xlI98++xj2/Jq7OyiEEELUJ217tOTL/e/xwPu3E9PMe0StZjLgbl7y1AqudjVwIaiB6lEB2LZsFw9fMIHEg8nV2kR1j+I5evQoo0aNom3btlx33XVERkayatUqoqIKP8u33nqLK664gmuuuYaBAwcSGxvLTz/9VK3vSQghxNmhcZtGfLZrKo9+di/NOnovrKnpFPK7l9w35/VsVCPTHrldbgD2bjjI05e9yJ51+6u1fn8ZeatoVVmavZJycnIIDQ0lOzvbaxJ/IYQoierWOLpYI2GxhqIHaxTs+7GCdSj5aHjQa4XfOSoOUsMm4jDtqFRMOoOBbpcNp9eIkViDgstdTr7/RH0j/2aFEBWhqiqr5qzn76+W4HZ5aHpBBz5duafkk90eAn9eW+Mx6Q16Lr99EDc+cw0NGkWUu9zp338nXj+18lIsQcZSy9nzXEzp+5d8b4paJ322EKKi1i/YzPyZi8jPtdHhsq68dGBXqefGvLsanbNmx6oqOoVBN57H6AnX0ahl+ed29/c+2+DrAIQQ4kx0BoUmgxXiLnCzdeE8dixeSnqISlTu4+i1si/C3LokMoI/wm7aCIqG2dme8Lw7sBu3VjpxC6C63ayb8yOHNq3jxhffwGiupTkehBBCiDpMp9PRb3hP+l7ZgwVfLGHu10tRwgLQlOIrnuhy7bUSk8ftYc6H81n95wbeX/syYVFVWxlVQ0EtY3VqrZ6sXC2EEEJ0v/gcul98Dv/+uIo/PvsHQ2szbmPx0ai6XAdKDSduATRV4++v/mXN3E28t3oKDVvEVK0+P+mz68f4YCHEWcnlsHN05zZSjxxC0zR+e/MlFn42jcQD23FadpIY+SAa7hLLWqMLR9cmh/8Pu3kDKIUPGThMO0mNfBZbwLJqiTEt4TDbFi8884lCCCGEH3M5XWxbvos96wsfd3zrzo94bez7bPt7K/oDKSWWMe4+DlDmhZmiU9Ab9dUSY8qRNH5+589qqUsIIYSorzweDztW7WHn6r2oqsrM575l8rVvsH7uRixrj5ZYJnD9cRSgUasyRsMqYDBVzxjRnPRcvn6hgo/b+jEZeSuEqJM2zf+Tf2fNxGkrAAoXDTt9lhdVl0NixHgC7RcRZLsIvRYGQKN+0ONJhdnPvo0nMa1Y3R4tn4DoLBzp1RPr5gV/cO6QodVTmRBCCFHPLP5uOe/e9yk56f8t/qnAqRPjmTYeArcHd4toMBpQChwYdxzDkJBOp/PaMemXx5l0zetsWVL8iRhN1YhuFkniobSi+eerYt5nixj7/Kgq1eHRdHjKmCOvrGNCCCGEL61fsJnX/u990o9lFu079eGYoJVHQYP8bg3RrEZ0+U4C1x0naH0izTs34ZUFz/H+A5+y5PuVxSvXoEFcBEkHS75pW1FLf1jJo5/dW6U6/KXPluStEKLOObxlEws//cBrX2nTc7sNCWQHfU5+2GzO6TiBdoM7ENMT9qxaTnLiilLbyE1PBWDetj0s2LG3xHNeGXkZet2Zv8wzjx/H43ahN5Q+l44QQgjhj/ZvPsRLN01FU0/pp0/rshVNw7zlCKZtCWhGA3q3hz6Xd+Pi18bQ78oebF+xu8TELYBiNJKaUoA+IKDwt4DHg2qv/HQLGUlZ5KTnEhJZ/vnqT6dqCqpW+mOWZR0TQgghfCUlIY3nhr+C0+7y2n/qpbYCBK86StCaY6gWA3q7m24XdeaSL0cyYGQfEg8k8++Pq0ttI/lQ9SRuAWx5do7uTaRx68ovcOovfbYkb4UQdc7mBRV/pNHtLuBg6rs0KBjJ30//TurhgxUqH2g2ERkY4LVPKef8N6rHjaOggICQqs2hJ4QQQtQ3f0z/2ztxWwZF1VAcLjRg5+q9DLp5AA9fMIG960teWVoxGNCZzUW5YEVRwGBAZ7FUOoGraRrpxzOqlLz1oMNTxuxzZR0TQgghfGX+zMXFErelUVQNfUHhuduX72LYPZfwzBVT2LFiN6pHLbVcFR+QKSb5UEqVkrf+0mdL8lYIUefkZVRuPoOM40eZN+3tSpVt3zCaG3qdU6myJqsVa1DlLwKFEEKI+ir9eEalymWlZPPC9W+VeY5iLPmJFsVgKHzGsxJXiAajntgqLn7iL6N4hBBCnF3SjlWuz3bYnEwe+UapT8PWFEVRaNapSZXq8Jc+u36kmIUQZ5XY1m1qvc2tRxN58se5TPrtbz79dy3HMrPLXbbTBRejlGN6BSGEEMLftOvVuuYqV8q4oKpkv9v7iu5YAy2VDKiQiu6MmxBCCFHXtOtd+T67thO3UBhvZMPwKtXhL312/YhSCHFW6X75CAwmU621p1MUgi1mIgKs5Nod7ExM4d2FK8qVwI1q0owLbrmjFqIUQggh6p6hdw4mKDywZipXS34sU9O0Uo+VJaZpFE/PGl/FoIQQQoj66aJR/YltHu3rMMoltEEIz8950tdh1BmSvBVC1Dmh0THc9OJbtZLA7da0EROvHMyTl1/I45ddwB0DewHgVlWW7ztcZtmQBtHcNOWtwjn4hBBCiLNQSGQwH216naCw6k/gqk5niSN9NJerwlMmBIUF8vHWNzCZq764qEdTzrgJIYQQdY3JYuK91VOIahzp61DKpNPr+Gzn24RWYX76E/ylz5bkrRCiTmrQpCnjPv2WyPhmNdpOVHAQAeaTSeK2sVEEmAov7LIKbGWWDW4Qhd5Q9YtAIYQQoj6Ljm/AT+kzOOfCjtVbsaqi2mxobjeapqGpKqrDgeZ0Vriq0KgQrEHW6gnrv/nzytqEEEKIuii0QQizjnzIhaP6+zqUUlkCzFVaWPRU/tJnS/JWCFFnGUwmbn39PfpeeyM6vb5G2vhn534y808mafckpVLgLFxVMzyw7Iu8dv0G1khMQgghRH2jKAqvL5zIAx/cgcFUuTWRFV0JF1Cqimq3o+bnoxYUFI66rYQLru9XqXIl0TQdahmbpskllhBCiLrt6a/H8+z3D2MOMFeqvMlac0/Jnn+d9Nmnq9wvq3rM6cqjwJZKgDUKkzEIVXXjdtswGgNRlPrxlybE2abfyBvpfvkIDm5ex/61q9m9chma6jljuQvH3sWelcs4tmt7qees3H+YuVt3ERpgxWTQk5qTB4DJoGdA6+allmvZozddBl9a8TcjhCi3/JwCju5JJDo+kvCYMDweD3mZ+QSFB6KvoRs6QoiqGXb3JVxyy/msm7eZVXPWsei75TgKzjxSdtRTV3NsXyJLZ6+s9pg6ndeO658YUW31eVDwUPpInbKOCeGv7AUOjuw8SkRsGA3iIlFVlbzMfAJCrBiMZ13aQYh6YeDIvvQb3pP18zez6o8NLP52GXlZBWcsN+TWC7EGWfjlvbnVHlPzzk24bcqN1Vafv/TZfv8tmpS6gf2Hfic3/zigoqpuABRFh05nxONxASpmUxjNmwyhRZMhPo1XCFEyc0AA7foOpF3fgVxy9wOsnP0N2xYvwFGQj+bxoGkaOoMB1e0mNCaW3iOuo/NFlxDXtgM/PP8/7Pl5JdY7qH0rNh9NJDk7l4w8B2GBVppHhjO4Q2uiQ4K8ztXp9HS99Apa9ehNfMcutfG2hTirbPxnK189/wO71+5HU1XcLg+qR0Wn12EOMOG0ufC4PYQ0CGbY3Zdw6+QbfB2yEKIEZquZ/iN60X9EL8Z/dBffvfoLc6bNJzs1B7fbg6ZqGEwG3E43kY3CueahYVz7yDASDySzffku0o9nVjkGRVG47PZB9Bvek56XdkWnq75BGqpGmY9ZqrW/ILcQtW7Xmr18MWk22/7dicfjwePy4HGroIA1yIrL4cTt9BAQEsCQWy/grjfGyI1XIeogg9FA76Hd6T20Ow+8fzu/ffAXP7z5O+nHM/C4VVSPitFswOVwExIZzJX3DuHm50aSk57HuvmbOLonsVriuOjG8zjvqt70G94TvaH6viv8pc9WtJJWAahhOTk5hIaGkp2dTUhISI21czRxOVt2flriseSkXP7vph9KLfvoY3fx7LNPExLcpKbCE0JUE9XjQfV40BsMOO12TFar1yJieRnpbJr/Jxv+/BWXw16pNowWK5ff/yitevSuUqy19f0nRHWprX+zq//cwLPDXi5xcaLSRDQM4/GZ9+F2eQhtEEy7Xq1rLD4hRPXweDy4nW5MFhMFuTasQRavxGpOei5/fvI33732K3kZ+ZVqw2gx8sD7d3Dp2AurFOvp338nXo9dfB2moNIfF3XmOZlxwffS14taV1t99s7Ve3lo4LN4XGd+Eu6EwNAAHv/8PgxGAwHBFjr0a1utN1WEENVPVVVcDhdmq5n8nAIsgWavmzAFuTbmf76Yb17+iYzjWZVqQ2/QM/aFUVz/+PAqxervfbbfjrzVNJWde78v9bjRqKdt+yivffl5To4mZANgd21h2dqJhIe2oluncZjNoTUarxCi8nR6fdGcuOaAgGLHgyIiOe+G0XS/YgRrf/2Btb/9WOE27nj3U6x1+MtciPruo0e/qFDiFiAjMYsnh7xQ9LpFl6Y898MjxLVqWN3hCSGqiV6vR28t7LMDQ4r32SGRwdzwxFUMv+8yfnh9Dl9M/h4qONTk/TUv07xTzQ3AODFPXlnHhfBnnz09q0KJW4D87AImjHi16HWjljE8+/0jtDq39GnKhBC+pdPpMFsL58Qtqc8OCLYy4r7LuOKui/l56p988tTXqB61Qm1M+uUxel/evVriLYm/9Nn1I8pKsDuycLlzSz0eERnAm+9d4bV17VZ4sRcUbOKCwS0ByMzex6Yd02slZiFEzbIGBTPwprFEN2tZ4bKVHbErhDgzj9tDwq5jVa7nwJbDTLzqtQongYUQdY810MLoCdfS69JzK1zWllezfbaKcsZNCH+2c/WeKtdxfH8y/7tiCk5H5RYhFELUHQajgWsfvZIht15Q4bJuZ8VuBFWUv/TZfpu8NRoDoQJ/CTnZdhbM2wfA5cPaYbUai46lZ+4kvyC5ukMUQvjIeaPGVOh8S1AwAWHhNRSNEEJv0GM0G898Yjkc2p7A9uW7qqUuIYTv3TL5egym8j8saDDqiWsVW4MRgUdTzrgJ4c9OjMSrqozETFbNWVctdQkhfG/UU1djDbJUqEzTDo1rKJpC/tJn+23y1qA3ExJU/sel/vhtFw67G6NRz7Cr2hc77nBmV2d4Qggfat61O1c+/DSRjcv3HdHt8isxGKsnsSSEKNm5F3WqtroykrKqrS4hhG+16d6SF39/ijY9yvfUzGW3Dya0Qc1Oc3TiEcyyNiH82YBrqrYGxKmkzxbCfzRsEcOrfz9H54HFc2olOe/q3jRu06hGY/KXPttv57wF6N55HEtWP4OqOss8z+X08PuvhaN0LhzcgoiIEubysEQV2yeEqL9a9+5H6979cNoKSD92lAMb1pCecJjEfXuw5WTjcbsJCAun++XD6XnlNb4OVwi/99SsBxndclylFyg6VVSTBtUQkRCirug2uAvdBnehINdG4sEUVvyyhiM7j7Jz1R4ykrJwOdwEhQcy7O5LuGXS9TUej4pS9srV9eQRTCEq6/73b2ftX5tIOZJW5bqCwwOrISIhRF3Rrldr3lw8GVu+nfRjGfz742oO7Uhg16o9pB3PxGlzYg2yMOTWC7nt5ZtqPB5/6bP9OnlrtTbg/D4vsWjFY5S10sHCBfvIyrShKHDVtSWP/Nl/+A86tr25hiIVQviKyRpAw1ZtaNiqTdE+TdNw2gowWawosgquELUiKDSQr/a/z/Vxd+IoKPum65l889LPTP7l8WqKTAhRVwQEW2nZpSktuzT12l/SCthCiJqj1+uZsWsqN7e4l8ykqj2h+vO7cxl008BqikwIUVdYAy00btOIUU9d5bU/P6cAs9WEwejX6chq5/dZCaslgviGA0o9rmkaP8/eDkDP3vE0aRpW4nlpmTtrIjwhRB2kKArmgEBJ3ApRywJDAxn50LAq17Nz1e5qiEYIUV8EhgTUauJWO8PCJ1o9GcUjRFWYLCZu+t/IKtdzcOuRaohGCFFfBIYE1Gri1l/67LMiM9Gu9XVEhrUr8djqlQkcTSi8W3jN9aXPt6dpHtIzd5Kde7hGYhRCCCEE3PjMNZx3ddXm0nM7PWxesp2dq/eiaaU/eSOEEJWhasoZNyHOBsPuuYTLbhuEolT+37zH7WHbsp1sW74Lj6dmV50XQpx9/KXPPivGKRsNAfTu9jhZOQfIzTuGwRBAwvElpGVs56fvtwHQpl0DOnUpfWXaAlsKqze+BoDVEkWPLg8QHBRXK/GfSYbLzbrsfMKNBvLcbjbm2mhmNXN5g1As+rMiPy+EEMJPmMxGJvzwKAe2HGbP+gMEhgWwaNYylv+yBtWjYg2yYMuzl1lHXlY+j144EYDIRhE8891DdOpf8k3c2laQa2PrvzuxBJhR9Apbl+ygQeNIBlzTh4Bgq6/DE0KUw5kWOKkvi58IUVU6nY6HP76b6x4fzvbluwgIsbJqznoWfbMMl9ONJdCMPd9RZh0el4eHBj4HQEhkMI/OuJe+V/SojfDPyGFzsGXpTnQ6hcCwANbP30JYVAgDrulDSGSwr8MTQpRDdffZU6ZM4aeffmLXrl1YrVb69evHK6+8Qtu2bUstM3PmTMaOHeu1z2w2Y7eXfU1zKkXzwZCUnJwcQkNDyc7OJiSkZleDLYvDmYPH4yAzex+bd3xcwdI6zulwG3GxfWskttOpmsZHCal8cTyNxOwcLN98Sv6/i0g+fgy33oA+piGWwUMJuH5M0Z3PCIOeH85tRYcguRgUoq6oK99/QpRXXfk3m5OeS15WPokHknlyyAsVKqsocMdrY7j24apPyVBef36ykJ/f+YOkAym0PLcZNz59DcmHU/nwkc9x2orP6WsNsvDC70/RZWCHWotRCFG207//TrwePv//MAaaSi3nynfy6yWf+fx7U5x96kqfnZ+dT1ZqDvZ8O3efW/E56K995ErufG10DURWsiXfr+C7134lYecx4ts14tpHh6PT63jzjmkU5NiKnW80GXjmu4fpN7xnrcUohChbbfXZl156KTfccAM9e/bE7Xbz9NNPs23bNnbs2EFgYMkLMM6cOZMHH3yQ3btPTu2mKAoxMTHlfn9nxcjb0phNhX8xVkskKWmbSExZW4HSKpt3fEJ+QTL5BUnodHoaRvciusE5NRLrxH3HmX40FYDst17CPm8OAPpmLVHy83Af2Eve9LdRTCYCrh4FQIbbwx3bDrG8T/saiUkIIYSoLSGRwYREBtOoZSxX3juE3z6YV+6ymgbTH/2C7NQc0o6lo3pU+g/vxXnX9K6ReTK/feUXPn3q66LX25fv5n9DXyqzjC3PzuSRr/PtsemygIMQddyJefLKOi7E2SwwNJDA0MIkxq2Tb2Dmc99WqPzsN37D5XRhy7VjL7DT89JzuejG8zCajNUe6/zPF/Pa2PeLXu/dcJCXbny7zDIup5sXb3yb749PL3qfQoi6qbx9dk5Ojtd+s9mM2Wwudv5ff/3l9XrmzJlER0ezfv16Bg4sffFFRVGIjS39af8zkasDQFF0nNvpHuIzLiAlfTPpmbvIzSvPxOka+w79VvTqWNJKmsRdSKe21XuXMNXpYsaxtKLXrq2bADD16kf4y++jOeykDL8AnA48yYleZffbHGzLLaBTcEC1xiSEEEL4yv3v3c6gmwey/Oc1bFu2kx0r95Sr3Hev/FL050XfLGfANb155ruH0VXj4oQOm4PvXvm5UmWz03LZuHArPS89t9riEUIIIXzppmeuoe+VPVjy/Qp2rt7LxoVby1Xul3fnFv15yfcrWfDFEl6a+z9M5upL4GqaxpeTZ1eqrNPmZNnPaxhy64XVFo8Qwnfi4+O9Xk+YMIGJEyeesVx2duEaWhEREWWel5eXR9OmTVFVlW7duvHSSy/RsWPHcscnydtTNIhoT4OI9hxMmM/OvZVb9fLIsUU0bngeYSHNqy2uLTkFuE6Z3cLY+Vw8xxNwrllB2v+NRMvPA6cDY+dzCbiueOK4wKNWWyxCCCFEXdChTxs69GnDom+Xlzt5e7p/f1zNmj830ueK7tUW1/H9yeRlFVS6/JnmBhRC+N6ZFjipL4ufCFFbWnRpSosuTdm8ZHu5k7en27x4O39/uZTLbx9UbXHlZuSRdDCl0uWlzxai7itvn52QkOA1bUJJo26LlVVVxo8fT//+/enUqVOp57Vt25bPPvuMLl26kJ2dzeuvv06/fv3Yvn07jRs3Ltf7kORtCaIju7Jz73dAOacDVjUse3Mw78/BujubxKRryGjcjIibb8LYsCFZP/+CJyuT4IsuIuyGG4qtxmnfs4e8fxahGI2EXDoEY5z3Qmjv7jzo9Trkof+BpmKf/zueQ/sLdxqNGFq0QRfsPUeHVafQNURG3QohhPBP3S/pgslixGl3Vayg0YguNISXHppF7NsLuPzGvjQd0JI5/24jJTGTPm3jGXV1P/Snjco9vj+Jf39cjaaq9L+qF/FtvfvsGRV8NPRUOr2OrheV/sNPCFE3SPJWiMrp0LcNYVEhZKXmnPnkErx738d898rPXHLLhZwzuh+/rNtOwpFUzo2LZexV52E2eo/KTTuWzpLvV2IvcNB7aDdadfUeYPXNy5V7UuaEXpfJkzJC1HXl7bNDQkIqPFf4uHHj2LZtG8uWLSvzvL59+9K378n1svr160f79u356KOPeP7558vVliRvSxAYEE2rZsNOTongUUFfwiOVmobxeAEx03djTPW+6+bYuZPE/z3jtS//32Wkvv8BLX79BUNkJADJr71GxqefFZ2T8sYbxD73LOE33EDB2rXM+fQLVl33f4Wrrfyn4IevsC/4A2OnroRNfgM1K5OM8bdj+/U7FL2e4PseKzr3mZaNMFXj46BCCCFEXRISEczdb9zCu/d9SrnXYDUa0cdEoeh0uFwqCfuSefOrf8hfuaHolJWHkpjx/b98PuVWGrdsCMDsN+bw8eNfFrXz6dOzGDPxOkY/dy37Nx/ixVFvkbDreKXfy8iHhxEcHlTp8kKI2iHJWyEqx2gycv/7t/PSjVPxuD0VLu92eji+P5n3f1hETvbhomvk5YlJzJy/is/uvY4O3VoCsOCLJbxx+7SidmY++y3D7r6EBz64g+P7k3hx1NvsWbe/0u/lolHn0bBF+RcbEkL4Rk312ffddx+///47S5cuLffo2ROMRiPnnnsu+/btK3cZSd6Wok2LEUToG5N4z0OQnIlq0ZN9SRwFXSNRA/RYd2Vj2p9D6JJk9LbydzyetDQS7ryL5j/+QP6qVV6JWwBUlaRJk8ldspSCZcv49P6nvBK3mt1G3oxpoGmYBwxCFxaBLiwCU6dzcKxYgmPDaoKBAJ3ClDaNub5hZDV9IkIIIUTdNOyeITTr3ITJ175BVnL2Gc/XhQSjnHJj0201kN8itNh5WaEWxt/5AT8sfJ5D2xOY/tgXxc75YuL37F67j82Lt1f68Um9Sc8dU27mmoeuqFR5IUTtkuStEJU3cGRf4tvG8dxVr5J0ILnkkxRKfQhWtRjIubCZ1zUyQF6EhQee/IS5v08mNyOPN++YVixBPOfD+Rzbn8SedfvJy8yvVPyKXuGGJ67i/14YVanyQojaVd19tqZp3H///fz8888sXryY5s0rPmWqx+Nh69atXH755eUuI8nbMtje+BLDgf8WCsv3EPXNQfjmYNmFysG+fTspb71Fzh9/lnyCppG/aBEAO5q39j7ksIPHDYB7787CfU4H7kMHAFAsVppbTazs06HKcQohhBD1xZLvVpQrcQugmE1er51R1mIXgSccVTSmP/YFq+duLLW+1X9sKPXYmYQ0COb74x+jN+grXYcQonZpUObK1eV8BkCIs9aK39aWnriFMv8jcjQLA2PJfWZ6lIWPHv2Snat243aVPMBqw4ItFYjUmznAxLfHphMUGljpOoQQtau6++xx48Yxa9Ysfv31V4KDg0lKSgIgNDQUq9UKwJgxY4iLi2PKlCkATJ48mT59+tCqVSuysrJ47bXXOHz4MLfffnu525XkbQncmZlkfv01efPn11gb6R9NL9d5kdmZJEdGFb3WhYZj7NIN15YN2P/+E8/ObXgK8lEz0wFoesVVzOrSskZiFkIIIeqaglwb82b8w5wPK9Bnuz1gOOUnUBnTLSiaxuw35lQhwpNMFhMetxuPu3Ah0ZimUUz86TFJ3AohhDgrOO1O/pm1jG9e+qnylZTVZ6sav743t/J1n0Jv1KMohVM1AIRFh/K/b8ZL4laIs9y0adMAuOCCC7z2z5gxg1tvvRWAI0eOoDvlKb/MzEzuuOMOkpKSCA8Pp3v37qxYsYIOHco/6FKSt6dxHjnC4ZtH406p/KqT1Wn40gVMv+pGr31hz79F/jcz0K9cgic1mQCzmZhu3bnhrnuYeMf/FVsQTQghhPBHmSnZPDzwWY7uSaxQOTUvD73l5Aqy5hQbeW000BXvP00nnsCpog792vL01w9iMBnYunQHwRFBdL2oE3q9JG6FqG9k2gQhKq4g18ZjgyZVaZ5ZAPOBTBSHG81cPJVh2Z1epbpPaN65CU/PGk9EwzA2LtyGJdBMt8GdMZqMZy4shKhTamLahDNZvHix1+u33nqLt956q0LtnE6StxSOtM1buBDV4SBv8ZI6k7gFuH7B76SFhjNnwGBcRiNoKr0axfDy9PfpFBzg6/CEEEKIWlWQa2P5z2vIzy5g5+o9FU7cAmgFNjyZ2ehCC+e+1Ts8BG5LIb9TtFcC17QnBePhjCrF26xjPPd/cDtdBpy8s37B9f2rVKcQwrckeStE+TjtTpb/spbM5CyO7kmscuIWQOdSCf1rH1lD24Dh5Mg28/4MAraUMRVDOcQ2j2bc1LH0uaJH0b7zr+1bRgkhRF3nL332WZ+8zZk7l+NPPY1mt9d4W7rwcDS3Gy03t/xlNI37f/yKMUvnk3zlCDqNuoGWjWRVSyGEEGefDX9vYfK1b5CfXVDlurTcXDx5eWAygkfFfMSNYddRnC0aoBn0hGXZCHWqpFawXkWnYDQb6di3Ddc/eRXdB3epcqxCiLrFXy4EhahJu9ft59lhU8gs53z0FWHdm4Hpkw3Y2jdAtRgITi6gsWLimFqx2SsVnYLeoKdtz5Zc89AVDLi6T7XHKoTwLX/ps8/q5K07LY3jjz+B5nLVSntqZibmzp1xbN1a7jL6Bg1o+uUXmCuxgp0QQgjhL+wFDl64/s1qSdwW0TRwOIte6nMdWDcfA8AJRPZpQ2pC+R/BDAix8uqC52jbs1X1xSiEqHP85UJQiJqiqiovXP9mjSRuT9DnOQlae7zotfmcZhUqbzQbmfzrE/S45JxqjkwIUZf4S599Vidvc/6cW67E7RcZGfyck81xlwu7phGh13OOxco9kZG0tVgq1GZFErfGuDia//or+iCZFF0IIcTZbfUfG8jNzK/VNnet2lPuc4PCA5m5+x1CG4TUYERCiLpA0xS0Mi72yjomxNlg+/LdJB2s3akID2w+VO5zjWYjn+14i9jm8kSrEP7OX/rsszp5qxaUb/TOWlsBGW4PjY1GHJrGIaeT+Xm5rC7IZ2HLVgScsopcdVFMJuLemSqJWyGEEAKw5VV8eiMN8MRF4IlvAHoduuQsDAdTUDxqtcam0+t49vtHJHErhBBCULk+G8AdG4a7WRSayYA+NQfjviQUl6eao4MnPr9PErdCiHrF75O3tu3byf1rHpqmYWwYC243xqZNCRowAHO7tuWq4/WGjTCfkqB9Jy2VD9PTyVZVDjqddKzg6NtiFAXFakUrKACDAes55xD3+msYGzasWr1CCCFEPXJ451EWfbMMp81JTNMoPG6VqPhI+gzrTosuTUChMCNbTq5zmuFpcfLiTI0OxdM4EvO/O6uUwLUGW7Dl2tHpdbQ4pylPfvkATds3rnR9Qoj6RUVBpYxHMMs4JoS/SDyYzMKv/iUvM4/YFjF43B7Co0PpN6IXzTvFozfq8VQg8epsH4erY3zRazU6FHeTBlgXbUdxuisdZ0CIlYIcG4pOIb5tHI9+dg/te7epdH1CiPrFX/psv03eerKyOP6/Z8hbuLDE4/qoKDzlXDjMrNPxd24un2Skk6eqHHIWzo8XodfTzGSserAWC23XrsGTkYEuNBSdyVT1OoUQQoh6wpZn470HPmP+54tLTM6GNAhGU9UKJW7VEKtX4vYELTwIT9MoDAcqvyL1d8c/xp5nJyDEitlqrnQ9Qoj6yV/mzxOiMpx2J58+PYtf3p2LWsKN0MDQAMxWU4USt6rFiKt9XLH9WrAVV+uGmLYnVDreDze8hiXIgiXAhDXIWul6hBD1k7/02X6XvNU8HlI/mEbG9OllzmfrSa3Y+tFpHjdb7Ccf/2hsNPJ+XGMCdfpKx3qCpXVrFL0eQ1RUlesSQggh6gtN0/j53bl8+vTXOAucpZ6Xk1a+m62n8kSHlnmsssnb2ObRWAMtWAOr+NSNEKLe8pf584SoqL+/XMJ7D3xW5uKh+dkFFV5cVI0KgVKmIvTEhEIlk7dBYYFEN2mA3lD1a3YhRP3kL3129U/W6kM58+az94ILSH///XItRFYRN4SFs71NW/5u0ZLLgoM56nLxyPFj5KtVn4MndPjwaohQCCGEqD82LNzKLa3vZ9r4GWUmbitLcZfeP5d17EwGXtOn0mWFEP7hxCiesrbKevnll1EUhfHjxxfts9vtjBs3jsjISIKCgrjmmmtITq780wNCVNTudfu585xHeOWW9yqcmC2XsvrlKvTZPS7tKolbIc5yNdln1ya/Sd7atm7l2EMP4UlNq7E2FEWhkdHInRGRAOxzOvkjp+KjgU6lCwkh7NqR1RGeEEIIUS8c35/Es1dMIbEKUxecif5YRqkXfPojlfutYDQZuP7JEVWISgghSrd27Vo++ugjunTp4rX/oYceYs6cOcyePZslS5Zw/Phxrr76ah9FKc42ORm5PHHxJA5uPVJjbeiTslHsJd/INRyu2BOzpxr7/A2VLiuEEHWJXyRvNY+HYw8/Amr1rh4NkOXx8Ft2Nk7t5ER7S/Pzi/5sq0qbej2xE56TOW6FEEKcVV4c9TZOR/U+IXM6xeXBtHYfnDrnnqZh2HUUfUp2peq887UxhEQEV1OEQoj66sQjmGVtADk5OV6bw+Eotc68vDxuuukmPv74Y8LDw4v2Z2dn8+mnn/Lmm29y0UUX0b17d2bMmMGKFStYtWpVjb9XIV695T3ys2012oaiaZhX7oFTfxtoGob9yRgOV+6G66inrqJRy9hqilAIUV+Vt8+u6/xiztvMWd/gSqj8JOZlyVc9PJmUyMTkJOKNJnJVD0nuwtUuA3U6BgdX7iIudMQIIm69BUu7dtUZrhBCCFGnLZm9kj3r9tdKW/qkLCx/bcTTMAz0enTJWehsFZ+iof+IXlz32JV06Nu2+oMUQtQ72hkeszxxIRgfH++1f8KECUycOLHEMuPGjWPo0KEMHjyYF154oWj/+vXrcblcDB48uGhfu3btaNKkCStXrqRPH5nKRdScrf/uZPUfG2qlLX16HgF/bMDTKALNZECfko0uz37mgqc5d1BnbnhiBN0GdznzyUIIv1fePruu84vkbfYvvxT92RLhxJlnQHVWz6DiYJ2ey4OD2Wq3k+By4tY0GpoMDIzX82RfC810eWTsCsKRbSx3nUGDBtHo5SnVEp8QQghRnyz4YnGttqe4PRgS0itdvl3v1kz86bFqjEgIUd9pwCkP5ZV4HCAhIYGQkJCi/WazucTzv/32WzZs2MDatWuLHUtKSsJkMhEWFua1PyYmhqSkpApGLkTFLPhiSa22p6gahqOV77Mbtohhyl//Q6+XeW6FEIXK22fXdX6RvFXz8jAFu4nrn4ElzI3qVsjYE0DWgQBUl46gODsB0Q5SNobicVTsizxEr+f1RnFFry2RDpoNTkdRoPCv2UZIvJ3DiyKxp5c9/YFiUAloHUeju4ZW/E0KIYQQfqAgp2YfvaxOkXERDLvnElRVRVfKKthCiLOPioJC6SN11P+OhYSEeCVvS5KQkMCDDz7IggULsFgs1RqnEFVVkFsDi5PVkNCoEEbcf1n9ycQIIWpFefvsuq7+Xom4HbD+c/j2Jhr1zaDpoFQsYYXTGegMGg065NPqilTaXJVMo17ZhDWz0+KyVILibFT+G12jcf/M/xK3J+kMGtGdc85YOrJ9Hk3OWY/+pxvg25vA465kHEIIIUT94fF4+GfWv0y+7g1yM/N8HU65pR/L4LVb3+fRiyZiLyh9rkohxNmlOufPW79+PSkpKXTr1g2DwYDBYGDJkiW88847GAwGYmJicDqdZGVleZVLTk4mNlbm8xQ1Y8Wva3npprc5vq/+jO7OTs1h2kMzubfnE+SkV21RcSGE/5A5b33J7YAvr4LDywGwApTjRrXBohI/IJP03QGkbAzFGukiol0u9gwTqkfBGOgma38gzpzCKRB0BpXYHtkERDvwOHXkJpowBpS8QFlA9Ik59E4khr3/ASh6lbAWp9y53PU7bPwCevxf+d+3EEIIUc+oqsrz173J8p/X+DqUStu6dCezX/uN0ROu9XUoQgg/M2jQILZu3eq1b+zYsbRr144nnniC+Ph4jEYjCxcu5JprrgFg9+7dHDlyhL59+/oiZOHn3r57On9MX+DrMCrtwObDzHjmGx6cdqevQxFCiGpTP5O3m78pStxWRmTbAqwRLizhLnQGCIk/OZomtKmd/b9Ho7oVGp+XQWBsYVLWGKAWjewticepQ2dQaXJRGq58A0nrTk7RoLd4aNQ7C6P1tMTvopehy/VgCqz0exFCCCHqstV/bKjXidsTZr/xG1c9eDlBYdJnC3G2UzUFpYyROmUtjHK64OBgOnXq5LUvMDCQyMjIov233XYbDz/8MBEREYSEhHD//ffTt29fWaxMVLvda/fV68TtCX9+spAxk64nPDrU16EIIXysOvtsX6qf0ybsmV/lKgKiChO3pzNYVEKa2ghvk1+UuC0PW7qRZhenYY1wExJvp/WVyTS5MI0mF6bR+spkghqW8LhlfjK81wOyj1XhnQghhBB1V22tUl3TbHl2xrS8j8M7j/o6FCGEj2nambfq9NZbb3HFFVdwzTXXMHDgQGJjY/npp5+qtxEhgDV/bvR1CNVC9aiMbjGOHSt3+zoUIYSP1XafXVPq18jbXX/Cqg/g2PoabSaqcw6aWrHse3Ccd3JW0UNgTDmSvznHYeEkuHp6hdoTQggh6rKN/2xl9htz2L5sl69DqTa5mXl88OBnvDL/OV+HIoTwoTPNkVfV+fMWL17s9dpisfD+++/z/vvvV6leIUqze+0+vn3lFzYt2ubrUKqNo8DBW3d9xMdb3vR1KEIIH6rpPru21J/k7boZ8Pv4WmnKYNGo1WUqd/xWmLx12WDjV7B3AZgCCqdUaHtZ7cUhhBBCVIPlv6xh8sjXUdV6ciu7AjYu3EZeVj6WQDMLv/6X5b+sQafXcf61/Tj/ur7odPXzoSYhRPn5y4WgEABb/93JExdPxuX0v8W0D21LIGH3MRq3acTSH1ax+LvluF1u+l3Zk8GjB2I0GX0dohCihvlLn10/krceFyx6qUabWHrYzcvLnKw97iGtoPBic9pQC3f3MNVou0VcNvh8GBxde3Lf9p+hzzi4tGbfuxBCCFGdZjzzjV8mbk9QVZWJV7/mNSXE8p/XsHLOWp7+erzvAhNC1Ap/mT9PCIAvJn7nl4nbExRF4fXbPmD+zMVF+1bNWc/i71fw0h9PozfofRecEKLG+UufXT+Gh6Tvh/yUGm1iQ6KHBQfcRFh98BfXYXjhiNtTE7cnrHofUvfUfkxCCCFEJeRk5HJ4Rx2eF1anA70elMr19+cO7syWJTtKnMt30TfL2bx4e1UjFEIIIWrNNj+a3uh0zTs3ISc91ytxe8KGBVtY+sOq2g9KCCEqoX4kbwMiCyeRrUGjuxjJeTKYeTcH1Gg7xYQ3h8ETYN/fpZ+zr/6v+CmEEOLsYA2yYAk0+zqMIorBgGI2o5hM6KxW9AEB6K1WdAEBKKaKPV0TERvGuKn/x9q5pS/osuZP/1igTQhROn9Z/EQIgPCYMF+HUCMCQwMY/9FdrCmrz54rfbYQ/s5f+uz6kbwNioIWF9RoE5EBOqzGWh512/NOuGc5hDQCYxlJ47KOCSGEEHWI0WTkguv7+zoMUBR0AQHoLBZ0RiM6kwlFrz/lsFK4z1C+GaQGjuzDpzvepkm7OCyBllLPK+uYEMI/FF7sKWVsvo5QiPIb8n8X+jqEatd5QHtm7nmHDn3alN1nB9Sdm81CiJrhL312/UjeAjRo5esIqo+ig/7jYehrYAos3Nfl+pLPNVgKp1UQQggh6on4to18HQKKyYRSjsXDFGPZi5UoClxyywU8+/0jBIUV9tkX3TSgxHN1OoULbzyv4sEKIeqVsi8Cy14YRYi6pnFr3/fZ1an35d14c8lkwqJCAbjwhv7o9CX/Hhh088DaDE0I4QP+0mfXn+StPcfXEVSf2xfBxZO897W9FPreB5zyD8dghas/hoCIWg1PCCGEqIqc9Fxfh1DuEbVnmvt24i+P89iMcV772vZoye0v3+x1MWgw6rn//Tto3LphhWMVQtQvWjk2IeqL7FT/uc6+/73beOH3p7z2xTSN4uGP78ZoOvm7QKdTuHXyDXTq3662QxRC1DJ/6bPLeWVTBzTpA5u/8XUU1cNkLXn/kBeh+1jYOx9MAdD+SkncCiGEqHc6ndee71791ddhlIvm8ZR5PCCo5D77+seHc+EN/Vjx6zp0eh39r+pFZMPwmghRCFHHnGmkTn0ZxSMEQMfz/CeBGRQeVOL+IbdeSK/LzmXZz2twO930vbIHsc2iazk6IYQv+EufXT+St448OLjU11FUD70RotqWfrxBK/+aIsJHjqQX8OmyAyzbm0auw02g2UDv5hFc2bURaXlOjDqFPIcbo17HRe2jCbEY8agah9PzCbEaCbUaWbE/HYfLQ5vYYPYk5RIRaKJ703CUSq5QLoQQZwOX08X6vzf7Ogw0t/uMUyJomobmcpV6XFEUOpVxURvdJIoR919W6RhFoYyUbH7+dAkrV+8jI0CPPsxKpw6NuXRAB3Lz7ShmHVmhTlDggrgWRFoC0DSNY0lZWMxGIsICWbf/KDk2O+1jGnBkbypWq4nOXZugL+VRWSGEEKCqKuvmlb6gV33T+/JzSz0WHhPGsLsvqcVo/FOWzc6XazYyf+9e0ihAb9LRKSaG69p2wmH3YNC7aRyVCIqT+IAeBBmjAEi2p6FTdESZI9idu49sVzaNLU3Yk52NQaejV1RTTLqaXaReiPqsfiRvZ10Ph5fVaBM/7XTx+AI7bvXkvucWOXh9hYPejfV8fXU1LRoW2rR66hEl2no0m69WHebnTcdwnvqXmevgYFo+365NKFYmwKTnup7xLNiezLEsGwB6nYJHLT6AvlGYhddHnkO/Vg1q7D0IIUR99uot77H4uxW+DgPN6UTT673mvdU0DVQVFAXN4ylM3KpqqXWERAZhMNaPn0r10Z7MNL7dspG/fl4De214QkPJa+3GHmNnA5l88/s2FJeCrakK/13PmXR6bo49h52/HiHheCYeI7gaGnHy39+jpmHKVglI9RAeZGb8w5dxYa8ybpoLURlnes6yvjyDKc56Hz78OT+/86evw6gWJouRwNBAX4fht44VZPDT4bV8t3kL6eke7DqFhjHZhIflkajbwqt7FhOmy6dPwwPsTS58qklBR/PgYSxIKeBQwVEUNKx6BZXCG+eqBim5wRzJCqeBxcDTXQcxLL6XL9+m8Ed+0mfX/SuSg0trPHELkOPQ2J/p/beWWqCRWgCNQ6rxb/OCJ6uvLuHltXm7eH/R/gqXK3B6mLn8kNe+khK3AMez7Nz4yWq6NwnjlZFdaBUdXJlQhRDCLx3ekVAnErcAaBpqQUHh3Lc6HZxhlG1JrnpwaA0FJ2ZsX8+kVQsLfy93NkJnAzpNRTOdPMcTpBUuEXzKAy9O1cNnxzYQZdOhMyvYI3WoOhWPGTQdKB4Fj1FPQUMdKWaNe+b+SbsNa3l+2GA6x8XW9tsU/upMC5zUk0cwxdktMzmLOdPm+TqMaiMLhtacv45v5rktP+DRVAiFgFBoYcklyOIsOifQ4qBP6AEMupM3xTVUDuT+So4rHgjCqHOjoBJicGBUPLg0PcZQNxHWAixGN98f38/i9F8Z2/wmOoV28ME7FX7JT/rsup+8PV47j3Hc2tXErV1NZz6xKoIbQscRNdvGWWrbsexKJW4ra/2RLAa/uZTW0UFMubozPZqVPjfx6gPp7EvNo1VUEL1bRNZajEIIUdv2rD/g6xCK0dzuSpULCA3g6gcvr+ZoBMCxvByeX/2P90AHBbSSZrko6fe0AtltwHJUj2oGd+DJ8zQjhSMoTpniaFdyKtd9NIu2uYE8ceMg+vYsfXqq3TuOs29PEjENQ+nWswU6Xf34QS9ql/bffYWyjgtR1x3Ychi3q+x53+sLo9nI2Odv8HUYfinPZefFbb8UJm7/YzW4ihK3QXo7saYsAvUOjDo3GsWnK4ozZ5LjtmLWuWhgyudE12rBTaDeQQYaLq0wNZXiSOalHW+RcKgTD3e7mKvalp7E3Z2Vwoa0Y0Rbgzi/YUsMOpkqSRTnL3123U/ehjb2dQTVJzcRlr0F5z/u60j8zh9bE33S7t6UPEZ+uJIhHWO4tX8z+rY4OZ1CZr6T//t8LRuPZBXtO7dJGJ/d0pPwwMIbBSv3p3MwLZ82MUFlJoCFEKI+iG7iP1PKFGQX8PETX/PA+7f7OhS/M+/QHjyn/1KuYI5UMwJ6cFtLKFtCXZpeYY+1gDu/+Y2Lf47h+iHd6X3+yfmM7XYXk5/+gbWrTt4Ibtq8AS+9MYro2FAAtu9NZN+RVOKiQ+neqYnMgX8W85fFT8TZLSref/psl8PFu/d9ysQfH/N1KH5nRdoeCjxOr30WfeGTTOcEHqFtQGLR/VJVg2w1EMdpd2NNOg8WvZswo53T74nqFAgx2kl3nlxsTq/T0Acc46GFf/LF1o3cfk4PLmvVBt1/DblVlUdW/cacwzuKysQFhvLp+dfRJrRwjt2dGSlsSkkkJjCI8+Oao5fE7lnLX/rsup+8bXcFhMRBzjFfR1I9Nnwhydtqpqoam05JkPrCvO3JzNueTNPIAH6+px8RQWYmzdnulbgF2HgkizGfrea1kefw6A+b2XYsp+hYx0YhzLqjD6HWshfYEUKIuqrLwA606NKUA1sO+zqUarHgi8WMmzoWvUEW0KhO23ZU4TedoqHTq+g1UHUa6L1/cGvgnbw9ZdYFTSlM4v6dn8zmZ3+kQWgAb390Kw0bhfP5x0u8Erdus8I2TxZ3v/EdrzxyNVM/+YflxxNwBisoHmhtCmX6o9cTFVHyyubCz2lK2Y9Z1pMLQXF2a9Iujm6DO7Ph762+DqVarPxtHdlpOYQ2CPF1KH5lw9HjxfZ5NB0NTZm0C/QeQKVTIFSXT6onxGsEbr7HjIKKSVfySG+TzoNBceHWDJzotS1mFyiwMSWRcQvmELbQzNdXXkfHRjHM3L3WK3Gr13mwcZyH1k7n3T638tqaNSxP3UxwoA1NUzCvj+bTC2+lRagMljor+UmfXfdvPxjMcPNP0PAcX0dSPQoyfB2BXziWZePN+bt56LtNXD99JSsPpPs6JAAOpxdwzYcrsTk9/Lk1qcRzth7L4dKp/3olbgG2H8/hjs/XYXd5mLX6CPd+vZ7Hf9jM6jry3oQQ4kwUReH5356g88D2vg6lWtjzHTjtzjOfKMqUkZXPFz+u4vmpf/LIc9+zZtbOwuE5p1Ih2GijdXgSMQHZJ/drJ08wB9gJDisgKNiGvmUBrm75aCaV8AY5dO27h/Mu3USvC3bQsFkqmp6T2381KP898amaFFwhetJzbNx/90xcLg8L5m4pajIn3kDiQCtZ7c3sCMrnio++ZL56lNymOhwRCvYoha2hOdz20Wzcbg9zV+zg6Q9+Z8L0uazYcrCmPkYhhKh2T88aT++h3fziSQLVo5KXle/rMOq9HJeNL/av4OkNPzFuyXd8tmQr2mlru+Y5TTS1pJVYXqeARXFhwINVcWBV7MSYsmlkzkLTIEyfT8+AAwwJ2coFQTtpYUpB0zTMeg9WvRPlvwVI8wosXvVmqQ5Gffsd2TY7Pxw82WeHWgto1yiZuIhsDJZEHto0hQMsIb5hBmEhNsJDCwiIOMR9qz5A0zTmHtrDA4vncP+iOfxxcDdqfXlmXtQZU6ZMoWfPngQHBxMdHc2IESPYvXv3GcvNnj2bdu3aYbFY6Ny5M3/+WbHFIuv+yFuA6HZw11JI3QOObFjyKuyd7+uoKqf5AF9HUO/N3ZrI+O824XCXvkK4Lx1My2fZ3lScnorHt+ZQBh0n/MWpRb9fd5THhrRl3IWlz9EnhBB1RXSTKN5cPJnj+5PIzcjj53f/ZOFX//o6rEpp06Ml1iCrr8Oo1zZuPMCTr88h3+4qXOxXBxgNhO1UyeoIoBFitnFB0910jjqGTld4EXUkJ4Ifd3cnP9cKOrBG2jBbTo7Y0elU9KFOAtoV0KXxUZT/hiMYguy065iAyeTm0P6GhXPp6kBxF24nnBhkkZVtY9nSXeQW2HEF6cmN15Mfb/CaM1fTg2osntjYQSY9n3oHh0HDYAdTNsxdsZPrLz6XR266sHo/SFGn+Mv8eUKENgjhhTlPkXIklczkbBZ9s4wf3/7D12FVSkzTKBq2iPF1GPXa3qTj3LH2C9I1W9E+U2twevSYFA+KAiadm1CTHQOlX+tqGlh1J29+hxsKCNHbSHaE0DPwEHql8EvSoHfS3pqIRediVV5LdAqY9W5y7FaS00OL1Zurd/Lrxh141CwuiN5FsNlGriGAdLf3IuIWkxuPqqCd8iiOzprMkN+fRDE5ybebSMoM5bcDOxnWvB3vXjjML25giJJVd5+9ZMkSxo0bR8+ePXG73Tz99NNccskl7Nixg8DAwBLLrFixglGjRjFlyhSuuOIKZs2axYgRI9iwYQOdOnUqV7v1I3l7QlSbwv+/4RtY+R78PcG38VSUKQgueMrXUdRb2TYn98/axNK9qb4O5Yy+W5dQ6bIl5XzfWrCHkd0bExNiKX5QCCHqoEYtY6ElPPH5/XQf3IXXxhaOeKirFEXxis9g1HPbSzf6MKL6zel08fK4z1h4IAlPWCBuI9iiFDwWBb0TzGkKnRNT6Dx4HzGhWQTpvUc4NwnJ4Mrmm/lmYz90igeTufijljoFmsekFSVuvco3T+ZQUgNUuwGdpqA4Tk6foHg0NIMeV4CCwebhr2U7yYsx4wiG/PjilWmlPaemU8htAKDgBOyRGsGH4LsFGxk2oBNtmkRV4BMT9YoGlPV1Vne/6oQoUXSTKKKbRNG2Zyu6XdyF54a/isdddxczU3QK2ilPcOh0CrdNuQmdzGtaKaqm8fKUr/nJswXHOVZAw6BT0elU3KoOj0eHETdtw1MJMBbOd5vlCSCOrGJ1aRoYFA+nT0CvVzTaBSQVJW5P1dSUxmZdPDbVhF7RyHKaiYjPRvUo5NtM2J0GNJcecg3kOBbyv46foVdcReU35jVhYZb3wmaKohXNZapXVEJNNnQxBUXH46My2XwgnjkHd3Ftm86c37h5JT89UedVc5/9119/eb2eOXMm0dHRrF+/noEDB5ZYZurUqVx66aU89ljhvNzPP/88CxYs4L333uPDDz8sV7v189tNb4DzxkPzkj+YOqHlRXDTj9DmUojpBOfeDHf8A426+jqyesfu8vDQdxs5d/KCepG4Bfh7Z0q11udWNeZvL3kaBiGEqMsUReHiMRdw4aj+vg6lVM07N+HNJZMYcE1vmnduwoWj+vP28hfpNriLr0Ord1weD6/8+A99nniHHy255DUPxGGF7DY6HA10uIMUHBEKhm659Bm+jZAQG1adq8S6WkSmEGopQG/2UNqAmCCTo8T9BoOKNdiBatVQjRraiWmLNQ1LpoamV/AE6CmINjA/7Qi5DRUKYqo26kY1KtiiC//896pdVapL1G0nFj8paxOivup1WTeueuByX4dRqoiGYXyw7hUG3TyA5p2b0P+qXry6cAIX3lB3f2fUVaqm8sOB37n1n3Fs7rWM0D4uFMVDkMVBoMWJ1eQmyOwkIjCPc6KPFyVuARJdYWS5iz+ddMDWoNRRrIG6kvtsnQLh+pNTXgSY3CgK6A0aIcEOgoPs6CwqQQ1y6dDEO3ELcG7QEVpbS79WDjI6ii2UZjJ4aNWo8Jr91707Sigl/EV5++ycnByvzeEo+d/r6bKzC6f7iogofU7llStXMnjwYK99Q4YMYeXKleV+H/Vr5O3prv0cfh0He/4CTYWASOh6M6yY6tu4LOEwbCqENYHWg898vihRWp6D1/7axQ8bjpY4GvVss+pABqP7NvN1GEIIUSnjP7oLRaew+NsVeNweAkKsXH77IH586w+fjsg1B5h56qsHaN65KZ3O84+5en0hx2Xjk81LmbVrPTaXhjvOjOb4L2Pq1ggwOmkSmUae3cLxrEjatT6KXl/4964rZciDokBggI18W3CJxwGcbj1Wo7vYflUDh9MICqgmDS3QjRakYcwEg12PpgNnAOQ2NVDsiu4UWtH/FNIHujBG2dBb3WhuHa50C64MMydGGDn/C3XrxiNwXanVCn8go2uFH7ttyo04bE7++uwfXA4XZquJS//vIv78ZCEuR8k33GqD3qjnf988RKuuzXnyiwd8Fkd951Ft7Ez4hk2ZX+Ow2GkVHU6eaqUpkGwL4lB+JOHGAlyajmxXALEBecW6ShUd6/Oa0diUSQNjLhoKKa4QUp1BxFpyS2zXpekxKiVf2Oer5sLYNAWdohJuyset6clzmbGYPNgcKt0ijmDUl/zvr31AInttsUWvNU0h1FBAm6BkIk35uDQ9qc5gjjvCONFnhwcVoNd52HpIBkn5vXL02fHx8V6vJ0yYwMSJE8sso6oq48ePp3///mVOf5CUlERMjPe0LjExMSQllf/fXv1O3gZEwKhvIDcJ8lOhQZvCBc72L4TkbTXbdnBDaHYeBDSA2M5wYBHkHIem/aHXHRAUXbPt+7H9qXlsPJLFpN+2keuou4/r1Lb9qXm+DkEIISrNGmjhyS8e4O43biH9eCaNWsViDbRwcFsC6+dvrtG2gyOC6HHJOQSFBdKud2u2LN3B8f1JtO3ZihH3XUZMU3m8vbISE7PYlXicF1N+J82TC+FgBAwN7DgPBeHJMXHpORu5rOMGzP8lWQ+kxLDfefJ3kksrnEvvdDa3gbxAPUaTA9WjoNN7//IOMtiJM2cSZrJhU02kuYJwaYU/bTOygnG6jICGLsqJZvXgAlyAvYtC8FIzedH6MyZui0bsqqALcmFplls0CljRq5gbFaAYVZzJAUWFPAaNxF3V+wSOqFvONLpWRt6K+s5gNPDA+7cz9oUbSD6cSsPm0QSGBlKQa2PBF0tqtG1roJleQ7sTGGKlQ9827Fl/gIPbjtCsUxOuuv8y4tvG1Wj7/iw5O4+DaQnYeBS37gAx/00rG2/KYIc9juOucJoGZtAuOAmjrjDJmuMyc9weVkqNCinuEOwUTneAAmEmO3luM0EG71GLQYqHdkaNhgbwaJDsgaOewr420x1AticATQO720iAwU2AofA3g0fNIzE/hDyDGbO++M3aE4ynTGzvVhWCDXb6RBzA8N80DXrFTbwlE7POxcETj8loGg3CcknZoUfTNJn31k+Vt89OSEggJCSkaL/ZbD5j3ePGjWPbtm0sW7as6oGeQf1O3p4QHFu4nXDt5zDrOsjYX7V6I1tD5iFQT7m7o+ih80i45AXvBO25N1WtLUFKrp0HvtnIqgMZvg6lTtqdlMtvm49z5TmNfB2KEEJUWlhUKGFRJxehePSze3l22Mvs23iwSvXGtYol9WgGTvvJuVMVRaHfiJ7c9doYr0VMLrnlgiq1JSAv387Lr/zOylX7yLvIgaOz9wWVogNj43y65h9mRNc1XsdaRCcT5shjr6fwt1u+asaoFBSbGmFVWgs8mh6dATyqAh7QmQovwmLM2QyP24JJfyLpW0C0MYfdBbFkuwPIdhbOEa8EudFZvRPDqlUjt58DDgeeXERNATRQ1JOz9J3Yr1C43xxlK3H6BmOkHXuWCcWhx+BQcAUrHAzUmLV0IzcOPLdCn6sQQtQlweFBBIcHFb2+9+2xpBxJY/Pi7VWqN6ZpFDmZedhybF77u1/chTtfG0OLLk2L9g0Ze1GV2hLgcLmZ+PPf/LF5F5d0W8OQbge8jisKtLUkkuEOIOS0pGuI0YFZl0aCIwLveWw19IpGqMHmdR9Ur2hkeIJQoaguCyoXB+QQ8N+ipCjQTAdBOtjk1JPhDqDAbcSp6jHqvG/U6nUasYE5pOQEcbCgAZpGiX3x7vwYHB49DreefJeZ9lGJRYnbU0UZ8zhuD8ODHkUHXeKP42mUyE9Hf+Ka+GvK+5EKPxQSEuKVvD2T++67j99//52lS5fSuHHjMs+NjY0lOTnZa19ycjKxsbGllCjOP5K3p2vQCu5bBwcXw9YfYe98yC/nCAiDBQZNgJYXQnR7sGXCzjngyIMWFxYumqbTn7keUWH3frWBdYczfR1GnaUBD323ifaxwbSOKf0RUiGEqE8aNIpg2vpX2bJ0B0tnr2TFb2tJTUgvV1mdXsfo566lzxXdaXVuc2x5Npb9tIac9Fw6D2xPy3OaoTdIn10TXnvjTxZvOYA9zoCrfQEljWfQmTTOb1zyBX6EOR+rzYFNM+NGT5YnAKvOiaJpJNtD2JzZmH25p9wkVxRcCRYUux5FrzHg/A2nJG4LGXUqsaZsjtkjMJkKj+kCSh6lowZraCEulHzjKSuZgaYA7v9G3ygaeK1UXfLTQIoOiHDhsakoZh0FkRr6PB1Tfl1Em0ZR9GhV9g96UQ/JgmXiLBUUFsjr/0xk15q9LPtpNUt/WEXigeQzF/zPyEeGMXBkX9r1aoXL4WLFr2tJO5ZB254tad+nDQajf6YnfO2Nv5byZ+IGLOcU0L3TnhLPMSgqMcYcbJw62lAjQHESZFTJcVvI9gR6lTHrXCU+wOJBz9a8eAw6D1adm6tCD59M3J6igR4UCkfU2j1G9KdNqWDWuQgz2NApKmoEtAhP56AzihZm73VwjtgiWJbRGodqRNUUQCHCVEBJFAWCDXav96LXa8xNnkOTwHh6RvQqsZyox6q5z9Y0jfvvv5+ff/6ZxYsX07z5mRe769u3LwsXLmT8+PFF+xYsWEDfvn3L3a7/fjvqdIWLhrW8CDxuyDkKBRmwfgYkrAWjFYJiQG+CvGRw5ELDLtDvfojpeLIeazh0G+O793EWcLg9fLRkv08St/k7lpCz5idc6QkoBhOWJl0Iu2AsxvCGtR5LeXhUja9XH2HilR3PfLIQQtQjXQZ2oMvADox75/9IOpSCLc/O3E8WsmnRNoxmIxExYZgDzWQmZ5GbkUeT9o255qEr6NCnTVEd1iArF48534fvwv95PCrzft/I/J0HcIfpyeyiEcjpa0qfFGop+eIJIMkWAnqFUKMdu9vIjgNN2JjRGFdkKXPgmjS0FCPWABuxwTklnhNpzEeHis1mQlE02jdIpGV4GhaDi3RbENszYkn9bw5dV4QHxaDDmHNKgl8pHGUbvlslo4OCdsovZc2lQzEUT+BqGmiqgmYBt7XwwtMd5sHVwMNXKzZK8tYvKZT+r54zHBOi/mvXqzXterXm9pdvJvlwKvYCB39/sYQ1czei6BQiYsMIDA0gMyWb7NQcYptHc9X9l3stBGqymLjgellkrCZpmsbi3Qf4MftfzG0dxAdnlroA6OlMios4Y2bRtEYNDdnss0exx1Z4nayhlDJnfeE+o87DUVs4wUY78cbSfwtE6DwkA+eF7+OgLZJ0V2EfHWHMp6E5G0XRMKASE52LXTOwzxlDpieQRsZM9KgkOCL4YNtAggJcGEwnk7/5HhNh2Eps06WVnAabl7RAkrd+qXr77HHjxjFr1ix+/fVXgoODi+atDQ0NxWotXMRvzJgxxMXFMWXKFAAefPBBzj//fN544w2GDh3Kt99+y7p165g+fXq52/Xf5O2p9AYIb1a4xXXzdTTiFOl5Dq7+YAWHM0r/Qq8puZvnk/HXOwAYQmPw2HIp2LMC+9EdNBr7Lvqg8FqPqTxW7EvzdQhCCFFjFEWhYfPCKQ7GTf0/H0cjTmUrcPDwPTPZnJOBK8KMO0DDEwBOuwFLUPEFRFS7jgMpscSUkGj1aAqH7A3I91hQ0MhLsXJse0OUBk702EtsX8sr/NnqNCulPjapoeBRFVIPhdOz2QFaRp9cCCIqII+B1v0sOtqadFsgikFDC3XhcYLefjKBq5oVDPkurMk6CuJO/lR2pZkxxxf/veJy6FEVXbF4VAusLjhS4nsR9ZyMvBWiyIl542+bchO3TZGpBOsKl8fDXbN+ZkX2XkJaOTDq3AQZnRxyNCDSmF/sfKeqJ90dRIDBBWheiVsAnaLRxppCtieARGcoVp3b66tOQUWvqChK4U1Ns86FU9WT7TSTpxkA5+lNAmBTFQpUCzpFo4U1jQLVjFvT09CcjUXnIlDnQPff9AcuTUeaK5h0T+F2gsXowZlhwRB7oo/W2J/XgDhrdrH2Mp1WCtSS5zI9kHu8zM9U1FPV3GdPmzYNgAsuuMBr/4wZM7j11lsBOHLkCDqdruhYv379mDVrFs888wxPP/00rVu35pdffilzkbPT6c58ijebzcayZcvYsWNHsWN2u50vvviiolWKs9DmhCwe+nYj/V7+xyeJW83jImvJTAAC2vQj7u5PaXT7NBSTFbUgi+xV39d6TOWVluc480lCCIH02aJ67E1P5+n5C+j/3nT+aZ5Pbgsz7gAFe7SKLsiJCx0ej3fmUvUoOPcHsWBTV2xOY7E6t+TFk+8pnJdWQ0E1KrgDwWU3ojqL/zzVnAqqXkM1qBRoZg5lNCgx1mO2MPZtikdx6mgeX3wFX52i0T48CY+7cEJbRQ9qg9Pm63VrGHM9mHJAf0oe2ZNuwXXMivbf6ZpWmLguyLGUOpIpy1RyIlrUc1o5NiEqSPpsUR2OFWTw+o45XDL/FdYGbcEYV9gPWQ0u9IqHA45oEp2hXmXcmo6/M9uzKr0ZBW4jAYqzxIVEAZqZ02hgzCfY4MCo86BqhYlbg049uaCnApHmAloGpqKhZ2leNFoJ34v5qsJOZwjO/0bCKgpEGvMIN+RjVlwE6exFiVsAo6ISZczl1C9ZVQM3OpxGDae78PeDWe8m3R3M5uw4nOrJm7Mp9iBWpLcoMRaA5PyzY2zjWaea+2xN00rcTiRuARYvXszMmTO9yl177bXs3r0bh8PBtm3buPzyyyvUboWSt3v27KF9+/YMHDiQzp07c/7555OYmFh0PDs7m7Fjx1YoAHH2mb89iWumreDnTcdxuNUzF6gBjsS9qLbCkUBBbQvnGTEER2Ju1A4A24H1PomrPMICTb4OQQhRD0ifLarDuqPHGPHV13y7dQuZZgfOcJX8hhr2UBVXCyeKoXBe2AKbiQKbEYdDT0Gmhfx9ISg5eo4rwby+dCjrDjcnu8BKijOYfzLbsTirrVc7OZn/zT2nKTiPBuDOMqG5FDS1cLEyt05Bi3eitbWBUeWvTd3Islm96ki3BfLr4t4UZAcQGpyPrpRfuZHWfAJNp4wAMmhopzxqacnQyG1iwh1swJoGlhQw5IIxG3Q7rWj/hpOXEEh2aiD5WdYyVzAONEqf7Zc05cybEBUgfbaoDvtzkxm9/D2+O7wChy6LkEA7FrMT0GgUmINOARUdS3PbsSi7HdsK4lib24yvk/uwNqspdtXE6tRmZDmspbYRoHNi1rkxKB4UNFya3ivBeqpYSw7BRjt/Z7ZgXn4D7OrJ78Y0j4Hf8qJwaN79ZJDeQaP/Rt2WdGPUoKhYdSf78AxHAMEBdgwBHuxOE3k2E8b/5sQ/UBDF3OSOLE5rzbyU9izLaEW+x1y0qOmpXB4didmRZX28or7ykz67QrcWnnjiCTp16sS6devIyspi/Pjx9O/fn8WLF9OkSZOailH4EU3TePHPnbhV3w5J8OScMsl5wMnpEfSBYQC4c1Kpq0b3aXrmk4QQZz3ps0V1eHnpUmyaC0+I5+Qtfw0cIRo6r1+RCh6PHo9Hj+rWo3doONvawaRxhBA+ShoARyE630FkO+8593IzA8hOO7maOaoOd5oFD6ALOm3kj1FDF2Un3WZi2sYLaB+WSKSlgLTsEPYciUNxFo6wKbAVvzA7waEaiAnPweYw4VEL35Sm11AAY46GNUXDGXIy86t34j3cwaNgOmLE1cwD+sLPo7RpHK5uKXPUCyHOTPpsUR0+3LsAm5pPhNWO4b8FwgpHmWrFFvlMcYeS4g7Foyok2wMJNjr/W3xMx56CWDoGHStxMbJ81fLffg29ouFSdSilDF1UFGhkyWJPXgzTk87ld1M2HQIycGoKJmMBJc01GqB3oqErNSEMYKDwhmuey8TB/AYEWZxYjG5ScoJRNe+yKjoyXYFF8QCk2QNxeAwEGpzoFZUsm5UjWeGcHyV9tqi7KpS8XbFiBX///TcNGjSgQYMGzJkzh3vvvZcBAwawaNEiAgMDz1yJOKslZNg4nF770ySUW2nPUNQR3ZuGc0vfZr4OQwhRD0ifLarK5nKxIfE4nlBP0fWVoqiFVz+m0p+cUTQNLdoNptP6VAOkGiw4FzYmMD4XNVglNzeAvGwrxS/gNJTAkh/ZVEwauDU8mp5tOY0hC5R0PXqPgmbUUDyQmxtAUkoYsdFZxcon24NRFAiy2snODwAN9PkQdMSDNQNsobri8ah4JXANNh1B+xWcER40Peg0cEVqXguctTFF8kjvgaV+TqL+0rSyfzLW8Z+Tog6SPltUhzVp+wg129HrTn4JKQoEGEqebxYKp0zQKcopCU8Nt6Jnny2GNgHJXue6ND0p7hCvfQZFRUVBX0ICV9PAoNMwK24cmpFjzjCOOcPQNI121kSaBGV5nW/zGDHoCn9feDQdBuXkbw09Ks2NBUTrnRzBw19ZEWzJjz0Zh17FanJS4DTj9Ogxl7DAqKqB26PgdhtIyQrC6Tk5rVO0MZgpvYeW+jmJ+stf+uwKJW9tNhsGw8kiiqIwbdo07rvvPs4//3xmzZpV7QEK/5JZUDfma9WHRBX9WS04OZG5578/G045XhcMaN2Asf2ac1H7aF+HIoSoJ6TPFlWVl2UDgwYKmEwuwsLzMVvcaBrYbCZy7WZUrfjcBIpTB6HuEmoELVDDYTOQfzAc54lpaxUoY4DNmelAC1ZxawroNdCroMGK3a0YELCLBoF5KAq4VB3HbaGnjMApvCjUOaFx2xQiuuWhJpg5viMWh/3kY5wKoHODetoMCHqnQthOHaZcDUXT8BjA1hg6d4/num7n8P/t3Xl8VPW9//HXObNPkpmQhGxsgoBssogKuC8oda1b61Lr3lXbKv211S7qtba217Z6e7XaW0W6WVtbtbVarUURVxQQFUUQ2SELEJLJJLOec35/DCSELCQQskzez8fj3Ouc9XumZD5zPvP9fr7njJ9wADclfZomLJNuppgtByqaTGLbyRaJ293Se9R93VvCcuPaIwibOBgGLG8Ywc50kJH+7fjNFHHHQ3U6TMppmULKxFcXLlfruJ9yXDgYxCw31q4esZZjEEt5eKVxNONS1QwL1gCZmvWFvihl/ih+I8Whnp3kueoBg6jlZ5IvSf6ua4zwxDgut5JF9cU8XD2VmnQmrnvdFo1JiMQCFOVGW42I2bE9RG1tLjhgJMBw2UwZVsJFYydz6dhpmO0VsJf+LUtidpeSt+PGjWPJkiWMHz++xfr77rsPgHPPPbf7WiZZ6ecvfNzbTQDAVzYGMxDCjkVoXPUaORNOJF2/g8TWjwAIjJreyy3MMICfXHg4Fx+l4VIi0jWK2XKgfnf3v/D6kyTzDQaXRDB3PRAaBgSDSTyeNNsjubTopZo2oMEFoXRboyEBSObbpHNbJn0dy8Gw9zzAgJgLgm30nLHaqE/mcQAbPE7Tw1rMa7B46whGl27D70qTctxYjoHPTJOw3eT4ksRrLE4dv5JwTixz0DA4dPpmlvxjErXbQk33YNpAEhzTwTEzyVxv1MHd6GCkHQwHPCn4yhHHcsXnju3U+yv92L5q5PWT+nnSdyhmy4G6f/Fiko0GBFtvS9pu4mkX/r16o1qOQTTtI1O9PsPGwHbANGBDYjAbEkUM9e/Eb6bbnZwzYbmJWR7CnswEY44DScdN3PYQt1z4zDT53nrStovKeAjLMQn5ktTaOdRGm3uV51oJTgxUM967qz6v47DTtrDdUfJdLb83GAYcn1fFYcF/8r8Vx/Fuw5CmUkjxlIdt9bnk+RJ43WlSKTe1O3OorwtiNpiYaQPDNrh06uHcccbs/X3Lpb/IkpjdpeTt+eefz5/+9Cc+//nPt9p23333Yds2Dz74YLc1TrJLbWOSVz7uG7VkDZeH/BOuoOb5+2hc/TpbHrwWK1aPk4xhBkKEZlzUK+0yDThhTKbXb1l+gG+eNpaiPF+vtEVE+jfFbDkQjuOw6Nl3GVxqUDfF3ZS43ZPHY+OzbPDYHFZUTWluBCttUpWbz9ZImHq/h9RevX2MBqNV4hYAV+aahgNmbhpXbgrMXTP4GrD7sdJxwEq21YPIaZG4hcyDXVGoAdM0SDqZoZGmAQF3GicN+OGI4RuaE7e778tncfipq3nlj0dmLrs7gZt08EWdpgdcA5g5tASf30soHOSqq46ntDR/H++sZAPD6bi3+AH1JJcBSTFbDtQzq1dRm8yjsLC+zSTrjkSQsJMg7Glkes4GRvi24zgGKxvKWLRjDJY70zsWDOKWh6A7BZCZnMx02i2NEDASFPjqMbGJ2n5cjsUIXw1+I822dC4fx4qZlF/ZtP/wnBo+2FmG7WoZy3NcCS4atJaJvkjTOtMwKHS5SDltl2oyDYOgCdeULObGj88jUpmL480kj+tTXur2LPHgwKzcoeQO8pPr8/KVY2YwtriozfNKdsmWmG04Ts9XeIhEIoTDYerq6giFQvs+QLLC9miCI+/8T283o4XoBy8ReetJUjs2Ybi9BEZMIf/EK/EUDOmxNnxqYgku0+TQwTlcOmM4ZeH2Z/eU/k+ff9Lf6N/swGTbNudM+C62ZVPzQxeeorbrzzZWBDny0HXkBFuWRUrZJnUpP5X1YdbuKAQMjCSYW72kA20kbwEs8A6O4cprOfTSSpmkEi4cy8R2DEzTIS8nhs+bJp12UdcQwEqZGJ6WX2k9rjSTh21t8yHWsg3q0z6mhTfhcbX9UPjqH48gWrOrxEI8Tf67NRz56al4CoOUFoU499TJDCsb1Oaxkh32/vzb/XrYvXdgBtqfFM+Oxdl046363JQep5g9cM36v19T1dDAsHEV5Oa0jMmOk+lcaOJwbdkrlHojLbbXpgM8WzuZ+nSAHYkcHAx8ZgqvaZHnSVDqjwAOLpwWk4EVuSMUuhuaXueZMYpcLcsVxGwPSxpHYGNiYxK3PSRtNyvrS7H3KCY/JljND0rewWu0jsnNP+S2VmcnSWFz8+JP88bGsQAYacjb4DBtxjCCpX4GB3O4+LDDGV+oEoTZLNtjdpd63gKsX7+eF154gWQyyYknnsikSZMORrskCxXl+pg6LJ/lm2p7uylNcieeTO7Ek3vseiUhHzXRJCnboTzs5+unjuGSo1USQUQODsVs2V+maTLj5PG88Z8PCL5rkTq17f1KQnWtErcAHtPGZ1qUhyKkYm42VhThXesm1UGuc3h5FYeOqMJtWmyL5bE2UkjKduPy2KSiXuyEG7cnxdAh2/G4mx/uBoUa2FKVT8L2tDif1221O8Szo1msdyuzPFRtT+DsiFNqmVz09U9x0XUn7vM4EZH9oZgtB+LUQw/l0ffeY8e2EHjq8LnSGIaD4xi4XDZgMCG4pVXiFiDfHeMQ/3bWJkpI2SbVsVy2xcJE4j68njRnDF+F12VhYWA5DgYOha56yjy1GDjYmKQdk0JXQ6u4GzBTjPVVsz6V6eUaMJLUG35y3Qki6eZOS0XuWJuJW8jUgHba6PVrOw4pMsfkR0OEo25S9WkKE14uOXUaXzx3Fi6znR+MRfqZLiVvX3rpJc4++2xisczwMrfbzbx587j88ssPSuMk+/zXuRO5/OHF1Mfbnsgkm/ncJi9/62QMAyKxNIU5Xkyzf9RXEZH+RzFbDtQ13z6Tj5ZvYMfLUVLHAXtV8UnF3eT569s93mXYpBwoG1THxnWlxMtsTAcc08BnpDli6DoOLaomnvIQsf0UFUWbji3yNzI0p5ZFFYdmEri+NFbCTVFRfYvELYBpOhQX1LNpe0GL9fGUB9s22iz5kHZMbBs27SxgVNH2VtuT1X7+79avEwx4qa9rJJQfxOVuf8IXGWCypH6e9B2K2XKgvj5jJq9t2MDGWptYYYyEO/ODZo430VTuoNxX2+7xRe4oaxMlhLxxquN55AWT5AaS1Cf8fNJQyGF51Zi7yhiVeGoZ7q1pOtaFRdBItPvDaNjVCJkqDBgG5JoJ9h4AviPto8F2k2O2zhMkHIeEY5O7RyLWdhzqncxJ6xqDfP/s71CcP4ja+hihHD8exWzZLUtidpd+hvjBD37AaaedxpYtW9ixYwdf+MIX+Pa3v32w2iZZaMqwfF646US+ceoYjh5ZgNfdP/5QusNfvjQLv8eFz+1icJ5PiVsROagUs+VADR05mAee+SbXXnkGxyw5FH99poeM40C8zkfNJ4OIJ7ztHr+7OqzXa0FeEorS2IPT+Aob+cKxL3L2xOWML9nKtKEbOHH4KorcLRPBed4EI/N27DpZ5lw5uS3r0+7m86XxGHuVW7BNqutyW7fLgVjSTUXVIN7ZPJxIvOVQulTS5MrxNxEKBXB7XAwqylPiVlpyOrGIdIFithyo4txc/vG5y/n+cacw2ZhC0MnEv7TVnPKJWu0PHU/YmX59ppH58RUyidYcb4KaRA7v1g5hc2OYmoSfIZ6dXWqbvVfayTQcvKQImnHchkWBp4Go5eO1xrbLGlRZDq/FBvFSYwE1lkPUTrHTSZDCxrINQq47KSsoxGWaFIZzlLiVlrIkZnep5+2KFSt4/fXXKSsrA+Duu+/m17/+NTt27KCwsPCgNFCyT2nYz02njeUmYEttA8f+ZGFvN+mg++mFhzNlWH5vN0NEBhDFbOkO4YIcPvulk/ksJ5NKW0y47xc4joGza0bnDVVFjBxW3aq3je00PwjWx3w4e+R4Z5WupSy3rtW1Bnsi7EwHsWh+6CoORFldV0I6tvsrq0F737LtlJnJ8Zq75hiLGWypHkyqzEtxSR1eb5pYzMuOrXnUbAmTKM6c55kPDmdkwQ4KcqI0JH18yn0qk448an/eLhko9vWw108eBKXvUMyW7pDn83H1EUdwNUcAMOGxH2H6E/jyI7hMh5UNJZyc/1GbMbvRyQyvSdmZEgi7uUyHlGWSsD1sjg1iZGAbrjZ62KZwkXZM3G2UPtiebv1D6vhwFWBkJpPa1adpRSJEwjaZGdxGgSvB9rSPf+8cyodWPpNyKwCoahxCiStKvitO0nHRuOnTnHnCpfv1fskAkSUxu0vJ20gkQlFR84x8wWCQQCBAXV2dgorsl6JcPy4DrH7yB7M/QgE3Fx+lurYi0rMUs6W7edwuwp4cdsbjTeuilo/lFeUcN+QTclxJYraXOstPNO0HDBwHNtS0LGcwNr+qzfObBuS6EtRZwaZ1KdskGfFipzMJ3fr6AOFwY6tj4wk3qfTumrcOJA28NZnXOz4pYMe6AnDA2NWD11eRIh0EKxcs28Wa7cWwvZi8iiRfuvnM/X2LZKDIkgdB6TsUs+VgGOqU88nOSrbaUJDXyJBQiq3pfErcdXh2JVlTjkl1OozfzJQg2BHPYddPoE1Slgu3aeMynaYfZlszqEkHCZlx/K7mkTA16SCV6XCLPXdXTDBpeSm/K80HiTD/t2UqCdtLcte1/A0xBo+OUuKrx8Gg0sqj0sojutHDNceqh7rsQ5bE7C5PWPb8888TDjf/8dm2zYIFC1ixYkXTunPPPbd7WidZz+d2cdbkcv7x7tbebsoBKwn5qIq0nLTFNOCXl0zrpRaJyECnmC3d7fyJE5i3dFnT66JghPNGvEuepzn+1SYDvJcYSl0swKZNRex0NRfLHZxbj9+Tavf8zl4PjB9vKiMV8TY93O3YnEfAn8Tra34wtCyDqur8PY4ycLxgux3MtIGBAXt1BAqVBCm+ZyU1pxQRHeMD2yL37R386OLz8bm7/PVYBposqZ8nfYtitnS3Cw6byE/eqCHa6CVaEcJd5BAr+IT1qcH4jUwsjjsewCDluNjSENr142sz2wHLMYkmfXhMi3XpQnamggzytP4hdV1yMGsaiyHtkOeOsymez5jQdsKe5h99HSfTS9eg7c/JkCeBz7SpT+8Ri315vPDtMsZclcehI3biTlpsfdbDzOHXkTu7g5lQRSBrYnaXv51eeeWVrdZ96UtfavpvwzCwLOvAWiVZL5FO8dePX2PBxnep8TeSWxwmWj2cvX/la2ImcYU2Y3hi2IkQdn05OH2rls3JhxUzc1Qhj7y+ns01jUwcEub6kw5lxij9Wi4ivUMxW7qDZds888ZK/rNkNdt3RMj3mtQGMtnQs8e1TNwC5HtjBK0Ebywbj5MyoCwFbhg7uIqRhTuotYOU0Hqis7RtUm9lEr22Ax9tHsKmbYMzG53M/3HybDbXhsn1J/C6LVKWSaQ2iJVq4yut6dDe94pJY4dy1f/O5PGf/YNNz3zCsPFDuPhbn2fGWdP3920SETkgitnSXZ7duJKn1q2guqaWQsdhx66JxjZbg6hL+gl748SdljXrl+wYRl0ygGuPST4dBxqSXnbH0pTtJmW7eXzrdC4pf5vQHknZqlSIzakC/B6LOsfPOzuHkeNJkW50U+6LEDCT5LriGLu625oddHf0uVr+yDuyoJRvfOdS/vSTJ3n7nXWUjirh/K+fySmXHneA75RI/9Gl5K1tt65fItJVKdviCy/cx1uVq5tXDgJv7mqsHWOw6g5psb/hq8U7bDHGHh/iduEakhtnQQdF13uS2zS4fOYIJg0Jc960Ib3dHBERxWzpNrf8+hkWLFvT9NrBIZyMUWBWU3ZObZvHTAxXsLDuSHCBsR08Qxs4pCAz+VhVKkR+MkaJN9K0fyLt5rH3ZuD52IsRslgbyKch0TLGGz4Lw5dJyEb3nGTM60Bqr0StA0aq/Z4Unz9qGrNGjWDWOUd2+n0Q2c1wMktH20W6QjFbusutbz/HHz5uHiFDAfi2xAm+FiH2+UE8uXEaF41YSq4n2bTL+3XlrIyW4TYdXLaNaTg4GKQtk5Ttwm04TXVpAaoSYf533SmcVf4+ue4EEStIbI9kcNgbx2+m8Lhs0rabiBUgYgWwMMk3MxOPtvfzqu1ATTKnxbqLhs9iSvlEppw0sTveIhlgsiVmm/vepfNs2+af//xnd55SstBz65a2TNzuYnoSmHmVrdZ7St9rkbgFML0NeAZ/dNDauC/DCwJN/10W9vO/l05j0pBwB0eIiPQtitnSGW+t3NgicQtgYGB4g+RtS7ZzFLix2T3vmCvqojgd3+PBz2BVrJRl9cNZGyti2Y7h/HThWVS8XkbtsiJqFpbgbPC1PKG965xtMAww3C2/ebsazKb6tsMHNcfngmCAH545m+NGjdj3zYu0p5tnrn7ggQeYPHkyoVCIUCjErFmz+Ne//tW0PR6Pc/3111NYWEhubi4XXnghVVVt146W7KSYLZ2xunZby8TtLqkhflxOCgyDiliYX606iSc2TOXfFeOYv24mC6rG4drVO9dyXKRsN2nbBYaBx7SxHKOpTu1u8bSHnVYuVen8Fonb3bxmJmbHLA/pXZOc1qaDpHfF5vY+KmPWOBK7RuHkuv18dcynOLNco2LkAHRzzO4t3VLUa82aNcybN4/58+ezbds2Uqn2a5mJPL++dUDZzQzUgJkEOxMADE8Dpj/S9r55FVA59WA0EYAcrwvHcQj63BTl+mhIphlRkMO1x43k5HHFrN/eQEMyzbjSEC6zf9RJERFRzJaueOHd1aQDmdBs7hqtO7pkB1efuIzJZVvBk2Y1ftakAi2O2/hxSYvXdmPrr5xR20806aemIkT4VS+eqIOBAy6Dgg8huNUhXmRgWOCJQPSI9tuZm/ZACoK4GewPEXPblA7L43NHTuGsCYexpS5CbWOMMcVFeF19q+ySyNChQ/nJT37CmDFjcByH3/72t3z605/mnXfeYeLEidx0000888wzPP7444TDYW644QYuuOACXnvttd5uuhxkitnSFc9tXNXGWgfD5VD3mRJccZtQbgIDWBUpxedNYRhgGnaLnrV7MgwH2zGxHSfTI9eBZMqN7Zhsa8xlSF5dq2PiaTf1cT/h3DhgsCORQ8gTx++CTbECBnssglg4to+CQAjD1UjAFWZS/qeYlH8GXxpTT02ynuHBwQTcrRPDIgPRfidvY7EYjz/+OA899BCvvfYaxx9/PLfeeivnn39+d7ZPskh9MsZ3Fs3nxU3vtbuPYdp4ij8gVTmVzECKjvq327Q/4GKvXTs+U5te/c7JDMrxtbv9kKKcdreJiPQlitnSVYl0mu8+9wL/2PARzhDAAU89TDa38/B1TxLcY8KwUpLkmRbvJHIBaKzzsfTf46G0+Xw7qsKkUy7cntb1GqNvDcITy4R1d10CM+QjaTr4txt4Iw62z0XaC94qF7FBrXvfGkl44uTLGDuqrN37GRIOMSQc2v83RGQPBvsYgrnr/0ciLTsg+Hw+fL7W3y3POeecFq9/9KMf8cADD/Dmm28ydOhQHn74YR599FFOOeUUAB555BHGjx/Pm2++ycyZMw/kVqQPUsyWrrJsmx++spA/rl4Gux5RDcPB5bYwXZmE69DcWr56yCtMDmcmCv+wvpT5W2ewLl6I08GETc3bDGzHIJU2sZ1MT9oP1w1hyOFxMBJ77A/vbR3KzkTOruQt2JjUpoKQyrTlvBFf5bTh7ZdAGOwPMdivmC3do7Mxu6/rctmEt99+my996UuUlpZy77338ulPfxrDMPjVr37Fl7/8ZUpKSvZ9EhmQfvDaHzpM3O7mCm/BO+IVXPnrMXz1OGlPm/sZhoO7+H06k5Y99ogIOeHWvwp2xDS6taqIiEiPU8yW/fXTha/w9w8+ao6wBqRCcOVpy1okbneb6I5R+XExb700gSceOIXGLTm4GzNH226HeNDg3ZWHkE41x1bHge1vFRFZlQ+WgyuWxkw7TBhRQmFDCu/WOgIrqwm++BG+DyrwbnHj3r5XbLah+HkLv9Etg8lEOmf3zNUdLcCwYcMIh8NNy1133bXPU1uWxWOPPUZDQwOzZs1i6dKlpFIpZs+e3bTPuHHjGD58OG+88cZBu0XpeYrZsr8eXPo2v33vHdKJTGw1DBuPL43LnalVm+tJcNeEfzQlbgEm5FVy2+h/MTa8jaAniWW3ncLaXfIAwLIM7F372Q0u6rbms/m9U0i8U070Xah5Hj6c62HnIylq6wLU1AVbnMtxDDZvLSTHlXcQ3gWRdnQyZvd1XfqmO3nyZCKRCJdddhmvv/46Eydmfi25+eabD0rjpH9zHIefvv03/rb6daKpOE4X+r6a/gimf8U+93MP2ogrVIEVLSG9fQyk2+gNa6R5u/41jFILt+dQ0tvH7/O8M0cVEA62nTQWEekPFLOlq377p9d5/Nll1MbiVB3haqpZu6dpQyvaPNbrstm6bATLPxmKYYGR5xBek2bH4SbJQsCAbclcXnp/AsPMWgatSRBdFyJZ68OTTGCkHfC5cYIePli1FQwXRiyB8eF6TCBYuRMzGAQrRLrIIV1gYSagYHGKw3IKGT5y8EF9b0T2x6ZNmwiFmnuPtdXrdrf333+fWbNmEY/Hyc3N5cknn2TChAksX74cr9dLfn5+i/1LSkqorGw9V4T0T4rZ0lV/+Hgp9334CjviDdiWgRMw8SRtvCmLVG7LCcZOKVzFIE+s1TlyXElmhtbxMmOpbsjBsk3MXZOT2Q4k0y7SKQ/JmItUygUGmDY4tW6MeOZZeenKCKY1Ev/afAY9+RGGnaaIrcSH5bFxUhE76uKEcmLYtkntzhxCZj5HlmqCb5Gu6lLydtWqVVx88cWcfPLJTJgw4WC1SbLEVc/fy+KK1hOTdTfDlcId3ozpryO58RiwWyZd3UWrMXYV6nMXfkK+PZLtNX7yfG7mTCrlhQ+rqIs1148qyPFy2zmayVJE+jfFbOmKH939T/71xkowDCw3OO2Uhd0eC1ISbGxz245opoeNY4BhGNgmJIoc9nyCtGwX6+1CGiriFH6UJEAM222QHhRodT6nrJBCw2bnBxvxBbycckgJ725LULmxEcNy8NRbBAIevv7Dsw/8DRDpin1NcLJr2+4JyDrjsMMOY/ny5dTV1fHXv/6VK6+8kpdffvmAmyr9g2K2dMXP3nuRX696PfPChIAnyQkjP2Z0eBsu06EmHuDNbSPZGC1ksC/Cp0o/bPdcRZ4oALmeJJXRPAx2J28NMtnaTJlCO5F5xrYB0zYxADPRXA8/PmoQObPH0vjvVbg9Lj6zs5BXAnmsjeykIRIAB3ymmx+dfppqz0vP6mTM7uu6lLxdu3Yt8+fP5ytf+QqxWIxLL72Uz33ucxjtVbeWAWdnPMrvP3yJtXWVPZK43ZPpq8c7/FWsnYdixwZhuOO48jfgymvZK+HyUz1cOuZUwgEPfo+LmoYkf1myiTXVUQ4dnMtnjxxKYW77PSNERPoDxWzZl2giwWOvvMP6TTt48e2VmLv+bRhpwMpMHLa3x1ZO5IfHtU4ovVNRyqpoAV6a64rFB4PjafvfW/0oN4XLkwDY3va/jh599alcfekscsJB/EEfDdEE/3lmOR9/VElpeZg55x7B4BLVxZMedhAeBL1eL6NHjwZg+vTpvP322/zP//wPF198Mclkktra2ha9b6uqqigtLW3nbNLfKGbLviSsJM9sfJM1O6v4w8YP2XN4zJkjVjAkp7lEYIE/xpyhK/n7+il8eui7pI32k6VVyUwM9bgsyKRucfb6DHP70ySizc/HjsfBvTOTvN3TkEuO5IePfJNArp9gXoBYKsXfP17J0qqtDA7k8NnxkzgkPGi/3wOR/TIQk7dDhgzhe9/7Ht/73vd48cUXmTdvHsceeyzpdJr58+dz3XXXMXbs2IPVVunDdsaj3PHmY/xr3dJebYfpa8As7biurss0KQn5m14X5Hj58omHHuymiYj0KMVsaU9DQ4JfPrGIP3y8Amv389wEF/7tDjkVDolB4EoZzdt2c+Cv705gZLCOy6e8h9eVmTzsg5pC7tl4OMVHV1G/MkTC48JIGBg7HdqbXiEZcrP1pDDByhTBzWlsE2yvgZlymnrxQCZmF5Y1P+jl5Pr49MUzuu/NENkPhrOPyU+64UHQtm0SiQTTp0/H4/GwYMECLrzwQiDTS3Pjxo3MmjXrwC8kfYJitrQnbqV4ZPWLPFv9r6ZuriMKoS7mp6IuRGkw0iJxu5vLcJhZ/AmD/VG2W0GG2HX4zJaThtZbPuptL4f6q/mgrvWknz5XipLcevK8CdJFLqprQ2yuKsCNQ2lOPQ2Gn8Z483O1aRgtYnbA4+GSCZO5ZMLk7no7RLqsJ2J2T9jv2R1OOeUUTjnlFOrq6vjjH//IvHnz+NnPfsakSZN47719T0ol2eO5dUv59qL5JO3WE5j0RZUNO4mnk/jd3t5uiohIj1DMlt2WLd/A9370BOvHOTjuPXp0GQbxwQbJkAM2GBaYMXC8u8ogOGAmM8+NP39jFo8sn8Kk0mpyR+9ku9+DqyxOWVmcknHbWVtbSG08iFEPrvoAltm6x4+70cAKmNSP9BEb7MH2GDguA8d2cFkOttcAGzbkJYjE4oQC/lbnEOk13dyL55ZbbuGMM85g+PDh1NfX8+ijj7Jw4UKef/55wuEw1157LXPnzqWgoIBQKMTXvvY1Zs2axcyZMw/kLqSPUsyW3T7YWcEXXnuUUO5GPLt+MHUZNn5XilAoRkmwnnxXlICRJOW4SO9VrD7flylzZGPyYbKEEZ6dDDIztW8jtpc6y8ux+WsAOC5vFfM2HsOGWDEAHtPi0IIduM1d1/WmGVZcw5HFaxmXW0nAm8KyDTbVFrEyVYLtuMj3GexM1jLIm98Tb49I5wzEnrdtCYfDfPWrX+WrX/0qy5cvZ968ed3RLuknauL1fPuV7knc1jz3LtF315OsrMNuiOMKBwkeVk7huUfiLe6+IZF/+/h1tsci/Pq067vtnCIi/YFi9sCWSKa5/a5/UONP47hbfwW0cUjngeMGcDIJ3AS40gZNaV4LMKEmFmStL0C5P9riHKYBI8I7qYsHcPIMAnYDDZEcHE/zA6UrBu765sRxOsfMnBfAY2B5d21zwQtr1lLx8BP86auXYpoaPix9RDc/CFZXV3PFFVdQUVFBOBxm8uTJPP/885x22mkA3HPPPZimyYUXXkgikWDOnDn86le/2u/mS/+gmD2wOY7DjYufoNGupXCPxG2uO4FhgIHN8QWrGeWvxjTAcaDe9lOdDuPsitrbE3kM9jTgM9MkHDerk4MxcPAZKUKuOOwRVnO8Ka4Z+jp3vn8mKb+XomC0KXG72xBvDYfnbml67TIdDinYhj+Z4P3oUKqNpfxoZRU/mnQ7PpfKEEofoeRta1OnTuWXv/xld55S+rh/r3+HpNU9PW53LnifdE0Ub2k+psdFans9kddX0/DBZkb++BJcge7rKfvy5hXM/2ABV008tdvOKSLSnyhmDzxvLVlHJBLDKWldysDBwQrQVELPnxfDDKRpiPswIi7c9S4Mx8AwwbAzE5oNKo20eR2PaZPnTRBJ+nHCJmW/XkMiL0D9rHIMw4srBpl0sIM7P4nhdkjv8OKkzRYPkrut2FLFPX9/mW+ef1J3vRUifcrDDz/c4Xa/38/999/P/fff30Mtkr5GMXvgebdmCxsaasj1NSdQ/a5U0xygk3M2MzpQDUDAsPEYDi4jhoNBdTqM7cDOdA4fNZQyOXdz03EOBh6j7ef3HH+Ko//zHu9VDyXnCgdyWm4/xL+91TEeUkwLbGBG8BN8pk3U8vFWlYfjy+848DdBRJp0KXl7yimn7HMfwzBYsGDBfjdI+pfGdLLbzpV/wnhCx4zFU5gHQPWfXmPnC+9j1TXS+OEW8qaP7LZrAfz3239jWvEopgzu3vOKiPQFitmyt3g8E7Pd9TY4JuwxEY7tAsflEAwkGTt+E3k5mVlIkmkXn2wrYntVGE+FB/wGpDJlFfae0GRPe25KHVJAwZ/W4uQPIjEy0xPHlZsi5/CduIKZLreOBbFP8khsym3zfL9fsIyjRw7n+KmjDuAdEOke2VI/T/oOxWzZW6OVAiCW9OA4mZDtMjKJXAOHQwNV+LCZ6E0yaFcP2ZQDa1JptiTzWd84GMtwsy0d5u16D4cEtuEz0jTaPrzeNAGz7QSu+zgfOVfVYEwNwpBM5ykTm5A7TsCVatrPxKbMXUeeK/N9Yfd3glxXHFK/p7p+IsV5Fx+U90akK7IlZncpebtw4UJGjBjBWWedhcfjOVhtkn7khKET+e+3/9Yt5yo8Z3qL14GxZex84X0ADE/LXkImYANe073fJRtsx+F3H77Iz0+8dr+OFxHpyxSzZW/Tpx2C222SNMEbgWQ4s97yOqRCNobbYeKE9fh9zXHV67YYV1rF8pSHep8fwzIw6114qlzUbgkRDG1rdZ20bRBN7houmQZ3KtN1J/ReNdtG5oPhkDulBtPf3JvIcEFwbD1Wo5v0jtb1bQ0LHn1uiZK30jc4RmbpaLtIFyhmy96mFQwl5PETd6LUJ32EfAkcx8Brpij01OMzLSZ7E4TM5syTx4Dx3hTvxdKs3eMH2ogVZFVjOTnuJAEzSZFdTz6Nra5pO7DGOwSoJPFUisAcL5gQcsdxmxYx29OUrC1xR5r+G/b4PdjJ/IBbEfmNkrfSN2RJzO5S8vanP/0pjzzyCI8//jif+9znuOaaa5g0adLBapv0A6Pzy/j8+JP5/cqXuvW8jm1T+/JKADyDQwTHD22x/aKxx3LVxNkU5+TzyIoXeODdf2F31AWoHevqqrqlvSIifY1ituytYFAOV3/+OO7/66v4a8AVd2gocUgNcsCAwnB9i8TtboYBZeE66hN+HLeD5UvjBGwqPi6kcEQtvmBzTxzHgbTjIteboD7px73Kz6zjx3HjHddSUJrP3z/8mHsXP98icbsn75DG1slb28Edgw2VO7v1/RDZb1lSP0/6DsVs2VvA7eHmyafxX+8+QWPKS9o28QcTjAhGMA3wOWaLxO2ejgjUsDxeTIPVvC5huxlkNJDnjrPdyqPQjuLfq/ft23WHMHjoJH688mZChXms4EMe+eS3uE0Lt+FQnQqTa1bjMmxCZrzD9sfT6w74PRDpFlkSs1sXPevAt771LT788EOeeuop6uvrOfbYYzn66KN58MEHiUTarnsm2e97Mz/L6SOmddv57ESKLfc9T+OKTbjCQYZ8/VOYnpYzZwY9fkbll5Lr8fO1aefwxcPn7Ne1RueXdUeTRUT6HMVsacvnPjuTS8/KjHTxNAJep6nOrM/b/kgWn6flNjvHxs61SbkMEpablG2StFzELQ+2Y1LkacSzOIh7nY/cbQmGjxtCbn4OnztmKrdceEK71zF9Nk5+CoqS4LHBAnc008SR5YUHevsi3WL3EMyOFpGuUMyWtlx0yFRum3IuAEnLTcCVZvfcnQ1WTrvHhcxUm+tzd/WUtXCxKlHG1lQ+9ZaPmnSQp6qm8FTFNIxKP8MOG0K4KMSxRTP5+mHXYJL5ITdqB1iXKCZpu/esvNRCZr1DwDN6P+9apHtlS8zuUvJ2t1mzZvGb3/yGiooKrr/+eubNm0d5ebkCywB266xLMNuaZaSL0nWNbPzJP2hYvgFPSZjht5yHb0hBq/2eWfs2abv5p8RB/rZr5HXEY7q5coImLBOR7KaYLXv70mXHE4xnvqmmg83r6xtblyto2hZvvc1dkMR0geWYpGw3acfVNMO135fGVe0muK6O9/+6jIZIrOm4IQ357V4n7QVjXCOMjuFMq8csi+NKZx4GP3/mkV28U5GDxOnEIrIfFLNlbxcdciTB2hAAeZ7m3q6rEwXt1p7fmg4Ss1qW33AZNm6zedSLhYvKdD4fJ8tYnypmyc5DsDb7ee/1GrbWNP97G8Rw9vxQq7cDfBgfSsxuu7xHpk0G5aEvd+1GRQ6WLInZ+5W83W3ZsmW8/PLLrFy5kkmTJqk+zwBWFAhxzaTZB3SOxJYaNtz5BIkN2wiMLWPE987HWxxqc99tsQjbGuuaXtclW9fs6YjHcHHd4acxsWj4AbVZRKS/UMyW3bweN18/9xiCmxpxNzR/Y400BKmJtO7Jk0y72FobbrU+HfXgWG3/cJuKuHE1JnHF08Qak2xe11wbd9uSGtKftP73Z9sGCbu5opdhgj0qQWq0xbnnTuXoiSO6dJ8iIv2VYrbs6bvjLyC5Io+GPX5k3Wn7eDcxqNW+CdvkxfohNFreFutTtknabjv9E0u5qV0bJl4ZxHIcVm9tjtmrNtdTvzPYYn8HWJ9sPRrGcSBme0l4riQcPLsrtygi+9Dl5O3WrVv58Y9/zNixY7nooosoKChg8eLFvPnmmwQCgYPRRuknvnnk+dx81IUY+9kDd8t9z5PeEQXAjifZfO+zbLjzCTbc+QS1i1a22NeF2aK3bb63/WEjbUk5Fg+8+y9+8tZf96utIiL9gWK2tOeCy2bx3a+dSdnSRIv1H6wdyobKQhIJN+mUybbteSzfNISktdc0CRa413upX5vX5vlrVxaQCvuom1xIZHIhg8vym7bl5Qfx/t4i+kmIVL2HdKOLxlofkbgPp41JI9LFCX7/3rt8+8nnDvi+RbrFvoZf9pNePNK3KGZLe+ZMHcvP51zEpuWjse3mOPl0/XAWRMvYnvbRYLl5p7GIX2ybwqp4IezxTO44EE35qEy03THq/S3DSbtNvMMbSY5IUl7Q/INtOODnw7dGtejl6zXSVFthVsTLqbP8JGwXO9MB3o+X83LDGH63dT3/u+ZebKft+vYiPSpLYnaXJiw788wzeemllzj99NO5++67Oeuss3C7u3QKyWKGYXDVpNm4TBc/WvyXLh/vpJvLICQ27mixLWfSsBavLWxuWvgQI0MlvLjpPQzAZZhYXQwQv/vwRS4ffxJD84q63F4Rkb5MMVv25fTzjsAf9PLth59k+wwf6RwT2zKoWjqYxk+KSeWaJPMN7HwLylOwu/x8CnybPLhSBtveKsaxDPJGRzDdDumYi50rCqhbH8b2AQZER+Xwrd8/w4Sjh/HCJ5+QtiysMXn4n01Sd3gBjsfEGhbHW5hou6GezLfqf77/EZ8/eipThqpevfSyLJn8RPoOxWzZl9mTxxDKu4Lv/ruRiYevJTc3juUY/LniML65cQRpn8HgcJQcT5LiQBT3rsnMbAdilgcbF5tjgzCAEl8Et2mTSLlZsXUoqbTJVTNeIceXJG2ZPFW5k0nODN6LvE3CSnDoOC8rl45gzOGb8fgsDBwMHKKWn4/s8qY2Og54jcwz/YeRD3i3djnTBh3RG2+XSLMsidmG47RXKaU10zQpKyujuLgYo70K1WSGeXQkEokQDoepq6sjFGr71x/p3x776BV+8/7zbInuwOfykLDaLpreEYOe+Tv64bGf4zNjj+uBK4no8096jmK2dNbLz7/PH36zkE+27cCTFySx67dU24T64S6S+SaO4WDn2OCA2WDiioMn3vzvynDb+N0WsYQb2zTYsxSejUOiGOyWIzjxb09R/PwOao4vpWGUhXdEQ5sToKSrfBgVmaGi3zj5GL5ywozufgtE2rT359/u16O+92Nc/vZrRFvxOGt/9F19bkqnKWZLZy1avY7/ffF11tVtwMr1sNO1+9+LQ7gkQl5eAnAIuDLP35ZjEHCncZsOjgPxtJua+hzMSJAdjW4OHVzN6eNWtLpONO1jR6p5dE28wcfyRYdy2MnrcLkspoU24m5jpqeUbfJB41AATh58Cp8bcUW3vwcibempmL1o0SLuvvtuli5dSkVFBU8++STnnXdeu/svXLiQk08+udX6iooKSktLO3Vv0MWet7feemuHwURkt0vGHc8l444naaVJ2Wl+tuRJnlrzJrF0slPHX3DoTK6aNJtz/37nQW4p/PaDF3l1y0ouGnssxw+ZcNCvJyLSExSzpbNOnHM4J845nFQqDRjM+90rPP2vd4k2JMj/JElkuJt4kQtX1AWOgzfiENxqk841sDwGJx45mm985kQunvsQTp7BHmVrsTyQDLdO3ALEizzsPL4YwwbbgFTcjTeQbrGPnTKguvngF55fweaFWzj9lImceNxhB+kdEenYvman7i8zV0vfoZgtnXXC2JGcMHYkybSF22Xyq2WLeWTFMnbEGqnflItZbhAMx4ntqnkba/SwdlsuPq+FYTpMLxzJL4+bw4V//ROWk2TqkA1tXifHlWBnKgd7V6VNf06CQ4/dTBI3YTOBq51uVh7TxmukSTpuXqpcw+uVv2dO+eGcPXQypnFAUy6J7JfujtkNDQ1MmTKFa665hgsuuKDTx61atapFcri4uLhL1+1S8vb222/v0slFvC43Xpeb22ZdyvdmfJaklcZtuvj3hnd4edMK3t22jo3121oc852jLuTqXZOfFQVCbI8d3NlV19RWsKa2gufXL2Pu9PP44uQ5B/V6IiI9QTFbusrjyXwt/NK1J3HdVSeQTKXxetwsfmU1r7z9Mcurqtm4roZdnXlwRRwuveAovnLViQCMzR/Ee1akqcyeY0AqBE4H3zbjYZNkkYWV42A1eLDTJm5fGsMAq9GFscWPae162HMcKlfWUBPbwWtvrOHCT0/na18+9WC9HSIiPUYxW7rK687UMrph+ky+esQMYukUQbeHVzdu4Pn1H7Eyuol3t20nncoE4UTczYVjJvKzE8/AMAymDy5nUdV6woFYm+c3DHCbFsk9JjkLBOJ4XWn8rhQOtDvTTZGnnqjl5/VtCRrTH/NK9ce8Wr2a/57+2e58C0R6xRlnnMEZZ5zR5eOKi4vJz8/f7+t2KXk7aNCgNn8RDIfDjB07lv/3//4fp5122n43RrKb23ThNjNB5uxRR3H2qKMAeLNiFQs3vY/f5eGsUUcxZlBz3ZzvzfgsNy18qMfaeN/yf/KZsce2mAxNRKQ/UsyWA+FymQRcmV47x5w0jmNOGgfAh6srePm1VWAYnHTsWMaPaa4/e9O1s/naj/9CXbkLTAPLzz5rIDkeBytn9w4G6YSbdCLz9dRsMPAldyduwd1gECs28WywMIAn/rGU886exrChBd178yIiPUwxWw6EaRjkeDIx+/gRh3D8iEMA+KhmG09/8hEJK83sEaOZWdY8j8wtx53I4sc3U9OQQ1m4rtU5HQfStqvFurRjEHJlRulELT8hd7zVcZZjkOtOkutOcmLJGl7YOo604+K5rSu4+JCjmV54SLfdt0h3ikRadhr0+Xz4fL5uO//UqVNJJBJMmjSJ22+/nWOPPbZLx3cpeXvvvfe2ub62tpalS5dy9tln89e//pVzzjmnS42QgW1m2WHMLGt76OMZI6eT6/Hz/Vd/T1WsdVDpbkkrzeKKVXxq5PSDfi0RkYNJMVsOhgljy5gwtu0Jw6ZOG8GvfnAJd85/npVOPc6uvKth7ep9u3dewgHDaj+zaxjgaszsZ6Yyh1sBk2TYxleXqd339tJ1St5Kz8uSyU+k71DMloNhXMFgxhUMbnPb+MGDefKiy/jN++CE/tWq5ny95W8qmQCZZG590k+BL9NTd3sqF7+ZwmtaLfZJ7jHcptDfyLhwFStqM52zXqteo+St9LxOxuxhw4a1WH3bbbd1y6iIsrIyHnzwQY488kgSiQQPPfQQJ510EosXL+aIIzo/oV+XkrdXXnllh9unTp3KXXfdpaAi3er4oRNZePFd3LTwIZ5b33GR/u6Q4w0c9GuIiBxsitnSGw6fPIw//+I6/vt3/+H3y9/FCpgYGJgpJzOJ2e6HQwuC1QAGbQ/YBDMOrjZK5SfzTHx1mYfFYLD7ekSIdJZq3kp3U8yW3jChuJh7Tr2RZ9cewps1vyWU00jaMalNBahP+zF3xWzbgY3RQURSPoY7OzOljXCxIVFIritOjpnEZ6axMNn7l9phOTubkre5HsVs6XmdjdmbNm1qUZO2u3rdHnbYYRx2WHNnxWOOOYZPPvmEe+65h9///vedPk+3Vow+++yz+eijj7rzlCIAGIbBvSd/gcfP+g4hz8FLrpYE85nVTi9gEZFsopgtB9O3r5jNv755NWXuIACGbWAmwExmetPmbQBPg4EnauDd0cYJbPDtcLWxgaYeEsGAl+OPGXNwbkBkX5wOFpFuppgtB9OZo87ja+P+QEXkSNY0FrM9FSLheEjYLhrSHt7bUcb2RC5J20NNIrjHkQZRK0BtOoCFi7aq4Fq7huG4DZMzh0zumRsS2VsnYnYoFGqxdGfJhL0dffTRrFmzpkvHdGvyNpFI4PW2MZ2wSDc5vPgQ/nXh7Xx+/MmU5xRQFhyEx+xSB/J2Ffhz+Z+Tv9BUl1dEJJspZsvBNnRIAc/ceBVfPOYohuWHKc/Lo2CHQe5GB9MCV9wmWGVR9grkrQUjnfkG7dvhUPyag6udLrm+iE1O0Mttt5xLTo568Ugv6OghUAlcOQgUs+VgG5ybw0+Oup3zyi+k1F9KvqeA+roS3q8YQsxq/re3tr6IisY80nYmUduQ8vJezRASVtvP0BuiBfhMNz+cej6lgXCP3ItIC30wZi9fvpyysrbLkLWne7Jeuzz88MNMnTq1O08p0kphIMT3Zn6W783MzFYZTcZ4Ys0brNi2gbRjAQZBt4/N0e0sqfyYtGN3eL4Cfx4/mHkxpwybjM/t6YE7EBHpfYrZ0hNyfT7+38nH8f9OPg6ARDzFi8++x4ql60kY4Iz14A942bG9nqV/X0fKtjHTDqRt3MUmDeUu9izEl590870vncbxs8aoZIKIDBiK2dITPKaHs8vP5ezycwFIWRbPr/iY19asxwrE8RUnCARc7EzEeHP7xyTtJJZlEt/h5z87fMwZ/yFus/nZu7J+MJ8dfg2nDzmcsDfY3mVF+pVoNNqi1+y6detYvnw5BQUFDB8+nFtuuYUtW7bwu9/9DsjUNB85ciQTJ04kHo/z0EMP8eKLL/Lvf/+7S9ftUvJ27ty5ba6vq6tj2bJlrF69mkWLFnWpASIHKtcb4IoJp7S5bWc8yrbGOuZ/sICnPnkT22n+WcXv8nDmyCO5/ZhL8bqUtBWR7KKYLX2Rz+/hjAumc8YFrScGjUbjbN8R5V9PLeMff1uCuTWJp95FosCN2+vixEmj+N6XP6XettLrVPNWuptitvRFHpeLs6eM4+wp41pta0gnqIrV8q8P1vLAqiVstnL5Q80gRpdXkuu1GJEzhVuO+ByDgzm90HKRZt0ds5csWcLJJ5/c9Hr35/eVV17J/PnzqaioYOPGjU3bk8kk3/zmN9myZQvBYJDJkyfzn//8p8U5OncfjtPpprZ38lAoxGGHHcZXvvIVRo4cuc/zRCIRwuEwdXV1LQoCixxMW6M1rNyxibLcAiYUDtv3ASIHgT7/pKcoZkt/trMmysoPtpI/KMiESUN7uzkyQO39+bf79Zhv/RiXz9/ucVYizsd3f1efm9JpitnSn9XF4yzZvIWQz8f0oUMwjda1b0UOtmyP2V3qefvSSy8drHaIHHTluQWU5xb0djNERHqEYrb0Z4MKcjnm+LG93QyRNqnnrXQ3xWzpz8J+P6eOPrS3myHSpmyJ2d1a81ZERERERCSr7WuCk37yICgiIpL1siRmK3krIiIiIiLSWVnyICgiIpL1siRmm73dABERERERERERERFpTT1vRUREREREOilb6ueJiIhku2yJ2UreioiIiIiIdFaWDMEUERHJelkSs5W8FRERERER6awseRAUERHJelkSs5W8FRERERER6aRsGYIpIiKS7bIlZit5KyIiIiIi0llZ0otHREQk62VJzDZ7uwEiIiIiIiIiIiIi0pp63oqIiIiIiHRStgzBFBERyXbZErOVvBUREREREemsLBmCKSIikvWyJGYreSsiIiIiItJZWfIgKCIikvWyJGYreSsiIiIiItJJxq6lo+0iIiLS+7IlZit5KyIiIiIi0llZ0otHREQk62VJzDZ7uwEiIiIiIiIiIiIi0pp63oqIiIiIiHRStsxcLSIiku2yJWYreSsiIiIiItJZWTIEU0REJOtlScxW8lZERERERKQr+snDnoiIyICXBTFbyVsREREREZFOypYhmCIiItkuW2K2krciIiIiIiKdlSVDMEVERLJelsRss7cbICIiIiIiMlDdddddHHXUUeTl5VFcXMx5553HqlWrWuwTj8e5/vrrKSwsJDc3lwsvvJCqqqpearGIiIj0JCVvRUREREREOmn3EMyOlq54+eWXuf7663nzzTd54YUXSKVSnH766TQ0NDTtc9NNN/H000/z+OOP8/LLL7N161YuuOCCbr4zERGR7NLdMbu3qGyCiIiIiIhIZ3XzEMznnnuuxev58+dTXFzM0qVLOeGEE6irq+Phhx/m0Ucf5ZRTTgHgkUceYfz48bz55pvMnDmzaxcUEREZKFQ2QUREREREZGDpbC+eSCTSYkkkEp06f11dHQAFBQUALF26lFQqxezZs5v2GTduHMOHD+eNN97o3psTERHJItnS81bJWxERERERkc5yOrEAw4YNIxwONy133XXXPk9t2zY33ngjxx57LJMmTQKgsrISr9dLfn5+i31LSkqorKzsrrsSERHJPp2M2X2dyiaIiIiIiIh0VieHYG7atIlQKNS02ufz7fPU119/PStWrODVV189sDaKiIhI1pRNUPJWRERERESkm4VCoRbJ23254YYb+Oc//8miRYsYOnRo0/rS0lKSySS1tbUtet9WVVVRWlranU0WERGRPkhlE0RERERERDqpu+vnOY7DDTfcwJNPPsmLL77IyJEjW2yfPn06Ho+HBQsWNK1btWoVGzduZNasWd1xSyIiIlkpW2requetiIiIiIhIZ3XzEMzrr7+eRx99lL///e/k5eU11bENh8MEAgHC4TDXXnstc+fOpaCggFAoxNe+9jVmzZrFzJkz9/s2REREsp7KJoiIiIiIiAwshuNgOO0/7XW0rS0PPPAAACeddFKL9Y888ghXXXUVAPfccw+maXLhhReSSCSYM2cOv/rVr7p0HRERkYGmu2N2b1HyVkREREREpLO6uReP04kHR7/fz/3338/999/ftZOLiIgMZOp5KyIiIiIiMrDsq0Zef6mfJyIiku2yJWZrwjIRERERERERERGRPkjJWxERERERkc5yOrGIiIhI7+vmmL1o0SLOOeccysvLMQyDp556ap/HLFy4kCOOOAKfz8fo0aOZP39+1y6KkrcywFj1SRreqSa2YjtOyu7t5oiIiEg7IpEYL/57BS8v+JDGhkRvN0ekye4hmB0tIiIDScxKsGjbcl6sWkok1dDbzRFp0t0xu6GhgSlTpnS6Bv26des466yzOPnkk1m+fDk33ngj1113Hc8//3yXrquatzJgRF7aROSF9bArZ2sG3RRcNg7/6EG92i4RERFp6dm/L+NX9/6bRCINQDDoZe53z+bEUyb0cstEyJrJT0REusNr29/jZx89SqMVB8BrevjSoedxdvmxvdwyEbo9Zp9xxhmcccYZnd7/wQcfZOTIkfz85z8HYPz48bz66qvcc889zJkzp9PnUc9bGRAallcTeX59U+IWwG5Ms+P3K7F3PRiKiIhI71vxznru+emzTYlbgMbGJD/9r7+zrTrSiy0TyVDPWxGRjC2127jzvUeaErcASTvFfR//lU+im3uxZSIZnY3ZkUikxZJIdM+orzfeeIPZs2e3WDdnzhzeeOONLp1HyVvJeo5lU/vkmra3JSyir27t4RaJiIhIe+68YT6Q6QhheU0sn4ltQszl8Ojf3sS2lRmTXqaatyIiAPzX/fdhu1p/6Dk4PLbhRdK21QutEtlDJ2P2sGHDCIfDTctdd93VLZevrKykpKSkxbqSkhIikQixWKzT51HZBMl60de34iTaDxqRFzYQfW0LBZePxz8qv+caJiIiIi289PQSqtIWTsiH7XOBaZD2Qe0Yk0QB/E/sfX73s9XcNed0Tp8ytrebKyIiMmAt/3ANq9w78LeTVvp35Xs8s3ktN447g8+OmNHDrRPpmk2bNhEKhZpe+3y+XmxNa0reSlar/dc6oi/ve7iG3Zhm+/+9j2dkiMBhBeTOKMMM6M9DRESkp/z5qbf51byFOCXNX5wtN1QfYZJuXsUOX4Ivv/Q0R79aznGjR3DpsVMpzA32QotlIFNpBBEZyJ7a8B63fPgPPOMKKCNT0sjEJuTJlE+oT/tJWSaG0cjv1v+eZbX/ZsqgcZxWcgpFvsLebLoMQJ2J2aFQqEXytruUlpZSVVXVYl1VVRWhUIhAINDp8yg7JVkrvmZnpxK3e0qti5BaF6FhaSUlX5mKGfQcpNaJiIjIbus2buf+eQtbrEuEDOqHGi0St7vZHlhWU8GK/1Ty5Nsf8McbLqE0P69nGiviOJmlo+0iIllqWzzK95Y8jYWDlfLQkPAwJLeO4Tk7ce3KklmOwXpzED63jWHAlvgGtlRs4OXqV/nBhO8wJFjey3chA0Yvx+xZs2bx7LPPtlj3wgsvMGvWrC6dRzVvJWtF36pstW5HYy0/eOFeZj7wGUbdfQqTf3kOlzx2IxtqW9a9tbbFiby4saeaKiIiMqD989l3WrxOeyFWZGIFjHaPsfyZL9tVdVF+8fiLB7V9InvShGUiMpA9tfodUk7zTOD1jT5G5NQ0JW4BXIbDqLwaPGbL8oUNVgMPL/9dj7VVpLtjdjQaZfny5SxfvhyAdevWsXz5cjZuzOSPbrnlFq644oqm/b/85S+zdu1avv3tb/PRRx/xq1/9ir/85S/cdNNNXbqukreSVRzLJl0TY9vD7xN/b3uLbTWNtZz7uy8xf9kTbGuoYWTBMAYHB7F0ywdURbe3Olf0ta3YsXSr9SIiInLgrLTF9q01/PDin/PEbxY0rU8FoHGwiWMYGB3Mc+KY4JiZb9zPv72SbZtbx3KRg0ITlonIAGNZFrXb6vjFFx/kwV882bS+wNfAtMGbcRt2q2MMA3LdcbxmCo+Zxtj14bgq9TGfvLu+p5ouA103x+wlS5Ywbdo0pk2bBsDcuXOZNm0at956KwAVFRVNiVyAkSNH8swzz/DCCy8wZcoUfv7zn/PQQw8xZ86cLl1XZRMkKzgpi9pn19G4tBon2faT3t2vPMTGugrGFo3k0Yt/TkluEQBJK4XTVld5B7betZi844aQe0w5rlzvwbwFERGRAcFKW/z2tj/zzwf/TX19AjMniGvMUJJeiBziIr27t60DZhIMCxzXXidxwPFArMzBVw0u4MoxX+Pcr36KC75xJsXDB/f0bckAYtiZpaPtIiLZwHEc/nL3P/jbPU+zs6oOgNQXJxJwJTl92ErKcyK79oOY7aHe9gPNo2bcpo1nVy9dj2ERtz2k0y6+euR3OPO6U7nwprMZOlYlFOTg6e6YfdJJJ7WdP9pl/vz5bR7zzjvvtN65C9TzVrJCzV9W0/BGRbuJW8dxePqjlwAozyvmsj9/k7G/OJ3T513Ns6texuduJzGbtKl/cRPV9y/Hqk8erOaLiIgMGA/cNJ8/3fUkEbcH6/CRWMOKMB2TyAizOXELYIDjMfBXGxipPU6we5gbgAmpQQ6eHQlSiTR/u+effPWom9mypqKH70pERCT7PPrjJ3jo5j9Qk45inxDEnuIjnZfDqcWrmxK3kOllG3SlCJotn5nTttliH5+ZovbDALZl889fv8BXj/wOq5Z80mP3I9JfKXkr/V56e4zYio6HSu5orKUuXg/AwnWLicSjhH15rNz2CV97+g6e+Whhh8dbOxPU/WdDdzVZRERkQIrU1PPMQwuIzTyU6OxxNI4vIjq5hB2TckgH2/5aatgGORvBSGd64Rp2y/pktg+82+M4rszxddsi/PqbqqcnB5HKJojIAJBKpnji3mewvpBP+g9DsH4wGOtnpQw7ZQvDCna2eUzAbP61NW2bpPYaOmMa0LDSj+3PDAKPReP88vrfHLybEMmSmK3krfR7qerGff7Bpe3mHrljCkfw2pcf47UvP8aYwhEAzF/2xD6vE1tWfUDtFBERGegqPqmiYdRg0qOKM11wAMsD6Zy96yI0c4xMAteApmVvnkAeHD0Jp7QQgMXPLsOyOiiYK3IANGGZiAwEO6vqqD3Cxv5sGNzN0TeYk2j3GAOH2qSfaNpD1PLSVtROluSz7brpRI8aAsDqtz+hsb6x29svAtkTs5W8lX7PyN136ebCYD5elweA8cWj8bo8eF0exhePBmBzXeU+z+GkbJJbowfWWBERkQGscMggUocWN722DbDdBmaKdn+IdSUcXJaBu50Q7Ntu44kBbheMHoaTF8S2bJb95/1ub78IkCnuuK9FRKSfyy8OY5+R22p9JO0nZbedSqqMhXh72wi2NOTTVuK2sdZPY10A3C4ajh1BfNQgAF594q3ubLpIsyyJ2UreSr9nVcf2uY/H5WbG0CkAfFT9CSkrTcpK81F1pr7OyEFDO3UtuyG1751ERESkTRtXbsXxeZpe7x5NadrgbqPTjZF2yNuQAsfBv93A2Kv8vKvRYfDS9B4HGLCr921tdV13N18EyJ5ePCIiHdm6phIn1DplZDkuVkeLW61P2wbv7iwnmXSx+v2h1EcCLbanEi7WvDGixbrYpBJAMVsOnmyJ2fvusijS17k69xvEt064jsWPvsvqHes55sGLAaiMbsNluLhh1uX7PN7wmHiH5h1QU0VERAY0x8FdHSE9rKDVJk89mGlIB8Exwd1gk786ibfRwVdrES904a01SQdtSKXwJtzkbrIx966O4PFgGDDpuHE9c08y8OyrRl4/eRAUEemI4ziY7yWwR7Se3Pvj+mJq4kEOC1UTdCfZkchhRW05NckcTMsgURtkyRvjKBpcS9gbId0YZNv6QViplikoO5j5QffwEyb0yD3JAJQlMVvJW+n3AocNYqfLAKvjv7pp5RN47JJ7ufuVh1hesRK/28fxI47kWydcx7TyfQeLvJOGYQb0JyMiIrK/Jh0/nsLNNVSVhsDjxrTA3tUR1wDcscwCkLMlhbcxE9s9DRa+mgSkbYxYAisaxRnTzqiZuihnf3kOZSNLDv4NiYiIZKlDJg5jyNt+Nh1nwaCWtekbkx7qk0WsbyhqdZyTzHSucjDYVjGI+JsxUocVNdW635N3Sz3HXTCD8TPGHJybEMkSykRJv2cGPQw6bzQ7n/h4n7+aHDX0cP5y6f+0v4NBq3MYXpP880aTc4QeAkVERA6E1+fh2z+/iju+8AANo0uwS/PBcreasMxXk8LbYDe9NhuSmKlMF1snmYKq7ZiDQ9j5oZbnd2yuu/kczv/qnIN+LzJw7WuYZX8Zgiki0hHDMPjWXV/iu1f/lMazfKROyCOd5yaW9BJPufH7UphtzCKarvE1/berIUn4lS3ELDcNE1uWWvAkLb509jF8/sZzDvatyACWLTFbyVvJCjlHleIdEaJxWRWJ9RGS6yNdO4EJhVdMxDcqTPSVLcTe34bjQGBiIXknDMX0609FRESkOxx73tHMnzKCf89fyMcrt/L26jrSAReJfA8Y4I7ZuHfXtnWclonbtIW9pYJvzfsKp3zuBJ762xL+88IKkok0M2aN5rOXzmTQoJzeuzkZGPY1wUk/mfxERGRfJp8wgd++ei/PP/ISqxdV8Kfjd2Ds+r01kXTj86abEriOA6lqP+naTPLWSNmUPPQ+X/jeZ7jkO+fxtzff58nFHxCJxZkxZjjXnHIU5QWhdq4s0k2yJGYrIyVZw1McJPypkQBU3buMVGVD5w40ofCy8QTGZervhU4dTujU4QermSIiIgNe2cgSrvyvTP35W67+DcuWb8Idz/S0dVJpSKbA7ca07KYeEY5tY1dU8Y37r+X0K08G4KKLZ3DRxTN65R5k4MqWXjwiIp1RWDaIy757AQD1L/6Tp6uXY3ptHMekMerFjJqYEQcrHcBJ7crsWg4FT63huhvP49KbzwfgolmTuWjW5N66DRmgsiVmK3krWan4+ilEXtpE49JqHMvGN3YQrjwv8Y9qsKMpMMDwu/AdEibvhKF4ioO93WQREZEB6c7fXMMTDy3k2T8vpjEa54gTxjC4wM/bz77DjgYLM+gj4HczdsxgLvraFxlzxKjebrIMdFky+YmISFf97OSzmPb+UH73wVJ2xBo5qng40+pM3n36LdYUpUkMyyXHNpmWDnLJL7/O4ceP7+0my0CXJTFbyVvJSobHRfj0QwiffkjLDWeM7JX2iIiISNtcbhef+fKpfObLp7ZYf83tn+2lFol0LFt68YiIdJVhGFw+eSqXT57acsMXVbdW+qZsidlmbzdARERERERERERERFpTz1sREREREZHOsp3M0tF2ERER6X1ZErOVvBUREREREemsLKmfJyIikvWyJGYreSsiIiIiItJJBvuon9djLREREZGOZEvMVvJWRERERESksxwns3S0XURERHpflsRsJW9FREREREQ6KVtmrhYREcl22RKzzd5ugIiIiIiIiIiIiIi0pp63IiIiIiIinZUlk5+IiIhkvSyJ2UreioiIiIiIdJLhOBgd1MjraJuIiIj0nGyJ2UreioiIiIiIdJa9a+lou4iIiPS+LInZSt6KiIiIiIh0Urb04hEREcl22RKzlbwVERERERHprCypnyciIpL1siRmm73dABERERERERERERFpTT1vRUREREREOstxMktH20VERKT3ZUnMVvJWRERERESkkwwns3S0XURERHpftsRslU0QERERERHprN29eDpaumDRokWcc845lJeXYxgGTz311F6Xc7j11lspKysjEAgwe/ZsPv744268IRERkSzVzTG7tyh5KyIiIiIi0kmGve+lKxoaGpgyZQr3339/m9v/+7//m1/+8pc8+OCDLF68mJycHObMmUM8Hu+GuxEREcle3R2ze4vKJoiIiIiIiPSSM844gzPOOKPNbY7jcO+99/L973+fT3/60wD87ne/o6SkhKeeeopLLrmkJ5sqIiIivUA9b0VERERERDqrk0MwI5FIiyWRSHT5UuvWraOyspLZs2c3rQuHw8yYMYM33nij225JREQkK6lsgoiIiIiIyADjdGIBhg0bRjgcblruuuuuLl+qsrISgJKSkhbrS0pKmraJiIhIOzoZs/s6lU0QERERERHpJMNxMDroqbN726ZNmwiFQk3rfT7fQW+biIiINOtszO7r1PNWRERERESkszo5BDMUCrVY9id5W1paCkBVVVWL9VVVVU3bREREpB0qmyAiIiIiIjLAOIDdwdKNz4EjR46ktLSUBQsWNK2LRCIsXryYWbNmdd+FREREstFBitn3338/hxxyCH6/nxkzZvDWW2+1u+/8+fMxDKPF4vf7u3Q9lU0QERERERHpJdFolDVr1jS9XrduHcuXL6egoIDhw4dz4403cueddzJmzBhGjhzJD37wA8rLyznvvPN6r9EiIiID1J///Gfmzp3Lgw8+yIwZM7j33nuZM2cOq1atori4uM1jQqEQq1atanptGEaXrqnkrYiIiIiISCd1d/28JUuWcPLJJze9njt3LgBXXnkl8+fP59vf/jYNDQ188YtfpLa2luOOO47nnnuuy712REREBpqDUfP2F7/4BV/4whe4+uqrAXjwwQd55plnmDdvHjfffHPb1zGMAyp3pOStiIiIiIhIZzl0XCOvi8+BJ510Ek5HD5aGwR133MEdd9zRtROLiIgMdJ2M2ZFIpMVqn8/XZq36ZDLJ0qVLueWWW5rWmabJ7NmzeeONN9q9TDQaZcSIEdi2zRFHHMGPf/xjJk6c2OnbUM1bERERERGRzsqSyU9ERESyXidj9rBhwwiHw03LXXfd1ebptm/fjmVZlJSUtFhfUlJCZWVlm8ccdthhzJs3j7///e/84Q9/wLZtjjnmGDZv3tzp21DPWxERERERkc6ygY5K1dk91RARERHpUCdj9qZNmwiFQk2r2+p1u79mzZrVYpLRY445hvHjx/PrX/+aH/7wh506h5K3IiIiIiIinXQw6ueJiIhI9+tszA6FQi2St+0pKirC5XJRVVXVYn1VVVWna9p6PB6mTZvWYrLSfVHZBBEREREREREREZEOeL1epk+fzoIFC5rW2bbNggULWvSu7YhlWbz//vuUlZV1+rrqeSsiIiIiItJZ+6prq563IiIifcNBiNlz587lyiuv5Mgjj+Too4/m3nvvpaGhgauvvhqAK664giFDhjTVzb3jjjuYOXMmo0ePpra2lrvvvpsNGzZw3XXXdfqaSt6KiIiIiIh0lpK3IiIi/cNBiNkXX3wx27Zt49Zbb6WyspKpU6fy3HPPNU1itnHjRkyzudDBzp07+cIXvkBlZSWDBg1i+vTpvP7660yYMKHT11TyVkREREREpLOUvBUREekfDlLMvuGGG7jhhhva3LZw4cIWr++55x7uueee/brObkreioiIiIiIdFYnZ64WERGRXpYlMVvJWxERERERkU7q7MzVIiIi0ruyJWab+95FRERERERERERERHqaet6KiIiIiIh0lmreioiI9A9ZErOVvBUREREREeks2wGjg4c9u388CIqIiGS9LInZSt6KiIiIiIh0Vpb04hEREcl6WRKzlbwVERERERHptH08CNI/HgRFRESyX3bEbCVvRUREREREOitLevGIiIhkvSyJ2WZvN0BEREREREREREREWlPPWxERERERkc6yHTocZtlPJj8RERHJelkSs5W8FRERERER6SzHziwdbRcREZHelyUxW8lbERERERGRzsqS+nkiIiJZL0titpK3IiIiIiIinZUlQzBFRESyXpbEbCVvRUREREREOitLevGIiIhkvSyJ2WZvN0BEREREREREREREWlPPWxERERERkc5y2Ecvnh5riYiIiHQkS2K2krciIiIiIiKdlSVDMEVERLJelsRslU3oJdFolNtuu41x48YRCAQoKy/niqu/yNqN1b3dNBEREelATeVOXnliMe8t+hCnn3zhE5FuZNv7XkSkT4jU1PPqk4tZ9p/3sCyrt5sjIj0tS2K2et72knPOOYeFCxficrkYNvIwqrZu5Pfzf8Pri9/igT88x3ETBxPwuXq7mSIiIrLLE/f9i9/f+STR7XVNX/SGji3nv576NsPHDenl1olIj8mSXjwi2eyZ55fwwI8eJ/pRJZ7tjQAUDS3ktr9+k3FHj+nl1olIj8mSmK2et73gww8/ZOHChQBcc+MP+cXvX+Tu+f8G4JOV7/Ls00/y9pq6XmyhiIiI7Pbkn1/huGO+xb2/e5EGw8QMhzBycgDYvHort5//3+qBKzKQ7H4Q7GgRkV6xaNEHTL/qLm5+/hU2HFnKjsunUvPp8dgek+2bd/D9c35CMp7s7WaKSE/Jkpit5G0Pi8bSvPRec2kEw8z8T2Aazf9TvPf2InbUp6iPpXu8fSIiIpKxs7qOT199N7e9uYTIjHJqjy2n6oLRJIbkYvq8GH4/AJtWbeX9V1b2cmtFREQGrsb6GNdcdg/f+ONzJAv8LbYlRw6i/vhDAKjbFuGNfyzphRaKiOw/JW97kG07vPxBDSXDxjJ81DgAHvrF95h7xan8v6tOa9qvZnslAMl0/6i9ISIiko2+cvU9rMtvWWHK8brYefwQbLeB4fM1ra/bFunp5olIb7GdfS8i0qNuvuKXrLRt0kGjze3x8YNxzMy2WsVskYEjS2K2krc9qGJngmTaweVy8f17HuWEORcSChdQtWUD46fMoHTIIQC4XG5cJoSDnt5tsIiIyAC14rWPWONue2ITx+siPiwPdj0Emm6TCccc1pPNE5Fe5Dj2PhcR6Tlb1lSwfNMOHAMw2k7eOh4XjjuT/ph84oQebJ2I9KZsidmasKwHNSSaHwKLisu58fb7m14nE3GuOXsyAOUjDmXCsFzcrrYDj4iIiBxcFZ9U4XjanzjUcZuQypQ3uuAbZ1FYNqinmiYivc3ZR0+dflI/TyRbVKytxnG7MGwwUw62p/VztLsqipm0OOniYxg5aXgvtFJEekWWxGz1vO1BoUBzrvyTVe8Ra4gCYFkWv73vDhqjmeEbV3/+MkaX5fRKG0VERARGHj4c34adbW+0HbxbopCIc/Wdl/Klu6/o2caJSO/KkslPRLLFiAlD8VTXYQDeSBt/g5ZN3usbuPCms7nlj9/olTaKSC/Jkpitnrc9qCTfSyjgIhKzePHpP7Hg6T9ROvQQamuqidTWAHDjjTdy/hkn9HJLRUREBrbR00ZyXEkJ/9nWQHpwyx9Uc96rZsroEn6+4X6MdoZnikgWs20wOhhm2U+GYIpki8FDC/nUiRP5x/oqPMVhzO0OyVyw3eCqjTO5Ms5DHz+I2630h8iAkyUxW59ePcgwDI6bUMC76yOMmTiNFcteo2rrBgxg+vTpfOUrX+Haa6/t7WaKiIgI8F9/nkvZLX/k74s/JFqaQygnwKdnTOQb99zU200TERGRPcz9vy8x6LY/87dnltEQDpCf6+eUo8fw7V9+Vj+0iki/p+RtD/N5TI4ek89Rt9/APbfdgGkqkIiIiPRFgRw/N/3yWm50HNKpNB6vJhIVEXYNsez/9fNEsonb4+YLP/4c1/3oMsVsEWmWJTFbydteYhhGexNhioiISB9iGIYeAkWkiWPbOB0MwewvM1eLZCPFbBHZU7bEbCVvRUREREREOitLevGIiIhkvSyJ2UreioiIiIiIdJbtgNH/HwRFRESyXpbEbCVvRUREREREOstxgI5mru4fD4IiIiJZL0tittnbDRARERERERERERGR1tTzVkREREREpJMc28HpYAim00968YiIiGS7bInZ6nkrIiIiIiLSWY6972U/3H///RxyyCH4/X5mzJjBW2+91c0NFxERGWD6SMx+/PHHGTduHH6/n8MPP5xnn322S9dT8lZERERERKSTHNvZ59JVf/7zn5k7dy633XYby5YtY8qUKcyZM4fq6uqDcAciIiIDQ1+I2a+//jqXXnop1157Le+88w7nnXce5513HitWrOj0NQ2nF/oI19XVkZ+fz6ZNmwiFQj19eRGRXhOJRBg2bBi1tbWEw+Hebo7IPilmi8hAtXfMjkQihMNhjuNM3HjaPS5Nild5ttXnps/nw+fztXnMjBkzOOqoo7jvvvsAsG2bYcOG8bWvfY2bb765e29MspZitogMVH05Zl988cU0NDTwz3/+s2ndzJkzmTp1Kg8++GCn7q9Xat7W19cDMGzYsN64vIhIr6uvr1fyVvoFxWwRGeh2x2yv10tpaSmvVu57qGNubm6rz83bbruN22+/vdW+yWSSpUuXcssttzStM02T2bNn88Ybbxxw+2XgUMwWkYGuL8bsN954g7lz57ZYN2fOHJ566ql939AuvZK8LS8vZ9OmTeTl5WEYRm80QUSkVziOQ319PeXl5b3dFJFOUcwWkYFq75jt9/tZt24dyWSyU8fu/ZnZXg+e7du3Y1kWJSUlLdaXlJTw0Ucf7WfrZSBSzBaRgaovx+zKyso296+srNxn23brleStaZoMHTq0Ny4tItLr1ONW+hPFbBEZyPaO2X6/H7/f30utEemYYraIDGTZHLM1YZmIiIiIiEgvKSoqwuVyUVVV1WJ9VVUVpaWlvdQqERER2dv+xOzS0tIDjvFK3oqIiIiIiPQSr9fL9OnTWbBgQdM627ZZsGABs2bN6sWWiYiIyJ72J2bPmjWrxf4AL7zwQpdifK+UTRAREREREZGMuXPncuWVV3LkkUdy9NFHc++999LQ0MDVV1/d200TERGRPewrZl9xxRUMGTKEu+66C4BvfOMbnHjiifz85z/nrLPO4rHHHmPJkiX83//9X6evqeStiIiIiIhIL7r44ovZtm0bt956K5WVlUydOpXnnnuu1QQnIiIi0rv2FbM3btyIaTYXOjjmmGN49NFH+f73v893v/tdxowZw1NPPcWkSZM6fU3DcRyn2+9ERERERERERERERA6Iat5Kn3XSSSdx4403tlo/f/588vPzAbj99tsxDINPfepTrfa7++67MQyDk046qdW2zZs34/V62/2lwzCMpiUcDnPsscfy4osvNm1ftGgR55xzDuXl5RiGwVNPPbU/tygiIpIVFLNFRET6B8Vskf5HyVvp98rKynjppZfYvHlzi/Xz5s1j+PDhbR4zf/58PvvZzxKJRFi8eHGb+zzyyCNUVFTw2muvUVRUxNlnn83atWsBaGhoYMqUKdx///3dezMiIiJZTDFbRESkf1DMFuk7lLyVfq+4uJjTTz+d3/72t03rXn/9dbZv385ZZ53Van/HcXjkkUf4/Oc/z2WXXcbDDz/c5nnz8/MpLS1l0qRJPPDAA8RiMV544QUAzjjjDO68807OP//8g3NTIiIiWUgxW0REpH9QzBbpO5S8laxwzTXXMH/+/KbX8+bN43Of+xxer7fVvi+99BKNjY3Mnj2byy+/nMcee4yGhoYOzx8IBABIJpPd2m4REZGBRjFbRESkf1DMFukblLyVrHD22WcTiURYtGgRDQ0N/OUvf+Gaa65pc9+HH36YSy65BJfLxaRJkxg1ahSPP/54u+dubGzk+9//Pi6XixNPPPFg3YKIiMiAoJgtIiLSPyhmi/QN7t5ugEh38Hg8XH755TzyyCOsXbuWsWPHMnny5Fb71dbW8sQTT/Dqq682rbv88st5+OGHueqqq1rse+mll+JyuYjFYgwePJiHH364zXOKiIhI5ylmi4iI9A+K2SJ9g5K30meFQiHq6upara+trSUcDrdaf8011zBjxgxWrFjR7q+Bjz76KPF4nBkzZjStcxwH27ZZvXo1Y8eObVp/zz33MHv2bMLhMIMHD+6GOxIREclOitkiIiL9g2K2SP+jsgnSZx122GEsW7as1fply5a1+PDfbeLEiUycOJEVK1Zw2WWXtXnOhx9+mG9+85ssX768aXn33Xc5/vjjmTdvXot9S0tLGT16tAKKiIjIPihmi4iI9A+K2SL9j3reSp/1la98hfvuu4+vf/3rXHfddfh8Pp555hn+9Kc/8fTTT7d5zIsvvkgqlSI/P7/VtuXLl7Ns2TL++Mc/Mm7cuBbbLr30Uu644w7uvPNO3O59/1lEo1HWrFnT9HrdunUsX76cgoIChg8f3rUbFRER6ecUs0VERPoHxWyR/kc9b6XPGjVqFIsWLeKjjz5i9uzZzJgxg7/85S88/vjjfOpTn2rzmJycnDYDCmR+DZwwYUKrgAJw/vnnU11dzbPPPtupti1ZsoRp06Yxbdo0AObOncu0adO49dZbO3dzIiIiWUQxW0REpH9QzBbpfwzHcZzeboSIiIiIiIiIiIiItKSetyIiIiIiIiIiIiJ9kJK3IiIiIiIiIiIiIn2QkrciIiIiIiIiIiIifZCStyIiIiIiIiIiIiJ9kJK3IiIiIiIiIiIiIn2QkrciIiIiIiIiIiIifZCStyIiIiIiIiIiIiJ9kJK3IiIiIiIiIiIiIn2QkrciIiIiIiIiIiIifZCStyIiIiIiIiIiIiJ9kJK3IiIiIiIiIiIiIn3Q/we3DXKGp48q+AAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 1680x400 with 5 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.rcParams[\"figure.figsize\"] = (4, 4)\n",
    "sc.pl.umap(adata, color=[\"leiden\", \"query_expression\", \"CD3E\"], wspace=0.4, ncols=3, legend_loc = \"on data\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "malva",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}