BY_COVID_update_August.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "faa7efb5",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pybel\n",
    "from pybel.io.jupyter import to_jupyter\n",
    "from utils import *\n",
    "import seaborn as sns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "6b327410",
   "metadata": {},
   "outputs": [],
   "source": [
    "itmp_chem = pd.read_csv('C:\\\\Users\\\\reagon.karki\\\\Documents\\\\ITMP\\\\EOSC Future\\\\ITMP ChEMBL data\\\\CHEMBL4495564.csv',sep=';',usecols=['ChEMBL ID'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "cfc7793f",
   "metadata": {},
   "outputs": [],
   "source": [
    "itmp_chem = list(itmp_chem['ChEMBL ID'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "de9ca28f",
   "metadata": {},
   "outputs": [],
   "source": [
    "len(itmp_chem)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "6b0c3d74",
   "metadata": {},
   "outputs": [],
   "source": [
    "itmp_dis = RetDrugInd(itmp_chem)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "8fdab82a",
   "metadata": {},
   "outputs": [],
   "source": [
    "# filename = ('data/normalized_data/bycovid_Aug.pkl')\n",
    "# outfile = open(filename,'wb')\n",
    "# pickle.dump(itmp_dis,outfile)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9f113052",
   "metadata": {},
   "outputs": [],
   "source": [
    "itmp_dis = RetMech(itmp_chem)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "448f3551",
   "metadata": {},
   "outputs": [],
   "source": [
    "# filename = ('data/normalized_data/bycovid_itmpMech_Aug.pkl')\n",
    "# outfile = open(filename,'wb')\n",
    "# pickle.dump(itmp_dis,outfile)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e3677996",
   "metadata": {},
   "outputs": [],
   "source": [
    "itmp_act = RetAct(itmp_chem)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "22838552",
   "metadata": {},
   "outputs": [],
   "source": [
    "len(itmp_act)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "11f4b7f1",
   "metadata": {},
   "outputs": [],
   "source": [
    "# filename = ('data/normalized_data/bycovid_itmpAct_Aug.pkl')\n",
    "# outfile = open(filename,'wb')\n",
    "# pickle.dump(itmp_act,outfile)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "4835429b",
   "metadata": {},
   "outputs": [],
   "source": [
    "infile = open('data/normalized_data/bycovid_itmpMech_Aug.pkl','rb')\n",
    "itmp_mech = pickle.load(infile)\n",
    "infile.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "24c359dd",
   "metadata": {},
   "outputs": [],
   "source": [
    "infile = open('data/normalized_data/bycovid_Aug.pkl','rb')\n",
    "itmp_dis = pickle.load(infile)\n",
    "infile.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "e1efd720",
   "metadata": {},
   "outputs": [],
   "source": [
    "infile = open('data/normalized_data/bycovid_itmpAct_Aug.pkl','rb')\n",
    "itmp_act = pickle.load(infile)\n",
    "infile.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "4a12ad1f",
   "metadata": {},
   "outputs": [],
   "source": [
    "prot_actMech = Ret_chembl_protein(itmp_act) + Ret_chembl_protein(itmp_mech)\n",
    "prot_actMech = set(prot_actMech)\n",
    "prot_actMech = list(prot_actMech)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "b1117409",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Filtering UniProt proteins from ChEMBL: 100%|██████████| 1470/1470 [00:04<00:00, 314.25it/s]\n",
      "Filtering human proteins from ChEMBL: 100%|██████████| 1436/1436 [00:04<00:00, 322.61it/s]\n",
      "Populating ChEMBL data for human proteins: 100%|██████████| 1419/1419 [00:04<00:00, 323.49it/s]\n"
     ]
    }
   ],
   "source": [
    "actMech2uprot = chembl2uniprot(prot_actMech)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "58ab1a12",
   "metadata": {},
   "outputs": [],
   "source": [
    "uprot = Ret_uprotid(actMech2uprot)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "47c8a207",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "reading uprot ids: 100%|██████████| 1419/1419 [05:09<00:00,  4.59it/s]\n"
     ]
    }
   ],
   "source": [
    "uprot_ext = ExtractFromUniProt(uprot)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "4234aa92",
   "metadata": {},
   "outputs": [],
   "source": [
    "# filename = ('data/normalized_data/bycovid_uprotExt_Aug.pkl')\n",
    "# outfile = open(filename,'wb')\n",
    "# pickle.dump(uprot_ext,outfile)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "bc2ec5ed",
   "metadata": {},
   "outputs": [],
   "source": [
    "infile = open('data/normalized_data/bycovid_uprotExt_Aug.pkl','rb')\n",
    "uprot_ext = pickle.load(infile)\n",
    "infile.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "7e70a959",
   "metadata": {},
   "outputs": [],
   "source": [
    "itmp_act"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "946b4744",
   "metadata": {},
   "outputs": [],
   "source": [
    "itmp_act = chembl2gene2path(actMech2uprot,itmp_act)\n",
    "itmp_mech = chembl2gene2path(actMech2uprot,itmp_mech)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "2b9cb155",
   "metadata": {},
   "outputs": [],
   "source": [
    "itmp_graph = pybel.BELGraph(name='ITMP Graph', version=\"0.0.1\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "f78da6ea",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Populating Chemical-MoA edges: 100%|██████████| 1787/1787 [00:00<00:00, 2317.76it/s]\n",
      "Populating Drug-Indication edges: 100%|██████████| 2451/2451 [00:04<00:00, 581.19it/s]\n",
      "Adding bioassay edges to BEL: 100%|██████████| 2494/2494 [00:09<00:00, 274.65it/s]\n"
     ]
    }
   ],
   "source": [
    "itmp_graph = chem2moa_rel(itmp_mech,'HGNC',itmp_graph)\n",
    "itmp_graph = chem2dis_rel(itmp_dis,itmp_graph)\n",
    "itmp_graph = chem2act_rel(itmp_act,'HGNC',itmp_graph)\n",
    "itmp_graph = gene2path_rel(actMech2uprot,'HGNC',itmp_graph)\n",
    "itmp_graph = uniprot_rel(uprot_ext,'HGNC',itmp_graph)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "id": "0cfd504a",
   "metadata": {},
   "outputs": [],
   "source": [
    "#Extract chembl ids from graph\n",
    "chembls = getNodeList('ChEMBL',itmp_graph)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "id": "7d99982e",
   "metadata": {},
   "outputs": [],
   "source": [
    "#add relation: chembl to side effects from opentargets platform\n",
    "itmp_graph = chembl2rxn_rel(ChEMBLids,itmp_graph)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "id": "cfd91368",
   "metadata": {},
   "outputs": [],
   "source": [
    "# filename = ('data/normalized_data/bycovid_itmpGraph.pkl')\n",
    "# outfile = open(filename,'wb')\n",
    "# pickle.dump(itmp_graph,outfile)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "a1dddfb2",
   "metadata": {},
   "outputs": [],
   "source": [
    "infile = open('data/normalized_data/bycovid_itmpGraph.pkl','rb')\n",
    "itmp_graph = pickle.load(infile)\n",
    "infile.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "id": "68ae19ce",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "            <h2>Metadata</h2>\n",
       "            <table>\n",
       "<tbody>\n",
       "<tr><td>Name   </td><td>ITMP Graph</td></tr>\n",
       "<tr><td>Version</td><td>0.0.1     </td></tr>\n",
       "</tbody>\n",
       "</table>\n",
       "            <h2>Statistics</h2>\n",
       "            <table>\n",
       "<tbody>\n",
       "<tr><td>Nodes          </td><td style=\"text-align: right;\"> 35952       </td></tr>\n",
       "<tr><td>Namespaces     </td><td style=\"text-align: right;\">     9       </td></tr>\n",
       "<tr><td>Edges          </td><td style=\"text-align: right;\">279462       </td></tr>\n",
       "<tr><td>Annotations    </td><td style=\"text-align: right;\">     0       </td></tr>\n",
       "<tr><td>Citations      </td><td style=\"text-align: right;\">     3       </td></tr>\n",
       "<tr><td>Authors        </td><td style=\"text-align: right;\">     0       </td></tr>\n",
       "<tr><td>Components     </td><td style=\"text-align: right;\">     1       </td></tr>\n",
       "<tr><td>Warnings       </td><td style=\"text-align: right;\">     0       </td></tr>\n",
       "<tr><td>Network Density</td><td style=\"text-align: right;\">     0.000216</td></tr>\n",
       "</tbody>\n",
       "</table>\n",
       "            <h2>Nodes</h2>\n",
       "            <table>\n",
       "<thead>\n",
       "<tr><th>Type             </th><th style=\"text-align: right;\">  Count</th><th>Example                       </th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "<tr><td>Abundance        </td><td style=\"text-align: right;\">  19931</td><td>a(ChEMBLAssay:CHEMBL4341379)  </td></tr>\n",
       "<tr><td>Pathology        </td><td style=\"text-align: right;\">   9005</td><td>path(SideEffect:paraesthesia) </td></tr>\n",
       "<tr><td>BiologicalProcess</td><td style=\"text-align: right;\">   5669</td><td>bp(GOBP:&quot;apical constriction&quot;)</td></tr>\n",
       "<tr><td>Protein          </td><td style=\"text-align: right;\">   1347</td><td>p(HGNC:CCR9)                  </td></tr>\n",
       "</tbody>\n",
       "</table>\n",
       "            <h2>Namespaces</h2>\n",
       "            <table>\n",
       "<thead>\n",
       "<tr><th>Prefix     </th><th>Name                            </th><th style=\"text-align: right;\">  Count</th><th>Example                                                      </th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "<tr><td>ChEMBLAssay</td><td>                                </td><td style=\"text-align: right;\">  15835</td><td>a(ChEMBLAssay:CHEMBL960125)                                  </td></tr>\n",
       "<tr><td>SideEffect </td><td>                                </td><td style=\"text-align: right;\">   7420</td><td>path(SideEffect:&quot;soft tissue haemorrhage&quot;)                   </td></tr>\n",
       "<tr><td>ChEMBL     </td><td>ChEMBL                          </td><td style=\"text-align: right;\">   4096</td><td>a(ChEMBL:CHEMBL254328)                                       </td></tr>\n",
       "<tr><td>GOBP       </td><td>Gene Ontology                   </td><td style=\"text-align: right;\">   2085</td><td>bp(GOBP:&quot;&#x27;de novo&#x27; NAD biosynthetic process from tryptophan&quot;)</td></tr>\n",
       "<tr><td>Disease    </td><td>                                </td><td style=\"text-align: right;\">   1585</td><td>path(Disease:Vulvovaginitis)                                 </td></tr>\n",
       "<tr><td>Reactome   </td><td>Reactome                        </td><td style=\"text-align: right;\">   1513</td><td>bp(Reactome:&quot;DCC mediated attractive signaling. [P06241-1]&quot;) </td></tr>\n",
       "<tr><td>HGNC       </td><td>HUGO Gene Nomenclature Committee</td><td style=\"text-align: right;\">   1347</td><td>p(HGNC:RAB2A)                                                </td></tr>\n",
       "<tr><td>GOMF       </td><td>Gene Ontology                   </td><td style=\"text-align: right;\">   1332</td><td>bp(GOMF:&quot;potassium channel activity&quot;)                        </td></tr>\n",
       "<tr><td>MOA        </td><td>                                </td><td style=\"text-align: right;\">    739</td><td>bp(MOA:&quot;Phosphodiesterase 10A inhibitor&quot;)                    </td></tr>\n",
       "</tbody>\n",
       "</table>\n",
       "            <h2>Edges</h2>\n",
       "            <table>\n",
       "<thead>\n",
       "<tr><th>Edge Type                              </th><th style=\"text-align: right;\">  Count</th><th>Example                                                                                                                    </th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "<tr><td>Pathology association Abundance        </td><td style=\"text-align: right;\">  73309</td><td>path(SideEffect:&quot;pulmonary arterial hypertension&quot;) association a(ChEMBL:CHEMBL779)                                         </td></tr>\n",
       "<tr><td>Protein association Abundance          </td><td style=\"text-align: right;\">  18083</td><td>p(HGNC:ERBB2) association a(ChEMBL:CHEMBL1201179)                                                                          </td></tr>\n",
       "<tr><td>Protein association BiologicalProcess  </td><td style=\"text-align: right;\">  17838</td><td>p(HGNC:VEGFA) association bp(Reactome:&quot;TFAP2 (AP-2) family regulates transcription of growth factors and their receptors.&quot;)</td></tr>\n",
       "<tr><td>BiologicalProcess association Abundance</td><td style=\"text-align: right;\">   2368</td><td>bp(MOA:&quot;C-C chemokine receptor type 2 antagonist&quot;) association a(ChEMBL:CHEMBL2110727)                                     </td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/plain": [
       "<pybel.struct.graph.SummarizeDispatch at 0x2a190a46530>"
      ]
     },
     "execution_count": 55,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "itmp_graph.summarize"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "id": "04ffd9c3",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[Text(0.5, 0, 'Node'),\n",
       " Text(0, 0.5, 'Number'),\n",
       " Text(0.5, 1.0, 'KG nodes in numbers (Total=35952)')]"
      ]
     },
     "execution_count": 56,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEWCAYAAACnlKo3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAmnklEQVR4nO3de7wVdb3/8ddbvGSpobnjIEhQYeWlSEk9p1TKVDRPWMdK8iiVR7S0n51TntBT6bEsO2aWZRopXrpIqJlkGCF56eKFjSEXlUDUBBFRNPCGYZ/fH9/vkmG59t6LYa+12Oz38/FYjzXrM9+Z+c6sy2fm+501o4jAzMysjM1aXQEzM+u5nETMzKw0JxEzMyvNScTMzEpzEjEzs9KcRMzMrDQnEWsZSQ9Jen83z3OepBHdOc9GkBSS3tyiZZ8g6TtNXN4ISYubtbxGkrSVpPsltbW6LhsLJ5FNWPWPtKSjJD0l6YD8ektJX5E0X9KzkpZIulHSwa2r9YaJiN0i4pZW12NjJWlL4EvAuZL2k/RMfjybE9szhcegDuZxpqSfNKGuP5G0VNJKSX+R9B+FcYNr1PfLhfEDJF0vaYWkxZJOrJp35HWuTHtJYdypkuZKWiXpQUmnVsZFxGpgAjCusWvfc2ze6gpYc0gaA3wb+EBE/CmHrwEGAMcCf86x9wEfAH7b9EraepO0eUSsWY9JRgH3R8QSYAmwTZ7PYOBBoO96zq+RvgEcFxGrJb0VuEXSnyNiZqFMR/X9CXAPcCSwK3CzpPkRcXOhzDsiYmGNaUX6TswG3gT8VtIjETExj/8ZMEvS6Tmp9Go+EukFJJ0AnAccUkkg+QjlIGBURNwZES/mx28i4pRO5hWSTpS0QNLTki6UpDxuM0lfkvSwpMclXSnptYVpj8njnpT0P1Xz3UzSOEkP5PGTJO2Qx70q75U+mZc5Q1K/Dur38tFX3mOelOuxKjd1DS+5buvsfRf2hDfPr2+R9DVJf8p7tr+S9DpJP8170jPyD3XRYZIWSXpC0rmSNivM/1OS7stHjlMlvaGqnidJWgAsUHJ+3uYrJc2RtHsHq3kocGtH26CwjJ0kTc578gslHZ/jI4HTgY/l9bwnxz+Z67sqr9MJXS2jKxExr/AjHfnxpjrqvg0wAjg7Iv4eEfeQdpg+Vedy/y8i7o6INRExH7geeHdh/GLgKWDf9VmfTZWTyKbv08BZwIER0V6Ivx+4M38h1tfhwLuAtwMfBQ7J8U/kx3uBN5L2cr8PIGlX4CLgGGAn4HXAwMI8PwscARyQxz8FXJjHjQFeC+ycpzsReL7Oun4QmAj0BSZX6lNi3epxFGn9BpB+7G4HLgN2AO4Dzqgq/yFgOLAn6QjhUwCSRpF+qD8MtAG/B66qmvYIYB/SXvbBwP7ALqTt9FHgyQ7quAcwv451mQgsJr0XRwJfl/S+iPgN8HXg5xGxTUS8I5d/nLTttgM+CZwvac9aM5Z0Q07StR43VJX9gaTngPuBpcCUqtk9rNRcdZmkHSuTVT1XhqsT622SHpP0ixoJvrJ8AfsB86pG3Qe845VT9D5OIpu+g4A7gDlV8R2BxyovJO2Qv8R/k/RCF/M8JyKejoi/AjcDw3L8aODbEbEoIp4BTgOOynvrRwI3RMRtee/yy8A/CvM8EfifiFicx58JHJmn/Tspebw5Il6KiJkRsbLO9f9DREyJiJeAH9P1F7+jdavHZRHxQET8DbgReCAibsrNLVcD76wq/82IWJGX9R1gdI6fCHwjIu7L034dGFY8GsnjV0TE86Ttsy3wVkB5uqUd1LEvsKqzlZC0M2nP+4sR8UJEzAIuITXx1BQRv87rHhFxK6k5dL8Oyh4eEX07eBxeVfYzed32A34BVI5MniAl+zcAe+UyP83TrAL+CHw5H8XuCfwb8OrCrA8ABpO22aPADZWjyipnkn4nL6uKryJty17PSWTT92nSHuollaaZ7Emgf+VF/kHqS/pCbtXFPB8rDD9Hblcn7bU+XBj3MKnfrV8e90hhec+y7t7yG4DrKnukpD29l/K0PwamAhMlPSrp/yRt0UUdO6rrqzr4sehq3eqxrDD8fI3X1fN6pDD8MGkbQdoW3y1sixWkPekBtaaNiN+RjrAuBB6XNF7Sdh3U8SnSD25ndgJW5B/jYv0GdFAeSYdKuiM3fz0NHEbaUdlgecfhD6Qj10/n2DMR0Z6bnJYBJwMHS6qs29HAENJ2uojUR7K4MM/bcvPt08ApuezbqtbpZFLi/ECNvo9tgae7Y/16OieRTd8y4EDSntwPCvHpwLskDaw5VTmPkn4AKwYBa3IdlpKaowCQ9GrS0UXFI8ChVXulr4qIJbld+38jYlfgX0jNJh3uFTfIs6y7J/tP3TDPnQvDg0jbD9K2OKFqW2xdOCECUv/A2hcRF0TEXqTmrV2AU6ltdh7fmUeBHQo/yJX6Lam1bElbAdcC3wL65Z2RKazbnFQsf6PWPauq+Lixk3ptTsd9IpU6bQYQEQ/nI562iNiHlNDu6mTeUayvpE+RzsA6sIMm37eROu57PSeRXiAiHiUlkpGSzs+x35Kaa34paR+l0323YMM6C68C/lPSkNy5WWk7X0Pq2Dxc0nuUTjM9i3U/fxcDZ1eabCS15b4BJL1X0h6S+gArSc03xaawZpgF7C9pkNLJAqd1wzxPlbR9bj46Bfh5jl8MnCZpNwBJr5X0kY5mIuld+T3cgpTsXqDj7TOF1JTToYh4BPgT8I3cHPR24DjS3jyknYLBWnsiwJako9flwBpJh5L6aTqa/6G5P6XW49C8Tq9XOiV9G0l9JB1Cau6bnsfvI+ktSidkvA64ALglNyUi6W2Sts2f63/P9fl2HrebpGF5vtuQTjpZQjr6RdLRpM/uQRGxqMb2HkDq57qjs+3YWziJ9BK53f19pH6Gb+Twh4AbSD8OT5NO8Tya9etMLppAanq6Lc/rBVKHORExDziJdHrkUlKzSnEP77ukju/fSlpF+oLuk8f9EykJrSR90W/Ny2maiJhG+pGfDcwkbbcNdX2e1yzg18CleVnXAd8kNd+tBOaSzqrqyHbAj0jb9GFSM+G5HZT9FfBWSTt1ML5iNKnP4FHgOuCMiLgpj7s6Pz8p6e7c7PX/gEm5Dh8nvZcbIkhNV5Uzob4FfC4iKvN9I/AbUt/EXFJfyejC9IcAi/K0JwIjI2J5HteP9F6uzGUGA4dHxN/z+K+RjpJnFI6QLi7M++PAFTWauHolhW9KZdarSBoL7BoRn2t1XXqa3HR3D7B/RDze6vpsDJxEzMysNDdnmZlZaU4iZmZWmpOImZmV1usuwLjjjjvG4MGDW10NM7MeZebMmU9ExCsugd/rksjgwYNpb2/vuqCZmb1M0sO14m7OMjOz0pxEzMysNCcRMzMrzUnEzMxKcxIxM7PSGpZEJO0s6WZJ9yrdlvSUHN9B0jSlW5BOk7R9jkvSBUq34pytwl3RJI3J5Rco3Su8Et9L6VagC/O0NS89bWZmjdHII5E1wOfzPSD2BU5SukXqOGB6RAwlXdZ5XC5/KDA0P8aSbiSD0n22zyBd0XVv4IxK4sllji9MN7KB62NmZlUalkQiYmlE3J2HV5Eu4T2AdC/pK3KxK0j3iibHr8y317wD6CupP+mSztPynfeeAqaR7ovRH9guIu6IdBXJKwvzMjOzJmhKn4ikwaT7S99JuvNZ5f7Pj5Gu7Q8pwRRvF7o4xzqLL64Rr7X8sZLaJbUvX768VhEzMyuh4f9Yz3cOu5Z0Q5mVxW6LiAhJDb8WfUSMB8YDDB8+3Ne+tx7r3d97d6ursNH442f/2OoqGA0+Esm367wW+GlE/CKHl+WmKPJz5cYuS1j3ntMDc6yz+MAacTMza5JGnp0l0u0+74uIbxdGTQYqZ1iNId0itBI/Np+ltS/wt9zsNRU4ON+LenvSvZKn5nErJe2bl3VsYV5mZtYEjWzOejdwDDBH0qwcOx04B5gk6TjS/aA/msdNAQ4DFgLPAZ8EiIgVkr4KzMjlzoqIFXn4M8DlwNbAjflhZmZN0rAkEhF/ADr638aBNcoHcFIH85oATKgRbwd234BqmpnZBvA/1s3MrDQnETMzK81JxMzMSnMSMTOz0pxEzMysNCcRMzMrzUnEzMxKcxIxM7PSnETMzKw0JxEzMyvNScTMzEpzEjEzs9KcRMzMrDQnETMzK81JxMzMSnMSMTOz0pxEzMystEbeY32CpMclzS3Efi5pVn48VLltrqTBkp4vjLu4MM1ekuZIWijpgnw/dSTtIGmapAX5eftGrYuZmdXWyCORy4GRxUBEfCwihkXEMOBa4BeF0Q9UxkXEiYX4RcDxwND8qMxzHDA9IoYC0/NrMzNrooYlkYi4DVhRa1w+mvgocFVn85DUH9guIu7I92C/Ejgijx4FXJGHryjEzcysSVrVJ7IfsCwiFhRiQyT9WdKtkvbLsQHA4kKZxTkG0C8ilubhx4B+HS1M0lhJ7ZLaly9f3k2rYGZmrUoio1n3KGQpMCgi3gn8F/AzSdvVO7N8lBKdjB8fEcMjYnhbW1vZOpuZWZXNm71ASZsDHwb2qsQiYjWwOg/PlPQAsAuwBBhYmHxgjgEsk9Q/IpbmZq/Hm1F/MzNbqxVHIu8H7o+Il5upJLVJ6pOH30jqQF+Um6tWSto396McC1yfJ5sMjMnDYwpxMzNrkkae4nsVcDvwFkmLJR2XRx3FKzvU9wdm51N+rwFOjIhKp/xngEuAhcADwI05fg5wkKQFpMR0TqPWxczMamtYc1ZEjO4g/okasWtJp/zWKt8O7F4j/iRw4IbV0szMNoT/sW5mZqU5iZiZWWlOImZmVpqTiJmZleYkYmZmpTmJmJlZaU4iZmZWmpOImZmV5iRiZmalOYmYmVlpTiJmZlaak4iZmZXmJGJmZqU5iZiZWWlOImZmVpqTiJmZleYkYmZmpTXy9rgTJD0uaW4hdqakJZJm5cdhhXGnSVooab6kQwrxkTm2UNK4QnyIpDtz/OeStmzUupiZWW2NPBK5HBhZI35+RAzLjykAknYl3Xt9tzzNDyT1kdQHuBA4FNgVGJ3LAnwzz+vNwFPAcdULMjOzxmpYEomI24AVdRYfBUyMiNUR8SCwENg7PxZGxKKIeBGYCIySJOB9wDV5+iuAI7qz/mZm1rVW9ImcLGl2bu7aPscGAI8UyizOsY7irwOejog1VfGaJI2V1C6pffny5d21HmZmvV6zk8hFwJuAYcBS4LxmLDQixkfE8IgY3tbW1oxFmpn1Cps3c2ERsawyLOlHwA355RJg50LRgTlGB/Engb6SNs9HI8XyZmbWJE09EpHUv/DyQ0DlzK3JwFGStpI0BBgK3AXMAIbmM7G2JHW+T46IAG4GjszTjwGub8Y6mJnZWg07EpF0FTAC2FHSYuAMYISkYUAADwEnAETEPEmTgHuBNcBJEfFSns/JwFSgDzAhIublRXwRmCjpa8CfgUsbtS5mZlZbw5JIRIyuEe7whz4izgbOrhGfAkypEV9EOnvLzMxaxP9YNzOz0pxEzMysNCcRMzMrzUnEzMxKcxIxM7PSnETMzKw0JxEzMyvNScTMzEpzEjEzs9KcRMzMrDQnETMzK81JxMzMSnMSMTOz0pxEzMysNCcRMzMrzUnEzMxKcxIxM7PSGpZEJE2Q9LikuYXYuZLulzRb0nWS+ub4YEnPS5qVHxcXptlL0hxJCyVdIEk5voOkaZIW5OftG7UuZmZWWyOPRC4HRlbFpgG7R8Tbgb8ApxXGPRARw/LjxEL8IuB4YGh+VOY5DpgeEUOB6fm1mZk1UcOSSETcBqyoiv02Itbkl3cAAzubh6T+wHYRcUdEBHAlcEQePQq4Ig9fUYibmVmTtLJP5FPAjYXXQyT9WdKtkvbLsQHA4kKZxTkG0C8ilubhx4B+HS1I0lhJ7ZLaly9f3k3VNzOzliQRSf8DrAF+mkNLgUER8U7gv4CfSdqu3vnlo5ToZPz4iBgeEcPb2to2oOZmZla0ebMXKOkTwOHAgfnHn4hYDazOwzMlPQDsAixh3SavgTkGsExS/4hYmpu9Hm/SKpiZWdbUIxFJI4H/Bj4YEc8V4m2S+uThN5I60Bfl5qqVkvbNZ2UdC1yfJ5sMjMnDYwpxMzNrki6TiKQ+ku5f3xlLugq4HXiLpMWSjgO+D2wLTKs6lXd/YLakWcA1wIkRUemU/wxwCbAQeIC1/SjnAAdJWgC8P782M7Mm6rI5KyJekjRf0qCI+Gu9M46I0TXCl3ZQ9lrg2g7GtQO714g/CRxYb33MzKz71dsnsj0wT9JdwLOVYER8sCG1MjOzHqHeJPLlhtbCzMx6pLqSSETcKukNwNCIuEnSq4E+ja2amZlt7Oo6O0vS8aQO7x/m0ADglw2qk5mZ9RD1nuJ7EvBuYCVARCwAXt+oSpmZWc9QbxJZHREvVl5I2pxO/iFuZma9Q71J5FZJpwNbSzoIuBr4VeOqZWZmPUG9SWQcsByYA5wATAG+1KhKmZlZz1Dv2Vn/kHQFcCepGWt+5bpXZmbWe9WVRCR9ALiYdNkRkS7bfkJE3Nj5lGZmtimr98+G5wHvjYiFAJLeBPyade8HYmZmvUy9fSKrKgkkWwSsakB9zMysB+n0SETSh/Ngu6QpwCRSn8hHgBkNrpuZmW3kumrO+tfC8DLggDy8HNi6ITUyM7Meo9MkEhGfbFZFzMys56n37KwhwGeBwcVpfCl4M7Perd6zs35JuqHUr4B/NKw2ZmbWo9R7dtYLEXFBRNwcEbdWHl1NJGmCpMclzS3EdpA0TdKC/Lx9jkvSBZIWSpotac/CNGNy+QWSxhTie0mak6e5IN+H3czMmqTeJPJdSWdI+mdJe1YedUx3OTCyKjYOmB4RQ4Hp+TXAocDQ/BgLXAQp6QBnAPsAewNnVBJPLnN8YbrqZZmZWQPV25y1B3AM8D7WNmdFft2hiLhN0uCq8ChgRB6+ArgF+GKOX5kvp3KHpL6S+uey0yJiBYCkacBISbcA20XEHTl+JXAE/gOkmVnT1JtEPgK8sXg5+A3QLyKW5uHHgH55eADwSKHc4hzrLL64RvwVJI0lHd0waNCgDay+mZlV1NucNRfo290Lz0cdDb+QY0SMj4jhETG8ra2t0YszM+s16j0S6QvcL2kGsLoSLHmK7zJJ/SNiaW6uejzHlwA7F8oNzLElrG3+qsRvyfGBNcqbmVmT1JtEzujGZU4GxgDn5OfrC/GTJU0kdaL/LSeaqcDXC53pBwOnRcQKSSsl7Uu6RP2xwPe6sZ5mZtaFeu8n0uXpvLVIuop0FLGjpMWkZHQOMEnSccDDwEdz8SnAYcBC4Dngk3nZKyR9lbXX6jqr0skOfIZ0BtjWpA51d6qbmTVRvf9YX8XavostgS2AZyNiu86mi4jRHYw6sEbZAE7qYD4TgAk14u3A7p3VwczMGqfeI5FtK8P5D32jgH0bVSkzM+sZ6j0762WR/BI4pPurY2ZmPUm9zVkfLrzcDBgOvNCQGpmZWY9R79lZxfuKrAEeIjVpmZlZL1Zvn0ivua/IXqde2eoqbDRmnntsq6tgZhu5rm6P+5VORkdEfLWb62NmZj1IV0ciz9aIvQY4Dngd4CRiZtaLdXV73PMqw5K2BU4h/QlwInBeR9OZmVnv0GWfSL6fx38BR5Mu3b5nRDzV6IqZmdnGr6s+kXOBDwPjgT0i4pmm1MrMzHqErv5s+HlgJ+BLwKP5gocrJa2StLLx1TMzs41ZV30i6/2PdjMz6z2cJMzMrDQnETMzK81JxMzMSnMSMTOz0pxEzMystKYnEUlvkTSr8Fgp6XOSzpS0pBA/rDDNaZIWSpov6ZBCfGSOLZQ0rtnrYmbW29V7KfhuExHzgWEAkvoAS4DrSJdTOT8ivlUsL2lX4ChgN9J/Vm6StEsefSFwELAYmCFpckTc24z1MDOzFiSRKgcCD0TEw+muuzWNAiZGxGrgQUkLgb3zuIURsQhA0sRc1knEzKxJWt0nchRwVeH1yZJmS5ogafscGwA8UiizOMc6ir+CpLGS2iW1L1++vPtqb2bWy7UsiUjaEvggcHUOXQS8idTUtZRuvEpwRIyPiOERMbytra27Zmtm1uu1sjnrUODuiFgGUHkGkPQj4Ib8cgmwc2G6gTlGJ3EzM2uCVjZnjabQlCWpf2Hch4C5eXgycJSkrSQNAYYCdwEzgKGShuSjmqNyWTMza5KWHIlIeg3prKoTCuH/kzQMCOChyriImCdpEqnDfA1wUkS8lOdzMjAV6ANMiIh5zVoHMzNrURKJiGdJt9ctxo7ppPzZwNk14lOAKd1eQTMzq0urz84yM7MezEnEzMxKcxIxM7PSnETMzKw0JxEzMyvNScTMzEpzEjEzs9KcRMzMrDQnETMzK81JxMzMSnMSMTOz0pxEzMysNCcRMzMrrdX3WLdN3F/P2qPVVdhoDPrKnFZXwazb+UjEzMxKcxIxM7PSnETMzKy0liURSQ9JmiNplqT2HNtB0jRJC/Lz9jkuSRdIWihptqQ9C/MZk8svkDSmVetjZtYbtfpI5L0RMSwihufX44DpETEUmJ5fAxwKDM2PscBFkJIOcAawD7A3cEYl8ZiZWeO1OolUGwVckYevAI4oxK+M5A6gr6T+wCHAtIhYERFPAdOAkU2us5lZr9XKJBLAbyXNlDQ2x/pFxNI8/BjQLw8PAB4pTLs4xzqKr0PSWEntktqXL1/enetgZtartfJ/Iu+JiCWSXg9Mk3R/cWREhKTojgVFxHhgPMDw4cO7ZZ5mZtbCI5GIWJKfHweuI/VpLMvNVOTnx3PxJcDOhckH5lhHcTMza4KWJBFJr5G0bWUYOBiYC0wGKmdYjQGuz8OTgWPzWVr7An/LzV5TgYMlbZ871A/OMTMza4JWNWf1A66TVKnDzyLiN5JmAJMkHQc8DHw0l58CHAYsBJ4DPgkQESskfRWYkcudFRErmrcaZma9W0uSSEQsAt5RI/4kcGCNeAAndTCvCcCE7q6jmZl1bWM7xdfMzHoQJxEzMyvNScTMzEpzEjEzs9KcRMzMrDQnETMzK81JxMzMSnMSMTOz0pxEzMysNCcRMzMrzUnEzMxKcxIxM7PSnETMzKw0JxEzMyvNScTMzEpzEjEzs9KcRMzMrLSmJxFJO0u6WdK9kuZJOiXHz5S0RNKs/DisMM1pkhZKmi/pkEJ8ZI4tlDSu2etiZtbbteL2uGuAz0fE3ZK2BWZKmpbHnR8R3yoWlrQrcBSwG7ATcJOkXfLoC4GDgMXADEmTI+LepqyFmZk1P4lExFJgaR5eJek+YEAnk4wCJkbEauBBSQuBvfO4hfl+7UiamMs6iZiZNUlL+0QkDQbeCdyZQydLmi1pgqTtc2wA8EhhssU51lG81nLGSmqX1L58+fLuXAUzs16tZUlE0jbAtcDnImIlcBHwJmAY6UjlvO5aVkSMj4jhETG8ra2tu2ZrZtbrtaJPBElbkBLITyPiFwARsaww/kfADfnlEmDnwuQDc4xO4mZm1gStODtLwKXAfRHx7UK8f6HYh4C5eXgycJSkrSQNAYYCdwEzgKGShkjaktT5PrkZ62BmZkkrjkTeDRwDzJE0K8dOB0ZLGgYE8BBwAkBEzJM0idRhvgY4KSJeApB0MjAV6ANMiIh5zVsNMzNrxdlZfwBUY9SUTqY5Gzi7RnxKZ9OZmVlj+R/rZmZWmpOImZmV5iRiZmalOYmYmVlpTiJmZlaak4iZmZXWkn+sm5m12q37H9DqKmw0Drjt1tLT+kjEzMxKcxIxM7PSnETMzKw0JxEzMyvNScTMzEpzEjEzs9KcRMzMrDQnETMzK81JxMzMSnMSMTOz0np8EpE0UtJ8SQsljWt1fczMepMenUQk9QEuBA4FdiXdp33X1tbKzKz36NFJBNgbWBgRiyLiRWAiMKrFdTIz6zUUEa2uQ2mSjgRGRsR/5NfHAPtExMlV5cYCY/PLtwDzm1rRcnYEnmh1JTYR3pbdy9uze/WU7fmGiGirDvaKS8FHxHhgfKvrsT4ktUfE8FbXY1Pgbdm9vD27V0/fnj29OWsJsHPh9cAcMzOzJujpSWQGMFTSEElbAkcBk1tcJzOzXqNHN2dFxBpJJwNTgT7AhIiY1+JqdZce1fy2kfO27F7ent2rR2/PHt2xbmZmrdXTm7PMzKyFnETMzKw0J5GSJL0kaZakuZKulvTqTsqOkPQvhdeX5/+41LuswZLmbmidu1thG9wj6e7KOkraSdI1XUw7QtING7DsP5Wc7kxJX8jDl0t6MK/D3ZL+uWx9NgaSjpAUkt6aX2/QNu5iWev1Gd5Urc/vQI1ph0k6rI5ywyVdsGE1bRwnkfKej4hhEbE78CJwYidlRwD/0sn4nqqyDd4BnAZ8AyAiHo2Ihv7ARER3bc9TI2IYMA74YfXIfGmdnmI08If8bM3R6e+ApM5OXhoGdJlEIqI9Iv7fBtWygZxEusfvgTdL+ldJd0r6s6SbJPWTNJj0wfrPvMeyX55mf0l/krSosken5Ny8VzNH0seqFyTpVZIuy+P/LOm9Of5qSZMk3SvpulyP4ZI+Jek7hemPl3R+A7bBdsBTeRkvHzl1VN+qddpB0i8lzZZ0h6S353ibpGmS5km6RNLDknbM454pTP/FPP97JJ1TWM8ZOXZtHXuItwFvztM+JOmbku4GPiJpdJ7/XEnfLCx3ZD6CuUfS9Bx7jaQJku7K6zsqx3fLsVl5PYfmsr/O08+t9X7XS9I2wHuA40inuldsl5cxX9LFkjarsf2OlHR5Hr5c0gUdfDa/n+dzE/D6wvRfydt6rqTxkpTjt+TteJekv1Q++5L6SPpWLj9b0mdzfC9Jt0qaKWmqpP5lt0eLVH4HRkj6vaTJwL21vgNKf0k4C/hY/kx8rJPPzstHlEpH0hPytl0kqfXJJSL8KPEAnsnPmwPXA58GtmftGW//AZyXh88EvlCY9nLgalIS35V0/S+AfwOmkU5X7gf8FegPDAbm5jKfJ53KDPDWXOZVwBeAH+b47sAaYDiwDfAAsEUe9ydgj27aBi8Bs4D7gb8Be+V4PfUdAdyQ498DzsjD7wNm5eHvA6fl4ZFAADtWbf9D8zq9Or/eIT+/rlDPrwGfrX4v8vtwZB7+CHBnHn4I+O88vFOuc1t+r38HHJFfPwIMqVru14F/z8N9gb8Ar8nreHSObwlsnd/vHxXq+doNeC+OBi4tvMd75W38AvDG/JmaVljfZwrTHglc3sVn88Os/WzuBDxdmNcOhXn9GPjXPHwLa78DhwE35eFPA9cAm1emB7bI9W7LsY+RPzcb84PavwMjgGcLn42OvgOfAL5fmFdHn50RrP2unJm301aky6U8Sf5ut+rhI5HytpY0C2gnfSguJf1jfqqkOcCpwG6dTP/LiPhHRNxLShiQ9iSvioiXImIZcCvwrqrp3gP8BCAi7gceBnbJ8Yk5PheYnYefIf3wHa7UVr5FRMzZkBUvqBzKv5X0I39lZS+0jvpWl/lxLvM74HWStqtap9+Qj3SqvB+4LCKey+VW5PjueW9wDukHtqP34tz8Po4l7cVX/Dw/vwu4JSKWR8Qa4KfA/sC+wG0R8WDVcg8GxuV53kL6sRgE3A6cLumLpGsQPQ/MAQ7Ke+v7RcTfOqhjPUaTt1V+rjRp3RXpAqUvAVeRtmlXan0292ftZ/NR0meq4r1KR75zSDsBxW39i/w8k7RzAek9+2HenpVt9xbSzs+0vO2+RPo+bexq/Q5A2u4P5uF6vgPQ8Wen2q8jYnVEPAE8ztr3qCV69J8NW+z5SG3pL5P0PeDbETFZ0gjSXkNHVhcn7e7KVbkEOJ10xHBZIxYQEbcrNTW94gJtLXI5cERE3CPpE6S9uVpOjYhaJwE8W3K5Av4tIqov8nmfpDuBDwBTJJ0QEb+TtCdpL/1rkqZHxFnrvUBpB9KP9x6SgnS0EMCv83NRVD1D+rEqqvuzKelVwA+A4RHxiKQzq+ZXmddLdP57I2BeRPS0kxtq/Q5Auc9Pzc+OpOokUXx/utquDecjke71WtZeu2tMIb4K2LaO6X9PaiPtI6mNtPd3V40yRwNI2oW0pzIf+CPw0RzfFdijMkFE3Em6xtjHSXuj3S4f5fQhHV7XU9+OyowAnoiIlay7TgeTmgurTQM+qdznkX9QIW3vpZK2qMy7pLuAAyTtqNTJPpp0hHgHqV9rSNVypwKfLfQLvDM/vxFYFBEXkJo93i5pJ+C5iPgJcC6wZ8k6Hgn8OCLeEBGDI2Jn4EFgP2BvpcsCbUZqIvpDnmaZpLfl+IfqWMZtrP1s9gcqfVuVhPGEUr9MPSdUTANOUO50zttuPtCmfIacpC0kdXYk35N09B2o/l2o+dnZ2DmJdK8zgaslzWTdSzv/CviQ1u1Yr+U6UjPUPaTmgv+OiMeqyvwA2Cw3Hfwc+ERErM7xNkn3kvoA5pH6KSomAX+MiFpNQmVtnddpVq7LmNxsUk99i84E9pI0GziHtQn4f4GDlTrpPwI8RvrivSw3c00G2nM9vpBHfRm4k5SI7i+7ghGxlHTm1s2k92VmRFwfEctJTWC/kHQPa5u/vkpq358taV5+DSkZzs113B24kpTo78qxM0jvWxmjSZ+domtzfAapb+k+UmKplBsH3EBqX19axzKuAxYA9+a63w4QEU8DPwLmkn4EZ9Qxr0tITT+z87b7eKT7AR0JfDPHZrHpnNHY0XfgZmDXSsc6HX92Nmq+7MkmIu8lbxERL0h6E3AT8Jb85SSf3XF+RExvZT3Xh6StgJciXSPtn4GLqpsOzKy13Cey6Xg1cHNuvhHwmYh4UVJfUpPMPT0pgWSDgEm5yeVF4PgW18fMqvhIxMzMSnOfiJmZleYkYmZmpTmJmJlZaU4iZg2idEXd8wqvv5D/jLc+83im61JmreMkYtY4q4EP53/ym22SnETMGmcN6f7Z/1k9QulKx79TuortdEmDcnyIpNuVrvj6tappTlW6Wu5sSf/bnFUw65yTiFljXQgcLem1VfHvAVdExNtJF3Ws3HTou6Q/Ve5B4Z/k+bIvQ4G9Sfeh2EvS/g2uu1mX/D8RswaR9ExEbCPpLODvwPPANhFxpqQngP4R8ff8B9GlEbGjpCeBf8rx7YBH8zy+RbosyNN59tsA34iIS1+5ZLPm8T/WzRrvO8Dd1H8F5Vp7diIljVfcfdGsldycZdZg+X4Zk1j3fiV/Yu0dCI8mXekV0gUji/GKqcCn8pVykTRA0usxazEnEbPmOI90J7qKz5IuYT8bOAY4JcdPAU7KV3wdUCkcEb8FfgbcnsddQ323FzBrKPeJmJlZaT4SMTOz0pxEzMysNCcRMzMrzUnEzMxKcxIxM7PSnETMzKw0JxEzMyvt/wNCxn4nu4eV4AAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "node_data = {'Node':['Pathology','BiologicalProcess','Abundance','Protein'],\n",
    "       'Number':[9005,5669,19931,1347]}\n",
    "node = pd.DataFrame(node_data)\n",
    "a = sns.barplot(x=\"Node\", y=\"Number\", data=node)\n",
    "a.set(xlabel='Node',ylabel='Number',title= 'KG nodes in numbers (Total=35952)')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "id": "3b0b95f3",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[Text(0.5, 0, 'Number'),\n",
       " Text(0, 0.5, 'Namespace'),\n",
       " Text(0.5, 1.0, 'KG Namespace in numbers')]"
      ]
     },
     "execution_count": 57,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbsAAAEWCAYAAAD/6zkuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAApE0lEQVR4nO3debxd0/3/8dcbISLEEENEiXmIIeSaihIUpUX7LeJLa2h/qVYpranV+kbnGlpFi1SVfhs1FVVfSgdBUxL3kuQKiSlBIkUIkZgiPr8/9jrsHHc4995z7rl35/18PM7j7r322nt/zkru+dy19j57KSIwMzMrsmXqHYCZmVmtOdmZmVnhOdmZmVnhOdmZmVnhOdmZmVnhOdmZmVnhOdmZWc1I2kPS9HrH0R5Jx0r6V73jsNpxsrPCkzRT0r659ZGS5knaM60vL+kcSdMlLZQ0W9KdkvZr45ghqVnSMrmyH0q6uqZvppeJiPsjYvN6x2HmZGdLFUnHAL8CDoqIe1PxTcAhwBeB1YANgV8CB7VzuHWBkTUK1XopScvVOwb7KCc7W2pI+gpwIbB/RPw7le0LfBI4JCImRMS76fXXiPhGO4c8Dzi3tQ83STdK+o+k1yXdJ2lobtvVkn6depALJI2XtI6ki1Kvc5qk7XP115X0J0kvS5oh6eTctp0kNUqaL+lFST9P5UNSD3SUpBckzZF0Wtl+D0h6LW27VNLyue1DJf1N0qvpuN9J5ctIOkvS05JekXSDpNVbaYO9JM3Krc+UdJqkKaldrpfUt5V9j5X0L0kXpDaZIelTZcfK99hHS/pD2Xs/TtLzaf8TJO2Yzv2apEs/ekpdmuKaJmmf3IYBkn6b2ml26sUvm4tzvKRfSHoFGC1pE0n3pmPNlXR9S+/Ruo+TnS0tvgp8H9gnIhpz5fsCEyJiVsu7telmYD5wbCvb7wQ2BdYCHgbGlm0/HPguMBB4B3gg1RtI1tssJa1lgL8Ak4HBwD7AKZL2T8f5JfDLiFgF2Bi4oew8I1Ic+wFn5hLEYuDUdL5d03G/ls65MvB34K9kPdhNgH+k/U4CDgX2TNvmkfWWK3U4cABZD3pbWm8/gJ2B6SnG84DfSlIHzrUz2Xs/ArgIOJvs33wocHhpKDtX9+l0rv8Bbs4l8auB98jaYXuytvxy2b7PAGsDPwJ+ANxNNlKwHnBJB2K2GnCys6XFJ4EHgeay8oHAf0orklZPf/W/Luntdo4ZwPeA7+V7RB9sjLgqIt6IiHeA0cB2kgbkqtwSEU0R8TZwC/B2RPw+IhYD15N9qALsCKwZEd9Pvc5ngN/w4RDqImATSQMjYkFEPFgWyrkRsTAimoHfAUem+Joi4sGIeC8iZgJXkCUwgE8D/4mICyPi7fQ+JqRtJwBnR8Ss3Hv7fAeG7y6OiBci4lWyJD6sjbrPRsRvUptcAwwiSyiV+kGK/25gIfDHiHgpImYD9/NhGwO8BFwUEYsi4nqyJHuQpLWBA4FTUju+BPyCJYewX4iIS1JbvkX2b7IBsG46v29+qTMnO1tafBXYDLiyrGfwCtkHKAAR8WpErAoMB1Zo76ARcQcwC/hKvlzSspJ+mob65gMz06aBuWov5pbfamG9f1reAFg3JeHXJL0GfIcPP/S/lN7bNEkPSfp0WZjP55afJeuNIWkzSbcrG2qdD/w4F9/HyHo5LdkAuCUXy+NkvcRKk9B/cstv5t5nm3Uj4s202Fb9cpW2McDsWPLJ+KW22gDoA8zJvecryHrsJfk2BjgDEDBR0lRJx3cgZqsBJztbWrxINky3B/DrXPk/gB0lrdeFY59Nlnz65cr+m+yml32BAcCQVN6RIbiS54EZEbFq7rVyRBwIEBFPRsSRZB++PwNukrRSbv+P5ZbXB15Iy5cB04BN0xDod3LxPQ9s1EY8nyqLp2/qLXWnhSzZ5ut08XiDy/4QKrXV82TDzANz73eViBiaq7vE9DER8Z+I+H8RsS7ZH0K/lrRJF+OzLnCys6VGRLxAlvAOkPSLVHY3cA9wq6SdlX0NoQ+wSweOOw54FDgmV7wy2QfkK2QfyD/uQugTgTcknSlpxdRr3FrSjgCSjpa0ZkS8D7yW9nk/t//3JPVTdoPMcWRDpKUY5wMLJG1B1vstuR0YJOkUSStIWlnSzmnb5cCPJG2Qzr+mpEO68P46axIwUlIfSQ3A57t4vLWAk9PxDgO2BO6IiDlk198ulLSKsht0Ni673rcESYfl/oCaR5YM32+tvtWek50tVSLiOWBvsmtMP0nFnyX7cP8DWbKYARwF7N/SMVrxXSB/R+LvyYbBZgOPkV0v7GzMi8muoQ1Lsc0FriTrMUJ2s8dUSQvIblYZma4bldwLPEXWi70gJXiA08h6oG+QXQP84I7BiHiD7DrnZ8iGEp8ku9GFdI7bgLslvZHeWykRdqfvkd2QMw84F7i2i8ebQHYzy1yym0w+HxGvpG1fBJYn+7ecR3YD0aCWDpLsCExI/ya3Ad9I11qtTuTJW82KSdIQsuTYJyLeq3M4ZnXlnp2ZmRWek52ZmRWehzHNzKzw3LMzM7PC8wNLe4iBAwfGkCFD6h2GmVmv0tTUNDci1myvnpNdDzFkyBAaGxvbr2hmZh+Q9Gwl9TyMaWZmheeeXQ/x+KxXGH767+sdhplZt2o6/4vdch737MzMrPCc7MzMrPCc7MzMrPCc7MzMrPCc7MzMrPBqmuwkrSPpujRbc5OkOySNknR7K/XHSZouaVJ63ZTKR0uK/OSHaZ6tSPNYIWmmpOa0X3N+fq00zUZrMV4kabYkJ34zs4Kq2Qd8mvH3FmBcRGwcEcOBbwNrt7PrURExLL3ykzE2AyNz64cBU8v2HRERw8gmcby4ghiXIZvL7Hmg1YkYzcysd6tlb2YEsCgiLi8VRMRk4H6gv6SbJE2TNDYlxvbcChwCIGlj4HWySRZbsgrZBIvt2YssYV4GHFkqlLRnrnf5SJqleZCk+1LZo5L2SHUvk9Qoaaqkc1PZ3pJuzR3vk5JuqSAeMzOrgVomu62Bpla2bQ+cAmwFbATslts2Npdozs+Vzweel7Q1WQ/vej7qHkmPks3M/N0KYjwS+CNZD/QgSX1S+WnAiamXuAfwFtmMznelsu2ASanu2RHRAGwL7ClpW+AeYAtJpee1HQdcVX7yNKTbKKnxvTffqCBcMzPrjHpdp5oYEbMi4n2ypDEkty0/jHl62X7XkSW6Q8kSVLkREbE1sA1wqaT+rQUgaXngQODWiJgPTAD2T5vHAz+XdDKwaprl+SHgOEmjgW0iopSdDpf0MPAIMBTYKrJ5k/4XOFrSqsCuwJ3lMUTEmIhoiIiG5fqt3FqoZmbWRbVMdlOB4a1seye3vJjKH1t2O/AF4LmUoFoUEU8DL5L1HFuzP7Aq0CxpJrA7aSgzIn4KfBlYERgvaYuIuA/4BDAbuFrSFyVtSNYL3CcitgX+D+ibjv874Oh0zBtTwjQzszqo5bMx/wn8WNKoiBgDkIb49ujsASPiTUlnAk+0VU/SWsCGQFtPwz4S+HJE/DHtsxIwQ1I/YFBENJMlwh3JhiTfAmZFxG8krQDsAEwGFgKvS1ob+BQwLsX6gqQXyIZT9+3sezYzs66rWbKLiJD0WeCilKDeBmaS3WjSlrEpsQDMjYglEkVEXNfGvvdIWgz0Ac6KiBdTeT9Js3L1fg0cAJyQO+5CSf8CPgPsLmkE8D5ZD/VOsuHT0yUtAhYAX4yIGZIeAaaR3dE5vvy9AGtGxOPtvGczM6shZZeXrBYkXQo8EhG/ba/uSutsGFt84dxuiMrMrOfo6qwHkprSTYJt8hQ/NSKpiWyI81v1jsXMbGnnZFcj6Uv0ZmbWA/gRWWZmVnhOdmZmVngexuwhtlxvDRq7aXp6M7OljXt2ZmZWeE52ZmZWeE52ZmZWeE52ZmZWeL5BpYd4d85Unvv+NvUOoxDWP6e53iGYWQ/jnp2ZmRWek52ZmRWek52ZmRWek52ZmRWek52ZmRVe4ZKdpHUkXSfpaUlNku6QNErS7a3UHydpuqRJ6XVTKh8tKSRtkqt7SiprSOszJTWn/ZolHZKru6DW79XMzCpTqK8eSBJwC3BNRIxMZdsBB7ez61ER0dhCeTPZDOU/TOuHkc1cnjciIuZK2hy4G/hzZ+M3M7PaKFrPbgSwKCIuLxVExGTgfqC/pJskTZM0NiXG9twKHAIgaWPgdWBuK3VXAeZ1JXgzM6uNQvXsgK2Bpla2bQ8MBV4AxgO7Af9K28ZKeist/y0iTk/L84HnJW1NlvSuB44rO+49KXFuBBxelXdhZmZVVbRk15aJETELQNIkYAgfJrvWhjEBriMbytwf2IePJrvSMObGwD8kjYuIiq7XSRoFjAIYPKBPB96KmZl1RNGGMacCw1vZ9k5ueTGVJ/rbgS8Az0XE/NYqRcTTwIvAVhUel4gYExENEdGw+krLVrqbmZl1UNGS3T+BFVKPCQBJ2wJ7dPaAEfEmcCbwo7bqSVoL2BB4trPnMjOz2ijUMGZEhKTPAhdJOhN4G5hJdqNJW/LX7OZGxL5lx72ujX3vkbQY6AOcFREvpvJ+kmbl6v08In5e4VsxM7MqUkTUOwYDth28Ytz+lU3ar2jt8qwHZksPSU0R0dBevaINY5qZmX2Ek52ZmRWek52ZmRWek52ZmRVeoe7G7M2WHzSU9c9p7XvtZmbWFe7ZmZlZ4TnZmZlZ4TnZmZlZ4TnZmZlZ4fkGlR5i2kvT2O2S3eodRrvGnzS+3iGYmXWYe3ZmZlZ4TnZmZlZ4TnZmZlZ4TnZmZlZ4TnZmZlZ4TnZmZlZ4S32yk7S2pGslPSOpSdIDabZzJO0uaaKkaek1KrffaEmzJU1K2y6TtEzadrWkGWnbw5J2rdf7MzOzpTzZSRJwK3BfRGwUEcOBkcB6ktYBrgVOiIgtgN2Br0g6KHeIX0TEMGArYBtgz9y209O2s4Arav1ezMysdUt1sgP2Bt6NiMtLBRHxbERcApwIXB0RD6fyucAZZMmr3PJAX2BeC9vuAzapduBmZla5pT3ZDQUebmNbU1lZYyovOVXSJGAO8ERETGrhOJ8Bmls6gaRRkholNS5asKgjcZuZWQcs7cluCZJ+JWmypIcq3KU0jLkWsJKkkblt56dEOAr4Uks7R8SYiGiIiIY+/ft0JXQzM2vD0p7spgI7lFYi4kRgH2BN4DFgeFn94WmfJUTEIuCvwCdyxadHxLCI+GREPFrtwM3MrHJLe7L7J9BX0ldzZf3Sz18Bx0oaBiBpDeBnwHnlB0k3uuwGPF3TaM3MrFOW6lkPIiIkHQr8QtIZwMvAQuDMiJgj6WjgN5JWBgRcFBF/yR3i1FSnDzAF+HX3vgMzM6vEUp3sACJiDtnXDVradh+wYyvbRgOjW9l2bHWiMzOzaljahzHNzGwp4GRnZmaF52RnZmaFt9Rfs+sptlhrC8afNL7eYZiZFZJ7dmZmVnhOdmZmVnhOdmZmVnhOdmZmVni+QaWHeGP6dO79xJ7tV6yCPe+7t1vOY2bWU7hnZ2ZmhedkZ2ZmhedkZ2ZmhedkZ2ZmhedkZ2ZmhedkZ2ZmhdehZCepX/u1eg9Ja0u6VtIzkpokPSDps2nb7pImSpqWXqNy+42WFJI2yZWdksoa0vpMSc2SJqXXx7v/HZqZGVSY7CR9XNJjwLS0vp2kXj0rtyQBtwL3RcRGETGcbBLX9SStA1wLnBARWwC7A1+RdFDuEM0sOenrYcDUstOMiIhh6fXvWr0XMzNrW6U9u18A+wOvAETEZOATtQqqm+wNvBsRl5cKIuLZiLgEOBG4OiIeTuVzgTOAs3L73wocAiBpY+B1YG73hG5mZh1R8TBmRDxfVrS4yrF0t6HAw21sayora0zlJfOB5yVtTdbDu76F49yThjAntHQSSaMkNUpqfH3Roo5Fb2ZmFas02T2frjmFpD6STgMer2Fc3U7SryRNlvRQB3a7jizRHQrc0sL20jDmzi3tHBFjIqIhIhoG9OnT8aDNzKwilSa7E8iG9gYDs4Fhab03mwrsUFqJiBOBfYA1gceA4WX1h/PRa3K3A18AnouI+bUL1czMuqKiZBcRcyPiqIhYOyLWioijI+KVWgdXY/8E+kr6aq6sdLfpr4BjJQ0DkLQG8DPgvPwBIuJN4EzgRzWP1szMOq3SuzGvkbRqbn01SVfVLKpuEBFBNvy4p6QZkiYC1wBnRsQc4GjgN5KmAf8GroqIv7RwnOtKN7KYmVnPVOkUP9tGxGullYiYJ2n72oTUfVJSG9nKtvuAHVvZNrqV8r1yy0O6HKCZmVVFpdfslpG0WmlF0up4LjwzM+slKk1YFwIPSLoREPB5fJ3KzMx6iYqSXUT8XlITMCIVfS4iHqtdWGZmZtVT8VBkREyV9DLQF0DS+hHxXM0iMzMzq5KKkp2kg8mGMtcFXgI2IPtS+dC29rPKrbz55ux53731DsPMrJAqvUHlB8AuwBMRsSHZl68frFlUZmZmVVRpsluUvkS+jKRlIuIeoKGGcZmZmVVNpdfsXpPUH7gfGCvpJWBh7cIyMzOrHmUPEmmnkrQS8DbZ1w6OAgYAYwvwyLAeY/11No0zjvp5t5zr6xd+plvOY2ZWa5KaIqLdkcZKv3qwME1ouhPwKnCXE52ZmfUWlT4b88vAROBzZF8of1DS8bUMzMzMrFoqvWZ3OrB9qTeXZgH4N9CrHwZtZmZLh0rvxnwFeCO3/kYqMzMz6/Eq7dk9BUyQ9GcggEOAKZK+CRAR3XNnhZmZWSdUmuyeTq+SP6efK1c3HDMzs+qr9G7Mc0vLkpYB+kfE/JpF1U0kLYiI/rn1Y4GGiPh6Wj8aOANYFngPeAg4LSJekzSOrB0aUt0G4ILSnHaSdgIuANYG3gSagJPT7OZmZtaNKr0b81pJq6Tv2z0KPCbp9NqGVl+SDgBOBT4VEUOBHchuylk7V20tSZ9qYd+1gRvJZj3fPCK2B/6Ke8JmZnVR6Q0qW6We3KHAncCGwBdqFVQPcTZZL242QEQsjoirImJ6rs75qV65E4FrIuKBUkFE3BQRL9Y0YjMza1Glya6PpD5kye62iFhEdqNKb7eipEmlF/D93LahwMPt7P8A8K6kEWXlW5MNW7ZJ0ihJjZIaF7z5ekfiNjOzDqg02V0BzARWAu6TtAHQ66/ZAW9FxLDSCzinpUqStkkJ8WlJR5Rt/iHw3c6cPCLGRERDRDT07zegM4cwM7MKVJTsIuLiiBgcEQdG5lk+nLW8qKaSXacjIppTMrwTWDFfKSL+mcp2Kdt3ePeEaWZm7an0BpW1Jf1W0p1pfSvgmJpGVn8/AS6QtF6ubMVW6v6Q7K7NkkuBYyTtXCqQ9Ll044qZmXWzSocxrwbuIpupHOAJ4JQaxNNjRMQdwMXAnZIek/RvYDFZO7RU9+Xc+ovASLJkOV3S48D+LPkUGjMz6yaVfql8YETcIOnbABHxnqTFNYyrW+S/Y5fWryZL7KX1a4BrWtl3r7L14WXrDwB7VCdSMzPrikp7dgvTw58DQNIugG8fNDOzXqHSnt03gduAjSWNB9Ykm+rHzMysx6v0cWEPS9oT2JxstvLp6bt2ZmZmPV5FyU7SssCBwJC0z36SPNuBmZn1CpUOY/4FeBtoBt6vXThLr7XWG8DXL/xMvcMwMyukSpPdehGxbU0jMTMzq5FK78a8U9J+NY3EzMysRirt2T0I3JLmsltEdpNKRMQqNYvMzMysSipNdj8HdgWaI6IIsx2YmdlSpNJk9zzwqBNd7cyZ8TQ/Orp+X108+w831e3cZma1VmmyewYYlx4E/U6p0F89MDOz3qDSZDcjvZZPLzMzs16j0ieonFvrQMzMzGql0ieorEk2X9tQoG+pPCL2rlFcZmZmVVPp9+zGAtOADYFzgZnAQzWKyczMrKoqTXZrRMRvgUURcW9EHA/02F6dpMWSJkmaKmmypG+l7wgiqUHSxfWO0czMuk+lN6iUZjiYI+kg4AVg9dqEVBVvRcQwAElrAdcCqwD/ExGNQGMdYzMzs25Wac/uh5IGAN8CTgOuBE6tWVRVFBEvAaOAryuzl6TbASTtmXqAkyQ9ImnlVH66pIckTZH0wc05km6V1JR6jKNS2bKSrpb0qKRmSaem8o0l/TXVv1/SFt3/7s3MDCq/G/P2tPg6MKJ24dRGRDyTpilaq2zTacCJETFeUn/g7fQM0E2Bncgei3abpE9ExH3A8RHxqqQVgYck/Yls2qPBEbE1gKRV07HHACdExJOSdgZ+TdnQb0qYowAG9Fux6u/bzMwybSY7See0sTki4gdVjqe7jQd+LmkscHNEzErJbj/gkVSnP1nyuw84WdJnU/nHUvl0YCNJlwD/B9ydEufHgRsllc61QvnJI2IMWVJk8Bqr+ek0ZmY10l7PbmELZSsBXwLWAHpFspO0EbAYeAnYslQeET+V9H9kE9OOl7Q/WW/uJxFxRdkx9gL2BXaNiDcljQP6RsQ8SdsB+wMnAIcDpwCvla4bmplZfbWZ7CLiwtJyup71DeA44Drgwtb260nSdwQvBy6NiMj1tJC0cUQ0A82SdgS2AO4CfiBpbEQskDSY7AadAcC8lOi2AHZJxxgIvBsRf5I0HfhDRMyXNEPSYRFxo7KTbhsRk7vzvZuZWabda3aSVge+CRwFXAPsEBHzah1YF60oaRLQB3gP+F+ymRvKnSJpBNns61OBOyPiHUlbAg+kxLgAOBr4K3CCpMfJhi4fTMcYDPyu9NUG4Nvp51HAZZK+m+K4DnCyMzOrg/au2Z0PfI7sutI2EbGgW6LqoohYto1t44BxafmkVur8EvhlC5s+1cphd2jhGDOAA9oJ1czMukF7Xz34FrAu8F3gBUnz0+sNSfNrH56ZmVnXtXfNrtLv4ZmZmfVYTmZmZlZ4TnZmZlZ4lT4b02ps0IYbc/Yfbqp3GGZmheSenZmZFZ6TnZmZFZ6TnZmZFZ6TnZmZFZ5vUOkh3p7zBo//6J/1DsMKZsuz926/ktlSwD07MzMrPCc7MzMrPCc7MzMrPCc7MzMrPCc7MzMrvF6T7CSdLWmqpCmSJknaWdKVkrZqoe6xki5t53jHSno5Hav02iptOz+d63xJa0qaIOkRSXt0MOZhkg7s2Ds1M7Nq6xVfPZC0K/BpslnS35E0EFg+Ir7cxUNfHxFfb6F8FLB6RCyWNBJo7uS5hgENwB1diNHMzLqot/TsBgFzI+IdgIiYGxEvSBonqQFA0nGSnpA0EdittGPqmf1J0kPptVvLp/ig/m1Af6BJ0pnAecAhqee3oqT9JD0g6WFJN0rqn/bbUdK/JU2WNFHSAOD7wBFp3yNq0TBmZta+XtGzA+4GzpH0BPB3sh7ZvaWNkgYB5wLDgdeBe4BH0uZfAr+IiH9JWh+4C9gybTtC0u658+waEQdLWhARw9KxXwQaIuLrqUf5XWDfiFiYkuE3Jf0UuB44IiIekrQK8CZwTmnflt6UpFFkvUgGDVirSw1kZmat6xXJLiIWSBoO7AGMAK6XdFauys7AuIh4GUDS9cBmadu+wFaSSnVXKfXGaH0YszW7AFsB49PxlgceADYH5kTEQyne+SmO9t7XGGAMwNaDN48OxGFmZh3QK5IdQEQsBsYB4yQ1A8dUuOsywC4R8Xa+sL1E1AoBf4uII8uOtU1nDmZmZt2jV1yzk7S5pE1zRcOAZ3PrE4A9Ja0hqQ9wWG7b3cBJuWMN60IoDwK7SdokHWslSZsB04FBknZM5StLWg54A1i5C+czM7Mq6BXJjuyGkWskPSZpCtlQ4ujSxoiYk9YfAMYDj+f2PRloSF9ZeAw4IbetdPNI6fXxtoJIw6THAn9McTwAbBER7wJHAJdImgz8DehLdu1wK9+gYmZWX4rwpaKeYOvBm8eNX7us3mFYwXjWAys6SU0R0dBevd7SszMzM+s0JzszMys8JzszMys8JzszMyu8XvM9u6LrO2hl30xgZlYj7tmZmVnhOdmZmVnhOdmZmVnhOdmZmVnh+QaVHuKFF15g9OjRdTl3vc5rZtZd3LMzM7PCc7IzM7PCc7IzM7PCc7IzM7PCc7IzM7PC61XJTtLiNBHqo5L+ImnVKh//O9U8npmZ9Qy9KtkBb0XEsIjYGngVOLHKx3eyMzMroN6W7PIeAAYDSNpY0l8lNUm6X9IWqfwzkiZIekTS3yWtncr7S/qdpGZJUyT9l6SfAiumnuPYVO+bqRf5qKRTUtkQSdMkXS3pCUljJe0rabykJyXtlOqtJOkqSRPT+Q+pQxuZmRm99EvlkpYF9gF+m4rGACdExJOSdgZ+DewN/AvYJSJC0peBM4BvAd8DXo+IbdLxVouIP0n6ekQMS2XDgeOAnQEBEyTdC8wDNgEOA44HHgL+G9gdOJisd3gocDbwz4g4Pg23TpT094hYWLuWMTOzlvS2ZLeipElkPbrHgb9J6g98HLhRUqneCunnesD1kgYBywMzUvm+wMhS5YiY18K5dgduKSUnSTcDewC3ATMiojmVTwX+kRJqMzAk7b8fcLCk09J6X2D9FDdp31HAKIABAwZ0tC3MzKxCvW0Y863U89qArLd1Itl7eC1dyyu9tkz1LwEuTT24r5AlnGp4J7f8fm79fT78A0LAf+ViWj8iHs8fJCLGRERDRDT069evSqGZmVm53pbsAIiIN4GTyYYk3wRmSDoMQJntUtUBwOy0fEzuEH8jd3OLpNXS4iJJfdLy/cChkvpJWgn4bCqr1F3ASUrdTUnbd2BfMzOrol6Z7AAi4hFgCnAkcBTwJUmTgalA6WaQ0WTDm03A3NzuPwRWSzeeTAZGpPIxwBRJYyPiYeBqYCIwAbgynbNSPwD6pONNTetmZlYHioh6x2DAuuuuG6NGjarLuT3rgZn1VpKaIqKhvXq9tmdnZmZWKSc7MzMrPCc7MzMrPCc7MzMrPN+g0kM0NDREY2NjvcMwM+tVfIOKmZlZ4mRnZmaF52RnZmaF52RnZmaF19tmPSisefMe54Ybd+rUvocfNrHK0ZiZFYt7dmZmVnhOdmZmVnhOdmZmVnhOdmZmVnhOdmZmVnhOdmZmVnhOdi2QFJL+kFtfTtLLkm7PlR0qaYqkxyU1Szq07BgDJS2SdEI3hm5mZi1wsmvZQmBrSSum9U8Cs0sbJW0HXAAcEhFbAgcDF0jaNneMw4AHgSO7J2QzM2uNk13r7gAOSstHAn/MbTsN+HFEzABIP38CnJ6rcyTwLWCwpPVqH66ZmbXGya511wEjJfUFtgUm5LYNBZrK6jemciR9DBgUEROBG4AjWjqBpFGSGiU1zp//XrXjNzOzxMmuFRExBRhC1kO7o4O7H0GW5CBLmi0OZUbEmIhoiIiGVVbxk9vMzGrFn7Btu43s2txewBq58seA4cDkXNlwYGpaPhJYR9JRaX1dSZtGxJO1DdfMzFrinl3brgLOjYjmsvILgG9LGgKQfn4HuFDSZkD/iBgcEUMiYgjZ9TzfqGJmVidOdm2IiFkRcXEL5ZOAM4G/SJoG/AU4I5UfCdxStsufcLIzM6sbD2O2ICL6t1A2DhiXW78ZuLmFeue2UDYF2LKqQZqZWcXcszMzs8JzsjMzs8JzsjMzs8LzNbseYrXVtuTwwybWOwwzs0Jyz87MzArPyc7MzApPEVHvGAyQ9AYwvd5xtGMgMLfeQbTB8XVNT48Pen6Mjq/rOhrjBhGxZnuVfM2u55geEQ31DqItkhp7coyOr2t6enzQ82N0fF1Xqxg9jGlmZoXnZGdmZoXnZNdzjKl3ABXo6TE6vq7p6fFBz4/R8XVdTWL0DSpmZlZ47tmZmVnhOdmZmVnhOdn1AJIOkDRd0lOSzurG835M0j2SHpM0VdI3Uvnqkv4m6cn0c7VULkkXpzinSNohd6xjUv0nJR1T5TiXlfSIpNvT+oaSJqQ4rpe0fCpfIa0/lbYPyR3j26l8uqT9qxjbqpJukjRN0uOSdu2B7Xdq+vd9VNIfJfWtZxtKukrSS5IezZVVrc0kDZfUnPa5WJKqEN/56d94iqRbJK2a29Ziu7T2e91a23c1xty2b0kKSQPTeo9ow1R+UmrHqZLOy5XXvg0jwq86voBlgaeBjYDlgcnAVt107kHADml5ZeAJYCvgPOCsVH4W8LO0fCBwJyBgF2BCKl8deCb9XC0tr1bFOL8JXAvcntZvAEam5cuBr6blrwGXp+WRwPVpeavUrisAG6b2XrZKsV0DfDktLw+s2pPaDxgMzABWzLXdsfVsQ+ATwA7Ao7myqrUZMDHVVdr3U1WIbz9gubT8s1x8LbYLbfxet9b2XY0xlX8MuAt4FhjYw9pwBPB3YIW0vlZ3tmHNP1D9avc/xa7AXbn1bwPfrlMsfwY+SfYkl0GpbBDZF94BrgCOzNWfnrYfCVyRK1+iXhdjWg/4B7A3cHv65Zub++D5oP3SL/muaXm5VE/lbZqv18XYBpAlEpWV96T2Gww8nz7QlkttuH+92xAYUvZBWJU2S9um5cqXqNfZ+Mq2fRYYm5ZbbBda+b1u6/9vNWIEbgK2A2byYbLrEW1IlqD2baFet7ShhzHrr/RhVDIrlXWrNFy1PTABWDsi5qRN/wHWTsutxVrL93ARcAbwflpfA3gtIt5r4VwfxJG2v57q1yq+DYGXgd8pG2a9UtJK9KD2i4jZwAXAc8AcsjZpoue0YUm12mxwWq5VnADHk/V2OhNfW/9/u0TSIcDsiJhctqmntOFmwB5p+PFeSTt2Mr5OtaGTnSGpP/An4JSImJ/fFtmfTnX5foqkTwMvRURTPc5fgeXIhmoui4jtgYVkQ3AfqGf7AaRrX4eQJeZ1gZWAA+oVTyXq3WZtkXQ28B4wtt6x5EnqB3wHOKfesbRhObIRhl2A04EbOnotsCuc7OpvNtk4e8l6qaxbSOpDlujGRsTNqfhFSYPS9kHAS+3EWqv3sBtwsKSZwHVkQ5m/BFaVVHqua/5cH8SRtg8AXqlhfLOAWRExIa3fRJb8ekr7AewLzIiIlyNiEXAzWbv2lDYsqVabzU7LVY9T0rHAp4GjUkLuTHyv0Hrbd8XGZH/QTE6/L+sBD0tapxMx1qoNZwE3R2Yi2WjNwE7E17k27MxYsV/Ve5H9tfMM2X/U0kXYod10bgG/By4qKz+fJW8WOC8tH8SSF7onpvLVya5drZZeM4DVqxzrXnx4g8qNLHlx+mtp+USWvLnihrQ8lCUvgD9D9W5QuR/YPC2PTm3XY9oP2BmYCvRL570GOKnebchHr+dUrc346M0VB1YhvgOAx4A1y+q12C608XvdWtt3NcaybTP58JpdT2nDE4Dvp+XNyIYo1V1tWLUPI786/yK7W+oJsjuPzu7G8+5ONlw0BZiUXgeSjYn/A3iS7O6p0i+AgF+lOJuBhtyxjgeeSq/jahDrXnyY7DZKv4xPpf/0pbu7+qb1p9L2jXL7n53ink4H7yxrJ65hQGNqw1vTh0aPaj/gXGAa8Cjwv+lDpW5tCPyR7PrhIrK/9r9UzTYDGtJ7fRq4lLIbiDoZ31NkH86l35PL22sXWvm9bq3tuxpj2faZfJjsekobLg/8IR33YWDv7mxDPy7MzMwKz9fszMys8JzszMys8JzszMys8JzszMys8JzszMys8JzszAoqPfn+wtz6aZJGV+nYV0v6fDWOZdYdnOzMiusd4HOlqV56ityTL8y6jZOdWXG9B4wBTi3fUN4zk7Qg/dwrPaT3z5KekfRTSUdJmpjmN9s4d5h9JTVKeiI9x7Q09+D5kh5Kc6d9JXfc+yXdRvYkErNu5b+wzIrtV8CU/ESZFdgO2BJ4lexxTVdGxE7KJvc9CTgl1RsC7ET2XMZ7JG0CfBF4PSJ2lLQCMF7S3an+DsDWETGji+/JrMOc7MwKLCLmS/o9cDLwVoW7PRRpuh1JTwOlZNVMNgFnyQ0R8T7wpKRngC3IJjndNtdrHABsCrxL9kxGJzqrCyc7s+K7iOxZhL/Llb1HuowhaRmy5xaWvJNbfj+3/j5LfmaUP2swyJ7DeFJE3JXfIGkvsimQzOrC1+zMCi4iXiWbJfpLueKZwPC0fDDQpxOHPkzSMuk63kZkD/G9C/hqmjoKSZulCW3N6srJzmzpcCHZ3GElvwH2lDQZ2JXO9bqeI3vy/J3ACRHxNnAl2Q0oD0t6FLgCjyBZD+BZD8zMrPDcszMzs8JzsjMzs8JzsjMzs8JzsjMzs8JzsjMzs8JzsjMzs8JzsjMzs8L7//Gh7/mbdZBJAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "nspace_data = {'Namespace':['ChEMBLAssay','ChEMBL','GOBP','GOMF','HGNC','Disease','SideEffect','Reactome','MOA'],\n",
    "       'Number':[15835,4096,2085,1332,1347,1585,7420,1513,739]}\n",
    "nspace = pd.DataFrame(nspace_data)\n",
    "a = sns.barplot(x=\"Number\", y=\"Namespace\", data=nspace_data)\n",
    "a.set(xlabel='Number',ylabel='Namespace',title= 'KG Namespace in numbers')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "id": "db746c14",
   "metadata": {},
   "outputs": [],
   "source": [
    "#to cytoscape compatible graphml \n",
    "pybel.to_graphml(itmp_graph,'data/export//bycovidgraph_2ndSep.graphml')\n",
    "\n",
    "#to regular BEL format\n",
    "pybel.dump(itmp_graph,'data/export//bycovidgraph_2ndSep.bel')\n",
    "\n",
    "#to neo4j\n",
    "pybel.to_csv(itmp_graph,'data/export//bycovidgraph_2ndSep.csv')\n",
    "\n",
    "#to sif\n",
    "pybel.to_sif(itmp_graph,'data/export//bycovidgraph_2ndSep.sif')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5d3b2344",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.10.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}