tools_demonstration.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "f1e2ec0f",
   "metadata": {},
   "outputs": [],
   "source": [
    "# %%capture\n",
    "# %run Preparation_for_future_project.ipynb import *\n",
    "\n",
    "from lookup_gene import *\n",
    "from differential_expression import *\n",
    "from find_targets import *\n",
    "from plotting_functions import *"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "42a2b40a",
   "metadata": {},
   "source": [
    "## differential_expression tools"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "df3d86b4",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.07144184318097363"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "diff_expression(gene='Cd44', cluster1='6_TREAT', cluster2='0_TREAT')\n",
    "# Finds difference in proportion of cells that express gene between two clusters"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "07ef5801",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The most differentially expressed genes are:\n",
      "\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "{'Mki67': 0.9637681159420289,\n",
       " 'Smc2': 0.9396135265700483,\n",
       " 'Top2a': 0.8952991452991453,\n",
       " 'Phgdh': 0.8292456335934597,\n",
       " 'Birc5': 0.761984392419175,\n",
       " 'Fasn': 0.7511148272017838,\n",
       " 'Cdk1': 0.7511148272017838,\n",
       " 'Cenpa': 0.7486993682645857,\n",
       " 'Dkc1': 0.7076365663322185,\n",
       " 'Nasp': 0.7012263099219621,\n",
       " 'Hmmr': 0.6923076923076923,\n",
       " 'Ryk': 0.6898922333704942,\n",
       " 'Hmgb3': 0.6874767744332961,\n",
       " 'Cenpf': 0.6862690449646971,\n",
       " 'Trmt10c': 0.6798587885544408,\n",
       " 'Cdc20': 0.679022668153103,\n",
       " 'Uck2': 0.6693608324043107,\n",
       " 'Mad2l1': 0.6621144555927164,\n",
       " 'Mpp6': 0.6569119286510592,\n",
       " 'Ptprs': 0.6524526198439242,\n",
       " 'Dut': 0.6512448903753252,\n",
       " 'Gart': 0.650873281308064,\n",
       " 'Tpx2': 0.646413972500929,\n",
       " 'Ide': 0.6424191750278707,\n",
       " 'Mrps6': 0.6387959866220736,\n",
       " 'Rrs1': 0.6351727982162765,\n",
       " 'Rrm1': 0.6307134894091415,\n",
       " 'Stmn1': 0.6270903010033445,\n",
       " 'Tex30': 0.6258825715347455,\n",
       " 'Ttc39c': 0.6218877740616872,\n",
       " 'Dtymk': 0.617056856187291,\n",
       " 'Aspm': 0.6141768859160164,\n",
       " 'Ccnb1': 0.6117614269788183,\n",
       " 'Ube2c': 0.6093459680416202,\n",
       " 'Naa15': 0.6086027499070978,\n",
       " 'Ncbp1': 0.6077666295057599,\n",
       " 'Nedd4': 0.6069305091044221,\n",
       " 'Ncapd2': 0.6057227796358231,\n",
       " 'Rad21': 0.6046079524340394,\n",
       " 'Pold2': 0.6041434410999629,\n",
       " 'Rcn1': 0.6020995912300261,\n",
       " 'Tcf3': 0.5953177257525084,\n",
       " 'Hells': 0.5948532144184319,\n",
       " 'Alad': 0.5944816053511706,\n",
       " 'Tpm1': 0.5908584169453734,\n",
       " 'Anp32e': 0.5904868078781123,\n",
       " 'Zc3h14': 0.5880713489409142,\n",
       " 'Cad': 0.5876068376068376,\n",
       " 'Mcm7': 0.5863991081382386,\n",
       " 'Cks2': 0.5856558900037161}"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "diff_exp_treat_6_0 = most_diff_exp(cluster1='6_TREAT', cluster2='0_TREAT', n_keep=50, df=data_treat)\n",
    "# Finds most differentially expressed genes between two clusters\n",
    "\n",
    "# 'cluster1' and 'cluster2' are the two clusters in which to find the most differentially\n",
    "# expressed genes\n",
    "# 'n_keep' specifies that you want to keep the top 'n' most differentially expressed genes\n",
    "# 'df' specifies the dataset you want to use - the expression of genes in cell clusters\n",
    "# differs in data that comes from different experiments, people, datasets, etc.\n",
    "diff_exp_treat_6_0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "3aaceb86",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "The list of genes that are differentialy expressed between cluster 6_TREAT and cluster 0_TREAT, that are not expressed in more than 0.5 of any other cluster, are:\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "{'Ptprs': '6_TREAT'}"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "compare_to_rest(cluster1='6_TREAT', cluster2='0_TREAT', thresh=0.50, n_keep=50)\n",
    "#     Returns dict - keys are genes that are differentially expressed between cluster1\n",
    "#         and cluster2, and are not expressed greater than 'thresh' in any other cluster,\n",
    "#         and keys are expression levels"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "03c2ad0b",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'Mki67': {'CTRL Median': 4.0,\n",
       "  'Median effect': -2.0,\n",
       "  'Pct of cells effect': 0.09399224806201545},\n",
       " 'Smc2': {'CTRL Median': 3.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': 0.11627906976744184},\n",
       " 'Top2a': {'CTRL Median': 7.0,\n",
       "  'Median effect': -5.5,\n",
       "  'Pct of cells effect': -0.02847346451997612},\n",
       " 'Phgdh': {'CTRL Median': 6.0,\n",
       "  'Median effect': -2.0,\n",
       "  'Pct of cells effect': -0.12380739415623132},\n",
       " 'Birc5': {'CTRL Median': 5.0,\n",
       "  'Median effect': -0.5,\n",
       "  'Pct of cells effect': -0.18135062611806796},\n",
       " 'Fasn': {'CTRL Median': 2.0,\n",
       "  'Median effect': 0.0,\n",
       "  'Pct of cells effect': 0.08512224209898633},\n",
       " 'Cdk1': {'CTRL Median': 4.0,\n",
       "  'Median effect': -2.0,\n",
       "  'Pct of cells effect': -0.12805605247465712},\n",
       " 'Cenpa': {'CTRL Median': 6.0,\n",
       "  'Median effect': -3.0,\n",
       "  'Pct of cells effect': -0.13774597495527718},\n",
       " 'Dkc1': {'CTRL Median': 2.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': 0.07252534287418011},\n",
       " 'Nasp': {'CTRL Median': 3.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': -0.04144305307096008},\n",
       " 'Hmmr': {'CTRL Median': 2.0,\n",
       "  'Median effect': 0.0,\n",
       "  'Pct of cells effect': 0.13901311866428145},\n",
       " 'Ryk': {'CTRL Median': 2.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': -0.0799791293977341},\n",
       " 'Hmgb3': {'CTRL Median': 2.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': -0.14683959451401318},\n",
       " 'Cenpf': {'CTRL Median': 3.0,\n",
       "  'Median effect': -2.0,\n",
       "  'Pct of cells effect': -0.11680083482409065},\n",
       " 'Trmt10c': {'CTRL Median': 1.0,\n",
       "  'Median effect': 0.0,\n",
       "  'Pct of cells effect': 0.15004472271914138},\n",
       " 'Cdc20': {'CTRL Median': 2.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': 0.04017590936195581},\n",
       " 'Uck2': {'CTRL Median': 3.0,\n",
       "  'Median effect': 0.0,\n",
       "  'Pct of cells effect': -0.18366129994036973},\n",
       " 'Mad2l1': {'CTRL Median': 2.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': -0.10032796660703636},\n",
       " 'Mpp6': {'CTRL Median': 2.0,\n",
       "  'Median effect': 0.0,\n",
       "  'Pct of cells effect': 0.0017889087656529634},\n",
       " 'Ptprs': {'CTRL Median': 1.0,\n",
       "  'Median effect': 1.0,\n",
       "  'Pct of cells effect': 0.24754025044722716},\n",
       " 'Dut': {'CTRL Median': 3.0,\n",
       "  'Median effect': 1.0,\n",
       "  'Pct of cells effect': -0.16137447823494333},\n",
       " 'Gart': {'CTRL Median': 2.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': 0.10450208706022668},\n",
       " 'Tpx2': {'CTRL Median': 3.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': -0.16234347048300535},\n",
       " 'Ide': {'CTRL Median': 2.0,\n",
       "  'Median effect': -0.5,\n",
       "  'Pct of cells effect': 0.17911449016100178},\n",
       " 'Mrps6': {'CTRL Median': 6.0,\n",
       "  'Median effect': -4.0,\n",
       "  'Pct of cells effect': -0.18038163387000594},\n",
       " 'Rrs1': {'CTRL Median': 2.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': -0.019528920691711305},\n",
       " 'Rrm1': {'CTRL Median': 3.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': -0.1545915324985092},\n",
       " 'Stmn1': {'CTRL Median': 13.0,\n",
       "  'Median effect': -10.0,\n",
       "  'Pct of cells effect': -0.2970333929636255},\n",
       " 'Tex30': {'CTRL Median': 2.0,\n",
       "  'Median effect': 0.0,\n",
       "  'Pct of cells effect': -0.01990161001788915},\n",
       " 'Ttc39c': {'CTRL Median': 1.0,\n",
       "  'Median effect': 0.0,\n",
       "  'Pct of cells effect': 0.505664877757901},\n",
       " 'Dtymk': {'CTRL Median': 3.0,\n",
       "  'Median effect': -1.5,\n",
       "  'Pct of cells effect': -0.0989862850327966},\n",
       " 'Aspm': {'CTRL Median': 1.0,\n",
       "  'Median effect': 0.0,\n",
       "  'Pct of cells effect': 0.15123732856291},\n",
       " 'Ccnb1': {'CTRL Median': 3.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': -0.1762820512820512},\n",
       " 'Ube2c': {'CTRL Median': 10.0,\n",
       "  'Median effect': -4.5,\n",
       "  'Pct of cells effect': -0.3216308884913536},\n",
       " 'Naa15': {'CTRL Median': 2.0,\n",
       "  'Median effect': 0.0,\n",
       "  'Pct of cells effect': 0.15488968395945146},\n",
       " 'Ncbp1': {'CTRL Median': 1.0,\n",
       "  'Median effect': 2.0,\n",
       "  'Pct of cells effect': 0.038237924865831774},\n",
       " 'Nedd4': {'CTRL Median': 4.0,\n",
       "  'Median effect': -1.5,\n",
       "  'Pct of cells effect': -0.31969290399522954},\n",
       " 'Ncapd2': {'CTRL Median': 2.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': -0.14430530709600475},\n",
       " 'Rad21': {'CTRL Median': 3.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': -0.04435002981514613},\n",
       " 'Pold2': {'CTRL Median': 2.0,\n",
       "  'Median effect': 0.0,\n",
       "  'Pct of cells effect': 0.047927847346451946},\n",
       " 'Rcn1': {'CTRL Median': 2.0,\n",
       "  'Median effect': -0.5,\n",
       "  'Pct of cells effect': -0.22957662492546205},\n",
       " 'Tcf3': {'CTRL Median': 4.0,\n",
       "  'Median effect': -2.0,\n",
       "  'Pct of cells effect': -0.183288610614192},\n",
       " 'Hells': {'CTRL Median': 2.0,\n",
       "  'Median effect': 0.0,\n",
       "  'Pct of cells effect': -0.17337507453786516},\n",
       " 'Alad': {'CTRL Median': 1.0,\n",
       "  'Median effect': 1.0,\n",
       "  'Pct of cells effect': 0.14288908765652952},\n",
       " 'Tpm1': {'CTRL Median': 11.0,\n",
       "  'Median effect': -9.0,\n",
       "  'Pct of cells effect': -0.29897137745974955},\n",
       " 'Anp32e': {'CTRL Median': 5.0,\n",
       "  'Median effect': -2.5,\n",
       "  'Pct of cells effect': -0.183288610614192},\n",
       " 'Zc3h14': {'CTRL Median': 2.0,\n",
       "  'Median effect': -1.0,\n",
       "  'Pct of cells effect': -0.057319618366129976},\n",
       " 'Cad': {'CTRL Median': 1.0,\n",
       "  'Median effect': 1.0,\n",
       "  'Pct of cells effect': 0.04949314251639836},\n",
       " 'Mcm7': {'CTRL Median': 2.0,\n",
       "  'Median effect': 0.0,\n",
       "  'Pct of cells effect': -0.19953786523553962},\n",
       " 'Cks2': {'CTRL Median': 7.0,\n",
       "  'Median effect': -2.5,\n",
       "  'Pct of cells effect': -0.19976147883124618}}"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "compare_treat_ctrl(diff_exp_treat_6_0, cluster_num=6)\n",
    "#     Compares expression of genes in treat and control of same cluster\n",
    "#         For each gene, returns several statistical features of the \n",
    "#         change in that gene's distribution"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b8aeb037",
   "metadata": {},
   "source": [
    "## find_target tools"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "b9d4e7c8",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['Mrpl15',\n",
       " 'Lypla1',\n",
       " 'Tcea1',\n",
       " 'Rgs20',\n",
       " 'Atp6v1h',\n",
       " 'Rb1cc1',\n",
       " '4732440D04Rik',\n",
       " 'Pcmtd1',\n",
       " 'Rrs1',\n",
       " '2610203C22Rik',\n",
       " 'Mybl1',\n",
       " 'Vcpip1',\n",
       " 'Sgk3',\n",
       " 'Snhg6',\n",
       " 'Cops5',\n",
       " 'Cspp1',\n",
       " 'Arfgef1',\n",
       " 'Slco5a1',\n",
       " 'Ncoa2',\n",
       " 'Tram1',\n",
       " 'Lactb2',\n",
       " 'Msc',\n",
       " 'Terf1',\n",
       " 'Rpl7',\n",
       " 'Rdh10',\n",
       " 'Stau2',\n",
       " 'Ube2w',\n",
       " 'Eloc',\n",
       " 'Tmem70',\n",
       " 'Ly96',\n",
       " 'Mcm3',\n",
       " 'Tram2',\n",
       " 'Tmem14a',\n",
       " 'Kcnq5',\n",
       " 'Ogfrl1',\n",
       " 'B3gat2',\n",
       " 'Smap1',\n",
       " 'Sdhaf4',\n",
       " 'Fam135a',\n",
       " 'Lmbrd1',\n",
       " 'Adgrb3',\n",
       " 'Phf3',\n",
       " 'Ptp4a1.1',\n",
       " 'Prim2',\n",
       " 'Rab23',\n",
       " 'Zfp451',\n",
       " 'Dst',\n",
       " 'Gm37233',\n",
       " 'Ccdc115',\n",
       " 'Imp4',\n",
       " 'Ptpn18',\n",
       " 'Gm28417',\n",
       " 'Fam168b',\n",
       " 'Plekhb2',\n",
       " 'Hs6st1',\n",
       " 'Uggt1',\n",
       " 'Neurl3',\n",
       " 'Arid5a',\n",
       " 'Kansl3',\n",
       " 'Lman2l',\n",
       " 'Cnnm4',\n",
       " 'Cnnm3',\n",
       " 'Ankrd23',\n",
       " 'Ankrd39',\n",
       " 'Sema4c',\n",
       " 'D430040D24Rik',\n",
       " 'Cox5b',\n",
       " 'Actr1b',\n",
       " 'Zap70',\n",
       " 'Tmem131',\n",
       " 'Gm38115',\n",
       " 'Inpp4a',\n",
       " 'Coa5',\n",
       " 'Unc50',\n",
       " 'Mgat4a',\n",
       " 'Tsga10',\n",
       " 'Lipt1',\n",
       " 'Mitd1',\n",
       " 'Mrpl30',\n",
       " 'Txndc9',\n",
       " 'Eif5b',\n",
       " 'Rev1',\n",
       " 'Chst10',\n",
       " 'Pdcl3',\n",
       " 'Rpl31',\n",
       " 'Cnot11',\n",
       " 'Rnf149',\n",
       " 'Map4k4',\n",
       " 'Il1r2',\n",
       " 'Il1r1',\n",
       " 'Il1rl1',\n",
       " 'Il18r1',\n",
       " 'Il18rap',\n",
       " 'Gm37915',\n",
       " 'Pou3f3',\n",
       " 'Mrps9',\n",
       " 'Tgfbrap1',\n",
       " 'AI597479',\n",
       " 'Fhl2',\n",
       " 'Nck2',\n",
       " 'Uxs1',\n",
       " 'Tpp2',\n",
       " 'Gm8251',\n",
       " 'Tex30',\n",
       " 'Poglut2',\n",
       " 'Bivm',\n",
       " 'Ercc5',\n",
       " 'Gulp1',\n",
       " 'Col3a1',\n",
       " 'Col5a2',\n",
       " 'Wdr75',\n",
       " 'Dnah7b',\n",
       " 'Slc39a10',\n",
       " 'Cavin2',\n",
       " 'Nabp1',\n",
       " 'Stat4',\n",
       " 'Stat1',\n",
       " 'Gls',\n",
       " 'Nab1',\n",
       " 'Nemp2',\n",
       " 'Mfsd6',\n",
       " 'Inpp1',\n",
       " 'Hibch',\n",
       " 'Pms1',\n",
       " 'Ormdl1',\n",
       " 'Osgepl1',\n",
       " 'Asnsd1',\n",
       " 'Stk17b',\n",
       " 'Gtf3c3',\n",
       " 'Ankrd44',\n",
       " 'Sf3b1',\n",
       " 'Coq10b',\n",
       " 'Hspd1',\n",
       " 'Hspe1',\n",
       " 'Mob4',\n",
       " 'Rftn2',\n",
       " 'Mars2',\n",
       " 'Satb2',\n",
       " '1700066M21Rik',\n",
       " 'Tyw5',\n",
       " 'Maip1',\n",
       " 'Kctd18',\n",
       " 'Sgo2a',\n",
       " 'Bzw1',\n",
       " 'Clk1',\n",
       " 'Ppil3',\n",
       " 'Nif3l1',\n",
       " 'Orc2',\n",
       " 'Gm15834',\n",
       " 'Fam126b',\n",
       " 'Ndufb3',\n",
       " 'Cflar',\n",
       " 'Casp8',\n",
       " 'Trak2',\n",
       " 'Tmem237',\n",
       " 'Als2',\n",
       " 'Fzd7',\n",
       " 'Sumo1',\n",
       " 'Nop58',\n",
       " 'Bmpr2',\n",
       " 'Fam117b',\n",
       " 'Ica1l',\n",
       " 'Wdr12',\n",
       " 'Carf',\n",
       " 'Nbeal1',\n",
       " 'Cyp20a1',\n",
       " 'Abi2',\n",
       " 'Raph1',\n",
       " 'Cd28',\n",
       " 'Ctla4',\n",
       " 'Icos',\n",
       " 'Nrp2',\n",
       " 'Ino80d',\n",
       " 'Ino80dos',\n",
       " 'Ndufs1',\n",
       " 'Eef1b2',\n",
       " 'Fastkd2',\n",
       " 'Klf7',\n",
       " 'Creb1',\n",
       " 'Mettl21a',\n",
       " '2810408I11Rik',\n",
       " 'Ccnyl1',\n",
       " 'Fzd5',\n",
       " 'Plekhm3',\n",
       " 'Idh1',\n",
       " 'Pikfyve',\n",
       " 'Map2',\n",
       " 'Unc80',\n",
       " 'Rpe',\n",
       " 'Kansl1l',\n",
       " 'Acadl',\n",
       " 'Lancl1',\n",
       " 'Ikzf2',\n",
       " 'Bard1',\n",
       " 'Atic',\n",
       " 'Fn1',\n",
       " 'Pecr',\n",
       " 'Xrcc5',\n",
       " 'Smarcal1',\n",
       " 'Rpl37a',\n",
       " 'Igfbp2',\n",
       " 'Arpc2',\n",
       " 'Aamp',\n",
       " 'Pnkd',\n",
       " 'Tmbim1',\n",
       " 'Slc11a1',\n",
       " 'Ctdsp1',\n",
       " 'Usp37',\n",
       " 'Cnot9',\n",
       " 'Zfp142',\n",
       " 'Bcs1l',\n",
       " 'Rnf25',\n",
       " 'Stk36',\n",
       " 'Ttll4',\n",
       " 'Wnt6',\n",
       " 'Cnppd1',\n",
       " 'Retreg2',\n",
       " 'Zfand2b',\n",
       " 'Atg9a',\n",
       " 'Ankzf1',\n",
       " 'Stk16',\n",
       " 'Tuba4a',\n",
       " 'Ptprn',\n",
       " 'Dnpep',\n",
       " 'Gmppa',\n",
       " 'Chpf',\n",
       " 'Obsl1',\n",
       " 'Stk11ip',\n",
       " 'Farsb',\n",
       " 'Acsl3',\n",
       " 'Ap1s3',\n",
       " 'Wdfy1',\n",
       " 'Mrpl44',\n",
       " 'Serpine2',\n",
       " 'Fam124b',\n",
       " 'Cul3',\n",
       " 'Dock10',\n",
       " 'Rhbdd1',\n",
       " 'Mff',\n",
       " 'Agfg1',\n",
       " 'Pid1',\n",
       " 'Trip12',\n",
       " 'A530032D15Rik',\n",
       " 'A530040E14Rik',\n",
       " 'Sp110',\n",
       " 'Sp140',\n",
       " 'Sp100',\n",
       " 'A630001G21Rik',\n",
       " 'Cab39',\n",
       " 'Itm2c',\n",
       " 'Gpr55',\n",
       " 'Psmd1',\n",
       " 'Armc9',\n",
       " 'Ncl',\n",
       " 'Ptma',\n",
       " 'Pde6d',\n",
       " 'Cops7b',\n",
       " 'Dis3l2',\n",
       " 'Eif4e2',\n",
       " 'Gigyf2',\n",
       " 'Inpp5d',\n",
       " 'Atg16l1',\n",
       " 'Dgkd',\n",
       " 'Usp40',\n",
       " 'Ugt1a7c',\n",
       " 'Mroh2a',\n",
       " 'Hjurp',\n",
       " 'Arl4c',\n",
       " 'Sh3bp4',\n",
       " 'Agap1',\n",
       " 'Ackr3',\n",
       " 'Cops8',\n",
       " 'Lrrfip1',\n",
       " 'Ramp1',\n",
       " 'Ube2f',\n",
       " 'Scly',\n",
       " 'Ilkap',\n",
       " 'Hes6',\n",
       " 'Per2',\n",
       " 'Traf3ip1',\n",
       " 'Twist2',\n",
       " 'Hdac4',\n",
       " 'Ndufa10',\n",
       " 'Cops9',\n",
       " 'Gpc1',\n",
       " 'Dusp28',\n",
       " 'Rnpepl1',\n",
       " '9430060I03Rik',\n",
       " 'Capn10',\n",
       " 'Gpr35',\n",
       " 'Kif1a',\n",
       " 'Mterf4',\n",
       " 'Pask',\n",
       " 'Ppp1r7',\n",
       " 'Hdlbp',\n",
       " 'Sept2',\n",
       " 'Stk25',\n",
       " 'Bok',\n",
       " 'Thap4',\n",
       " 'Atg4b',\n",
       " 'Dtymk',\n",
       " 'Ing5',\n",
       " 'D2hgdh',\n",
       " 'Pdcd1',\n",
       " 'Fam174a',\n",
       " 'St8sia4',\n",
       " 'D1Ertd622e',\n",
       " 'Ppip5k2',\n",
       " 'Gin1',\n",
       " 'Pam',\n",
       " 'B230216N24Rik',\n",
       " 'Cdh20',\n",
       " 'Pign',\n",
       " 'Relch',\n",
       " 'Gm7160',\n",
       " 'Zcchc2',\n",
       " 'Phlpp1',\n",
       " 'Bcl2',\n",
       " 'Kdsr',\n",
       " 'Vps4b',\n",
       " 'Serpinb8',\n",
       " 'Dsel',\n",
       " 'Tsn',\n",
       " 'Nifk',\n",
       " 'Clasp1',\n",
       " 'Ralb',\n",
       " 'Tmem185b',\n",
       " 'Epb41l5',\n",
       " 'Ptpn4',\n",
       " 'Tmem177',\n",
       " 'Tmem37',\n",
       " 'Dbi',\n",
       " '3110009E18Rik',\n",
       " 'En1',\n",
       " 'Insig2',\n",
       " 'Ccdc93',\n",
       " 'Ddx18',\n",
       " 'Actr3',\n",
       " 'Slc35f5',\n",
       " 'Lypd1',\n",
       " 'Nckap5',\n",
       " 'Mgat5',\n",
       " 'Tmem163',\n",
       " '2900009J06Rik',\n",
       " 'Ccnt2',\n",
       " 'Rab3gap1',\n",
       " 'Zranb3',\n",
       " 'R3hdm1',\n",
       " 'Ubxn4',\n",
       " 'Mcm6',\n",
       " 'Dars',\n",
       " 'Pfkfb2',\n",
       " 'Yod1',\n",
       " 'Fcmr',\n",
       " 'Il10',\n",
       " 'Mapkapk2',\n",
       " 'Dyrk3',\n",
       " 'Eif2d',\n",
       " 'Rassf5',\n",
       " 'Ikbke',\n",
       " 'Srgap2',\n",
       " 'Fam72a',\n",
       " 'Rab7b',\n",
       " 'Slc41a1',\n",
       " 'Rab29',\n",
       " 'Nucks1',\n",
       " 'Elk4',\n",
       " 'Klhdc8a',\n",
       " 'Nuak2',\n",
       " 'Tmcc2',\n",
       " 'Dstyk',\n",
       " 'Rbbp5',\n",
       " 'Nfasc',\n",
       " 'Mdm4',\n",
       " 'Ppp1r15b',\n",
       " 'Snrpe',\n",
       " 'Gm38394',\n",
       " 'Lax1',\n",
       " 'Atp2b4',\n",
       " 'Btg2',\n",
       " 'Tmem183a',\n",
       " 'Cyb5r1',\n",
       " 'Adipor1',\n",
       " 'Klhl12',\n",
       " 'Rabif',\n",
       " 'Kdm5b',\n",
       " 'Ppp1r12b',\n",
       " 'Ube2t',\n",
       " 'Lgr6',\n",
       " 'Ptpn7',\n",
       " 'Arl8a',\n",
       " 'Gpr37l1',\n",
       " 'Rnpep',\n",
       " 'Timm17a',\n",
       " 'Shisa4',\n",
       " 'Ipo9',\n",
       " 'Nav1',\n",
       " 'Csrp1',\n",
       " 'Phlda3',\n",
       " 'Tnnt2',\n",
       " 'Tmem9',\n",
       " 'Kif21b',\n",
       " 'Inava',\n",
       " 'Gpr25',\n",
       " 'Camsap2',\n",
       " 'Ddx59',\n",
       " 'Kif14',\n",
       " 'Zfp281',\n",
       " 'Gm19705',\n",
       " 'Gm4258',\n",
       " 'Ptprc',\n",
       " 'Nek7',\n",
       " '2310009B15Rik',\n",
       " 'Dennd1b',\n",
       " 'Zbtb41',\n",
       " 'Aspm',\n",
       " 'Cdc73',\n",
       " 'Glrx2',\n",
       " 'Ro60',\n",
       " 'Uchl5',\n",
       " 'Gm29170',\n",
       " 'Rgs2',\n",
       " 'Rgs1',\n",
       " 'Pla2g4a',\n",
       " 'Ptgs2',\n",
       " 'Ptgs2os2',\n",
       " 'Odr4',\n",
       " 'Tpr',\n",
       " 'Ivns1abp',\n",
       " 'Swt1',\n",
       " 'Trmt1l',\n",
       " 'Rnf2',\n",
       " 'Fam129a',\n",
       " 'Edem3',\n",
       " '1700025G04Rik',\n",
       " 'Tsen15',\n",
       " 'Rgl1',\n",
       " 'Arpc5',\n",
       " 'Ncf2',\n",
       " 'Smg7',\n",
       " 'Gm28791',\n",
       " 'Nmnat2',\n",
       " 'Lamc2',\n",
       " 'Lamc1',\n",
       " 'Dhx9',\n",
       " 'Rgs16',\n",
       " 'Glul',\n",
       " 'Ier5',\n",
       " 'Mr1',\n",
       " 'Stx6',\n",
       " 'Xpr1',\n",
       " 'Acbd6',\n",
       " 'Qsox1',\n",
       " 'Cep350',\n",
       " 'Tor1aip1',\n",
       " 'Tor1aip2',\n",
       " 'Gm2000',\n",
       " 'Soat1',\n",
       " 'Abl2',\n",
       " 'Tor3a',\n",
       " 'Fam20b',\n",
       " 'Ralgps2',\n",
       " 'Rasal2',\n",
       " 'Cryzl2',\n",
       " 'Cop1',\n",
       " 'Tnr',\n",
       " '4930523C07Rik',\n",
       " 'Mrps14',\n",
       " 'Cacybp',\n",
       " 'Rabgap1l',\n",
       " 'Rc3h1',\n",
       " 'Zbtb37',\n",
       " 'Gas5',\n",
       " 'Dars2',\n",
       " 'Cenpl',\n",
       " 'Klhl20',\n",
       " 'Prdx6',\n",
       " 'Fasl',\n",
       " 'Suco',\n",
       " 'Pigc',\n",
       " 'Dnm3',\n",
       " 'Dnm3os',\n",
       " 'Eef1aknmt',\n",
       " 'Vamp4',\n",
       " 'Prrc2c',\n",
       " 'Prrx1',\n",
       " 'Gorab',\n",
       " 'Scyl3',\n",
       " 'Sell',\n",
       " 'Slc19a2',\n",
       " 'Blzf1',\n",
       " 'Sft2d2',\n",
       " 'Tiprl',\n",
       " 'Gm38190',\n",
       " 'Mpc2',\n",
       " 'Mpzl1',\n",
       " 'Rcsd1',\n",
       " 'Creg1',\n",
       " 'Gm16565',\n",
       " 'Cd247',\n",
       " 'Pou2f1',\n",
       " 'Gm15853',\n",
       " 'Tada1',\n",
       " 'Pogk',\n",
       " 'Uck2',\n",
       " 'Tmco1',\n",
       " 'Aldh9a1',\n",
       " 'Mgst3',\n",
       " 'Pbx1',\n",
       " 'Nuf2',\n",
       " 'Hsd17b7',\n",
       " 'Ddr2',\n",
       " 'Uap1',\n",
       " 'Uhmk1',\n",
       " 'Olfml2b',\n",
       " 'Atf6',\n",
       " 'Dusp12',\n",
       " 'Fcgr2b',\n",
       " 'Fcgr4',\n",
       " 'Fcgr3',\n",
       " 'Sdhc',\n",
       " 'Tomm40l',\n",
       " 'Fcer1g',\n",
       " 'Ndufs2',\n",
       " 'Adamts4',\n",
       " 'B4galt3',\n",
       " 'Ppox',\n",
       " 'Usp21',\n",
       " 'Ufc1',\n",
       " 'Gm20045',\n",
       " 'Dedd',\n",
       " 'Nit1',\n",
       " 'Pfdn2',\n",
       " 'Nectin4',\n",
       " 'Arhgap30',\n",
       " 'Usf1',\n",
       " 'B930036N10Rik',\n",
       " 'Ly9',\n",
       " 'Slamf7',\n",
       " 'Cd48',\n",
       " 'Slamf1',\n",
       " 'Cd84',\n",
       " 'Vangl2',\n",
       " 'Ncstn',\n",
       " 'Copa',\n",
       " 'Pex19',\n",
       " 'Dcaf8',\n",
       " 'Pea15a',\n",
       " 'Igsf8',\n",
       " 'Atp1a2',\n",
       " 'Pigm',\n",
       " 'Slamf9',\n",
       " 'Igsf9',\n",
       " 'Tagln2',\n",
       " 'Slamf8',\n",
       " 'Aim2',\n",
       " 'Ifi206',\n",
       " 'Ifi214',\n",
       " 'Ifi213',\n",
       " 'Ifi209',\n",
       " 'Ifi208',\n",
       " 'Ifi207',\n",
       " 'Ifi204',\n",
       " 'Mndal',\n",
       " 'Ifi211',\n",
       " 'Ifi203',\n",
       " 'Ifi205',\n",
       " 'Rgs7',\n",
       " 'Fh1',\n",
       " 'Kmo',\n",
       " 'Chml',\n",
       " 'Exo1',\n",
       " 'Cep170',\n",
       " 'Sdccag8',\n",
       " 'Akt3',\n",
       " 'Zbtb18',\n",
       " 'Adss',\n",
       " 'Desi2',\n",
       " 'Cox20',\n",
       " 'Hnrnpu',\n",
       " 'Kif26b',\n",
       " 'Smyd3',\n",
       " 'Tfb2m',\n",
       " 'Sccpdh',\n",
       " 'Ahctf1',\n",
       " 'Cdc42bpa',\n",
       " 'Psen2',\n",
       " 'Itpkb',\n",
       " 'Parp1',\n",
       " 'Acbd3',\n",
       " 'H3f3a',\n",
       " 'Sde2',\n",
       " 'Pycr2',\n",
       " 'Nvl',\n",
       " 'Cnih4',\n",
       " 'Wdr26',\n",
       " 'Lbr',\n",
       " 'Enah',\n",
       " 'Srp9',\n",
       " 'Degs1',\n",
       " 'Fbxo28',\n",
       " 'Trp53bp2',\n",
       " 'Capn2',\n",
       " 'Brox',\n",
       " 'Aida',\n",
       " 'Mia3',\n",
       " 'Taf1a',\n",
       " 'Dusp10',\n",
       " '1700056E22Rik',\n",
       " 'Gm37800',\n",
       " 'Hlx',\n",
       " 'Marc2',\n",
       " 'Rab3gap2',\n",
       " 'Iars2',\n",
       " 'Bpnt1',\n",
       " 'Eprs',\n",
       " 'Rrp15',\n",
       " 'Gpatch2',\n",
       " 'Kctd3',\n",
       " 'Kcnk2',\n",
       " 'Cenpf',\n",
       " 'Ptpn14',\n",
       " 'Smyd2',\n",
       " 'Rps6kc1',\n",
       " 'Angel2',\n",
       " 'Vash2',\n",
       " 'Flvcr1',\n",
       " 'Tatdn3',\n",
       " 'Nsl1',\n",
       " 'Batf3',\n",
       " 'Atf3',\n",
       " 'Nenf',\n",
       " 'Pacc1',\n",
       " 'Ppp2r5a',\n",
       " 'Dtl',\n",
       " 'Ints7',\n",
       " 'Lpgat1',\n",
       " 'Nek2',\n",
       " 'Slc30a1',\n",
       " 'Gm26670',\n",
       " 'Traf5',\n",
       " 'Rcor3',\n",
       " 'Hhat',\n",
       " 'Utp25',\n",
       " 'Traf3ip3',\n",
       " 'Lamb3',\n",
       " 'Hsd11b1',\n",
       " 'Cr1l',\n",
       " 'Fam171a1',\n",
       " 'Nmt2',\n",
       " 'Rpp38',\n",
       " 'Acbd7',\n",
       " 'Dclre1c',\n",
       " 'Suv39h2',\n",
       " 'Hspa14',\n",
       " 'Gm45902',\n",
       " 'Fam107b',\n",
       " 'Frmd4a',\n",
       " 'Gm36988',\n",
       " 'Prpf18',\n",
       " 'Bend7',\n",
       " 'Sephs1',\n",
       " 'Phyh',\n",
       " 'Mcm10',\n",
       " 'Camk1d',\n",
       " 'Cdc123',\n",
       " 'Nudt5',\n",
       " 'Sec61a2',\n",
       " 'Upf2',\n",
       " 'Echdc3',\n",
       " 'Usp6nl',\n",
       " 'Celf2',\n",
       " 'Gata3',\n",
       " 'Taf3',\n",
       " 'Atp5c1',\n",
       " 'Kin',\n",
       " 'Prkcq',\n",
       " 'Pfkfb3',\n",
       " 'Rbm17',\n",
       " 'Il2ra',\n",
       " 'Il15ra',\n",
       " 'Fbh1',\n",
       " 'Ankrd16',\n",
       " 'Mindy3',\n",
       " 'C1ql3',\n",
       " 'Rsu1',\n",
       " 'Gm13270',\n",
       " 'Vim',\n",
       " 'Hacd1',\n",
       " 'Stam',\n",
       " 'Nsun6',\n",
       " 'Arl5b',\n",
       " 'Nebl',\n",
       " 'Mllt10',\n",
       " 'Dnajc1',\n",
       " 'Commd3',\n",
       " 'Bmi1',\n",
       " 'Pip4k2a',\n",
       " 'Armc3',\n",
       " 'Otud1',\n",
       " 'Arhgap21',\n",
       " 'Thnsl1',\n",
       " 'Apbb1ip',\n",
       " 'Pdss1',\n",
       " 'Abi1',\n",
       " 'Acbd5',\n",
       " 'Mastl',\n",
       " 'Yme1l1',\n",
       " 'Spopl',\n",
       " 'Il1rn',\n",
       " 'Psd4',\n",
       " 'Ehmt1',\n",
       " 'Arrdc1',\n",
       " 'Zmynd19',\n",
       " 'Dph7',\n",
       " 'Mrpl41',\n",
       " 'Pnpla7',\n",
       " 'Nsmf',\n",
       " 'Nrarp',\n",
       " 'Tor4a',\n",
       " 'Nelfb',\n",
       " 'Tubb4b',\n",
       " 'Ndor1',\n",
       " 'Tmem203',\n",
       " 'Tprn',\n",
       " 'Ssna1',\n",
       " 'Anapc2',\n",
       " 'Man1b1',\n",
       " 'Dpp7',\n",
       " 'Uap1l1',\n",
       " 'Sapcd2',\n",
       " 'Npdc1',\n",
       " 'Fut7',\n",
       " 'Abca2',\n",
       " 'Paxx',\n",
       " 'Fbxw5',\n",
       " 'Traf2',\n",
       " 'Edf1',\n",
       " 'Phpt1',\n",
       " 'Ajm1',\n",
       " 'Rabl6',\n",
       " 'Camsap1',\n",
       " 'Ubac1',\n",
       " 'Nacc2',\n",
       " 'Tmem250-ps',\n",
       " 'Qsox2',\n",
       " 'Gpsm1',\n",
       " 'Dnlz',\n",
       " 'Snapc4',\n",
       " 'Entr1',\n",
       " 'Pmpca',\n",
       " 'Inpp5e',\n",
       " 'Sec16a',\n",
       " 'Notch1',\n",
       " 'Egfl7',\n",
       " 'Dipk1b',\n",
       " 'Med22',\n",
       " 'Rpl7a',\n",
       " 'Surf1',\n",
       " 'Surf2',\n",
       " 'Surf4',\n",
       " 'Rexo4',\n",
       " 'Cacfd1',\n",
       " 'Slc2a6',\n",
       " 'Brd3os',\n",
       " 'Brd3',\n",
       " 'Wdr5',\n",
       " 'Rxra',\n",
       " 'Col5a1',\n",
       " 'Olfm1',\n",
       " 'Mrps2',\n",
       " 'Ralgds',\n",
       " 'Gtf3c5',\n",
       " 'Tsc1',\n",
       " 'Gtf3c4',\n",
       " 'Ddx31',\n",
       " 'Ttf1',\n",
       " 'Setx',\n",
       " 'Med27',\n",
       " 'Rapgef1',\n",
       " 'Trub2',\n",
       " 'Coq4',\n",
       " 'Urm1',\n",
       " 'Cercam',\n",
       " 'Odf2',\n",
       " 'Gle1',\n",
       " 'Sptan1',\n",
       " 'Wdr34',\n",
       " 'Set',\n",
       " 'Zdhhc12',\n",
       " 'Zer1',\n",
       " 'Tbc1d13',\n",
       " 'Endog',\n",
       " 'Spout1',\n",
       " 'Kyat1',\n",
       " 'Lrrc8a',\n",
       " 'Nup188',\n",
       " 'Sh3glb2',\n",
       " 'Dolpp1',\n",
       " 'Ptpa',\n",
       " 'Ier5l',\n",
       " '1700001O22Rik',\n",
       " 'Ntmt1',\n",
       " 'Asb6',\n",
       " 'Ptges',\n",
       " 'Tor1b',\n",
       " 'Tor1a',\n",
       " 'BC005624',\n",
       " 'Fnbp1',\n",
       " 'D330023K18Rik',\n",
       " 'Gpr107',\n",
       " 'Ncs1',\n",
       " 'Ass1',\n",
       " 'Fubp3',\n",
       " 'Exosc2',\n",
       " 'Abl1',\n",
       " 'Nup214',\n",
       " 'Fam78a',\n",
       " 'Prrc2b',\n",
       " 'Pomt1',\n",
       " 'Uck1',\n",
       " 'Swi5',\n",
       " 'Golga2',\n",
       " 'Dnm1',\n",
       " 'Ciz1',\n",
       " '1110008P14Rik',\n",
       " 'Lcn2',\n",
       " 'Slc25a25',\n",
       " 'Fam102a',\n",
       " 'Dpm2',\n",
       " 'St6galnac4',\n",
       " 'St6galnac6',\n",
       " 'Ak1',\n",
       " 'Eng',\n",
       " 'Fpgs',\n",
       " 'Cdk9',\n",
       " 'Sh2d3c',\n",
       " 'Tor2a',\n",
       " 'Ptrh1',\n",
       " 'Stxbp1',\n",
       " 'Fam129b',\n",
       " 'Lrsam1',\n",
       " 'Rpl12',\n",
       " 'Slc2a8',\n",
       " 'Ralgps1',\n",
       " 'Angptl2',\n",
       " 'Mvb12b',\n",
       " 'Pbx3',\n",
       " 'Mapkap1',\n",
       " 'Gapvd1',\n",
       " 'Hspa5',\n",
       " 'Rabepk',\n",
       " 'Fbxw2',\n",
       " 'Psmd5',\n",
       " 'Traf1',\n",
       " 'Cntrl',\n",
       " 'Rab14',\n",
       " 'Gsn',\n",
       " 'Ggta1',\n",
       " 'Dab2ip',\n",
       " 'Ndufa8',\n",
       " 'Rbm18',\n",
       " 'Mrrf',\n",
       " 'Pdcl',\n",
       " 'Rc3h2',\n",
       " 'Zbtb26',\n",
       " 'Rabgap1',\n",
       " 'Strbp',\n",
       " 'Crb2',\n",
       " 'Dennd1a',\n",
       " 'Lhx2',\n",
       " 'Nek6',\n",
       " 'Psmb7',\n",
       " 'Rpl35',\n",
       " 'Arpc5l',\n",
       " 'Scai',\n",
       " 'Ppp6c',\n",
       " 'Arhgap15',\n",
       " 'Gtdc1',\n",
       " 'Zeb2',\n",
       " 'Zeb2os',\n",
       " 'Acvr2a',\n",
       " 'Orc4',\n",
       " 'Mbd5',\n",
       " 'Epc2',\n",
       " 'Mmadhc',\n",
       " 'Rnd3',\n",
       " 'Rbm43',\n",
       " 'Nmi',\n",
       " 'Tnfaip6',\n",
       " 'Rif1',\n",
       " 'Neb',\n",
       " 'Arl5a',\n",
       " 'Stam2',\n",
       " 'Fmnl2',\n",
       " 'Prpf40a',\n",
       " 'Arl6ip6',\n",
       " 'Nr4a2',\n",
       " 'Gpd2',\n",
       " 'Ermn',\n",
       " 'Cytip',\n",
       " 'Gm13546',\n",
       " 'Acvr1',\n",
       " 'Pkp4',\n",
       " 'Tanc1',\n",
       " 'Wdsub1',\n",
       " 'Baz2b',\n",
       " 'March7',\n",
       " 'Ly75',\n",
       " 'Rbms1',\n",
       " 'Tank',\n",
       " 'Psmd14',\n",
       " 'Dpp4',\n",
       " 'Ifih1',\n",
       " 'Cobll1',\n",
       " 'Scn2a',\n",
       " 'Ttc21b',\n",
       " 'B3galt1',\n",
       " 'Stk39',\n",
       " 'Cers6',\n",
       " 'Spc25',\n",
       " 'Fastkd1',\n",
       " 'Ppig',\n",
       " 'Phospho2',\n",
       " 'Ssb',\n",
       " 'Mettl5',\n",
       " 'Ubr3',\n",
       " 'Gorasp2',\n",
       " 'Tlk1',\n",
       " 'Mettl8',\n",
       " 'Dcaf17',\n",
       " 'Dync1i2',\n",
       " 'Slc25a12',\n",
       " 'Hat1',\n",
       " 'Metap1d',\n",
       " 'Dlx2',\n",
       " 'Itga6',\n",
       " 'Pdk1',\n",
       " 'Rapgef4',\n",
       " 'Map3k20',\n",
       " 'Cdca7',\n",
       " 'Sp3',\n",
       " 'Sp3os',\n",
       " 'Ola1',\n",
       " 'Cir1',\n",
       " 'Scrn3',\n",
       " 'Gpr155',\n",
       " 'Wipf1',\n",
       " 'Chrna1',\n",
       " 'Atf2',\n",
       " 'Atp5g3',\n",
       " 'Lnpk',\n",
       " 'Hoxd8',\n",
       " 'Mtx2',\n",
       " 'Hnrnpa3',\n",
       " 'Nfe2l2',\n",
       " 'Agps',\n",
       " 'Rbm45',\n",
       " 'Osbpl6',\n",
       " 'Prkra',\n",
       " 'Fkbp7',\n",
       " 'Plekha3',\n",
       " 'Sestd1',\n",
       " 'Cwc22',\n",
       " 'Ube2e3',\n",
       " 'Itga4',\n",
       " 'Itprid2',\n",
       " 'Pde1a',\n",
       " 'Dnajc10',\n",
       " 'Nckap1',\n",
       " 'Dusp19',\n",
       " 'Nup35',\n",
       " 'Zc3h15',\n",
       " 'Itgav',\n",
       " 'Fam171b',\n",
       " 'Calcrl',\n",
       " 'Ctnnd1',\n",
       " 'Selenoh',\n",
       " 'Tmx2',\n",
       " 'Med19',\n",
       " 'Zdhhc5',\n",
       " 'Clp1',\n",
       " 'Serping1',\n",
       " 'Ube2l6',\n",
       " 'Timm10',\n",
       " 'Ssrp1',\n",
       " 'Tnks1bp1',\n",
       " 'Ptprj',\n",
       " 'Nup160',\n",
       " 'Fnbp4',\n",
       " 'Mtch2',\n",
       " 'C1qtnf4',\n",
       " 'Ndufs3',\n",
       " 'Ptpmt1',\n",
       " 'Celf1',\n",
       " 'Rapsn',\n",
       " 'Psmc3',\n",
       " 'Slc39a13',\n",
       " 'Spi1',\n",
       " 'Madd',\n",
       " ...]"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gbm_genes = genes_expressed(cluster='6_TREAT')\n",
    "# Returns a list of all the genes expressed in a given cluster\n",
    "gbm_genes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "5c2095e9",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['H13', 'Rpn2', 'F3', 'Cad', 'Bst1', 'Met', 'Rpn1', 'Spn', 'Kitl']"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Surface genes expressed in GBM treatment group\n",
    "cell_surface_genes = match_surf(synonyms)\n",
    "[gene for gene in gbm_genes if gene in cell_surface_genes]"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d508ae0d",
   "metadata": {},
   "source": [
    "## lookup_gene tools"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "71548039",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'6_CTRL': 0.748062015503876,\n",
       " '11_TREAT': 0.6666666666666666,\n",
       " '8_CTRL': 0.6448362720403022,\n",
       " '8_TREAT': 0.6397849462365591,\n",
       " '9_CTRL': 0.6392914653784219,\n",
       " '12_CTRL': 0.6203703703703703,\n",
       " '24_CTRL': 0.6078431372549019,\n",
       " '18_CTRL': 0.6039603960396039,\n",
       " '2_CTRL': 0.5915119363395226,\n",
       " '3_CTRL': 0.5892439594699922,\n",
       " '26_TREAT': 0.5853658536585366,\n",
       " '1_TREAT': 0.5683975432719152,\n",
       " '2_TREAT': 0.5658852061438965,\n",
       " '12_TREAT': 0.5444444444444444,\n",
       " '21_TREAT': 0.5384615384615384,\n",
       " '1_CTRL': 0.522331566994701,\n",
       " '9_TREAT': 0.5221238938053098,\n",
       " '18_TREAT': 0.5066079295154186,\n",
       " '10_CTRL': 0.5,\n",
       " '24_TREAT': 0.5,\n",
       " '15_TREAT': 0.479020979020979,\n",
       " '11_CTRL': 0.4725,\n",
       " '21_CTRL': 0.46551724137931033,\n",
       " '6_TREAT': 0.46153846153846156,\n",
       " '0_CTRL': 0.4528149048197651,\n",
       " '3_TREAT': 0.43859649122807015,\n",
       " '10_TREAT': 0.41594202898550725,\n",
       " '5_CTRL': 0.4020172910662824,\n",
       " '17_TREAT': 0.4,\n",
       " '0_TREAT': 0.39009661835748793,\n",
       " '13_TREAT': 0.38016528925619836,\n",
       " '14_TREAT': 0.3746397694524496,\n",
       " '15_CTRL': 0.3644859813084112,\n",
       " '17_CTRL': 0.3463855421686747,\n",
       " '4_TREAT': 0.3433734939759036,\n",
       " '5_TREAT': 0.33076923076923076,\n",
       " '13_CTRL': 0.30446194225721784,\n",
       " '16_CTRL': 0.28703703703703703,\n",
       " '22_TREAT': 0.27485380116959063,\n",
       " '26_CTRL': 0.25,\n",
       " '4_CTRL': 0.23529411764705882,\n",
       " '16_TREAT': 0.22388059701492538,\n",
       " '14_CTRL': 0.21551724137931033,\n",
       " '19_TREAT': 0.1794871794871795,\n",
       " '20_CTRL': 0.16326530612244897,\n",
       " '22_CTRL': 0.14285714285714285,\n",
       " '20_TREAT': 0.14173228346456693,\n",
       " '25_TREAT': 0.11363636363636363,\n",
       " '25_CTRL': 0.1111111111111111,\n",
       " '19_CTRL': 0.10638297872340426,\n",
       " '7_CTRL': 0.07630522088353414,\n",
       " '7_TREAT': 0.07091172214182344,\n",
       " '23_TREAT': 0.040983606557377046}"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sorted_counter_dict(gene='Cd44', data=data, proportion=True)\n",
    "#     Returns a list (technically, a \"dict\") of clusters and proportion of cells \n",
    "#        in each cluster that express 'gene'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "f0665507",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'25_TREAT': 0.6818181818181818,\n",
       " '26_CTRL': 0.25,\n",
       " '25_CTRL': 0.2222222222222222,\n",
       " '3_TREAT': 0.03508771929824561,\n",
       " '10_TREAT': 0.013043478260869565,\n",
       " '16_TREAT': 0.007462686567164179,\n",
       " '12_TREAT': 0.0022222222222222222,\n",
       " '1_TREAT': 0.0011166945840312675,\n",
       " '4_TREAT': 0.0008605851979345956}"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "lookup_gene(gene='Prom1', sure=True)\n",
    "#     By default, will see whether input gene is unique,\n",
    "#         If the input matches more than one gene in list, returns list of genes that match\n",
    "#                 (set 'sure' equal to 'True' if you are sure this is the gene you want to look up\n",
    "#                  - apply this function to gene='Cd5' with 'sure' equal to 'True', and then separately\n",
    "#                    'sure' equal to 'False', to see the effect)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "20749d6a",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The genes and their most common clusters are:\n",
      "{'Prom1': '25_TREAT', 'Cd44': '6_CTRL', 'Hmgn1': '26_CTRL', 'Actb': '1_CTRL', 'Cd9': '22_CTRL', 'Sox2': '26_TREAT'}\n",
      "\n",
      "The most common cluster is 25_TREAT, followed by 22_CTRL, followed by 26_CTRL\n"
     ]
    }
   ],
   "source": [
    "lookup_gene(gene_list=['Prom1', 'Cd44', 'Hmgn1', 'Actb', 'Cd9', 'Sox2'], counter=True)\n",
    "#         In order to input a list of genes, you need to set 'counter' equal to 'True'"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "81d06c62",
   "metadata": {},
   "source": [
    "## plotting_tools tools"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "1412908e",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The expression distribution of Sox2 in cluster 6_CTRL:\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "{'Median': 1.0,\n",
       " 'Mean': 1.4061135371179039,\n",
       " 'Std': 0.8126495590884275,\n",
       " 'High': 6,\n",
       " 'Low': 1,\n",
       " 'Pct of cells': 0.22189922480620156}"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "expression_distribution(gene='Sox2', cluster='6_CTRL')\n",
    "#     Returns a list (technically, a \"dict\") of statistical features and those features \n",
    "#         values for the distribution of a particular gene's expression in a cluster"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "e06dcef7",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA48AAAGbCAYAAABksaWrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAABPe0lEQVR4nO3debit53z/8fcnk8wJEpH5SBBECSKoKTXGrFVEUzORtn5oo6ZWpYqqGqpVQ4iGiogxNChqViKCGCJSGSUyTzKIkPj+/rjvnaysrLXXOufsYZ1z3q/rWtde+xnu53s/072+z5iqQpIkSZKk+ay33AFIkiRJkmafyaMkSZIkaSKTR0mSJEnSRCaPkiRJkqSJTB4lSZIkSROZPEqSJEmSJtpguQNYKNtss02tWLFiucOQJC2B7373uxdV1bbLHceawjZSktYNi90+rjXJ44oVKzj++OOXOwxJ0hJIcuZyx7AmsY2UpHXDYrePXrYqSZIkSZrI5FGSJEmSNJHJoyRJkiRpIpNHSZIkSdJEJo+SJEmSpIlMHiVJkiRJE5k8SpIkSZImWpLkMcnOSb6c5KQkJyZ5Ye/+z0l+muSHST6RZOvefUWSq5Oc0D/vXIo4JUmadePaVEmSFttSnXm8Fji4qu4I3Bv4iyR3Ar4A3Lmq7gL8H/DygXFOraq9+uegJYpTkqRZN65NlSRpUS1J8lhV51bV9/r3K4CTgB2r6vNVdW0f7Fhgp6WIR5KkNdW4NnV5o5IkrQuW/J7HJCuAuwHfHur1LOCzA//fJsn3k3w1yf2XKj5JktYU87SpkiQtuA2WcmJJNgc+Bryoqi4f6P43tMtwjuidzgV2qaqLk9wDODrJnoPj9PEOBA4E2GWXXZaiCmuEFS/79IKXecbrH7XgZUqSVt24NnWg/4K2kYvRtmjV2CZLWi5LduYxyYa0Ru6Iqvr4QPenA48GDqiqAqiqa6rq4v79u8CpwO2Hy6yqQ6tq76rae9ttt12KakiStOzGtamDbCMlSQttqZ62GuAw4KSqevNA9/2AlwKPrapfDXTfNsn6/ftuwO2A05YiVkmSZtm4NlWSpMW2VJet3hd4KvCjJCf0bq8A/hW4GfCF1hZybH+y6gOAVye5FrgOOKiqLlmiWCVJmmUj29Sq+szyhSRJWhcsSfJYVd8AMqLXyIauqj5GuxxHkiQNmKdNlSRpUS3501YlSZIkSWsek0dJkiRJ0kQmj5IkSZKkiUweJUmSJEkTmTxKkiRJkiYyeZQkSZIkTWTyKEmSJEmayORRkiRJkjSRyaMkSZIkaSKTR0mSJEnSRCaPkiRJkqSJTB4lSZIkSROZPEqSJEmSJjJ5lCRJkiRNZPIoSZIkSZrI5FGSJEmSNJHJoyRJkiRpIpNHSZIkSdJEJo+SJEmSpIlMHiVJkiRJEy1J8phk5yRfTnJSkhOTvLB3v0WSLyT5Wf9784FxXp7klCQnJ3n4UsQpSZIkSRptqc48XgscXFV3BO4N/EWSOwEvA75YVbcDvtj/p/fbH9gT2A94e5L1lyhWSZIkSdKQJUkeq+rcqvpe/34FcBKwI/A44H19sPcBj+/fHwd8qKquqarTgVOAfZYiVkmSJEnSTS35PY9JVgB3A74NbFdV50JLMIFb9cF2BM4aGO3s3m24rAOTHJ/k+AsvvHBR45YkSZKkddmSJo9JNgc+Bryoqi6fb9AR3eomHaoOraq9q2rvbbfddqHClCRJkiQNWbLkMcmGtMTxiKr6eO98fpLte//tgQt697OBnQdG3wk4Z6lilSRJkiTd2FI9bTXAYcBJVfXmgV6fAp7evz8d+ORA9/2T3CzJbYDbAcctRaySJEmSpJvaYImmc1/gqcCPkpzQu70CeD3w4STPBn4OPBGgqk5M8mHgJ7Qntf5FVV23RLFKkiRJkoYsSfJYVd9g9H2MAA8eM85rgdcuWlCSJEmSpKkt+dNWJUmSJElrHpNHSZIkSdJEJo+SJEmSpIlMHiVJkiRJE5k8SpIkSZImMnmUJEmSJE1k8ihJkiRJmsjkUZIkSZI0kcmjJEmSJGkik0dJkiRJ0kQmj5IkSZKkiUweJUmSJEkTmTxKkiRJkiYyeZQkSZIkTWTyKEmSJEmayORRkiRJkjSRyaMkSZIkaSKTR0mSJEnSRCaPkiRJkqSJTB4lSZIkSRMtSfKY5L1JLkjy44FuRyU5oX/OSHJC774iydUD/d65FDFKkiRJksbbYImmczjwNuD9cx2q6slz35O8CfjlwPCnVtVeSxSbJEmSJGmCJUkeq+prSVaM6pckwJOABy1FLJIkSZKklTcL9zzeHzi/qn420O02Sb6f5KtJ7j9uxCQHJjk+yfEXXnjh4kcqSZIkSeuoWUgenwIcOfD/ucAuVXU34K+ADybZctSIVXVoVe1dVXtvu+22SxCqJEmSJK2bljV5TLIB8EfAUXPdquqaqrq4f/8ucCpw++WJUJIkSZIEy3/m8SHAT6vq7LkOSbZNsn7/vhtwO+C0ZYpPkiRJksTSvarjSOBbwB5Jzk7y7N5rf258ySrAA4AfJvkB8FHgoKq6ZCnilCRJkiSNtlRPW33KmO7PGNHtY8DHFjsmSZIkSdL0lvuyVUmSJEnSGsDkUZIkSZI0kcmjJEmSJGkik0dJkiRJ0kQmj5IkSZKkiUweJUmSJEkTmTxKkiRJkiYyeZQkSZIkTWTyKEmSJEmayORRkiRJkjSRyaMkSZIkaSKTR0mSJEnSRCaPkiRJkqSJTB4lSZIkSRNNlTwmeUqSO/bveyT5WpIvJbnD4oYnSZIkSZoF0555fA1wSf/+RuA44GvA2xcjKEmSJEnSbNlgyuG2rarzk2wM3A/4Y+C3wEWLFpkkSZIkaWZMmzxemOS2wO8B36mqa5JsCmTxQpMkSZIkzYppk8d/AL4LXAc8uXd7MPCDxQhKkiRJkjRbprrnsaoOB7YHdqqqL/TO3wb2n2b8JO9NckGSHw90OyTJL5Kc0D+PHOj38iSnJDk5ycOnro0kSZIkaVGMPfOYZDix/PVQ95W53/Fw4G3A+4e6v6Wq3jg03TvRktI9gR2A/0ly+6q6biWmJ0mSJElaQPOdebyW9lCccZ+5/hNV1de44WmtkzwO+FBVXVNVpwOnAPtMOa4kSZIkaRHMd8/jbZZg+s9P8jTgeODgqroU2BE4dmCYs3s3SZIkSdIyGZs8VtWZizztd9AexFP975uAZzH6Ca41qoAkBwIHAuyyyy6LE6UkSZIkad57Hv+TMUnboKp62qpMuKrOH5jWu4Fj+r9nAzsPDLoTcM6YMg4FDgXYe++9J8YqSZIkSVo18122espiTjjJ9lV1bv/3D4G5J7F+CvhgkjfTHphzO+C4xYxFkiRJkjS/+S5b/fuFmkiSI4F9gW2SnA28Ctg3yV60s5tnAM/r0z0xyYeBn9AeyvMXPmlVkiRJkpbXfGcebyTJQ2mv0LhVVT0myd7AllX1pUnjVtVTRnQ+bJ7hXwu8dtrYJEmSJEmLa75XdVwvyf+jPeDmZ8ADeuergdcsUlySJEmSpBkyVfIIvAh4SFW9Hvhd7/ZTYI/FCEqSJEmSNFumTR63AM7q3+eearoh8JsFj0iSJEmSNHOmTR6/BrxsqNsLgC8vbDiSJEmSpFk07QNz/h/wX0meC2yR5GTgcuAxixaZJEmSJGlmTJU8VtW5Se4J3BPYlXYJ63FV9bv5x5QkSZIkrQ2mSh77+xgvrqrjgON6t52T3KKqfrCI8UmSJEmSZsC09zx+gPaAnEEbAf+5sOFIkiRJkmbRtMnjLlV12mCHqjoVWLHgEUmSJEmSZs60yePZSe4+2KH/f87ChyRJkiRJmjXTPm31LcAnk7wBOBXYHXgx8NrFCkySJEmSNDumfdrqu5NcBjwb2Jn2tNWDq+qjixibJEmSJGlGTHvmkar6CPCRRYxFkiRJkjSjpr3nUZIkSZK0DjN5lCRJkiRNZPIoSZIkSZpobPKY5NiB769amnAkSZIkSbNovjOPt0+ycf9+8FIEI0mSJEmaTfM9bfWTwP8lOQPYJMnXRg1UVQ9YjMAkSZIkSbNjbPJYVc9Mcj9gBXBP4LClCkqSJEmSNFvmfc9jVX0D+EaSjarqfas6kSTvBR4NXFBVd+7d/hl4DPAb4FTgmVV1WZIVwEnAyX30Y6vqoFWdtiRJkiRp9U31tNWqem+SP0jy3iSf638ftBLTORzYb6jbF4A7V9VdgP8DXj7Q79Sq2qt/TBwlSZIkaZlNlTwmeQ5wFHAe8HHgXOCDSZ47zfhV9TXgkqFun6+qa/u/xwI7TRu0JEmSJGlpzXvZ6oCXAA+tqh/MdUhyFPAx4N0LEMezaMnpnNsk+T5wOfC3VfX1BZiGJEmSJGkVTZs83hL4yVC3k4FbrG4ASf4GuBY4onc6F9ilqi5Ocg/g6CR7VtXlI8Y9EDgQYJdddlndUCRJkiRJY0x12SrwDeDNSTYFSLIZ8M/AN1dn4kmeTnuQzgFVVQBVdU1VXdy/f5f2MJ3bjxq/qg6tqr2rau9tt912dUKRJEmSJM1j2uTxIOAuwC+TnA9cBtwVeN6qTjjJfsBLgcdW1a8Gum+bZP3+fTfgdsBpqzodSZIkSdLqm+qy1ao6F3hgkp2AHYBzqursaSeS5EhgX2CbJGcDr6I9XfVmwBeSwA2v5HgA8Ook1wLXAQdV1SUjC5YkSZIkLYlp73kEoCeMUyeNA+M9ZUTnw8YM+zHag3gkSZIkSTNi2stWJUmSJEnrMJNHSZIkSdJEE5PHJOsleVCSjZYiIEmSJEnS7JmYPFbV74BPVtVvliAeSZIkSdIMmvay1a8lufeiRiJJkiRJmlnTPm31TOCzST4JnAXUXI+q+rvFCEySJEmSNDumTR43AY7u33danFAkSdI0krwXeDRwQVXdebnjkSStG6ZKHqvqmYsdiCRJmtrhwNuA9y9zHJKkdcjUr+pIcsckr0zytv7/HknusnihSZKkUarqa8Alyx2HJGndMlXymOSJwNeAHYGn9c5bAG9epLgkSdJqSHJgkuOTHH/hhRcudziSpLXAtGceXw08tKoOAq7r3X4A3HVRopIkSaulqg6tqr2rau9tt912ucORJK0Fpk0eb0VLFuGGJ63WwHdJkiRJ0lps2uTxu8BTh7rtDxy3sOFIkiRJkmbRtK/qeAHw+STPBjZL8jng9sDDFi0ySZI0UpIjgX2BbZKcDbyqqg5b3qgkSWu7aV/V8dMkd6C9U+oY4CzgmKq6cjGDkyRJN1VVT1nuGCRJ655pzzxSVb9K8r/A6cA5Jo6SJEmStO6Y9lUduyT5OnAG8GngjCTfSLLrYgYnSZIkSZoN0z4w5320h+ZsXVW3Am4OfKd3lyRJkiSt5aa9bPUewMOq6rcAVXVlkpcCFy9aZJIkSZKkmTHtmcdjgX2Guu0NfGthw5EkSZIkzaKxZx6TvHrg31OBzyT5NO1JqzsDjwQ+OM1EkryX9qTWC6rqzr3bLYCjgBW0eymfVFWX9n4vB54NXAe8oKo+t1K1kiRJkiQtqPnOPO488NkY+DhwDXCr/vcTvfs0Dgf2G+r2MuCLVXU74Iv9f5LcCdgf2LOP8/Yk6085HUmSJEnSIhh75rGqnrlQE6mqryVZMdT5cbQXHEN78M5XgJf27h+qqmuA05OcQrtk1ktkJUmSJGmZTP2exySbArcFNh/sXlXfXMVpb1dV5/Yyzk1yq959R9o9lnPO7t1GxXQgcCDALrvssophSJIkSZImmSp5TPI04G3Ab4CrB3oVsNBZW0Z0q1EDVtWhwKEAe++998hhJEmSJEmrb9ozj28AnlBVX1jAaZ+fZPt+1nF74ILe/WzafZZzdgLOWcDpSpIkSZJW0rSv6vgN7Z7EhfQp4On9+9OBTw503z/JzZLcBrgdcNwCT1uSJEmStBKmTR5fCbw5yTarMpEkR9IeeLNHkrOTPBt4PfDQJD8DHtr/p6pOBD4M/AT4b+Avquq6VZmuJEmSJGlhTHvZ6v8Brwb+PLn+lsQAVVUTX6NRVU8Z0+vBY4Z/LfDaKWOTJEmSJC2yaZPH/wTeDxzFjR+YI0mSJElaB0ybPN4S+Luq8ommkiRJkrQOmvaex/8AnrqYgUiSJEmSZte0Zx73AZ6f5G+A8wd7VNUDFjwqSZIkSdJMmTZ5fHf/SJIkSZLWQVMlj1X1vsUORJIkSZI0u6ZKHpM8a1y/qnrvwoUjSZIkSZpF0162OvywnFsDuwP/C5g8SpIkSdJabtrLVv9guFs/G3nHBY9IkiRJkjRzpn1VxyiHA89eoDgkSZIkSTNs2nseh5PMTYE/BS5b6IAkSZIkSbNn2nserwVqqNsvgOcubDiSJEmSpFk0bfJ4m6H/r6qqixY6GEmSJEnSbJr2gTlnLnYgkiRJkqTZNW/ymOTL3PRy1UFVVQ9e2JAkSZIkSbNm0pnHD4zpviPwAtqDcyRJkiRJa7l5k8eqOmzw/yS3BF5Oe1DOUcCrFy80SZIkSdKsmOo9j0m2TPIPwCnAdsDdq+rAqjp7UaOTJEmSJM2EeZPHJJskeTlwGnBH4H5V9dSqOnVJopMkSZIkzYRJ9zyeDqwPvAE4HtguyXaDA1TVlxYpNkmSJEnSjJiUPP6a9rTVPxvTv4DdVnXiSfag3Ts5Zzfg74CtafdVXti7v6KqPrOq05EkSZIkrZ5JD8xZsZgTr6qTgb0AkqwP/AL4BPBM4C1V9cbFnL4kSZIkaTpTPTBniTwYOLWqzlzuQCRJkiRJNzZLyeP+wJED/z8/yQ+TvDfJzUeNkOTAJMcnOf7CCy8cNYgkSZIkaQHMRPKYZCPgscBHeqd3ALvTLmk9F3jTqPGq6tCq2ruq9t52222XIlRJkiRJWifNRPIIPAL4XlWdD1BV51fVdVX1O+DdwD7LGp0kSZIkreNmJXl8CgOXrCbZfqDfHwI/XvKIJEmSJEnXm/SqjkWXZFPgocDzBjq/IcletFeBnDHUT5IkSZK0xJY9eayqXwG3HOr21GUKR5IkSZI0wqxctipJkiRJmmEmj5IkSZKkiUweJUmSJEkTmTxKkiRJkiYyeZQkSZIkTWTyKEmSJEmayORRkiRJkjSRyaMkSZIkaSKTR0mSJEnSRCaPkiRJkqSJTB4lSZIkSROZPEqSJEmSJjJ5lCRJkiRNZPIoSZIkSZrI5FGSJEmSNJHJoyRJkiRpIpNHSZIkSdJEJo+SJEmSpIlMHiVJkiRJE22w3AEkOQO4ArgOuLaq9k5yC+AoYAVwBvCkqrp0uWKUJEmSpHXdsieP3R9U1UUD/78M+GJVvT7Jy/r/L12e0CRY8bJPL3iZZ7z+UQtepiRJkrRYZvWy1ccB7+vf3wc8fvlCkSRJkiTNQvJYwOeTfDfJgb3bdlV1LkD/e6tli06SJEmSNBOXrd63qs5JcivgC0l+Ou2IPdk8EGCXXXZZrPgkSZIkaZ237Gceq+qc/vcC4BPAPsD5SbYH6H8vGDPuoVW1d1Xtve222y5VyJIkSZK0zlnW5DHJZkm2mPsOPAz4MfAp4Ol9sKcDn1yeCCVJkiRJsPyXrW4HfCLJXCwfrKr/TvId4MNJng38HHjiMsYoSZIkSeu8ZU0eq+o04K4jul8MPHjpI5IkSZIkjbLs9zxKkiRJkmafyaMkSZIkaSKTR0mSJEnSRCaPkiRJkqSJTB4lSZIkSROZPEqSJEmSJjJ5lCRJkiRNZPIoSZIkSZrI5FGSJEmSNJHJoyRJkiRpIpNHSZIkSdJEJo+SJEmSpIlMHiVJkiRJE5k8SpIkSZImMnmUJEmSJE1k8ihJkiRJmsjkUZIkSZI0kcmjJEmSJGkik0dJkiRJ0kQmj5IkSZKkiZY1eUyyc5IvJzkpyYlJXti7H5LkF0lO6J9HLmeckiRJkrSu22CZp38tcHBVfS/JFsB3k3yh93tLVb1xGWOTJEmSJHXLmjxW1bnAuf37FUlOAnZczpgkSZIkSTc1M/c8JlkB3A34du/0/CQ/TPLeJDdfvsgkSZIkSTORPCbZHPgY8KKquhx4B7A7sBftzOSbxox3YJLjkxx/4YUXLlW4kiRJkrTOWfbkMcmGtMTxiKr6OEBVnV9V11XV74B3A/uMGreqDq2qvatq72233XbpgpYkSZKkdcxyP201wGHASVX15oHu2w8M9ofAj5c6NkmSJEnSDZb7aav3BZ4K/CjJCb3bK4CnJNkLKOAM4HnLEZwkSZIkqVnup61+A8iIXp9Z6lgkSZIkSeMt+z2PkiRJkqTZZ/IoSZIkSZrI5FGSJEmSNJHJoyRJkiRpIpNHSZIkSdJEJo+SJEmSpIlMHiVJkiRJE5k8SpIkSZImMnmUJEmSJE1k8ihJkiRJmsjkUZIkSZI00QbLHYAkLZQVL/v0gpd5xusfteBlSpIkrYk88yhJkiRJmsjkUZIkSZI0kcmjJEmSJGkik0dJkiRJ0kQmj5IkSZKkiUweJUmSJEkTmTxKkiRJkiYyeZQkSZIkTTTTyWOS/ZKcnOSUJC9b7ngkSZoFto+SpOUws8ljkvWBfwceAdwJeEqSOy1vVJIkLS/bR0nSctlguQOYxz7AKVV1GkCSDwGPA36yrFFJ0jpixcs+veBlnvH6Ry14mesg20dJ0rKY5eRxR+Csgf/PBu41OECSA4ED+79XJjl5Aaa7DXDRApSznBa8DvmnhSxtKmvDcoB56rEM83RVrfXLYj4ztpzW+GWRf1qwOuy6AGWsqSa2j7AobeQav/51a3w9FnA7Wk5rQx1g7ajH2lAHWDvqsRB1WNT2cZaTx4zoVjf6p+pQ4NAFnWhyfFXtvZBlLjXrMDvWhnqsDXWAtaMe1kHdxPYRFr6NXFuW3dpQD+swO9aGeqwNdYC1ox5rQh1m9p5H2pHUnQf+3wk4Z5likSRpVtg+SpKWxSwnj98BbpfkNkk2AvYHPrXMMUmStNxsHyVJy2JmL1utqmuTPB/4HLA+8N6qOnEJJr2gl8EuE+swO9aGeqwNdYC1ox7WQbaPq29tqId1mB1rQz3WhjrA2lGPma9Dqm5ym4QkSZIkSTcyy5etSpIkSZJmhMmjJEmSJGmidTJ5TLJfkpOTnJLkZSP675vkl0lO6J+/W44455PkvUkuSPLjMf2T5F97HX+Y5O5LHeM0pqjHmrAsdk7y5SQnJTkxyQtHDDPTy2PKOsz0skiycZLjkvyg1+HvRwwz08sBpq7HTC+LOUnWT/L9JMeM6Dfzy2JdZRs5G2wfZ4dt5GywfZwRVbVOfWgPFzgV2A3YCPgBcKehYfYFjlnuWCfU4wHA3YEfj+n/SOCztPeB3Rv49nLHvIr1WBOWxfbA3fv3LYD/G7FOzfTymLIOM70s+rzdvH/fEPg2cO81aTmsRD1melkMxPlXwAdHxbomLIt18WMbOTsf28fZ+dhGzsbH9nE2Puvimcd9gFOq6rSq+g3wIeBxyxzTSquqrwGXzDPI44D3V3MssHWS7ZcmuulNUY+ZV1XnVtX3+vcrgJOAHYcGm+nlMWUdZlqft1f2fzfsn+Engs30coCp6zHzkuwEPAp4z5hBZn5ZrKNsI2eE7ePssI2cDbaPs2FdTB53BM4a+P9sRu8A7tNPi382yZ5LE9qCmraea4I1ZlkkWQHcjXY0bNAaszzmqQPM+LLol4GcAFwAfKGq1sjlMEU9YMaXBfAvwEuA343pv0Ysi3WQbeSaZY1ZDmtD+wi2kcvN9nH5rYvJY0Z0Gz5q8T1g16q6K/BvwNGLHdQimKaea4I1Zlkk2Rz4GPCiqrp8uPeIUWZueUyow8wvi6q6rqr2AnYC9kly56FB1ojlMEU9ZnpZJHk0cEFVfXe+wUZ0m7llsQ6yjVxzrDHLYW1oH8E2chbYPi6/dTF5PBvYeeD/nYBzBgeoqsvnTotX1WeADZNss3QhLoiJ9VwTrCnLIsmGtAbliKr6+IhBZn55TKrDmrIsAKrqMuArwH5DvWZ+OQwaV481YFncF3hskjNolz0+KMkHhoZZo5bFOsQ2cg2xpiyHtaF9BNvIWWP7uHzWxeTxO8DtktwmyUbA/sCnBgdIcusk6d/3oc2ni5c80tXzKeBp/YlN9wZ+WVXnLndQK2tNWBY9vsOAk6rqzWMGm+nlMU0dZn1ZJNk2ydb9+ybAQ4CfDg0208sBpqvHrC+Lqnp5Ve1UVSto+9gvVdWfDg0288tiHWUbuYZYE5bD2tA+gm3k0kY6nu3jbNhguQNYalV1bZLnA5+jPVXuvVV1YpKDev93An8M/FmSa4Grgf2ramZOFwMkOZL2RKltkpwNvIp24/BcHT5De1rTKcCvgGcuT6Tzm6IeM78saEeRngr8KO06fIBXALvAGrM8pqnDrC+L7YH3JVmf1lh8uKqOGdq2Z305wHT1mPVlMdIauCzWObaRs8P2cabYRs4G28cZkDVgfkqSJEmSltm6eNmqJEmSJGklmTxKkiRJkiYyeZQkSZIkTWTyKEmSJEmayORRkiRJkjSRyaOklZLknUleudxxSJI0a2wjtbYzedQ6Lcn+Sb6d5KokF/Tvfz73gtkljOMZSb6xlNNcVVV1UFX9wzTDJjk8yWsWOyZJ0sKzjVx5tpFa25k8ap2V5GDgrcA/A7cGtgMOor0MeKNlDG2k/lLcdT4GSdLis41cM2OQFpvJo9ZJSbYCXg38eVV9tKquqOb7VXVAVV3Th7tZkjcm+XmS8/vlKJv0fvsmOTvJwf2I7LlJnjkwjbHjDsVyR+CdwH2SXJnkst798CTvSPKZJFcBf5DkUUm+n+TyJGclOWSgnBVJKsnT+zQvSvI3A/33SXJ8H/f8JG8e6He/JN9Mclkv9xnzxHD9kdKBefCKPr0zkhzQ+x0IHAC8pNfrv3r3lyb5RZIrkpyc5MELsEglSQvENtI2UhrH5FHrqvsANwM+OWG4fwJuD+wF3BbYEfi7gf63Brbq3Z8N/HuSm085LgBVdRLtaO63qmrzqtp6oPefAK8FtgC+AVwFPA3YGngU8GdJHj9U5P2APYAHA3/XG15oR5DfWlVbArsDHwZIsgvwWeDfgG17vCfME8OwWwPb9Po9HTg0yR5VdShwBPCGXq/HJNkDeD5wz6raAng4cMaIMiVJy8c20jZSGsnkUeuqbYCLqurauQ4DRxWvTvKAJAGeC/xlVV1SVVcArwP2Hyjnt8Crq+q3VfUZ4EpgjynHncYnq+p/q+p3VfXrqvpKVf2o//9D4EjggUPj/H1VXV1VPwB+ANx1INbbJtmmqq6sqmN79wOA/6mqI3s9Lq6qE8bFMCbOV1bVNVX1VeDTwJPGDHcd7QfJnZJsWFVnVNWpKzVHJEmLzTbSNlIayeRR66qLgW2SbDDXoap+vx/RvJi2bWwLbAp8tzeYlwH/3btfX85g4wr8Cth8ynGncdbgP0nuleTLSS5M8kva0dhthsY5b0Q80I763h74aZLvJHl0774zMF/jdNY8/QAuraqrBv4/E9hh1IBVdQrwIuAQ4IIkH0oyclhJ0rKxjbSNlEYyedS66lvANcDj5hnmIuBqYM+q2rp/tqqqzecZZ1XHrSm7fxD4FLBzVW1Fuw9kqqfeVdXPquopwK1olwt9NMlmtIZv9/lGnVD0zXs5c3YBzhk3blV9sKruB+za+//TNPFLkpaMbaRtpDSSyaPWSVV1GfD3wNuT/HGSzZOsl2QvYLM+zO+AdwNvSXIrgCQ7Jnn4FOWv7LjnAzslmfQEuy2AS6rq10n2od1rMZUkf5pk2x7bZb3zdbR7Lh6S5ElJNkhyyz4fVsbfJ9koyf2BRwMf6d3PB3YbiGGPJA9KcjPg17QfD9et5LQkSYvINtI2UhrH5FHrrKp6A/BXwEuAC2g78XcBLwW+2Qd7KXAKcGySy4H/od1oP42VGfdLwInAeUkumqfMPwdeneQK2oMFPjxlLAD7AScmuZL2YID9+z0iPwceCRwMXEJ7EMBdx5ZyU+cBl9KOpB4BHFRVP+39DqPdu3FZkqNp93K8nnbU+TzaEd5XrMS0JElLwDbSNlIaJVWTzrZL0mhJ9gU+UFU7LXMokiTNFNtIrY088yhJkiRJmsjkUZIkSZI0kZetSpIkSZIm8syjJEmSJGkik0dJkiRJ0kQmj5IkSZKkiUweJUmSJEkTmTxKkiRJkiYyeZQkSZIkTWTyKEmSJEmayORRkiRJkjSRyaMkSZIkaSKTR0mSJEnSRCaPkiRJkqSJTB4lSZIkSROZPEqSJEmSJjJ5lCRJkiRNZPIoSZIkSZrI5FGSJEmSNJHJoyRJkiRpojU2eUxSSW67BNN5Z5JXLvZ01lRJDk/ymjH9Dkjy+aWOaSGl+Y8klyY5bspxxs6TdU2SQ5J8oH9f0bfbDVahnDslOX7hI1xYSY5Lsuc8/Rdlm1iq/aHWHrahs2FNaEOTbJfka0muSPKmKcc5I8lDFju2pTK4vaxOG5/kyCSPX9DglliSjyfZb57+909y8iJMd61ap9ZkS5o8Jtk/ybeTXJXkgv79z5NkKeNYGVV1UFX9w3LHsSaqqiOq6mGrMm6Spyf5bpLLk5yd5A2DSUeSOyb5UpJfJjklyR8Ojf/gJD9N8qskX06y6ypW437AQ4GdqmqfEXE+I8k3VrHsiZI8Kck3ez2+MqL/Xn0+/ar/3Wug3/5JTu7z6IIk70uyZe93sySHJTmz/yD4fpJHLFY9FsA/AG+E6WLv8+2k3v8n8zXWSZ6f5Pgk1yQ5fKjfzkmOTXLJ8I+mJP+dZO+h4t4IvHrctFZnm1hsi3XQI8nd+w/PK5Ocn+SFU4yTJC9I8uPeXpyd5CNJfi/JZ3tZVyb5bZLfDPz/ziT7Jvld//+Kvg08c6j8NTLZtg1dtyxmG7qSDgQuArasqoNHTGtRD5gm+YckP0pybZJDhvr9Qe93WZKLk3wiyY6LFcvqSHIX4K7AJ/v/2yf5VJJz+j5pxdDwhw/t365Msv6YsuedD0lukeSoJBf1zxEDvwe2SvK5Pu4Rg9NI8u4M/b4CXg+8dlw9q+rrVbXHSs6eJZGBA9oLXO5uSY7pbc5FSd4wYfjBZfq7JFcP/H9Aj/O3/f/L0n4H3mdg/MF2bvBzn6HpHN63mx36/68YGPbXSa4b+P/Eaeu7ZMljkoOBtwL/DNwa2A44CLgvsNFSxbE2WY2GYE2wKfAiYBvgXsCDgRfD9fX+JHAMcAtaw/aBJLfv/bcBPg68svc/HjhqFePYFTijqq5a1YqspkuAf6HtrG8kyUa0+fAB4ObA+4BP9u4A/wvct6q2AnYDNgDmGvgNgLOABwJb0ebVh4cbr1mQZHvgD4Cje6d5Y+8N5geAvwK2BP4a+GCSW42ZxDm0+fLeEf1eTpuvtwEen54sJnkycFpVDZ8N/RTwBz3mdcqo/VHfFv8beBdwS+C2wDRnUt4KvBB4AW0bvj1t+T+qqh5RVZtX1ebAEcAb5v6vqoP6+Of0/lsCfwm8O8lM/piZlm3owltX29BVsCvwk6qqhQltpZ0CvAT49Ih+PwEeXlVbAzsAPwPesXShrZTnAUcMzMff0faPT5hnnMH92+ZVdd2Y4SbNh9fQfifsBuxO238cMhDX93u3FcAfAvREZPuq+sTghKrqOGDL3PTg6VpvTDu3EfAF4Eu0ffNOtN8gYw0uU+DnwGMGuh3RBzuq998G+DLwkaFizhlaNzavqm8NxLUZbd36JXBAn+7rBqZ7EPCtgXHHXjU1qgKL/qH9wLsKeMKE4R5J2wCuAH4BvHig318D59J+6D0LKOC2vd/hwNuBzwJX0n4035r2o/tS4KfA3QbK2gH4GHAhcDrwgnliOhx4Tf++DS1huYz2o/7rwHpjxrsDbWW6BDgZeFLvvnvvdveBWC4C9u3/fwX4R+C4vsA/Cdyi91vR6/1s2sr2td79WcBJva6fA3bt3QO8Bbigl/VD4M5TzOtHAyf0en4TuMtAv7sB3+vjHQV8aG7+jJgHzwC+MfB/0VbWn/VY/x3IlOvQXwH/1b/fuS/nDPT/PPAP/fuBwDcH+m0GXA3cYUzZO9B+9F9Ca6Se27s/G/g1cF2f3t8PjXfHof6XDawz/05r6K4Avg3sPmndmFD/5wBfGer2sL7sBufDz4H9Roy/OfB+4DPzTOOHzLONAs/t69kVfd0ZXIdHbk+0xukDQ+vvBgPrx2m9vNOBA8ZM92nA/0yYP9fHTvuhdMFQ/wuB+0wo4zXA4UPdPgvs0b9/CHgSLSH5PrD1mHK+ADx9obcJYH3gFcCpfZ59F9h5oJy5/eFXgOeMmiZj9gm0bea3wG/6ujy3rU1ath+lNZKXD05zYJjXAf85zTY+MM7taNvUPlMMezhD+x9gX+DsoW4XAE8cmu+3XZm4lvODbaht6EruL0aUdX0bOqb/7wPf6fX8DvD7A8tvcN/wkKHxxu07zqAlqz/sZR4FbDzNPJonxg8Ah8zT/2Z92f9knmF2ph1cvhC4GHjbQL+R68HAvB/cXlZlnT4NuN+I7hv08leM23ZWcn9xk/lA27b/fOD/vwA+17+/g5Z4QjtQ/RJae3MsA79dhqbxbuBVY/rty8A+eNK6MGL8cb81zphb/4bnzYhpvpS2XV5B2388GNivr6e/7evqD/qwWwGH0faPv6D9Flh/YBv8X9p+4JJRy4O2DXx9ZZfT0PwZ3q4Oof926v/fqa8j246q75hyn0Y70P5C4Mcj+j+Dgf3LSsW8qpVdyRmzH3At/UfjPMOdC9y/f7/5wAqzH3A+7UfOZsAHuemGfBFwD2BjWvZ/ep9x6/cV4ct92PVoP7r+jna0djfaBv3wMTFdv4LSNsZ3Ahv2z/0ZsePuMZ4FPJO2U7h7j2/PoQ1jU9oO6o0D436lr7xzdf0YN/3x/f7ebxPg8bSE5459Wn9LT5yAh/e6bk1rBO9IO4o037y+O62hvFefd0+nrdg36/PrTNqR/A2BP6ZthCvT8B3T49mFtvO+SaIzpqyjgdf377/HTZPHLwCf6N/fCrxjaPwfM+aHF/BV2g+njYG9elwPnmbjGtW/rzOXAPv0ZXIE8KFp1o15pjMqefxL4LND3Y4BDh74/360nXXRfnw+bEz529ES4XEJ9hNp6+U9+7p0W9rR6Hm3J8Ykj30+XM4Nidn24+YB7UzLv88zb24UO229/Srw2P798cDZwGYT5vGo5PGfgef3dfYU2nb5VsYkh32cfwXevNDbBO3H/4+APfoyuCtwy4Fypkke59snHM6NG+Nplu1v+/xdD9hkRMxf6vPrm7T9yn8Bu0xYDgcBZ065X7hRzL3bvvRGtcf1WNoR/sHk5/r5tSZ8sA21DV3J/cWIso6mt6Ej+t2CljA9tc+Dp/T/5/Yvh4+LcZ7t8AxaAr9DL/8k4KBJ82hCHUYmj31eXEbbzn8LPGPM+OsDP6AlApvR1vX79X5j14PhfQarvk5f/+N/qN98yeMl/fNdJh88GjsfaMn6Z2jr6s1p2/iLer+/oLV1m9ASpUf1dfRV80zrr4CPj+m3LzdNHkeuCyPGHflbY6CcickjrY08C9hhYLvfvX8/hIGkbGDbeFdfRrfqsT5vYBu8Fvh/fTmNaufeC/wnLUG/iLYP+r1ptsvheg10uz5O2n7j9b3sDYbrO0+5XwTeQPuNdC19HzXQ/xmsYvK4VJetbgNcVFXXznXo1+9e1q/zfUDv/FvgTkm2rKpLq+p7vfuTgP+oqh9Xu3zwkBHT+ERVfbeqfg18Avh1Vb2/2in+o2hH+6CtkNtW1aur6jdVdRrtCMr+U9Tjt7QfubtW1W+rXdddI4Z7NO1Sx/+oqmt7PT5GayioqnfTjhx+u5f3N0Pj/+dAXV8JPGnoOvdDquqqqrqadrnBP1bVSX3+vg7YK+0ev98CW9CO4KYPc+5AXUbN6+cC76qqb1fVdVX1PuAa4N79syHwL73+H6UdoVwZr6+qy6rq57TT8HtNGiHtfqW96fe80Y6CXwD8dZINkzyMdgnjpr3/5rSEadAvafNiuOydaQnWS6vq11V1AvAeWiO6Oj5eVcf1ZXIEN9Rz3nVjJU2sZ1V9o9plqzvRGoczhgtJsmGP8X1V9dMx03oO7fKZ71RzSlWdyeptT78D7pxkk6o6t6rGXW+/Ne3o4U2Mir1v8++n/UC+pv99Xq3apcf/SPsx8FXaUf4NgbsA/5Xkg2n38j1/aJwreszTmnabeA7wt1V1cl8GP6iqi1diOjD/PmHYNMv2W1V1dFX9ru+Phu1E+2H4QtoPm9OBIyfEeEvaD/PVsUOSy2hXHHwC+Kuq+v5qlrmcbENtQ+csRBs67FHAz6rqP/v8PpLWzj5mJWMb9q9VdU5VXUI7cDQX63zzaKVV1c+rXa65DS3pG9eO7UNLYP66L/9fV9XcMwvmWw/mM+06vXX/O7ItG+NfaVdi3Iq2Hh+e5L7jBp4wH75HS0Iu7p/raAfNoZ1124q2PX2dlmA/FfiXJO/o7dzwPa0r286NWxeGjfutsTKuox2suVOSDavqjKo6ddSASbYDHkFLpK+qqgtoBxcG92fnVNW/9W1jXDu3P2157UC76mzwFqJV9aSBduy5wB8PtgH0dm7os1mv1y60230+WFXn0xLJp69mPNdbquTxYmCbwWuFq+r3+0p+8UAcT6BdCnJmkq8O3Pi5A+0owpxRK9L5A9+vHvH/5v37rgzNcNqlYNtNUY9/ph2Z+nyS05K8bMxwuwL3GprGAbTLgOa8m3Zk9N+q6pqh8YfruiFtZzCq/67AWwemcwntaM2OVfUl4G20H73nJzk0/QZpxs/rXYGDh2LfmbYMdgB+MbRjXNmN+ryB77/ihuUyUtqDTl4PPKKqLgKoqrmzHY/q5R0MfJh2dgnaWckth4raktE77R2AS6pqsN+ZwOrecD+untOsG9Oaup5V9QvavRUfGuyeZD3aEbPf0M6wjbMz7XLJYau0PfUfdU+mnWU6N8mnk9xhzOCXMjrxHxl72tPY3kA7MrcR7cDCezLwMKFpVdUlVfXkqror7Qzav9GOQL6Mdjb7IcBBSe40MNoWtKO/05p2mxi3DKY2YZ8wbJple9bIMW9wNS0p+U5PSv4e+P0kW80zzsW0H2Or45zevmxJa9AftJrlLTfbUNvQOavdho6ww4g4FrsdHDePVllPTObu/R91P+vOtKsarh3Rb+x6MGGy067Tl/W/N2nLxqmq71XVxT1p+QztQOkfTTHeqPnwEeD/+vS3pLUlH+jD/7qqDqyqu1TVy2jJ0yto29z6tDb0XrnxE1ZnuZ07hXa/7yHABUk+lP7AmBF2pe0fzh1Y9u+iJexzpmnnvlFVn62q39AO0tySdhZ7dXy47+O3o/3euMdQ/3Oqauuhz9xB8qcCJ1U7IQJt3fmTfsB9tS1V8vgt2lGlx803UP+B8TjaQjualhBAOwq988Cgu6xGLGcBpw/N7C2q6pGTRqyqK6rq4KrajXZE7q+SPHjMNL46NI3Nq+rPAJJsTruX5DDgkCS3GBp/uK6/pZ2uvj6UoWk9b2ham1TVN3vM/1pV9wD2pD104q9793Hz+izgtUPlbdqPRJ4L7Jjc6Ml+q7Ms5tV3VO+m3Uj8o8F+VfXDqnpgVd2yqh5Ou3Rq7lUaJ9Iu6ZsrZzPafTKjzmydA9wiyeAOfRfaZRPTGHWEcT7zrhsr6UTgLkPL4y6Mrie0Sy52n/unj3cYbcf0hJ6Uzxf37mO6r+r29LmqeigtUfgpbVmP8kPaunu9CbHvRbuX6fhqZ8S+QzuiurqP+D4QOLaqfky7dPr43lD8iPYjds4daUduF9q4ZTDsKm44Cw9DBybG7RO46bo8zbKdtP7/cGiYue/zPR30i8BOWYCHMfSk4qXA72XNfjy+baht6Eqbrw0dcg7tR/SgxW4Hx82j1bUBbZmMOih2FrDLmMRy3vVgnGnX6f6j/lSG2rKVVMy/7xw0PB/uSjvbe1VVXUm71PYm22xfZ1JV/80N7VzRHjx4l4FBZ72d+2BV3Y+2XhfwT3O9RkzvGmCbgeW+Zd344TEr284tqGoHfZ5H29dNe2D1acBuSc5Lch7wZtoBtAV5qv6SJI9VdRntiPPbk/xxks2TrNfPBMydYt0o7fG0W/UfgpfTTj1D2yk/I+1db5sCr1qNcI4DLk/y0iSbJFk/yZ2T3HPSiEkeneS2fcc/F9+oJ18dA9w+yVPTLqvcMMk9k8wdhXgr8N2qeg7t9PY7h8b/04G6vhr4aI1/wtY7gZenv1su7ZHLT+zf75nkXv1Iw1X0h7tMmNfvpp1JuVeazZI8qidX36JdN/2CJBsk+SPaZSALLsmDaEdKnlDtyV7D/e+SZOMkmyZ5MS0BObz3/gTtcsgnJNmYdm/OD2vEJZlVdRbtfqx/7OXdhfYwhSOGhx3jfNoP3WkvT5i0bgzXc/1ehw2A9XqMc0eOvkJbbi9Ie33F3Nm3L/VxD0iyS1+Ou9Ierf3FgeLfQWsAHlOjL8UY9B7gxUnu0cu7bS9zlbantPeGPbYn9tfQzqKOW8e/ANy9z4dpYv8OcP++fyHJ3WiXnv5wTCwb9LLXB9bv83iDoWFuRbsv5JDe6XTaU1U3p10Odlof7ma0o4NfmK/+q+g9wD8kuV1fBndJcssRw50A/FHfNm5LW5/n6jFyn9B7n087CDNnlfeVA/4D+MO0V8psSLv06hu9TRipqn5Gu5zqyLTHkW/Ul8n+GX9Uf6ye4L+Jth8YNFfu3Gf9EaPPBNtQ29CVNakNHfIZ2vz+kx7Xk2kP6DhmyskN7zsmmW8e3URf/hvTfrNuMLi9JvmjJHv07WFb2o/k71c7+zbsOFoC//o+zY1zw2WgY9eD+azEOg1tPj9waPyNaZdYAtxssJ0b2tYfBvwp7eF+o+KYNB++Azynb7Ob0A6G/mCojI1pZ6n/snc6Hdg37ffNfentXPdA2j1+C23cb41hJwCPTHsFya1pZxoB6PPhQWnt8a9pZwYH27kVaVcuUe0y9M8Db0qyZZ9/uye50XKa4APAvZM8pK+XL6IdrDppJcqYV//9+jnaw4zmlXYVxO60fcte/XNn2i08C3Ppaq3CjZKr+qGdAj+Odsr6QtrZgANpl5ZtRLus7lLaBvgdBp5KRbtM7DzGPylu8MbZGz1chHbD7bUD/+9Au+/mvD69Yxm6WXVg2OvLpm1QZ9AakbOBV85T1z1ojdrcE72+1Bfg42hH8+ae/rY57ZKHA/r/X+GGJ8VdTrs2fJu64abfYuihCbTT0z/qw58FvLd3fzDtB/OVtBX5iD69SfN6v97tMtqO9iPAFr3f3rQnTc49Ke4oVu5m/9uOmrcjxv0yrZG9cuDz2YH+/9zjv5K2A7vt0PgPoZ3NurrP0xXzLKudaI3kJbQjgweNq8OIcTfqy/kS2j1JN6kXN715fOS6Mc88rKHP4QP970a7kf5q2j0Ngw8FeS1tPZ1bXw/lhgcgzB2N+/XQPB75xNM+zkG0p5ZdSbuE4m6TtifGPzBne9p9hL/s69lXgDvNM+2PAE+eNnbaZayn0NbT07jxQ4ReMbQuHTJiHh8yNP33c+Onde5M239dCrxpoPsTGfMQgQXYJtan3cdyeq/Xd2jvH71RObSji5/vw/xvr9/cA3NG7hN6v9txw9MPj16ZZTvfB/gz2j7vUtr+bOcpxgntPskTae3FL2j7mj2HhrvJ/GL001Y37fV9zMD8Gv7c5Gmxs/bBNnQvbEMXpA0dMfz9aG3JL/vf+00znXn2HWcMrhPc9OmRY+fRmHVoeHt9Ru/3/2j7xKto6+OHGHhK6oiydqGdKb64L9N/nbQeDM97Vn2dvjNtn5ahcm/0Gej39b48LqclevsPlXclNzy0ad75QHvd1H/1el9CW39vN1Teq2n3g879vxWtLfklLfGYewLpPWmJ6bh67stNH5gzdl0YMf643xrXl0N72NFRfd78sC+HuQfm3IW2/V/R63oMNzw855bAN2jb7vcG6vmOvvx+Sds+9x+1Dc4T8x/R9kOX0/ZB8z4EcWjcG82fcfOI9oCpq2hnlPelPTfiyqHPE2gHQj42Yjr70A7W32Jl6jbqk16AZkTai+A/UFXvWe5YpFmRdk/h+2ivb5jZnVaSbwPPrnZpq6QlZhuqWZbkg7R72Y5e7lhWVZKPAYdVuw9T66BR131L0kypqp/QjnbOtKq613LHIEmaTVX1J8sdw+qqqicsdwxaXiaPkqQlleT+jLlfpqrmfXqkJEmzLu11GT8Z0/tO1V63s0byslVJkiRJ0kRL9aoOSZIkSdIazORxmSQ5I+1l5uuEJF9J8pz+/YAkn1/umGZJkvsnOXkZprtHku8nuSLJC5Z6+muCtNedXJkZfpXDYkryjCTfGPi/0l4BIi0a20jbyEG2kbPLNnLdayNNHhfI4MqS5JAkH1jumGZVVR1RVQ9b7jhmSVV9var2WIZJv4T2SP4tqupfh3smeVKSbyb5VX+K4bz6e8LOTHJVkqNz05d3Dw67Rvw4rKqfV3tB+bh3d10vyYq+L/B+cmmAbeT0bCNvyjZydtlGrntMHqWVtJbt9HalvXdqnEuAf6G9OHhe/eXK76K9K2s72rvo3r6qgc3CfJ6FGFbFmhq3pDXfWrb/sY2c8RhWxZoa98xYlZdD+hn5ks+ivUh5P+A3wG9pL+z8wZjhzwBeTHu56S9pLzvduPe7Oe2lphfSXmR6DP1l4HXDiz1Po70A9XTGvNid9kLQb3HDS3jfBmw0ZtgVvQ7PpL0c91Lai1rv2WO8DHjb0DjPAk7qw36OG7+Q9qHAT3vd3kZ7IfxzBuIffPHxW/s0L6e9nPj+A/0OAT5Me0n7FbSd+N7zLIc7AF+g7dBPBp7Uu+/eu929/78D7QXB+/b/v8INL5b+JfBJbniR6ty8eTbwc+Br89Wf9oLztwAX9LJ+CNy593sk7elbV9BedP3i3n1fbvxS3Tv2mC7rdX7sQL/DgX+nvUD7CtqLwnefZ548tpdxWS/zjr37l4DrgF/T1tXbz1PGjV4aPmaY1wEfHPh/d9q2cJMXPwP/SXvB7dV92i+ZZz5/hPbS418CX2Pg5bvzzYsJy2ET4E3Amb3fN3q3m8TA0IvFJ6wrP+/Dzr2w9z60/cJX+7AXAUdN2AYPpL3I/Vzg4IH+69Fe9H4q7UXPH2bCOjpiGo+jvcz78l7Ofr37VsBhfZq/AF7DDS+EfgbzvKjcj59pPthG2kbaRtpG2kauFZ9lD2Bt+QyuLLSd+QcmDH9G37B2AG5B28Ee1PvdEngCsCmwRd8xHN37bdZX6j36/9sP7iiGpnEP4N60V7Ks6NN40Zhh5zasdwIbAw+j7TCPBm4F7EjbwTywD/944BTaDnwD4G+Bb/Z+2/QY/xjYEPhL4FrGN4x/2uu8AXAwbSe48cC8/DWtQVmftkM6dkwdNqM1sM/sZd2dtiPas/d/bp8Hm9IasjcOjPuVvkO4cy/nY3PLcGDevL/322RC/R9Oa+C3pu2c7whs3/udS2/4aT+A5hrqfekNY59npwCvADYCHkTb6c8t88Npjfw+fdpHAB8aM09uD1xF+6GyIa0BOoX+A6nX+zlTrN/TNIyfBF461O1K4B7zbAMPGbEOXj+fe/dn0baDm9GO8J4wMM7YeTFhOfx7r/uOtPXq93v5o5b1XLfBhnHSurLBQIxHAn9Da9g2Bu43YRs8spf7e7Qfxw/p/V8EHAvs1GN9F3DkfPNuqPx9aI3zQ3ssOwJ36P2O7uVtRtvejwOeN2Z7XesbRj8L/8E20jbSNtI20jZyrfgsewBry4dVaxj/dOD/NwDvHDPsXsCl/ftmtKNjTxi18k+Y5ouAT4zpN7dh7TjQ7WLgyQP/f4zesNLe0fbsgX7r0S7B2BV4GgONF22ndDZjGsYRsVwK3HVgXv7PQL87AVePGe/JwNeHur0LeNXA/58CfkQ7wnazge5fAV4/NJ3f0Haac/Nmt4H+89X/QcD/0X6UrDcUz8+B5wFbDnXflxsaxvvTfhysN9D/SOCQ/v1w4D0D/R4J/HTMPHkl8OGhOH/BjY8mL1TD+EX6j7uBbtdPa8w2MKph3G2eaWzdh9lq0rwYtxz6PLh6bh0bsx3sNqLbYMM4aV0ZbBjfDxzKwJmRCdvgHQa6vQE4rH8/CXjwQL/taWdvNmC6efcu4C0jum8HXMPA/gR4CvDlUdsr60DD6GfhP9hG2kbaRtpG2kauFR/veVxe5w18/xWwOUCSTZO8q99UfTntsoCtk6xfVVfRGoCDgHOTfDrJHUYVnuT2SY5Jcl4v53W0I57zOX/g+9Uj/p97gfeuwFuTXJbkMtqRrdCO1OxAO7oJQLWt6SzGSHJwkpOS/LKXtdVQnMPzaeMx16vvCtxrLqZe1gHArQeGeTftaNi/VdU1Q+MPxngm7SjkNmP6j61/VX2JdhnSvwPnJzk0yZZ9vCfQdt5nJvlqkvuMqMcOwFlV9buheHYc+H/kujOmrDPn/ullnjVU1kK5EthyqNuWtCPCK+P6+Zxk/SSvT3JqX4fP6L3mWz82B5hnOWxDO7p56jQxTNF/1Loy6CW0deO4JCcmedZKlr1D/74r8ImBde4k2iVV200Z986MrvOuPf5zB8p+F+3oqrScbCOxjRxRD9vIzjbSNnI5mDwujlrN8Q8G9gDuVVVbAg/o3QNQVZ+rqofSjqr8lLazH+Udvf/tejmvmCtjAZxFO2W/9cBnk6r6Ju2yk53nBkySwf8HJbk/8FLgScDNq2pr2mUDqxLnWcBXh2LavKr+rE9rc9olHYcBh4x4ytlgjLvQjlhdNNBtcLnOV3+q6l+r6h7AnrTLYv66d/9OVT2OttM5mnZN/rBzgJ2TDG6fu9COUK6sc2g7PuBGy2JVyprkROCuA9PajXbpyP+NGX7cdjLY/U9o9yA8hPaDacVc8dMENGY5XES7zGv3+UadUPS4deUm41XVeVX13KragXZE/e2Z/zHew2Wf07+fBTxiaJ3buKoGl+V8cZ/F6DqfRTuqus1AuVtW1Z7zlCWtDttI20jbSNvIuRhsI9cwJo+L43xgxdCObWVsQTuCeVnfeb9qrkeS7ZI8NslmtJX5StqRlXHlXA5c2Y+8/tkqxjPKO4GX96eHkWSrJE/s/T4N7Jnkj/rRzxdw4yObwzFeS7tufYMkf8dNj8xN6xjg9kmemmTD/rlnkjv2/m8FvltVz+kxvnNo/D9NcqckmwKvBj5a4x89Pbb+fZr3SrIh7V6KXwPXJdko7f1dW1XVb2nLZlT53+7jvaTXYV/gMcCHVmGefBh4VJIH93gOpq0335xm5H5Uc2PaZR/rJdm4lzPKEcBj0t7HtRltHn68qsYdVT0f2G1CCFv0eC+m3Yfzumni7rGPXA79yPJ7gTcn2aHX8T5JbjZt2YxfVy6kPeTg+noleWKSnfq/l9Iar/keaf7KfmZlT9q9SUf17u8EXptk117utkketxIxHwY8s68L6yXZMckdqupc4PPAm5Js2fvtnuSBK1G2tDJsI20jbSNtI+fisI1cw5g8Lo6P9L8XJ/neKoz/L7SbkC+i3fz73wP91qPt3M6hXQbyQODPx5TzYtpRqStoR16PGjPcSquqTwD/BHwo7VKJHwOP6P0uAp5Ie3T1xcDtgP8dU9TnaPdG/B/t8oNfM/lyiHExXUF7iMH+tPlzXo/xZn0Hsh/tUiaAvwLunuSAgSL+k3Z/wHm0SzbGvhB4vvrTGvZ303aCZ9LmwRt7v6cCZ/RxDqI9CGG47N/Qnv72CNo68HbgaVX10ylnxWBZJ/dp/Fsv6zHAY/o0pvFU2o+0d9DuM7magaP4aS8Gvn+f1om9TkfQHhyxBePXTWgPdvjbtMtAXjxmmPfT5uEvaE/gO3bKuGH+5fBi2n0936FtR//Eyu0PR64rVfUr4LXA//Z63Zv2NMZvJ7mSdj/RC6vq9HnK/irtgQ1fpD2wYu5l4W/t438+yRW0eXGvaQOuquNoDe1baGcuvsoNR9yfRnvwxE9o8+ujtLM20mKwjbSNtI20jbSNXEOlXWovrdvSXu77gap6z3LHotm2WOtKkhW01wpsWFXXLmTZkrQ6bCM1LdvItZ9nHiVJkiRJE5k8SpIkSZIm8rJVSZIkSdJEnnmUJEmSJE1k8rgOSHJ4ktcsdxyzKEll/vcJaRUleUWSRX+4wlJNZ1pJ9k1y9hJM54wkD+nfD0nygcWeprS2sX0cz/Zx8dg+Lvp0bB8XkcnjGi7JLv1R0HOfSnLVwP/3X+4Yl8ss/ChI8owk35gwzJOSfDPJr/pTyuYbdkl2vAuhql7X3xc20ers3FdmOpLWHbaP49k+Li/bR63JTB7XcFX186rafO7TO991oNvXF3J6SdZfyPKGyt5gscqecZfQ3lv2+oUobBbm4yzEIGndZvu4VrB9lGaMyeO64+ZJPp3kiiTfTrL7XI8kd0jyhSSXJDk5yZMG+h2e5B1JPpPkKuAPkuyQ5GNJLkxyepKxLwpO8qgk309yeZKzkhwy0G9FPxL87CQ/B77Uuz8ryUlJLk3yuSS7zlP+R5Kcl+SXSb6WZM/e/UDgAOAl/Qjzf02aQUluluSNSX6e5Pwk70yySe938yTH9Dpf2r/vNDDuM5Kc1ufv6UkOSHJH4J3AfXoMl42ablX9T1V9mPbS5vni24z2sugdBo6c79CPSn40yQfSXqz8jCT7JPlW2kt4z03ytiQbDZRVSQ5K8rNen39Pkt7vtkm+2ufpRUmOGhhvz4F15fwkr+jdR8Vw/dHSgWV9YJJzekwH9377Aa8Antzr9INx83TMfBk1naf35XhRkr9ZpGV+iyT/0etzaZKjh8o+OMkFva7PnCeGseUkeXSSE/py/GaSu4wrR9Iqu3lsH+e1mvtK20fbx6OHyrZ9XJNVlZ+16AMUcNuhbofTjt7tA2wAHAF8qPfbDDgLeGbvd3fgImDPgXF/CdyXdrBhU+C7wN8BGwG7AacBDx8Tz77A7/Vx7wKcDzy+91vR431/j2MT4PHAKcAdezx/C3xznvo+C9gCuBnt6OQJQ/V+zbTzq4//KeAWvcz/Av6x97sl8IRe/y2AjwBHD8zDy4E9+v/bD8y/ZwDfmHLZPQf4yoRh9gXOHup2CPDbPu/W6/PxHsC9+zxcAZwEvGio3scAWwO7ABcC+/V+RwJ/08vaGLhf774FcC5wcO++BXCveWI4hPai4MFlfWSfX7/Xp/mQgfE/MBDf2Hk6Yp6Mms67ewx3Ba4B7jhm3FVa5r3/p4GjgJsDGwIPHFhG1wKv7t0fCfwKuPmYGMaVc3fgAuBewPrA04EzgJv1/meMm39+/Pi56QfbxxOG6m37aPto++hnpT/LHoCfBV6g4xvH9wz8/0jgp/37k4GvDw3/LuBVA+O+f6DfvYCfDw3/cuA/pozvX4C39O9zO7LdBvp/Fnj2wP/r9R3LrlOUvXUvb6uB2KdqHIEAVwG7D/S7D3D6mPH2Ai7t3zcDLus70k2GhnsGS9M4fm3CeC8CPjFU7/sN/P9h4GX9+/uBQ4Gdhsp4CvD9MeXfJAZGN1p3GOj/BuCw4WEnzdMx0x6ezk4D/Y8D9h8x3uos8+2B3zGiwevL6Gpgg4FuFwD3HjHsfOW8A/iHoW4nc0PjeQY2jn78TP3B9tH2cfR4L8L2cXg820c/Yz9etrruOG/g+6+Aufs/dgXu1U/7X9YvHTkAuPXA8GcNfN+VdlnI4PCvALYbNdEk90ry5X5pwy+Bg4BthgYbLv+tA2VfQtuJ7Tii7PWTvD7Jqf1SkDN6r+Hyp7Et/ajxwLT/u3cnyaZJ3pXkzD6trwFbJ1m/qq6i/cg4CDg37fKnO6xCDKtjcB6S5Pb9MpLzeryv46bzZdw68RLaPD8uyYlJntW77wycOm0MUwxzJrDDqIEWYJ6Oq9ugVV7mtHlxSVVdOmb6F1fVtVPEMF85uwIHD21rOzNmnklaZbaP87N9tH28DNtHdSaPOgv4alVtPfDZvKr+bGCYGhr+9KHht6iqR44p/4O0yx52rqqtaPc4ZGiY4fKfN1T+JlX1zRFl/wnwOOAhwFa0o2oMlF8jxhnnItrRsD0HprtV3fCQhYOBPWiXoWwJPGBwWlX1uap6KO1I2U9pl4WsbAzTGFfecPd39Dhu1+N9BTed76MLqjqvqp5bVTsAzwPenva49rOA3ecbdYridx74vgs33Mdyk3HnmacLZXWW+VnALZJsvZoxzFfOWcBrh7aFTavqyNWcpqTp2D42to9zBdk+2j7K5FEcA9w+yVOTbNg/90y7mX2U44DLk7w0ySb96Oadk9xzzPBb0I4c/TrJPrQGbT7vBF6eG27s3yrJE+cp+xrgYtoRstcN9T+fds/JRFX1O9rO9y1JbtWnvWOShw9M62rgsiS3AF41N26S7ZI8Nu2G/WuAK4HrBmLYKQM34w/r83Bj2v0X6yXZOMmGYwY/H7hlkq0mVGkL2j0RV/Yjkn82YfjBeJ6YG258v5TWcF1HW1duneRFaTfSb5HkXtOW272yH7Hck3Yf0dzDBs4HViRZr8cw3zxdEKuzzKvqXNolZG9Pe3DAhkkewEqaUM67gYP62Ykk2SztARtbrGqdJa0U20dsH4fisX1sbB/XYSaP67iqugJ4GLA/7SjXecA/0W6wHzX8dcBjaNe3n047OvUe2pHNUf4ceHWSK2gPEfjwhHg+0af/obRLIX4MPGLM4O+nXdrxC+AnwLFD/Q8D7pR2ScPR8023eyntYQTH9mn/D+3IGrR7UTah1fdY2uUbc9ajHYU7h3YZ0QNp9Yb2hLwTgfOSXDRmuk+l7YTfAdy/fx95FLGqfkq7qf60Xq9xl2i8mPZD5Ipe1lFjhhvlnsC3k1xJOyr+wqo6va8rD6Ut//OAnwF/sBLlAnyVNo+/CLyxqj7fu3+k/704yfeYf54upFVd5tCW229pR30voN03sypGllNVxwPPBd5G+5FyCu0eIUlLwPbxRmwfG9vH5l+wfVxnpWqhrxqQpBtLsoL2Y2rDoXsdJElaZ9k+ak3jmUdJkiRJ0kQmj5IkSZKkibxsVZIkSZI0kWceJUmSJEkTmTxKkiRJkiYyeZQkSZIkTWTyKEmSJEmayORRkiRJkjSRyaMkSZIkaSKTR0mSJEnSRCaPkiRJkqSJTB4lSZIkSROZPEqSJEmSJjJ5lCRJkiRNZPIoSZIkSZrI5FGSJEmSNJHJoyRJkiRpIpNHSZIkSdJEJo+SJEmSpIlMHiVJkiRJE5k8SpIkSZImMnmUJEmSJE1k8ihJkiRJmsjkUZIkSZI0kcmjJEmSJGkik0dJkiRJ0kQmj5IkSZKkiUweJUmSJEkTmTxKkiRJkiYyeZQkSZIkTWTyKEmSJEmayORRkiRJkjSRyaMkSZIkaSKTR0mSJEnSRCaPkiRJkqSJTB4lSZIkSROZPEqSJEmSJjJ5lCRJkiRNZPIoSZIkSZro/wOtQZCaSWIN9QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1080x360 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_expression_distribution(gene='Gsdme', cluster_num=6)\n",
    "#     Returns plot of distribution of gene transcripts levels for 'gene' in 'cluster' \n",
    "#         Note - this function takes a cluster number (so it can plot the CTRL and TREAT)\n",
    "#         as opposed to a particular cluster (e.g., '6_CTRL')'''"
   ]
  }
 ],
 "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.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}