B_cells_versus_T_cells_FP.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## B cells versus T cells\n",
    "The data used in this notebook is lymphocyte data for one patient's B cells and T cells. We use this data to show the proteome variation between the cell types. Here, we note presence/absense changes and quantify fold changes using a T test.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "First, we import our loader module. This brings the functions defined there (in our repository at ~/load_data.py) into scope so we can use them in this script. Then we can load our data and store it as <code>data</code>.\n",
    "\n",
    "Calling <code>head</code> shows the first several lines of the dataframe, which provides an idea of the type of data present and the structure of the dataframe."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Subject1_FACS145_B_022920_C10</th>\n",
       "      <th>Subject1_FACS145_B_022920_C11</th>\n",
       "      <th>Subject1_FACS145_B_022920_C12</th>\n",
       "      <th>Subject1_FACS145_B_022920_C13</th>\n",
       "      <th>Subject1_FACS145_B_022920_C9</th>\n",
       "      <th>Subject1_FACS145_T_022920_D10</th>\n",
       "      <th>Subject1_FACS145_T_022920_D12</th>\n",
       "      <th>Subject1_FACS145_T_022920_D13</th>\n",
       "      <th>Subject1_FACS145_T_022920_D9</th>\n",
       "      <th>Subject1_FACS145_T_062920_D11</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Protein ID</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>A0A0B4J2D5</th>\n",
       "      <td>819194.25</td>\n",
       "      <td>1203040.90</td>\n",
       "      <td>439446.10</td>\n",
       "      <td>1906699.40</td>\n",
       "      <td>1330848.20</td>\n",
       "      <td>283948.94</td>\n",
       "      <td>78067.72</td>\n",
       "      <td>0.00</td>\n",
       "      <td>123710.484</td>\n",
       "      <td>181092.27</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>A0AVT1</th>\n",
       "      <td>501930.78</td>\n",
       "      <td>718435.75</td>\n",
       "      <td>824289.10</td>\n",
       "      <td>495488.97</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>401786.880</td>\n",
       "      <td>559705.25</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>A6NHR9</th>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.000</td>\n",
       "      <td>0.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>A8MWD9</th>\n",
       "      <td>738986.20</td>\n",
       "      <td>971444.40</td>\n",
       "      <td>789952.70</td>\n",
       "      <td>960253.10</td>\n",
       "      <td>319829.47</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>344603.060</td>\n",
       "      <td>394223.78</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>B5ME19</th>\n",
       "      <td>785554.06</td>\n",
       "      <td>1040235.90</td>\n",
       "      <td>632014.44</td>\n",
       "      <td>709235.70</td>\n",
       "      <td>362616.84</td>\n",
       "      <td>403412.84</td>\n",
       "      <td>674793.25</td>\n",
       "      <td>974256.94</td>\n",
       "      <td>720808.750</td>\n",
       "      <td>736337.70</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q8WWY7</th>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.000</td>\n",
       "      <td>0.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q92597</th>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.000</td>\n",
       "      <td>0.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q92953</th>\n",
       "      <td>7166249.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>6292137.50</td>\n",
       "      <td>0.00</td>\n",
       "      <td>9075473.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>5203961.000</td>\n",
       "      <td>1345615.80</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9BWT6</th>\n",
       "      <td>0.00</td>\n",
       "      <td>2075055.20</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>2112621.80</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.000</td>\n",
       "      <td>0.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q15388</th>\n",
       "      <td>0.00</td>\n",
       "      <td>788927.20</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>409172.50</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.000</td>\n",
       "      <td>535433.44</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1564 rows × 10 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "            Subject1_FACS145_B_022920_C10  Subject1_FACS145_B_022920_C11  \\\n",
       "Protein ID                                                                 \n",
       "A0A0B4J2D5                      819194.25                     1203040.90   \n",
       "A0AVT1                          501930.78                      718435.75   \n",
       "A6NHR9                               0.00                           0.00   \n",
       "A8MWD9                          738986.20                      971444.40   \n",
       "B5ME19                          785554.06                     1040235.90   \n",
       "...                                   ...                            ...   \n",
       "Q8WWY7                               0.00                           0.00   \n",
       "Q92597                               0.00                           0.00   \n",
       "Q92953                         7166249.00                           0.00   \n",
       "Q9BWT6                               0.00                     2075055.20   \n",
       "Q15388                               0.00                      788927.20   \n",
       "\n",
       "            Subject1_FACS145_B_022920_C12  Subject1_FACS145_B_022920_C13  \\\n",
       "Protein ID                                                                 \n",
       "A0A0B4J2D5                      439446.10                     1906699.40   \n",
       "A0AVT1                          824289.10                      495488.97   \n",
       "A6NHR9                               0.00                           0.00   \n",
       "A8MWD9                          789952.70                      960253.10   \n",
       "B5ME19                          632014.44                      709235.70   \n",
       "...                                   ...                            ...   \n",
       "Q8WWY7                               0.00                           0.00   \n",
       "Q92597                               0.00                           0.00   \n",
       "Q92953                         6292137.50                           0.00   \n",
       "Q9BWT6                               0.00                           0.00   \n",
       "Q15388                               0.00                           0.00   \n",
       "\n",
       "            Subject1_FACS145_B_022920_C9  Subject1_FACS145_T_022920_D10  \\\n",
       "Protein ID                                                                \n",
       "A0A0B4J2D5                    1330848.20                      283948.94   \n",
       "A0AVT1                              0.00                           0.00   \n",
       "A6NHR9                              0.00                           0.00   \n",
       "A8MWD9                         319829.47                           0.00   \n",
       "B5ME19                         362616.84                      403412.84   \n",
       "...                                  ...                            ...   \n",
       "Q8WWY7                              0.00                           0.00   \n",
       "Q92597                              0.00                           0.00   \n",
       "Q92953                        9075473.00                           0.00   \n",
       "Q9BWT6                              0.00                           0.00   \n",
       "Q15388                              0.00                           0.00   \n",
       "\n",
       "            Subject1_FACS145_T_022920_D12  Subject1_FACS145_T_022920_D13  \\\n",
       "Protein ID                                                                 \n",
       "A0A0B4J2D5                       78067.72                           0.00   \n",
       "A0AVT1                               0.00                           0.00   \n",
       "A6NHR9                               0.00                           0.00   \n",
       "A8MWD9                               0.00                           0.00   \n",
       "B5ME19                          674793.25                      974256.94   \n",
       "...                                   ...                            ...   \n",
       "Q8WWY7                               0.00                           0.00   \n",
       "Q92597                               0.00                           0.00   \n",
       "Q92953                               0.00                           0.00   \n",
       "Q9BWT6                         2112621.80                           0.00   \n",
       "Q15388                          409172.50                           0.00   \n",
       "\n",
       "            Subject1_FACS145_T_022920_D9  Subject1_FACS145_T_062920_D11  \n",
       "Protein ID                                                               \n",
       "A0A0B4J2D5                    123710.484                      181092.27  \n",
       "A0AVT1                        401786.880                      559705.25  \n",
       "A6NHR9                             0.000                           0.00  \n",
       "A8MWD9                        344603.060                      394223.78  \n",
       "B5ME19                        720808.750                      736337.70  \n",
       "...                                  ...                            ...  \n",
       "Q8WWY7                             0.000                           0.00  \n",
       "Q92597                             0.000                           0.00  \n",
       "Q92953                       5203961.000                     1345615.80  \n",
       "Q9BWT6                             0.000                           0.00  \n",
       "Q15388                             0.000                      535433.44  \n",
       "\n",
       "[1564 rows x 10 columns]"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import load_data\n",
    "data_raw = load_data.load_FragPipe(version=\"June_selected_FP\")\n",
    "#data_raw=load_data.load_max_quant()\n",
    "\n",
    "data_raw"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "These are some methods that assist with the data analysis, such as checking if a value was measured or not and if the protein was measured in each type."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import data_utils"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Before any analysis, we log normalize and median normalize across runs."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "data_log2_normalized = data_utils.normalize(data_raw)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Our next steps will require separating the runs by cell type, to compare B cells to T cells. As you see above, the pattern for naming the columns is \"Intensity \", B or T indicating cell type, and the replicate number. Here we define the list of cell types that the functions should look for based on the naming pattern."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "cell_types=['B_', 'T_']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here we record which proteins are only identified in one cell type. We'll set these aside for now, in the <code>data_in_one_type</code> dataframe. Note that <code>min_reps</code> is the minimum number of replicates required to list a protein as present in one type. This can be changed, i.e., <code>min_reps=5</code> to determine the proteins identified in all 5 replicates of a single type; do not set it to a number greater than the replicates run or no proteins will be returned."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Proteins identified only in:\n",
      "B_    69\n",
      "T_    56\n",
      "dtype: int64\n"
     ]
    }
   ],
   "source": [
    "present_only_in = data_log2_normalized.apply(data_utils.check_presence_absence, axis=1, \n",
    "                                             min_reps=3, cell_types=cell_types).dropna()\n",
    "print(\"Proteins identified only in:\")\n",
    "print(present_only_in.value_counts())\n",
    "data_in_one_type =data_log2_normalized.loc[present_only_in.index]\n",
    "data_in_one_type['Type'] = present_only_in"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, we select the proteins that are measured in at least three samples from each group, allowing the calculations to proceed without imputed zero-handling."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "indecies = data_log2_normalized.apply(data_utils.check_three_of_each_type, axis=1, cell_types=cell_types)\n",
    "data = data_log2_normalized[indecies]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Total groups identified:\n",
      "1564\n",
      "Proteins identified in at least 3 of each cell type:\t\n",
      "904\n"
     ]
    }
   ],
   "source": [
    "print (\"Total groups identified:\")\n",
    "print(data_raw.shape[0])\n",
    "\n",
    "print (\"Proteins identified in at least 3 of each cell type:\\t\")\n",
    "print (data.shape[0])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we calulate p values using a ttest. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "D:\\Program_Files\\Anaconda\\lib\\site-packages\\scipy\\stats\\_distn_infrastructure.py:1932: RuntimeWarning: invalid value encountered in less_equal\n",
      "  cond2 = cond0 & (x <= _a)\n"
     ]
    }
   ],
   "source": [
    "#Call ttest function\n",
    "tstats = data.apply(data_utils.ttest_wrapper, axis=1, cell_types=cell_types)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here we calculate fold change."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "tstats['log2(B)-log2(T)'] = data.apply(data_utils.get_fold_changes, axis=1, cell_types=cell_types)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Fold change and p values allow us to create a volcano plot."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAEWCAYAAABsY4yMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO29eXxU9fX//zyzJZMFEkNYFFxqlRptFIKy2F+r1VoXKkUQq4LbRwFta9uPW/vph27Uz0+K1i4WAWtRcakK+mlLa7WL1k9FbQkoVZBaRAuoEDGBrDOZmff3j1mY5c5kQubOTIbzfDzyIDNz595zZ8jrnnvOeZ8jxhgURVGU0sNRaAMURVEUe1CBVxRFKVFU4BVFUUoUFXhFUZQSRQVeURSlRFGBVxRFKVFU4JWSQkTeFpEzI79/R0QeLLRNdiEiV4jIX+MeGxH5aCFtUooLFXglL0SEt1tEOkSkVUR+KyJjCm2XHYjI0sh5doiIX0R64x4/VWj7lIMHFXgln3zOGFMFjAJ2AT8tsD22YIyZb4ypipzr/wCPRh8bY84ptH3KwYMKvJJ3jDE9wCqgId02InKIiKwQkXcjHv//xr02VUReEZE2EVkrIo19HVNEykXkQRHZE3nf30VkhMV2XxeRVUnP/VhEfhL5/QoReUtE2kVkm4hc2p9zz2DftMg57RORrSJyduT5oSJyr4i8JyI7ReT7IuLMYn/nisimiJ07ReTGXNipDC5chTZAOfgQkQrgIuClDJutBDqA4yP/Tom8dzzwC+BzwDpgNvBrERlrjPFl2N/lwFBgDOADTgK6LbZ7BPiWiAwxxuyLiOksYLqIVAI/AU42xmwRkVHAIVmedlpE5BTgAWAm8CfCdzjVkZfvJ3y381GgElgDbAeW9bHbe4FZxpj/E5Fa4KiB2qkMPtSDV/LJ/4pIG7AP+Ayw2GqjiHCeA8w3xrQaY3qNMX+JvHwNsMwY87IxJmiMuZ+wYE/q49i9QB3w0cj7mo0x+5I3Msa8A6wHPh956tNAlzEmejEKASeIiNcY854x5vVsTz4D/wH8whjzB2NMyBiz0xjzRuQO4xzgq8aYTmPMbuBO4AtZ7LMXaIhcqFqNMetzYKcyyFCBV/LJ540xNUAZ8CXgLyIy0mK7McCHxphWi9eOAG6IhFnaIheMMcChfRx7JfA08MtI2OcHIuJOs+3DwMWR3y+JPMYY00n4zmM+8F4kUfyxPo6bDWOArRbPHwG4I8eKnusyYHgW+5wBnAu8IyJ/EZHJObBTGWSowCt5J+JBPwEEgU9YbLIdOEREatK8dqsxpibup8IY80gfx+w1xnzXGNNAONwzFbgszeaPA6eJyGhgOhGBj+znaWPMZwiHUd4A7sl8tlmxHTg6zfM+YFjcuQ4xxhzf1w6NMX83xkwjfDH4X+CxHNipDDJU4JW8I2GmAbXA5uTXjTHvAU8BS0SkVkTcIvLJyMv3APNFZGJkP5Uicp6IVCfvJ+mYp4vIxyMx9X2EQxhBq22NMS3Ac8AKYJsxZnNkHyNE5PxILN5HODdguY9+ci9wpYicISIOETlMRD4W+RyeAe4QkSGR144WkU/1ca4eEblURIYaY3oj55sLO5VBhgq8kk9+IyIdhAXnVuDyDDHsOYRF+A1gN/BVAGPMOsJx+LuAVuBfwBVZHHsk4cqdfYQvKn8BMi2Cehg4kzjvnfDfyw3Au8CHwKeA67I4dkaMMX8DriQcX98bse2IyMuXAR5gE+HzXUX47qEv5gBvi8g+wiGl2QO1Uxl8iA78UBRFKU3Ug1cURSlRVOAVRVFKFBV4RVGUEkUFXlEUpUQpqlYFw4YNM0ceeWShzVAURRk0NDc3f2CMqbd6ragE/sgjj2TdunWFNkNRFGXQICLvpHtNQzSKoigligq8oihKiaICryiKUqKowCuKopQoKvCKoiglSlFV0SiKMrgJhQx7Ov34A0E8Lid1lR4cDim0WQctKvCKouSEUMiwZVc71zywjh2t3Yyu9XLPZRMYO6JaRb5AaIhGUZScsKfTHxN3gB2t3VzzwDr2dPoLbNnBiwq8oig5wR8IxsQ9yo7WbvwBnTVSKFTgFUXJCR6Xk9G13oTnRtd68bicBbJIUYFXFCUn1FV6uOeyCTGRj8bg6yo9KduGQoaWdh87W7toafcRCungITvQJKuiKDnB4RDGjqjmyetOzVhFo8nY/KEevKIoOcPhEOqryzistoL66jJLwdZkbP5QgVcUJa9oMjZ/qMAripJXNBmbP2wVeBGpEZFVIvKGiGwWkcl2Hk9RlOKnP8lYZWDYnWT9MfB7Y8xMEfEAFTYfT1GUIifbZKwycGwTeBEZAnwSuALAGOMHNIuiKEosGavYi50hmo8ALcAKEdkgIj8XkcrkjURkroisE5F1LS0tNpqjKIpycGGnwLuA8cDdxphxQCfw9eSNjDHLjTETjDET6ust58YqiqIoB4CdAr8D2GGMeTnyeBVhwVcURVHygG0Cb4x5H9guImMjT50BbLLreIqiKEoidlfRfBl4KFJB8xZwpc3HUxRFUSLYKvDGmFeACXYeQ1EURbFGV7IqiqKUKCrwiqIoJYoKvKIoSomiAq8oilKiqMAriqKUKDrRSVGUGKGQYU+nX5uAlQgq8IqiADpKrxTREI2iKICO0itF1INX+oXewpcuOkqv9FCBV7JGb+FLm+govXiR11F6gxsN0ShZo7fwpY2O0is91INXskZv4UsbHaVXeqjAK1mjt/Clj47SKy00RKNkjd7CK8rgQj14JWv0Fl5RBhcq8Eq/0Ft4RRk8aIhGURSlRFGBVxRFKVE0RKMcMLqqVVGKGxV45YDQVa2KUvxoiEY5IHRVq6IUP+rBH+QcaJhFV7UqSvFjq8CLyNtAOxAEAsaYCXYeT+kfAwmz6KpWRSl+8hGiOd0Yc5KKe/ExkDCLrmpVlOJHQzQHMQMJs+iqVkUpfuwWeAM8IyIGWGaMWZ68gYjMBeYCHH744f0+gJbqHTgDDbPoqlZFKW7sDtGcaowZD5wDfFFEPpm8gTFmuTFmgjFmQn19fb92Ho0hT1/yAqcuepbpS15gy652QiGTI/NLGw2zKEppI8bkRwxF5DtAhzHm9nTbTJgwwaxbty7rfba0+5i+5IUUD/TJ605VzzJL9A5IUQY3ItKcLsdpW4hGRCoBhzGmPfL7WcD3cnkMLdUbOHaEWfSioSjFgZ0x+BHAkyISPc7Dxpjf5/IAWqpXfOgKV0UpHmyLwRtj3jLGnBj5Od4Yc2uuj6Ex5OJDV7gqSvEwqMsktVSv+NCwmaIUD4Na4EFL9YoNDZspSvGgzcaUnKJhs8FJKGRoafexs7WLlnaflhqXCIPeg1eKCw2bDT40MV66qAev5Jxo2Oyw2grqq8tUJIocTYyXLirwinKQo4nx0kUFXilKNCacP6KJ8Xg0MV4aqMArRYf2GMovmhgvXfLWiyYb+tuLRilNtMdQ/tH2EoOXgvSiUZQDRWPC+UfXk5QmGqJRig6NCStKblCBV4oOjQkrSm7QEI2ScwYaz9XFUoqSG1TglZySq1WRGhNWlIGjIRolp+iqSEUpHlTglZyiFTCKUjyowCs5RStgFKV4UIFXcopWwChK8aBJViWnaAWMohQPKvBKztEKGEUpDlTgFWWQo31klHSowCsFQUUpN+g0JiUTtidZRcQpIhtEZI3dx1IGB9oOOHfougMlE/moovkKsDkPx1EKwIEM5lBRyh267kDJhK0CLyKjgfOAn9t5HKUwHKgnrqKUO3TdgZIJuz34HwE3A6F0G4jIXBFZJyLrWlpabDZHySUH6omrKOUOXXegZMK2JKuITAV2G2OaReS0dNsZY5YDyyE80ckue5TcEU2QdvkDB+SJR0UpOTGootR/dN2Bkgk7q2hOBc4XkXOBcmCIiDxojJlt4zEVm4mv2lgwtYHRtd6U0Xp9eeIqSrlF1x0o6bAtRGOM+YYxZrQx5kjgC8CfVdwHP/FhmaXPbWXRjMYDCg9ERemw2grqq8tU3BXFBrQOXukX8QnSDdvbuP3pLSyY2sBxI6vxelzqiStKEZGXZmPGmOeMMVPzcSzFXpITpBu2t7FwzSa8Hpd64opSZGg3SaVfaNXGwcmBrHdQCo+GaJR+oQnS0iRT6whthzB4UQ9e6TeaIC0t+lqwpiuPBy8q8ErRomGB/NCXgOvK48GLhmiUvJNNJ0kNC+SPvgQ8mljv73oHpfCoB6/klWz712hYIH/01TpCE+uDF/Xglbyyp9PPnX8I187XeN20dfdy5x+2cOv0xoTVmBoWyB99tY7QxPrgRQVeyZpcDOkIhUJcPuUoblm9MSYmi2Y0Egol9qPTsED+yEbAtR3C4ERDNEpW5GpIR9AQE3cIe+W3rN5IMGk3GhbIL1oZVZr06cGLSDlwHfAJwAB/Be42xvTYbJtSRKSLiT953an98uyMMZahF2MSFV7DAooycLIJ0TwAtAM/jTy+GFgJXGiXUUrxkYuYeChkEJGsQy8aFsgenXGrWJGNwI81xpwY9/hZEXnVLoOU4mSgMfFoiOfOP2xh0YzGhBi8hl4GhpaUKunIRuA3iMgkY8xLACIyEXjBXrOUYmOgQzriQzwt7X4WTG2grtLDoTVeRg4pVyEaALkKnymlRzYCPxG4TET+HXl8OLBZRP4BGGNMo23WKUXDQGPiyW2G561sBuCFW05XcR8gWlKqpCMbgT/bdiuUQcFAYuJa9mgf+tkq6eizTNIY806mn3wYqQx+tOzRPgr92WrPoOJFksvTCsmECRPMunXrCm2GYhNa6WEfhfpsNcFbeESk2Rgzweo1Xcl6kJNPYTiQEI9eFLKjUCWlmuAtblTgD2KK3fsqdvsUTfAWO33G4EWkXUT2Jf1sF5EnReQj+TBSsYdi79hY7PYpfXeiVApLNr1ofgjcBBwGjAZuBO4Bfgn8wj7TFLspdu+r2O1TCp/gVTKTVZmkMWZi3OPlIvKSMeZ7IvJfdhmm2E+xl9ela2sgouGZYkF7BhU32XjwIRGZJSKOyM+suNeKpwRH6TfF7n05BRbNaEywb9GMRpyqHUWFdqIsXrLx4C8FfgwsiTx+EZgtIl7gS+neFOlC+TxQFjnOKmPMtwdmrpJLit37cjgc3L92W8JwkPvXbuPW6aWzeDoUMnzQ6aOnN4hTBK/HSY23eL4DZXDTp8AbY94CPpfm5b9meKsP+LQxpkNE3MBfReSpaE8bpTgo5o6NdZUevvaZsVn3vxlsJZVWVUKLZzYyYkg5R9ZVFrXtB8Jg+35KgWz6wY8m3Cr4VPb3g/+KMWZHpveZ8AqqjshDd+RHQzpKn8QLQV2Vh19/6VS6/ZlFYTCWVFpVCd20aiMLp51Adbm7aC+8B8Jg/H5KgWxi8CuAXwOHEq6k+U3kuT4REaeIvALsBv5gjHnZYpu5IrJORNa1tLRkb7lSkiRPjrpgyVp27fMxaqg3Y3x3MJZUpqsSqvA4S65SaDB+P6VANgJfb4xZYYwJRH7uA+qz2bkxJmiMOYlweeUpInKCxTbLjTETjDET6uuz2q1SwhyoEAzGksp0NeRdkbuVUmIwfj+lQDYC/4GIzI54404RmQ3s6c9BjDFtwHNoZ0qlD/oSgnSNrQbjghurKqbFMxs5oq6iaCqZcsVg/H5KgWyqaK4C7gLuJBxDXxt5LiMiUg/0GmPaIhU3ZwKLBmCrUsTkKoGWqTY/Uxx3oANJCkG0iumJ66bQ0xvCKZRsFc1g/H5KAdu6SYpII3A/4CR8p/CYMeZ7md6j3STtwe7qhVwm0DLta0+nn+lLXkgR/0fnTsLjclLrddPa3atVGkWKVtHYwwF1kxSRn5Kh6sUYc32mgxpjNgLjsjVSsYd8VC/ksqNgptr8dOGbHa3d3PD4q1qVUeQUc0luqZIpBr8OaM7wowwC8lG9kOsEWrqVkeniuG3dvVqVoSgWpPXgjTH3Jz8nIiONMe/ba5KSS/JRvZAcNx83pobrzziGoAknRA/0Vjz5lr7W6+aBq07hnT1dVHicdPmD1Fa6+e6vN9lyXooy2OlvP/jfAePtMESxh3w0FItPoNVXlXHz2WO5adXGAYWErEJLD1x1Cr7eEAt+9VrsuTsuPBHI3UVFUUqJfiVZRWSDMca2uLomWXNPvlYQxnvbFy1/KeWCki4eny7x1tLuS0morrji5Ji4x+/7tgs+jtMhA76oKMpgJJcj++7JgT1KHslXQ7Fo3Hxna1fWIaFMFx+r0FKFx2m578MPqeCSn7+sY+MUJYl+CbwxZknfWynFRj6rF/oTEspUfWO1ny5/kLMahjOjaUysu+Tq5u04HUJ9VVlC18mlz2096OLxWoaoJKMzWZUErBKb/akt78+ClkwJ4FFDvSn7+ejwSq4/41jmP9gce27p7Caqyp0pcf/FMxvxekpzlaSVkAPazEtJwbaFTgeCxuALS3LI5KyG4SmCmo1oZOtJWsXZ4+P1yfsxGC5YsjZl+8fmTWbWshdTnn/iuikMry7P0adTHKQLa40YUsb5d6X/LJXSJVMMPpteNMpBQnLIZEbTmJi4Q/Y19NlO+OlrolTyfnoDIUuPPxC0fr43EOr/h1DkpAtrdfu1mZeSioZolBjJIZMar9tW0ehvAjhdfN/pEFvq8IuRdGGtoKGo5+sqhUE9eCVG8krRtu5e2zsA9meep5XHf8eFJ/LA2m0snd3E6Fov48bUcPPZY1nwq9f45A+eY/qSF9iyqz3WdTLfpOt+eaCkW81b7nYU9XxdpTBoDF6JkasYfC7ssIrhh0KGt/d0JqxkPbSmHH8gRKc/yEeGVdIbDPWrDt9OMpWBApaJ0r5yF/3dZ/z7tcqmNMkUg1eBVxIYaBVNLo7f326SC6Y2sHDNJpbNbmJohZtPLHo2Zb8v3HI6h9VW2Ga3FemSyL/+0qns2udLXaUbCGVVBXMgQq0j80oXTbIqWZMcMnG5HFmHUA6E5BBGW3f65mjp4s/RXMG8B8M98IplsEQ6e7v9wZRzfGdPV9ZN4foT1oqiI/MOTjTJqhQMK6/yvitP5rYLPo7b6YgtWNqwvS3mrVolEtu6e4GwaO3t7mXRjEZuWb2/Jn7p7CZqve68n1/U3vqqMuafdjQ1Xjdd/iAOB1mv0s1VQltH5h2cqMArBSPZq6yvKqOl3cfXn/hHTJwXzWjk/rXbYqGIZbObmBeXE1g0o5Hbn94C7K8iWfrc1oRVrSLQ2t2b9xh8XaWHB646hV37ehIWYS2b08RZDcN5ZtPu2LZd/qCtVTD5aDqnFB8aolEKRrJXOf+0o2NCCGEP85bVG/nv8xpiceZRNeUsnHYCT143hQeuOoX7125jw/a2mHCubt7Ohu1tzFvZzEXLX2Lhmk0F81QdDqGq3JVyTvNWNvPf5zUkVLwcUVdhaxVMX2sOlNJEPXilYCR7lenq7p0OicWZa7weRg4tj7Umvv6MY/jGuQ1UepwcUuHhK2cey6b32lPuAMYf3pj38wPSLs5yOiRW/+92OfA4hS5fkIeunojLITmfzZqvpnNKcaECrxSM5L412YQp+hKqscOrefjqiexu97Gn08/9a7fxtc+MLZinmik0Em3H8PaezpQwzj2XTaDGm1ubdWTewYeWSSq2k6msLxQyvL+vh3fbuukNhnCIcMPjrw54WEix1Hv3VZ7Y0u7jtZ17Lfvcax8ZJRty2Q9eUfpFXwLncAjGGGYufTG2CnXlVacQNIYKj4tyt4P39nanLHjKJOBRTzUQCLG7w8f21i7cTgfDq8Jln/mkrzsOfyBoewWNcvCiAq/YSqae71Hv1Otx8vi8yZS7HVz70PrYheDu2U389E//pKXdz/VnHMNRwyqp8Djp8AW47Bd/y+jlBwIh3tjVntJa+GMjqgsi8nWVnthFaU+nPybyHpfT9goa5eDFtv/pIjJGRJ4Vkc0i8rqIfMWuYynFS1/116GQYdc+Hx2+QEzco9tc+2Az3zyvga+f8zEW/Oo1Trv9OS64ey279vVQX1UW285qwU5Lhy+lE+b8B5tp6fDZfcopBAIhNr+/j+lLXuDURc8m9Mepq/RwRF0Fi2c2aoWLknPsdGUCwA3GmOOAScAXRaTBxuMpRYiIWK4sFQl721EPP12YImRMLCYffe6mVRuZf9rRCdslhzP86VoIB/PbQjgUMry7N1waabWK1OEQjqyrZOzIah6dO4nnbz6dJ66bUpAWArlujKYUHtsE3hjznjFmfeT3dmAzcJhdx1OKE6fAohmJ3umiGY04I9oV9fDTda4Mhkza9gTx24kIgUAoJlAuh/WFJd+iuafTz+52X8a7GIdDOKQy3Hrg8EMqGF5dXhBx37KrPeEuY/P7+wiUYE/9g4m8BCNF5EhgHPCyxWtzRWSdiKxraWnJhzlKHnE4HNy/dhsLpjbw6NxJLJjawP1rt+FwhP/rRcsIlz631fJC8P7eHkuh7vIHE7b7zq9f441d7bz9QSdfengD3/3N6yy5dHzC/sJj/PITf496w93+AMOHlHNWw/CUcyimGLtVrmTeymbe3dutnvwgxvYySRGpAv4C3GqMeSLTtlomWXr0VUUT/3p9VRmLLzyRti4/ezr9LH1uK0DKvNVlc5qoryrj3x92xbaLrmZdOO0EPC4Htz+9hfpqD9849ziCIYPTIZS7HNRXldueZLU65yWXjueuP7/JM5t2x87huJFDcuapD7Q0dGdrF6dadOFcNX8yR9RVarlmEVOwMkkRcQOrgYf6EnelNOmrTDD5da/HiS8QZOGjm2LiWFXm4sH/mMgHHeHFSz/+4z/5r3MbmLn0xYRj7WjtpsLj5IbHX2XB1AaWPreVLl8woXdNPlrkWnnD1z20nhVXnMzcTx7NyKHljOwjDNMfwe5vK2CrfadbkLWn08+ooaU11/ZgwjaBl3AW7V5gszHmh3YdRyl++lpBGV9G2O0PUut1s3hmIw4R2rp7UypsAC4+5Yi0nSWjMfrrzzgmJu6wP7lp9zDudJVDH3b6uWj5S7xwy+kZ7yKsBHvZ7CZG1ZRbti/IphQ1077vuWwCx9RXsWxOUywZXAxtHpSBY+e96qnAHODTIvJK5OdcG4+nDFKSE3yzlr8EwG1PvcG8lc04RFIE8yd/epNlkTF9sD8Wv/S5rbEY/ZHDKiyFtqfX3sRhurF6vcEQK644OTYrNl1s2zIe/mAzr27fazl+sD+tgNNdDFq7e2NtHlbNnxzLlRSyzcPBgN2VS7Z58MaYvwLayUjpEyvRuWnVRhZOO4Er7/u75UKglg4fo2rKeeK6KXT5gmz7oJPbn95CS4cv5u329IasFxA5hd3tPfQGQra0MkjusROdHVvudvDFhzckeM4jhpTR7U8Mw8QL9rgxNbFe8sOry/j/n9rMrdMbEzzz/rQCznQxcLnKGF1bgdfjYtTQcsYf3qgNyWwkH1O2tBeNkjfSxZXTJfiev+k0nJHOiskj7pITtW3d4fBOIGRwO4SKMicfdPhpafclJGjvu3ICbqeTlkgzstXN2/naZ8bmPC4ff64iglNgxtIXU0Q4ehGzGk1YX1XGjZ8dmzC85GeXjOewmnKGxYWY+iMU6cYIat+b/JOr70J70SgFJ5MIpfNAt7Z0MnJoOaOGeqnxejK2uk2+AESrVlra/bHhH1VlrhRxv3zKUdz5hy0pXvFAic87hEKG9/Z2c8eFJyZMqYomhSExbl7rdbN0dhMt7b6YuEe3+eLD63no6okcEjJpE9WZ7kqs7i501WxhyMeULRV4JS9kSgRaiU50UlNLhy/m0aQT4HRVKwumNjBvZTPzVjYzbkwN3zn/eC79+cspScQZTWNy9kdlNbT8zZaOtOcWHTcYtdsfCNLaDT/50z+55ezjLAWgpd1HhceV8Hlk2wpY+8IXD/mYsqUCr+SFTN5KVHQenTsptqr19qe3sGF7W+y9B7LveK90/mlH88WHE3vd3LJ6IwumNsTKBAeK1V3Kw1dPTLn43LI6nF8odzv4we+3xN4f/eP2B4I8s2k3M5rG2FK6qH3hi4N83E2pwCt5oS9vJdpZMb7vTPI2kLmGO3m49YghZYyu9XJR02iOHzXEMkRSV+lheHVZTv6orO4kWrt6LS8+Rw+vJBA0seZn8X/cezr9sdW9P7tkfOzCFH/XcfKRjbS0+9QLH8Tk425Kk6xKXsgmERjd5s4/bGFG05iY+B461IvL5UhZ9RptIVztdbK3M8Cu9h5WvLAt9t766jLKXEJLRy/XPtgce8/hdRW819bNAy++zTfPa+CwyP4HSnKyeNyYGu6YdWKstXGUaCItvoVwcr/76Hle1DSaz48fTW8wRDBkWLXu31x0yhH4AiFbqy+UwUOmJKsKvJI3slmdGQiE2LK7PWHBTXQhzu4OH++2ddPTG6Sq3MWX4koOH7lmEgvXvM7lU45KqDp56OqJXPrzly0rUu6e3cTHhlfhducm5plcFREdAp5s07I5TYwaWp5SHpn8WbV1+3mvrSdhJW60TcMFd6/VShgF0CoapUjIJvbb2t2b0lr3zj9s4StnHpsg+otnNlJfVcaO1m52tHaza18PM5rGpFSdtEQ6OS6Y2pDy2rUPNodFMUcCnxxTrav08Mym3QmVPG3dvdR43Zx/1wsZvW+HQwiGSFmJO29lMw9fM9H26gulNFCBV/JCtr1VrBKmM5rGpIj+Tas2snhmI8GQYeTQcspcToYPKUt5755OP2c1DOfY4VWWMfhcimJyTDXaC3/D9jbmrWwGwp72nbNOShD8dGWa6ZLHzsh+7ay+UEoDFXjFdvqzEMcqGVtX6UkRuvqqMqrK3VwbF75Ycul4zmoYzjObdse2W//2Hr58xrHMiRvxF1+mmGtRTK5/j/foz2oYzoKpx9MbDNHhC3DbU2/Q0uFj0YxGQqHU9gnpEtNej1Nr2fNMMQ1y7w8ag1dspz8r9tKVGl4SqV+PsuKKk3nkb+8wo2lMzBNe3bydr59zXMK81geuOsUyyblw2gmMHFpue2IyKgyhUIgPOv0pzbyiF5rH5k3m0BpvynvTXRiBrARnsApTMZGPlgIDQWPwSkHpz4o9q9KxWq87xWM9dkRVSvJy0YxG/IEgj1wziQ87/VSXu2jvCVge+4i6ClvPOZlAyKSEmaJ1+PNWNqHd5Z8AACAASURBVGPlaPVVRpcunxEv6sGQ4fu/3RTrQ19MwjRY6E+3zmJDBV6xnf6u2LNKxkaFLhQKETTQGwzhD4QSEq23rN7IfVeegi8Qoqc3SFdkmpLVsd1OBzvbuhla4aKu0p7WwfGe3x0Xnmh5oanxumOfRTpvuz8iYuVtLprRSEu7nw3b2waNMBUT+WgpYBf5mV+mHNREq0viW/v2N2Yc7Rn/YVcvs5a9yKcWP8eCX73GjZ8dy7gxNUD4j87lFH7w+81ctPwlVrywDY9TuNtidN/1j2zgxsdf5b02n21zR+M9v5AxrLjiZB6dO4llc5oYN6Ym1tb4nssmUOt1p8xEtWoN3J9jwv47heiQ8sEiTMVEuvbPgyGprR68Yju5WrGXTryiYY7RtV7e39vDM5t2M25MDZdPOYqL7wnXwC+cdgIfqa/krZbOWHuABVMb6OkNsqu9hxHVuR/lF/X8xo2pwSHCgl+9lljmWV3GIZUearyenIUB0nmb0SHlg0WYionB3KCtuJKsx1abdUuaCm2GUqT4AiE2/Ls15fmGUUPY2tLJsSOq2PZBJx2+AMeOqOadPV344rzV4w8dyuvv7qWqzMWYQyp4q6UTXyBImcvJsSOqqPC4kByGpv3BEK/v3McRdRUptpS5nBx/2BA8TkfGcxt3eC1l/bjwRI8ZPVZVmYvDaisodznoCYQodzsodzlzep4HA8ZAbyiEMSACboejaD5D+cxf0iZZNUSjDCrKkrzPMpcTt8vBkcMqcTsd9AbDDovLIQmCCuG4fZnLyaE13pi4A/gCQf65qwN/MLehGrfDwbEjq3E7HSm2+AJB4n0rEetzi4qIMWHx9gVC+IMh0vll0WOWuZyxC9nbH3Ty6o423v6gk2jEJ5t9KfsRAY/TQZnLgcdZPOLeF8UVohkyFs58rtBWKEVK295uulo6UypnvhIZLVdV5qQnsghqwakNLNy0KRauGDemhv9qPA6vx8lbvcGUgd0Az087jcPrKnNmrwDekOH9fT3cuCx12MejF06KhatcQPmudr6cFAZwjagmBFmX6UWPOTqSrD1v+UuJx21JLDu1q7JGyzPzSfrPtbhCNFoHr2Tgw04fW3d3MqzKQ4jwf+vu3iC1FR7K3Q7Ov+uFWEfJQ4eWY4DrHlpPfVUZN589lptWbaS+qozFF57IFStSa+OtatH7S7yweT1OAiGDCRk+6PAn9JRZPLORH/w+XAOfrra91uumtbsXfyDIRclCHbHXGJNWQNNNylo1f3LCBS7XfWyKvW681NA6eGXQEwoZdu3z8bXHXomJxp2zTmR0rZeQMXT7g9RXlbFhextLn9saE/lHrplEMGSYfe/LsXLKmx5/lcUzGxNG+S2d3cTwqoEJXHK3y+hFJbqK9YGrTqG9J0B1uSvWLgFISKbGr4Ltq8Ty3bZuZi59Ma2ApitP3dPpT9lXLitrBnPdeKmhMXhlUGAlGl977FXeeL+DUxc9y0XLX+Lms8cyq2k03z6/AY/TQZc/yL92d+B0kCByG7a38YPfb+GXcyfx/E2n8di8yYwdXjXgKpqojfVVZfwg7gIC8Mym3Vz2i7/R5Q/wzp4u5n7qI7FySSuBjT/ftu5eyzK9qFBHBTRZuK3KU5dcOp7VzdtT9pXLyprBXDdeaqgHrwwK0olG/EzTm1Zt5OFrJrKztTuhJNGqR01Lh4+W9h7augIcUVfB++0+Rg0ZWKmkPxCMtSXe25066KO+qoyqMhc3rUod4JEssPHnu/S5rSya0ZiSe7j96f3ToKwENLk8FeDHf3yTy6ccxab32ve3IJ7dlNOSv3yMolOywzaBF5FfAFOB3caYE+w6jnJwkE40ojNNx42pYf5pR+MQSfCcd7SG57M+dPXEBFH72SXjUmrTl81p4riRQ/oVJ46PuYsIt5zzMW58/FUWTG1Isff6M47h2odSxwY+fPXEFIGNP98N29u4/ektLJx2AkcPr8LlEL7z69diIZ7oZ2EloPErYVvafax9aw9v7u6IdbPs8gcZVVOe09j4YK4bLzXsDNHcB5xt4/6VgwircMPimY0sfW4r48bUcONnx7JwzSZ2RuLs8UQfL5x2Ao/OncTimY3hTpRJYjtvZXNKmCMT0Th5dPXprGUvMnJIOTtau2Ned7y9h9dVWLf/dUiKwCafb0uHj5FDy6kqc+IU+MqZx/Z7ZXB0ny0dPuatbOaGx19l5NByarzW7wuFDC3tPna2dvFuWze79ob768evro3fJvpa/J3DC7eczpPXnZoyuSv5PYo92ObBG2OeF5Ej7dq/cnDhcAjH1Ffx6NxJ9PSGaOnwUe520NLhSxjmEY1XJ3v6TofgD4Z4qvk9po07jNZO/4DjxFZ5gUDIxHrA3/70lthQ71GRqh4r20T2C198Fc0x9VWx8Irb5aCjJxAbFHJWw3Aevnoizsgs22zKEPuzojhdT5v7IyWp0aqfTNUyVglVrbDJL5pkVfJKJu+tr9febOngu795nXZfL8YYljz7LxZMbeCY4VUp8epkT/8nf3wTj9PBdZ/+KP5AiN5gyDJxKf1YwWKVF+j2B2LH37C9jYVrNuELhHAIOB2k2LZoRiNOSb0bmL7kBd5s6aCu0sNhtRUIktD2+JlNu7nk5y/jcTmpry7LWhyjwntYbUVC6Cb5M0/XFmJG05hYQjddtUxbtz/t95juPf25c1Kyp+BJVhGZC8wFOPzwwwtsjWInffU3z+TZxQtDS7ufm88eyzfOPQ6nCE6HWMaro71nnly/k2njDktIUt5x4Yncdcm4hLmuP7tkPMaYWJihL6Jx8mjtfY3XTWWZizWv7kyY2HT/2m18f/rHEYT7125Lee3W6Y19lhbaUZmS6fvI1NMm/rhWieTkObLx36NW2OSXgnvwxpjlxpgJxpgJ9fX1hTZHsZFM3ltfnl28MGzY3sbF97zM6bf/BadD8HqcLJ7ZmBCvrqvy4HYKV973d85oGJEyj/WGx1+lxuthxRUn89vrP8GKK06m3O3gn7s6svYm6yo9PHDVKXwnUpZZ4XESMnDxxCMZ6nVz21NvsP7tPXzrc8fT7Q/SGwyxYGoDC9ds4qLlL7FwzSa+9pmx1FV6CIXCr8V3m4wXPjs6Gmb6zNMdLxoC87iclttcf8YxKXNk47/HwdyZcTBScA9eGfwMZN5qJm/QSuCsSu9qvB5GDCln4bQTqKlwR8R1MzOaxjC61hvzOpP3vWtfDzc8/ip3X9qEuKDS4+K1HW0cO6Iqq/N2OIShFS527evhkb+9w+VTjoolbkfXernvypPp6Q3xhcgq1OiCqjVfPpVO3/7PKhRZ6bpwzaa05ZN2VKZk+j5GDfWmHC9qU/xxk7c5alhlxu+x1uvm4asnsrvdx55OP6ubt8cuckrusbNM8hHgNGCYiOwAvm2Mudeu4ymFYaDzVkfXenG7HAiZB0lnEjiHQziyrpLqcnfCsv6Wdj+LZjTS5Q+mLbHc0drNtQ81xzzrJZeOx+vJ3pvs8YfCvW/iEr0QFrXtH+6vx48+N//BZh6bN5nDaitin9+7e7tTvN7k8slctVyOJ9NF02qAuFPg1umNCcdNtslg0u4zmkeJ/w6XzWnimPoqTbDahG0hGmPMxcaYUcYYtzFmtIp7adKfpFm6UseOnkBsLF+60r++Su/iqzbiQzm3P72FQyrd3D27KSW5ufS5rbHtx46s5rYLPs5df36TQD/K9oLGJMSm46nwOC292d5gKCGZubvdl1X5ZHToicflxB8IRma9HniJYV+DWOITsofWeBkx1JuS0E1O2g6rLEu7T6v/K/NWNtMaWcug5B4N0SgDor/zVkcMKePOWScxrMpD0Bje39vDbU9t5tbpjZYeKoSrPKy81mjVTfxrIol3Ahu2t/FhZy+P/O0dFkxt4Oj6SrZ/2M2vNuyMJUa7/EHea+vm60/8gx9ddBKmH6JZ7nYmxKbjP4t0dw69wXBIK5o83dPpz2rlZ65LDNPNvx1IF8hMdxqaYM0/BU+yKoOb/ibNeiM9yOf84m+c+cPn+foT/+DyKUcRCoVi3uCooZFkaXsPm9/fZznGzqqscMuudjxO4Y4LT0zwII8cVsEzm3Yzb2UzNz2+kdpKN9PHHxZLdi741Ws4RKivKuOrj75CMLL/bBhWWcY9cyawunl7SgnkYbXlLE26c7j70vE8tfFd/IEgO1u7cIgwusabMlZw2ZzU9gF2lxgaDP9u7Rrw2MBkrz56gdAEa/7RdsHKgOivV/luWzezrHqjz50UE/bo/qJx8eRtn7zuVACmL3kh5bXH502mJxDA5XASNAanCG6nMHNp+JjjxtRw+6wTufwXqe2Co6P/nr/pNAAOHerNqjfNh50+Xt2+l2FVHqrKXThECIYMtz21mRqvhy9++qMAOB2C2yG0+3q5YsX+z2vxzEaeXL+Tcz4+iiOHVeB2OhhpMUIwXfvfF245PRbT7w9W3120jXG0DUJyK+GB9HnXRU72oO2CFdvob/LPRGLW8exo7ea9vT3s6wkwYkhZTADSVb9kqrpxOqDTF2L+g3+Picjdl45nxZUns/j3b3D5lKPSrmKt8bpjHuYlP38569403f4gP/nTm8w/7WiO9lRiBG57ajPPbNrNrKbR7O3u5bq46poll45nykfqeKx5Bztau2NJ2ivv+zuja8MDOayOmesmXsl3BPVVZfT0hlh8YSNbWzpjLY2jn/dABdqORLGSGQ3RKAMm3S25Felu06Ni0+3fH6dN1ybX63EiIqyaPzlWMw5wVsNwuntDzE+qSLn2ofVUuJ38d6TSJRrzTt5vlz/I4pmNeFwO6qvKsu5NU1nm5I5ZJ1JX6WFrSye3PbWZy6ccxaym0XzpjGNi4h6157qH1nPNJz8Se3/04hL9PVpCmExfSdH+Eh8Tj/bzWfCr1zjzh8+zcM0mbvzsWM5qGB67gOQiRNSf/yvKwFEPXskrVuWO0da3O1q7CZr9/Vqs2uQ+cNUp7Nrns6zP/vo5x9GSpiLlvb09lEWE22q/S2c30d7Tyw9+v4UffeEkvnN+A9N+trbPBGAoZNjZ1sO8lc2W9uze15O2QiZKfFfM6MVu1NDylGPl0gMOhQwOEf74n5/CKWAgoRWCValmtB1y/Ercpc9t1SRpEaMCr+SVqEg9Nm8y77Z1s6fTz+1Ph2O+o2u9lLsdsQvAhu1t3L92G4/Pm0TIhEsSXQ4Htz21OUWIVlxxMnu7e9NWpEQX1Sy+8ET2dPgIGcNtF3yccreTEUPKcQjsbOumvtpDb9BwSGVZVuGPDzp9MXGPt2fB1AY+jKzQtbInmreMj3vHXxzGH96Y9vMb6FQkq1DL/Ved0mepptfjTJhSFbW9P+sGlPyiAq/kHYdDGDmknL3dvXz10f0j+O65bALDKsO11FEv1etxWnrsLe3+WCJwR2s3H3b6aevujVWzJA/H+NWGnVw+5ajYLNaoOA3xuli45nWe2bQ7Fq83JkTIiGUlSzyhkKHLZ136F637trpbuHt2E1VlDv70n5/C7RKcInzzvOPY0+mPdWu0c2WnVajl33u6+ozvB0Impdf+Tas28sR1U2yzVRkYKvBKQegr3BDf6TAqRtGhHmUuBz+Y2cjNqzbGPP9ouODGz47l/rXbEjz625/ewvzTjk5ZaXrTqo3cfuGJsUlP0Xj9yv84BYcIZRkqaEIhw/v7etj2QaelMNZXl7H8+a2WbYMfWLuNZf/3duwiU1nmYnh1GaOGljP+8EbbE49W9eg/+dOb3H3p+IRWC8nx/d5AyHrhViBkm63KwFCBVwpGNuGGqBhFk4DxnvBdl4yjpzfEiCFl7NrnA4jFvsMXDUeszLKu0mMpTskyuqO1G6cIrV1+7vrzm9w6vTHFxmiIo9MX4Cd/ejPBQz+rYTjfPK8BEfjW547nmv/vaP7nd5tZuGYTd186HodDaO0KxI5106qN/HLuJA4d6k0R9YGUJGbCqhqnpcNHT2+Ix+ZNxhhjeTwdxTf4UIFXbCMXAhUVlWQPvL6qjG5/MCEevGx2E8OqPMx/cD0A3z6/gYXTTqDC46SmwpMgTuPG1HD9GcdQV1XGsjlNsZLAaHy8yx/kG+ceF2srEG93NMSxYGoDLR0+bn96Cz+66CRGDi3jw85eLv35yzGb7px1IosvbKS1s5dv/ep1Wjp8/OyS8by5u4MN29vY0dpNyGAp7nbVjFsluhfPbKSuysPIIenH9+kovsGHLnRSbCFXAhXvLc9c+mLs+WVzmiwXQT02bzKzlr2Yskhq3JiaWIKwvqosJVkYTW7+xyc+QrnbwRfj+sTfM2cCY0fut3vX3m7e3dtNfXU5vcEQb3/Qhdsp9AZNQnOxqE0Lp52APxhi3srm2HPRRVWja708ce0URCThQrin02+5kCt+0dFACIUMH3SGvXanhBOoNd6+L8B23VUoB44udFLyTl8DLLIlGqt/f19PggeebhFUbzDE4pmNuJ2OhLj98OoyDqn08OjcSRiItfCNvu+W1eFQiT8Q4sNOPwumNsS8+mtW7rc7FDK0+wKEDAltgH85dyK9AcMdF54YywdEPfSaCjf+uDh1/KKqe+ZMoMMXiJUoRi+EQ8pdtvZtcTiE4dWppZjZvC8XFxglP6jAK7aQy8ZS0aqb+PBAukZe7+3tIRgyjBpaxuPzJhMyhhsefzXBU6+tSLdCNpQgtNH6/PjVnPt6/HjdTq5Y8feEcNGHnYmrVaPvbenwMdTrpr2nl2VzmmLNzY6oq+CJ66bgckhszmrUjmseWMdj8yZbnp/X40zbfE1RklGBV2wh1wm55Kobr8eZEg9eNrsJl1P4j/vDz6244uSUfuy3rN7Iw9dMsrTtnT1dlvXsC9dswuNyEgiE2N7aw5ByV8KCn/rqMstFQgunnYDHFa7b/8a5xyWEfe6+dDy1lW4ChJucxYeS5p92NIFgiIevnsj3f7spVsJptcjrnssmMGJIOB+hgq8kozF4xRby0VgqGkfu8gXZ9kEn5W5HQp32o3MncdHyl1Le99vrP0FbV2/KStYF//tarLY+yp9v+BROh3DokHJaOv3MWvYiq+ZPZvuHXXztsVeZ8pE6vnzGMbzb1p0QmgH4439+kpseD5dyrpo/OSGHEI3DL1yzKbbQCUipFFo2p4lhlR4cDgcGwwVL1lrG+KN9bLR518GHxuCVvJOPxlIOhyAIs+8NV608OndSgvhZ9WiPb4OwYGoDI4eUM9TrpqXdR0uHL2H/o2u9vLm7g4VrNrF8ThM1FW7qq8oIGaitLOOX10ykrTvAxfe8ZBma2drSGavMSe7XEo3DR0slo4nY5Fr9eSubY/H/na1dCecS3aYispL0QPMcSumizcYU28hHY6n4WH9yc7Klz21NGMYdFeColz1vZTPv7+th9r0v8/i67Txw1SmxBmZnNQyPbbujtZu5K5sB4fufP54t77ezp8NHT8CkNDa7ZfVGrj/jGBbPDL832j1ydfP2BLvj+8/saO3m6OFVfGxkdca8RaZB2FbbK4p68MqgJj7Wn9wWoKXDh9fj5EcXnUR9dRlOh/C937yeEIapq/RQX1XGtHGHJSRYl81u4vWdbbGpT23dvbgcsK8nEIvrr5o/2VKQj6irwOkItx84pNJDe08vX/70MWx6rz3F04dI8tTtjP2eLm+Rrn49Gt5J3l5RVOCVQU286EWbkz189UScDkFE+M6vX+OZTbsZN6aGb5/fwJeShHZYVRnXn3FMamjkwWYevmYi77b1AOBxOtjXE2DFC9ti26VrJPbm7g7qKj2Uu53c8NirbNjexrgxNSyY2hAr17z1t5ti4Zv4xUKZFhIlh73cLgcdPYFYaEkXHinJaJJVGfSkW3yTnOg9q2E43zn/eP6xcx/Dq8uoKnMBBpfTwem3/yVlv3++4VMJXv3imY1UeJy8v89HjddNyBjcTkdCw7Q7LjyRe//6FjeffRyVZU4uXJo6veqJa6fQGwwRNFDudjCscn/4qr8LiXThkZIpyaoCr5Q08QLodjnw9YZYuOZ1Lp9yVMxr/+XcSdwYqZWH/W0MxhxSwdaWjoQ2Bg9dPTHWiuCshuHc+Nmx7GztocLjpMsf5JDK8OCOCo+TDl+ALn8ooSpGq1yUXKMCrygRAoEQO/d2x0QawoL+3WnHc91D6y3bGPzskvE4BN7d28PR9ZVsbemkxuvmkEoPV9739xQP/dG5kwgZw/WPvEJ9tYf/Pq8Bp0PUw1ZsIZPA21pFIyJni8gWEfmXiHzdzmMpSja4XA5cDkkQ5Q3b23jwxXd46OqJ/PCiE1N6nn/x4fW8u7eHhWs20e0Psbp5Oxctf4m93b2WSdZAyFDhcXLXJeO4dXojo2srdESdUhBsE3gRcQI/A84BGoCLRaTBruMpSrZYlRuufWsPHqfgELEU7WjN+rUPNTOjaQxA2tmuHpeDuqpyFXWl4NjpwZ8C/MsY85Yxxg/8Ephm4/EUJSushlffPbuJ+17YhtMhGWvNo2IPxMoy4/fz04vHZRwUoij5xM4yycOA+NUdO4CJyRuJyFxgLsDhhx9uozmKEsZqlW1NuYu6Uz9CmUu4e3YT1z7YnLZmPSr20bLMR66ZxAcdPtq6ehnqdVPj1TJFpTiwU+Ct7ktTMrrGmOXAcggnWW20R1FiWLW9PbQm7IkPLffw2LzJBEMhQoaEmvUll47nrj+/CYTF/mtnjqWyzIlDyhhdW6FJVKWosFPgdwBj4h6PBt618XiKkhNcLgeH1ngJBEL4An6+eV4Dcz95NHs6/fz21Z1887wGFkxtoMzt3F/DXlloqxUlFTsF/u/AMSJyFLAT+AJwiY3HU5Sc4nI5cLnKcRLAFfH4Tzh0CEO8Dio9mjxVih/bBN4YExCRLwFPA07gF8aY1+06nqLYRXm5i8PKtauHMviw9X+tMeZ3wO/sPIaiKIpijdZzKYqilCgq8IqiKCWKCryiKEqJogKvKIpSohRVN0kRaQHeKbQdaRgGfFBoI7JEbbUHtdUe1NaBcYQxpt7qhaIS+GJGRNala8lZbKit9qC22oPaah8aolEURSlRVOAVRVFKFBX47FleaAP6gdpqD2qrPaitNqExeEVRlBJFPXhFUZQSRQVeURSlRFGB7wcislBENorIKyLyjIgcWmib0iEii0XkjYi9T4pITaFtSoeIXCgir4tISESKrgRtMA2PF5FfiMhuEXmt0LZkQkTGiMizIrI58t1/pdA2pUNEykXkbyLyasTW7xbapmzRGHw/EJEhxph9kd+vBxqMMfMLbJYlInIW8OdI2+ZFAMaYWwpsliUichwQApYBNxpj1hXYpBiR4fH/BD5DeIjN34GLjTGbCmpYGkTkk0AH8IAx5oRC25MOERkFjDLGrBeRaqAZ+Hwxfq4iIkClMaZDRNzAX4GvGGNeKrBpfaIefD+IinuESixGEBYLxphnjDGByMOXCE/UKkqMMZuNMVsKbUcaBtXweGPM88CHhbajL4wx7xlj1kd+bwc2E57jXHSYMB2Rh+7IT9H+7cejAt9PRORWEdkOXAp8q9D2ZMlVwFOFNmKQYjU8viiFaLAiIkcC44CXC2tJekTEKSKvALuBPxhjitbWeFTgkxCRP4rIaxY/0wCMMd80xowBHgK+VMy2Rrb5JhAgbG/ByMbWIiWr4fHKgSEiVcBq4KtJd8hFhTEmaIw5ifCd8CkiUrThr3h0DlkSxpgzs9z0YeC3wLdtNCcjfdkqIpcDU4EzTIGTLf34XIsNHR5vE5F49mrgIWPME4W2JxuMMW0i8hxwNlDUiWxQD75fiMgxcQ/PB94olC19ISJnA7cA5xtjugptzyAmNjxeRDyEh8f/usA2DXoiict7gc3GmB8W2p5MiEh9tApNRLzAmRTx3348WkXTD0RkNTCWcMXHO8B8Y8zOwlpljYj8CygD9kSeeqmIK36mAz8F6oE24BVjzGcLa9V+RORc4EfsHx5/a4FNSouIPAKcRrit7S7g28aYewtqlAUi8gng/4B/EP57AvivyBznokJEGoH7CX//DuAxY8z3CmtVdqjAK4qilCgaolEURSlRVOAVRVFKFBV4RVGUEkUFXlEUpURRgVcURSlRVOCVokNEOvreKu17H4p0fnwt0lnRHffa50XkW5HfvyMiOyOdQd8QkbtFxBF57XYR+XSa/R95oJ0aRaRCRH4bOd7rInJb0utfFZHLRORnEbs2iUh35PdXRGRmJtsUJRkVeKXUeAj4GPBxwAtcHffazcCSuMd3RpafN0S2/1Tk+Z8CdrUFvt0Y8zHCvVdOFZFzAETERbhn0MPGmC9G7DoX2GqMOSnys8pm25QSQwVeKVokzOKIN/4PEbko8rxDRJZEvOA1IvI7EZkJYIz5XaT7nwH+RqSLpogcC/iMMR9YHMoDlAOtkX28A9SJyMg+7CsXkRUR2zaIyOmR5ytE5DEJ9+J/VEReFpEJxpguY8yzkWP4gfXs7/L5aWB9XAdQS7K1TVFABV4pbi4ATgJOJLw8fHGkj/gFwJGEve6rgcnJb4yEZuYAv488dSphQY3na5EOge8B/zTGvBL32vrIezLxRQBjzMeBi4H7RaQcuA5oNcY0AguBJgv7aoDPAX+Ks6+5j+P1xzZFUYFXippPAI9EOvntAv4CnBx5/nFjTMgY8z7wrMV7lwDPG2P+L/J4FNCStE00RDMcqBSRL8S9thvoa2LXJ4CVAMaYNwi3rzg28vwvI8+/BmyMf1MkHPMI8BNjzFsZ7EtHNrYpigq8UtRYterN9Hz4RZFvE+5r859xT3cTDsOkYIzpJezpfzLu6XKgW0QmxiU5z8+FfcBy4E1jzI+ysc+C8sj2ipIRFXilmHkeuCgybKGesAD/jfDItBmRWPwIws21ABCRq4HPEh6rF4rb12bgo1YHiXQ2nAJsjXv6WOA1Y8zLcUnO5C6SzxMe/BKN8R8ObInYNyvyfDSBGz3W94GhwFeT9pXWPguOZRC0qlUKjwq8Usw8t7VRvAAAARNJREFUSTi88SrwZ+DmSEhmNeE+7a8RnuP6MrA38p6lwAjgxYjXHZ269TwwLiLmUaIx+NcIz0ZYArH4/UeBvmbDLgGcIvIP4FHgCmOML/J8vYhsJNyyeSOwV0RGA98kXLWzPmJftMrnKRLvICzph22Kot0klcGJiFRFhiDXEfbqT42If6b3/Bj4jTHmj31sNx0Yb4xZcIC2OQG3MaZHRI4mnEg9NlI5k+l9TxK+iL1pl23KwYVOdFIGK2silSgeYGFf4h7hf4CJWWznAu4YgG0VwLMRb1uAa/sS9whfJ5xsTSvwObBNOYhQD15RFKVE0Ri8oihKiaICryiKUqKowCuKopQoKvCKoigligq8oihKifL/AP6E06Zsjoy+AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#graph volcano\n",
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns\n",
    "from numpy import log10\n",
    "tstats['-log p'] = -log10(tstats['pvalue'])\n",
    "\n",
    "figure = sns.scatterplot(x='log2(B)-log2(T)',y='-log p',data=tstats)\n",
    "figure.set_title(\"B cells vs T cells\")\n",
    "#plt.axvline(1, color ='orange')\n",
    "#plt.axvline(-1, color ='orange')\n",
    "plt.axhline(-log10(.05), color ='orange')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next we will identify the altered proteins."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "altered = tstats.apply(data_utils.is_altered, axis=1, change_factor=False).dropna()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "These are significantly altered (p<0.01) without restricting by the magnitude of change:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "T cells    56\n",
       "B cells    23\n",
       "dtype: int64"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "altered.value_counts()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Here are counts for other p and change cutoffs:\n",
      "P<0.05, >150% of the other type\n",
      "T cells    63\n",
      "B cells    30\n",
      "dtype: int64\n",
      "\n",
      "P<0.05, >200% of the other type\n",
      "T cells    23\n",
      "B cells    16\n",
      "dtype: int64\n",
      "\n",
      "P<0.05, >500% of the other type\n",
      "B cells    2\n",
      "T cells    1\n",
      "dtype: int64\n",
      "\n",
      "P<0.01, >150% of the other type\n",
      "T cells    47\n",
      "B cells    14\n",
      "dtype: int64\n",
      "\n",
      "P<0.01, >200% of the other type\n",
      "T cells    21\n",
      "B cells     9\n",
      "dtype: int64\n",
      "\n",
      "P<0.01, >500% of the other type\n",
      "T cells    1\n",
      "dtype: int64\n"
     ]
    }
   ],
   "source": [
    "print(\"Here are counts for other p and change cutoffs:\")\n",
    "print(\"P<0.05, >150% of the other type\")\n",
    "print(tstats.apply(data_utils.is_altered, axis=1, pvalue=.05,change_factor=1.5).value_counts()) #50% increase from either\n",
    "print(\"\\nP<0.05, >200% of the other type\")\n",
    "print(tstats.apply(data_utils.is_altered, axis=1, pvalue=.05,change_factor=2).value_counts()) #Double\n",
    "print(\"\\nP<0.05, >500% of the other type\")\n",
    "print(tstats.apply(data_utils.is_altered, axis=1, pvalue=.05,change_factor=5).value_counts()) #5 times larger\n",
    "\n",
    "print(\"\\nP<0.01, >150% of the other type\")\n",
    "print(tstats.apply(data_utils.is_altered, axis=1, pvalue=.01,change_factor=1.5).value_counts()) #50% increase from either\n",
    "print(\"\\nP<0.01, >200% of the other type\")\n",
    "print(tstats.apply(data_utils.is_altered, axis=1, pvalue=.01,change_factor=2).value_counts()) #Double\n",
    "print(\"\\nP<0.01, >500% of the other type\")\n",
    "print(tstats.apply(data_utils.is_altered, axis=1, pvalue=.01,change_factor=5).value_counts()) #5 times larger"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Differentially expressed proteins"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEKCAYAAADzQPVvAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO2deZhU5ZX/P6e27uou6IamG1QgLkGUGFRAWXTiFjVGRmNcSBQMuACiIXHc4kxInDCTcY1jYhAwCSjBFU3cjTEumbhgQBQVg/xcEAxLA93QTXfX+v7+qIVabvVat5bu83meerRuVd176jb1veeec95zxBiDoiiK0rtxFNoARVEUxX5U7BVFUfoAKvaKoih9ABV7RVGUPoCKvaIoSh/AVWgDrBg0aJA58MADC22GoihKSbF69eodxphaq9eKUuwPPPBAVq1aVWgzFEVRSgoR2ZjtNQ3jKIqi9AFU7BVFUfoAKvaKoih9ABV7RVGUPoCKvaIoSh+gKKtxFEUpbSLhMOHmeiTsxzjLcPpqcTidhTarT6NiryhKTomEw0S2rcP9yIXQ+DlUDyd0wQMweJQKfgHRMI6iKDkl3FyPKy70AI2f43rkQsLN9YU1rI+jYq8oSk6RsH+f0Mdp/BwJBwpjkAKo2CuKkmOMswyqh6durB6OcXoKY5ACqNgripJjnL7aaIw+LvixmL3Tl9myJRIOE9y9ldCujQR3byUSDufZ2r6DJmgVRckpDqcTBo8iOOPPSDiAcXosq3E0kZtf1LNXFCXnOJxO3FVDcA0cjrtqiKV4ayI3v6jYK4pSEDSRm19sFXsR2U9E7hORehFpE5F1InKCncdUFKU00ERufrFN7EWkGngNEOBM4HDg+8B2u46pKErp0JVErtJz7EzQXg9sMcZcnLTtUxuPpyhKCdHZRK6SG+wM43wLWCkiD4vIdhF5R0SuEhGx8ZiKopQQnUnkKrnBTrE/GJgDfAKcDtwF3AxcafVmEZkpIqtEZFV9vWbjFUVRcokYY+zZsUgAWGWMmZS07efAOcaYw9v77Lhx44zOoFUURekaIrLaGDPO6jU7PfstwLq0bR8Cwy3eqyiKotiInWL/GjAybduhQNbp54qiKIo92Cn2dwITROQ/ROTLInI+MBf4tY3HVBRFUSywTeyNMX8nWpFzAfA+8N/APGCBXcdUFEVRrLG1EZox5hngGTuPoSiKonSM9sZRFEXpA6jYK4qi9AFU7BVFUfoAKvaKoih9AJ1UpSiKJZFwmHBzPRL2Y5xl2qSsxFGxVxQlAx0Z2PvQMI6iKBnoyMDeh3r2SrfR2/zei44M7H2o2CvdQm/zezeJkYHJgq8jA0saDeMo3UJv83s3OjKw96GevdIt9Da/d6MjA3sfKvZKt9Db/N6Pw+nEUTWk0GYoOULDOEq30Nt8RSkt1LNXuoXe5itKaaFir3Qbvc1XlNJBwziKoih9ABV7RVGUPoCGcZScoKtpFaW4UbFXeoyuplWU4kfDOEqP0dW0ilL82Cb2InKTiJi0x1a7jqf0jEg4THD3VkK7NhLcvZVIONzpz+pqWkUpfuwO46wHTkx63nkFUfJGT8MwuppWUYofu8M4IWPM1qSH3tcXIT0Nw+hqWkUpfuz27A8WkS+AALAS+HdjzCc2H1PpIj0Nw+hqWkUpfuwU+5XAdOAfQB3wY+B1EfmKMWZn+ptFZCYwE2D48OFdPpiW/nWfXIRhdDWtohQ3toVxjDHPGWMeMcasNca8CEyOHe97Wd6/2Bgzzhgzrra2a7f/iZjzklNx/XI07iWnEtm2rktJxr6MhmEUpfeTtzp7Y0yziHwAjMj1vsPN9fuSi5CIOQdn/Fm9zU6gYRhF6f3kTexFpBw4DHg55/vW0r8eY1cYRsNrilIc2Flnf7uInCAiB4nIeGAFUAncl+tjJWLOyWjpX8HR8JqiFA92ll4OBR4kWmv/OOAHJhhjNub6QBpzLk50Za2iFA+2hXGMMd+xa9/paMy5ONHwmqIUD72mEZqW/hUfurJWUYoHbYSm2IaG10qPnvRIUoqbXuPZK8WHhtdKC21V3btRz16xFYfTibtqCK6Bw3FXDVHRKGI0od67UbFXFAXQhHpvR8M4SlGji7LyhybUezfq2StFiy7Kyi+aUO/dqGevFC3a8yi/aEK9d6NirxQtGkPOP7pepfeiYRylaNGeR4qSO1TslaJFY8iKkjs0jKPYRk8raTSGrCi5Q8VesYVcrcbUGLKi5AYN4yi2oKsxFaW4ULFXbEEraRSluFCxV2xBK2kUpbhQsVdsQStpFKW40AStYgtaSaMoxYWKvWIbWkmjKMWDir2i9CK0S6iSDRV7paCoOOUOnTSltEfeErQi8u8iYkTk7nwdUylutIVxbtG1DUp75EXsRWQCcDmwNh/HUwpDV4dVqzjlFl3boLSH7WIvIlXAcuBSoMHu4ymFoTteuopTbtG1DUp75MOzXwysMMa81N6bRGSmiKwSkVX19erZlRrd8dJVnHKLrm1Q2sPWBK2IXA58GZjW0XuNMYuJXhgYN26csdMuJTckJ1ddmC576XFxcqUlFFWcuoeubVDawzaxF5GRwM+BfzHG6H15LyOj8uOiR7s8rFrFKffo2gYlG3aGcSYCg4D3RSQkIiHgBGBO7HmZjcdWbCYjbPPqLXD2gi6HEBxOJ+6qIbgGDsddNUSFXlFsws4wzh+BVWnblgAbiHr86u2XMBnJ1c2r4C83YaY/QziCeumKUmTYJvbGmEagMXmbiOwFdhlj3rfruEp+SCRXkwW/eTsh8eAeqGEERSk2tOul0i208qPv0tX1FEpxkNd2CcaYE/N5PMU+NLnaO+mofYW2ZChd1LNXuo0mV3sXnVkYp6ueSxcVe6Xo0bBBfuiMkOuq59JFu14qBaMzHS81bJA/OiPklol5XfVcEqhnrxSEzvbS0bBB/uhM+wpNzJcu6tkrBSHcXB/11n11cPrPwTsA194thFqG4Oi3Tzg0bJA/OtO+QhPzpYuKvdJlcjFwRML+qNCf/BN48qqEuDin/J5IxcDE/jRskD86K+TakqE00TCO0iVyNXDEOMvghBv2CT1EPfaHp6aEaDRskF+0wqr30qFnLyLlwBzgeMAAfwPuMca02WybUoQkwi9pMfTgjD93ydtz+moxoUOQDkI0GjZQlNzQmTDO/UAT8KvY8+8Cy4Dz7TJKKV5yEUOPh4FcTk+nQjQaNugaOtdXsaIzYj/SGHNk0vOXReRduwxSipuextBTSil9ddFOmU/M0X72OUJLVZVsdCZmvyY2QxYAERkPvGafSUox09MYekopZaxTJmfegZn7bjQUpKLUI7RUVclGZzz78cDFIhJ35YYDH4rIe4Axxoy2zTql6OhpDN2yNfLy8wnPfQ+3hmp6jJaqKtnojNh/w3YrlJKiJzF0LaW0Fz2/SjY6DOMYYza298iHkUrvQUsp7aXQ51f7GBUvuqhKyStaSmkvhTy/mhwublTslQT5KtnrbhhISwo7R6FKVXO1BkOxBxV7BSh+r6zY7VM0OVzsdBizF5EmEdmT9tgkIn8QkYPzYaRiP8Veslfs9imd65qpFI7O1Nn/ArgOOAAYClwL3As8BPzOPtOUfFLMXlkkHMYV2lu09ilRCp0cVtqnU6WXxpjxSc8Xi8ibxpifici/22WYkl+KuWQv3FyPe9fHRWufEkWT78VNZzz7iIhcICKO2OOCpNdMtg+JyJUisjYp9POGiJzZc5MVOyhmr0zCfnj1Fjjr7hT7zJTfF4V9yj60a2bx0hnP/iLgLmBB7PkbwFQR8QJXtfO5zcANwAaiF5XvAX8UkbHGmLXdN1mxg2L2yoyzDJq3w0s/Sww6IdhC2Lc/riKwL1ckGsSF28DhIOz04kjq7a8oPUGMyeqc5/5gIruAG40xi9p737hx48yqVavyZJVS7MQrcdInKLXXR6fUyjQj4TBm2zqcSd+RsxcQrhyC1Bxc1LZ3h1L7+5QKIrLaGDPO6rXO9LMfSrS98XHs62f/A2PM5i4Y4CTaEtkHvN7Zzyl9mxRBqKwldOlLEGzt8K6jFMs0rWrUeWIOzjPvIFjWr1fVqZfi36c30JmY/RLgSWB/ohU5T8W2dYiIfFVEmgE/sBA4xxjzXpb3zhSRVSKyqr5ey+n6OlYTsWjaiqPqgA5jwaVYppmtGgp3Ra+rOCrFv09voDNiX2uMWWKMCcUeS4HOZsXWA0cBE4B7gPtE5AirNxpjFhtjxhljxtXWatKtr9MTQSjmMtJsZKtRJ9jS6yqOSvHv0xvojNjvEJGpIuKMPaYCOzuzc2NMwBjz/4wxq4wxNwLvAFf3xGClb9CRILTXcKsUF/c4fbWE06qhOHsB4aqDel3FUSn+fXoDnanGuQS4G7iTaMz+9di27uAAyrr5WaVEyEXyrb26/45ivvEy0vSEbjGLpsPpJBKrhurt1Til+PfpDdhWjSMiNwPPAJuAfsCFREsxzzTGPNfeZ7Uaxx7yUQHRncqZru4n3FwfjeGnXQgi058jEgljnGVIZQ1m786iKyNVouz7t6h/n1zSrWocEfkV7SyaMsbM7eC4Q4Dfx/67G1gLnGGM+VOHFis5J18VELnqfNhe3X8kS4jHsWczjt+dvs9T1OqOokWHyOef9sI4PXKtjTHTe/J5Jbfkq/1sLpNv2QQhW4iHvfWJ42lrXUVJJavYG2PuS98mIkOMMVvtNUmxg3xVQKQI8dBxcNwPobIWh8NBJBzutqedHILCWUb44qdx7lwP7goItoC3Bp6/ft8HtLpDUVLoaj/7Z4Exdhii2Eu+Gp0lkm+v/g+MnwVPXhUNsfQgbJQRghp5JuaE6+GZa/atNj1ncfTNObzAKEpvoksJWhFZY4w52kZ7AE3Q2kGuEqedPRZNW3AsPSPj4hKc8WfcFqGV9pLHwd1bUxOyU34Pf/r3zDDOt+6BSDhxgbHzOypKMdKjdglp3JsDe5QCkM9GZw6nk1AkjKOTYaOOkscZISjvAMuQlOl/AHL/WToWT1Es6JLYG2MWdPwupVjJZwVEV8JGHSWPM/bV2gAjz4SjvhsV/tYGeOdBcDjBV7evM2ZrA7z2v30ydq+NxpR0dAatkpWUpKi7AhMJd1o8urJwpqPkcca+Pn8Lc8L1yCPTEvs2Fywj7OmP65Sb4Ik5KZ0jcXtzeVqKhmyCro3GFCvy2uK4s2jMvvCkxPh9dZAmop2JhXd24UxGTB4y4vvJ+3I4HJb5gMj057qUJyhlurPorDeeByWV9mL2nemNo/RBUhqRHffDfUIPnW5K1tmpRZ2ZkpW8r0gkbH0nEAn2mQZb7TWK00ZjihUaxlEsSRGMLAnRXIlHV5PHWRdVOZyZ20eeicPhILRrY6+KXbcn6MU8T1gpHOrZK5akdCZsbbC9S2FXZpdadog8ZzHy6i3RGH18e6we37H0jERP/Mi2dSkdMvNJe506u0p7nSOLeZ6wUjg0Zq9YkouYfS5syJaANDs/wbn7030raCvrwL8nmowtq8J4KjHGFE0Mv70YO5DyPfc1ccueDO9o3URH+RKt1umdtBezV7FXspIsGLi9sWqc/HQp7E4CkjPvgOXnRz3cKcuhrD/yyyMz9h2a+x6ugcMztttJtiR06NKXoGnrvu8ZuxtJrjTKdmHtbufIfC6wU/KLJmiVbpEcWnH1q+10mKU7ZIQ4WnZ1OQGJuyLx//LwRdGLVJEMychmszPcmvo9j/ruPqGPvSdbMrwroa9kdCxg30QTtErBsaoLNxetgH/9JTjdicVRbF7VbgKS1oZ9zxs/h+BeOOvulPYJ4QLFrhM2++qi1U3eAdHwkzg7tTo4l5U0Wq3TN1GxVwpOxgpaXx3SvB2emrtvcdRZd8PKRSkJyOQwBGcvgL/ctG+n1cNh9+boRSK+orbffhhPf5wFCFU4fbWEL3wc596tqYu+piyPrgZe/0z0jfFkuI2VNFqt0zfRMI5ScDI8TYu6fp68CnP6fyfi0o7Bo4hMfw4u+VP0DsBdAc3bo++Pi3/sboCHp8Ifr4Dt6yDYmv8vSGzsYFm/jO8lD1+EOf2/94Wb3nkQc8EyWytptFqnb6KevVJwMjzNLKGMCM6EV+5wOon0249QSwOuxy+PhkfOvAMz8BAirkocgT3RuwNIvTP45jF5/GapZAufRHASia0xwO3FGQlgLn4SHE5b5tDmsymeUjyo2CsFJyMsE2yxDDNEnB6S5SibaDmdTkKhGpwXP4k0b4tOsFq5iNAJNxbUe80WPok4PbirhhAKBXFu/yClEsd5wTLC3moc5FaIdSxg30NLL5W80V5td/w1V7gNad4CkQj8cXZC9CJTlkPdV7o8tLyYhlp3WBvfuDl7z5/qoQWzWykdctnPXlG6RUedGOOeZmjXRlwv/BhO+SlM+0N0GElrA1JZS3j3ZsJpi6vaWxiU7L2GQkFo2oKJBDEONxHfYFwud17PQUfhk6y9fSLBvNqp9E40Qavkhc7UdkfCYRyuMvjm7fDElfCrsbD8PHCVIc9cE2158Oy/Ibs/J9ywEdn9Oe5n/63DVgihYABn/TocS89AfnkUjqVn4Nz+QfQCkGccTidOXy3G6UHCfsLN9QmbjcNtvS7Akd+LktI7sU3sReRGEfm7iOwRkXoReUpEjrDreEpx48iSnHTEarvjnr9j21p45OLUSpxHLobTfgbfewomXoncfxbOu0ZHp1KNnxWdO5tlYVAkHMbZ9AXy8NTUKphHpuFo3mb3184gFArC9g9wLzk14yIV8Q3OqMQxFywj4hucdzuV3oednv2JwAJgEnAyEAJeFJGBNh5TKVIk3pEymerhiCP6TzDh+bsrrFfHNm2FUNu+OH58+5NXRUs1Y8/TFwaFm+ujSdoiCI9EwmGcezbjePgiyzscl8tNuO4rRKY/h5n7DpHpzxGu+0rew00AkYihvsnPFw0t1Df5iUSKL7endA3bYvbGmNOTn4vINGA3cBzwlF3HVYoTI86M1aycdXd0O0lliVkWFdHakLUkE++AxPsiDjehUISG1iCBUJjBxh+txrFaRORwIzZ/72TCzfW4s114Yhcpl8sNsWSsUJg4ayRiWL+ticvvX8XmhlaGDvCyaNpYRtb1w+XSyG+pks+/XL/Y8RqsXhSRmSKySkRW1ddrj47eRgSBlYuiq1mnPxP978pF0e0klSW+9r/Ri0Jy++Kz7o5uz9JqOb49eP5yfvHaLv6xrYnPduzlqgfW8GF9IDqfNm2fZsrv8xIeSfaQJZx04Un7DsW0enXn3kBC6AE2N7Qya9lq/rm7VT38EiZvpZci8ggwAhhnjGm3kbeWXvY+OtOSN6Wl8gk3YGpGIDs3wKu3RFfCDh2X0Wo5MmU5Ye8gmgMRbnl1Bw+t/oKhA7zMP/sIPC4HT63ZzI/GRah+8/bogPLKWoxvMGHf/rg8ZfZ+5zQP+akZh/LVt38azTMk3eF0p6y0o+Pu3BsgEArjcTmpqfTgcHT+HuaLhhaOu+XljO0rZk/kSzWV1Paz97wp3afgpZci8gvgeOD4joRe6Z10VHZo9brfU01ZRRPu+ErY5u2EKofgv/hPOCNB/MaJ3zOQY/8nVZg2N7RS4XFyzaPvMm/yKG5etZWfnH47jU3NbGmKsPiVen54ag0jB3dNBLtKuof8kxe3suSb10YvPKf/PHHhifQfiiuL0HdVuK1CMPdePI6Rg/tZfs5q/x6Xk6EDvAm7AYYO8LJzb4D9qsp7eFaUQmG72IvIncB3gJOMMZ/YfTyleOlo1Wb89bgAtbSEaAwP5eCpz+MmhHG4ueuNBhb+30eJzyyZfoylMDW2Btnc0Eq1183pX92f0+59P+U9H2xp5vE5k6jrZ594BULhlGOu2bSHGc/C/VNuo7Gpme1NhmFVg6nLkoDNFjsfVOnB4XBYCr9VCOby+1fxhznHZXjk2S4MI2p9LJo2llnLVie233LuaO57/VPGDB+d47Ok5AtbY/YichdwIXCyMeYfdh5L6R3EBeicBa9xwm2vcOWD77C2sYzvPryZtY1lLPy/z1Le/8u/bGDR1LEMHeAFSAjTwlc+ZugALy2BMAcOqkgRXYiKYFswYut3iXvIydQ3B3l9m5O5z+5gj3MgbWGTtdolW+z8nc27OWfBa6zf1pTxufQLTPxzgVDmDXW2C0NDa5CRdf144LLxrJg9kXmTR3Hf659y9akjqaksntxCb8PuCijbPHsR+TUwDfgW0CAicZeu2RjTbNdxldLGSoBueGwt8yaPoiUQzvDi65v97FddzuNzJtHiD/Ppjr3c/qf11Df7WTR1LPtVl9MWjFh6/x6nsL2pjWAo0q3YdkfUVHq49+JxKZ7zbeeN5g9vf8H13xjJdSvWpnjUg/uX0RrYF05JFu6jh1Uz+8RDqPa6qetXRq2vzNJjzxaC8bgyw0TtXRhcrjKGDqjA63GxX1U5Y4aPzvn5UfbR1fBbd7AtQSsi2Xb8n8aYm9r7rCZoez/ZYtHZkoN/ve5EfOUutu3xZ/1BRCKGxtYArYEwYWNwOxy4nUJja5D6Jn+KuC6dMQ6300l9k5+dewM8tnoTV586Mqc/rvTv6XY5cDmE1kCYKYvfzBDk+WcfwYylf098rxqfh28veJ1aXxnXnj6SGx7bZ/+vLxzD8jc38oOvj+CAARUpx+usaNQ3+TlnwWsZdliFfBR7ydXfoiAJWmOMugCKJe0JUjbPFKDa66Ha6+EPc47LmrBMvxj8+sIx/PrlDdQ3BZg3eRTVXjcOEUD4pH4vFR4nHqeDOSd9mTv/vJ7/Pmd0ToXO4ZCU/UUihtZAK3ecfySNrUEWvvIxazY1JpLKsC+c8uRVx3HvxePYurstIfTx16984G3uv+RYvB5nxvFGDu7X7jmKY3Xnce/F4zRUUwC6En7rLtoITck77SURayo9lsnB/3pmXUKIs4mx1X6vfOBt5k0exaxlq5m1bDUAD14+nvomP/OeeD8lvDLjuINy9uOyunMBMi5yt5w7OhF2amzdt6J3c0MrrYEwg/uX4XY6LIVg194AvvLMn3D6BSYbXbkwKPbSlfBbd1GxV/JOe16MwyEMqvQkvPDG1iC3/2k9azY18tN/bV+Is+033VMdUlXOtN++lXJRuG7FWpZdemxOflzZ7lxqfB7LfMT8s4+g3O3g1ufXJ/YR/6G3BsJ8XN9sWylkZy8Mir3k4y5LxV7JOx15MQ6Hg/lPr+vQy0n3nt0uB0MHeKn1lSWSmS0x73joAC+TDq7hipMOwSHCvMmjEiEUiAqvQyQnPy6rO4w7/7ye/zhzlOXF6JC6SkJhQ32zP/Fd4z/0eD7h1xeO4coH3rYshezpIiql8OTjLkuHlyh5p6MkYvLrtb4y5p4ygoMGVVJR5mRQZVnW9xw62EdrMEx9k58lr33KuWOHUVPpocbnob/XxZZGP7N/vzrx/uE1FWxpbOWOFz6ivtnP41dMoq5/z+vu05PMRw+r5trTRxIIRRKhozjxJFxc2NN/6PHvufLjek4ZtR/BcIRwxLBi1ed8a8wwRtT62FDfbGsVh1I6tJegVbFXCkJH3mi8smZLYxuzfr86RchG1PrY3uznn42ttAXDVHndXLE86vU+NHMCv/vbJ3xv0kEp1SvLLxvPRb9ZaVnZctt5oxncv5wDaypzIpDplRWLpo1l/tPrLI/d0SIpgFAowvrtTSl5jHhjsobWoFbUKAnaE3ttYacUhHis+IABFdT2K8sQOYdDCEdICD3sC4es397EBYve4LyFb/Cjx9+j2R+i1hcVNgHOHTsso3qlvsnP5oZWZp94SMZr161Yi6/clTNPOB5/jVcR1VR62NzQyppNjdz+p/XMmzyKh2dO4KHLJ3DXix8x/n9eyrpICqChNZgQ+rjNs5atTnT2tLuKQ+kdaMxeyStdiS9bCdm5Y4dlCN91K9Zy23mjCUcMg/uXs391NG6f/NmdewOcNqqOQ+t8lmWPwVDuVtOmx19FJJGjWLOpkVnLVjN0gJc7LziKc8cO49LjD6axNZi19LM9Qc9HFYfSO1CxV/JGV1cJWglZ3EtOptZXhq/czRVJ4Z7bzhvNrc+vTyRg3/5sJ98/5VCm/e4ty7LHXItjcpVLJGIyKi0euHw84bBBhEQS9nuTDiISybzotCfoWiufX0o5Ga4xeyVvdHWVoNXF4YHLxnPhb1am7GPJ9GN48K2NnDt2WKJc87HVm/jusV9KrEi9/5Jjufh3b2Uce/7ZR1DXv8z2wRzJIuH1ODNyEfHqmpvOOoL9q70Zn+0ood2RAJWySBUL+Whp0FM0QasUBdlaIbx2w0kpS/6TSRepAV53RvXJQzPHs3Fna0ri85ZzRzNsoBdjYMvuNgb3L+Ok21/N2P8r157Igys/41tjhtn+o03+LlbtEuZNHsUR+/e3PBfdEevkz4Qjhv96Zh0vrNtelCJVCpRCe4mC97NXFOjeKkGrRT8jB/fjyauOS/TAEYT7Xv80Y7HS0hnHAtGkrTFYHjsUMVxw7Jeob/LT2BpgYKU9P9pkr/CO84/Muvgrfi6sxL0rgmLlhd5y7mjqmwKs2dSYte2xkp1ST4ZrNY6SN9KrVHoSX962x8+UxW/ytVtfYcriN/nepIM4elh14vXNDa24nMKtz3/Izc/9g53NARamtUJecNEYbn3+Q06541WuffRdtjS22TZ2L77QqtZXxsBKDytmT2TRtLEJm4cO8FLXr4yaSk9Km+fjbnm53Uqdjo6XfgGcfeIhieelIlLFglXL6lJKhqtnr+SNXK0SbK8Ncrz/zdABXrbubqO+KcC1p4/k6kfeodZXxvyzj+DAQRVEDNzy3Ie8sG57on1wazDM1j1tDOlfnvPwRiAUTtTZz1j694xY/Q++fij7V3lxOIT6Jn/W3kGd9cSzeaHVXnfi/JSKSBULpZ4MV7FX8kouerF01ANn6AAvC6eOZd4f30+pq9/c0JpI2C6dcUxC6NMXOtkRz/a4nMw9ZURGjf8Nj63lkVkTUy4wuQgXpIfMjh5WzdxTRlDj87Bk+jF8qaaiZESqWCj1xnEq9krJkS32X+V18/DMCbQEwtT6PNQ3+6n2ui2F0+WI9tGxWmRlRzy7ptLDQYMqLW0xxqQIRke5jc4ka5O90FpfmeWwFK3Q6Tql3KP7mMMAAB/VSURBVDhOY/ZKyVFT6ckYRXjbeaO5fsVarnn0XWr7lfG7v33CHecfmZhulcxpo+oodzu456IxlnX7dsSzHQ6hosw65gukjKFrL7fR2Xh+shd694VHJ4Q+/v0uv38V/9zd2qO8QGewe9Se0nm09FIpSRpb/OxoDuIQMEBzW5Bg2LBfVTn/+dQHiRDN9d8Yia/Mleidc9qoOr5/8qFcsTzaEO22849k+pLM+vtcePZxzzkSiRA24BTYsTeQ0uPmzguO5OfP/oP6Zn/W2vnkCVciwgWL3siw9+GZE7J659lKXlfMnsh5C9/I+fdO/v7FXpfe29DSS6VXEYkYtuxOnUi14KIx7FcVFan6pkDivXvaQvjKXDw0cwI7mwNUeJyJBOnmhlaue/RdbjtvdEaIo6fx7LjQ3fnn9SlN2U4bVcf9lxzL7tYgjS1B3LGFXOnho3i4IF0wV8yeaHknsrmhlWsefddSTLOFhXbuDWTsJ5d3NO0NqSnVUEgpo2EcpeSwEpE5y9/m/X/uYcriN7n+GyO5YOxQfnrWKDxOBy2BMBu2NWMwNPtDKaK3ZlMjtz6/ngcvn8BfrzuRR2ZNZEStr8eeZ9zGiycemJITeGHddi7+3Vu0BcMEwhFCYcOt543m6GHVlmKb/l137g1YhoIaW4MJMU0Xcauw0KKpY3ls9aaM/YjkzuMu9br03oZ69krJ0V5ZYbwx2gOXj+eLhtaM0YP7V3szvNzafh6C4WhPmn/ubqEtGO5xu+N4qeV+1d4MW2t9ZfjKXCl3E/ESzPRyyPTvuvCVj7nl3NEZq4Vv/9P6xHlIF1OrKhLB8L1JB7FuS1PKfpw5jK5ok7biwlbPXkS+JiJPisgXImJEZLqdx1P6BtkWt8RnuNb6ynCKZCQlr1uxFqdDUpK7p42q46qTR3Dx797i5Njiqm172mhsTfWOO0NyMlJE+Om/juLznS0Zts49ZUQihxC37YbH1vLjM0dlhI/Sv+uaTY3c9/qnPDJrIn+9/iTmn31EYmxj/DxYiWl6S2kTW3Ucb7c8b/Io7nv9UxyO3ElCLhfRKT3Hbs/eB7wP3B97KEqPiVfjpDcSu/1P6xNJ2S272yy9/7ZgmIE+D/PPPoJBPg9VFR4uvPfNjIvCwzMnQGXnbUqPrZ82qo4bv3k4//nUugxP/MBBFdYlmFm+a/pCnqtPHUm520EwFKGuf5nlOMPOnMOrTx3Z6QVCkYhhx14/bcEwThG8HifV3syBM+mlnO3VpWvpZ36xVeyNMc8CzwKIyFI7j6X0LWr7lbF0xrHs9YfoV+7i5uc+ZM2mRpZMP4brVkRX01qFEMpdTgKhCOXuaCzfH7K+KIQihkjEdFp80mPr544dxmc7Wqhv9icGlsRn4nqcDkvbPqnfS4XHlUjMxoWwxudJ9AJyuxw0t4U46+7XEheVBy4bj9MhXRLMriwQsqqqSZ/u1V7lTWc7mmqljr1oglYpCB3VX2d7PS4SP/7jezT7g7QEQtz83IecO3YYK2ZP5ODa6MKleGw7vRfOz57+gB8+9A61/cpxOoQqr5vTRtWlHHvoAC9bdrdlJDrbIz22Xu1188u/bIg2H2v2M2vZaq559F0G+TyUuRwZfXpuOXc0v/zLBgKhcEYt/bcXvM62PX72q/IiSEqr5hfWbefC36zE43JaTvxqj/TQDmB5zq0S4tetWMvGnS2Jc5RtyPrWPW2Wf+NslTpdOedK1yiaBK2IzARmAgwfPrzA1ih20pWB4+mvJ4tEfVOA678xkhu/eThOEVwO4ZMdexk6wJsyArCm0sPg/uXMf/qDRK+ceG19/CIAJNr/LrhoDD994gPuvvDoTn+neGy91lfG7BMPocbnYe4pI3hizRcpXv2gfmUMjA1Nn3/2EVR4nDS2BlOGqLRXsmhXhUt75zzbMSs8zsRx099z9LBqvjfpoMSagPS/sVbq5J+i8eyNMYuNMeOMMeNqa2sLbY5iIx15de29niwSazY18t17V3LS7a/idAhul5P73/gs4dGv2dTI/KfXUe52YjC8sG67ZXuEOcvf5seTv8Izc49nyfRj8DiF+mZ/l8oQayo9PHD5eO644EiGDoj20R9Y6WH2iYfwl3XbeOnDrYwY7CMQirB1Txv9PC6GVJVzzaPvMmvZ6sSiqgFeN4FQmDvOPzKlK2ZcCO3qvNjeOc92zJZAOHHc9Pdka0MR/xuXegfJUqRoPHul9Olswq0jr667M1fjScc7/7yem7/9Vfav9rJxZwvz/vg+c08ZwdAB3qy9crbvaaPM5URcUFnm5jffG9vlMsSmthB3vfhRyiKqeNO1tmCE78QGlsQbtY2s86XEzK0Gs6SPTrSr82J753y/Km/GMeMx+/hx0+3qqA1FTaWH+y85lo07W6jwOGkJhLU5m82o2Cs5oSsJt2yCLRIN4XR35mo86fjf54zOmAb1y79s4LbzRid65VitJp3/9DrmTR7F/KfXseCiMYnVrZ1hZ6wNwrzJozI82k279tX7x7fN/v1qHp45gf1ibY0Btje1WbZunn/2EQypKk/5jrnuvNjeOY8f8/E5k2gLRnAKGdU47Q1ZT99fHH8okrIO4t6LLVf5KznC7jp7n4gcJSJHxY41PPZcg/K9jK4k3Kzqr285dzQ3Pfk+67c1McDrzlqfnSwqr91wEn+Yc1zKBSW5K6HVStkjDujPPRbJ0YWvfMzmhlYOG9KPpTOOoSUQxlgWQ1oT94yt7hwqPE5LL9cfiqQ0H2sLWnvXh9RWWn7H/aqi32HL7tYeNxnrqCbe4RDq+pUzfGAFBwyoSOQdkklO+A7pX97u/jRBm3/s9uzHAckdmP4z9rgPmG7zsZU80pWEW8JTvGISewMhwhGTGDQST0Sme68DvO52Q0TpISSHhWdZ3+zHH4rwq798xJLpx7C7NUgwHCEcMfzojMNoCYTZvsfP+YveiLYUmDaWAd7OVbjEPePG1mDGcbPdTYQjJqVXjDOLN+x0iOUA8VyWLlrdMXR0zru6v+TPa4I2/9jq2RtjXjHGiMVjup3HVfJPdxJu9U1+pv32Lb7+i7/yo8ff49rTR1LrKyMQCqd4r04H7bb1tWr7uzcQ4o7zj0zxLH994RgiJpqovX7FWlyxoPyPHn+PKYvfZN4T7xMxJtGnZtay1eyILVjqiLhn/NjqTRklnwcMKM8otbznorGs2biLeZNH0RIIUd/kx1fuZMFFYzJaN3s9mefQbs/YYPi8oaVHLZDTSzutmrMlowlae9EWx0pO6KqnWd/k55wFr2V4sUtnHMvASjcDK/d1fNy6uy0l5h1/b9wjzrav300fR4XHRShscDqEPa1B+nvdfDe2YnbF7In88OF3Mj6XPN7w1etO5IAqL65OxO937fXz7qbdDPJ5qChz4RQhYgw3P/ch1V4PV578ZYCopy4QCEeY+pt9JaALp47lqXc2M+bAGmoqPdT4PLgcwgHVFZ1uW/zaDSdxwICKDm1NJ9vCqVufT23FED/nPV39qouq7EFbHCu209XEYbbb+MaWAMFwhGqvJ+G93nH+kV2u3qn1lREMGb6zdF8FzJ0XHInLCYumjeWuFz9iQJaKkeQ5rSKwfnsThw/p36EItQbCzFj698RM25FD+rFxRwv1TQGqvR52twaZE+uJE6/nn3RwDY+s3pxI2qbP0X38ikmWx811k7HkO4W4/W6ng1tjQ2HWbGpMnPNcCHWpj/grRYqmzl4pfdq7bU8n2218XHSSa+rjcfD093o9Tuqb/ISNYcn0YxI16QA3nHFYoncOREX86kfexe10MrhfGfMmf8WySVly3P2280bjcTq468WPOhUe8XqcPDRzArdfcCQep4OrH3qHeU+8z4/OOIy5Xx+REPq4PXOWv83lXzs48fnkC038eVsoYnmsXDcZi5/r+Eze+U+v47yFbzBj6d+59vSRHD2sOnExyVUIqSv/XpSeo569UhCsSijjNeXpNfVWbX3vv+RYtu3xW4Yd6pv9DKkqt/Tad+0N0BIIU+Z2JNoZJO934dSxVHldzJs8ilufX8+dU45ixnEHdZg4jEQM2/b4ufbRdzO+zzWPvsvSGcda2uNMErjkzp3x59lq/XPpGYdiF5SXrjkBl9OR0Rguvfxzy+5Wan1liZXBja1BFr7ysSZXixwVe6UgxMXqkVkT+WdjKzv3BhKteq1q6m//03rmn30Ehw72ISKEIhE27myj1lfG5obY1KkVa1l26bF8tK0ZDO3W0z86eyJzTxlB/3IXS6YfQ7M/RGNLkOoKF62BCAtf+Zj6Zj+hiGFIVXmHNfc79vota+TjYRmnWNsTz3dGk7Zj+NVLGxLPsyVnk89hTyc+hUIR/rGtidmxu6Bsk7AOqfMxtDq6JsDrcWYMMO/IVqXwqNgrBcPhEIb0L2d3azCRKM1WUx8IhfF6nBnefNx7jseUt++JNh178PLxGeMG4++t9ZWxszmQMdhkYKWb7z/wDvXNfm47bzS1/cpoDYTo73XjasdjjkQMLf7sA1WGDvCydU9bxl3EPVPHUlnm4KVrTiBiwOsRLjn+YC49/mBaAmEG9y+n2mvvitLtzf6E0MO+SVjpFyWv25m4awhFjOWsgMfnTLLVVqVnqNgrBaWjcESy91rfFBX65BBCSyDM9d8YyXfvXZno1wJw6/Pr+elZo1g641gaWwIpdw5Lph+TInBxsbr5219NVJ5Ep11NoGFvkKoKCGaJnUcihq172vg01oAtXSRbAuFEeAlINGbbr6ocg+H8hW+mXHAOGlRJOLaKOB8Jy2A4kmKzVcgsPRcQDEUsL2zZzpFSHGiCVik4nU3UxUf9xROI8dp4X5mLWf9yIA9cNp6DaitZMv0YABa8/P8ocwllLgfzn16XCBFlGx7idjpSnocjhv2ry3lp3VbCFhXK8aqUfza2JuL/yQnTpTOO4dDBPg6o9vKLKUdy/TdG8tjqTZS5HPyzsQ1jolVD8eNdt2It4YixPA8dtYTuLu5Yb/048UlYD8+cYLlCGbRGvlRRz17JC7mYSuRxOZl7yoiM3jO/emkDPzjlUC78zcqURGtTW5C5D77DiDof919yLLv2Rj38bXv8CS88XmZYU+mhyuvm6GHViYuCAJsaWjnhsMEI0fh2cr19vCpl3uRRiSEl/zvlKGp8HnY0BdjZHOCapITtbeeN5trTR/Kbv37K65/s5J6LxiTuSuLfpb2Lih016XW+MhZOHZu40xk6wMvcUw5lcL/yrGsL7GrGptiLLqpSbCdXYhWJGD7f1cKJt7+Ssn3RtLHMf3pdRghl/tlHMGPp3zNej48uXPLapxkdKuODvy89/mDK3Q6ufGBN4rXF08ZyWFK9/RcNLdz14gbmnHQIBvhsRwtup/Cjx99LNFSzsikQjjBr2WqGDvCy7JJjOemOVxOvP37FJEQk5aK4c2/ActFYfIFTTwmFImxv9hMKR3A5HdT5yjpcRKYjBYsTXVSlFJT2hnF0RawcDqGiLHMxUbZ2ul+qqUhpa3z0sGquOe1QhlSVU+5y8NN//UpKZ8x4Bc1DMycQCEXYtTfAvMmjWPjKx6zZ1MjMZat5fM4k6vqVA1BZ5mTapAOZ+tt9q2CXXXostb4yDq3zccf5RybKEuMJ5OoKN4FYbHtzQ2tKNc6908bR7A8lJlHFL4r9y1229pFxuRzsX+3t+I1J5KISSMkvKvaK7eSy6dWgyrKMEEJdvzLL5OiuvYHoMBKXg0dnTUSElKqfZZda174HQpEUwU2u+GkLRoU6EjG0BCJckZbo3RGbnjXN4vP1zX6qvG6a2oIsmjaWmkoP5W4HK288GYfDgdNBYrZsfH+X37+KR2ZNzPh+p42qQ0T4oqFFPWulU6jYK7aTy6X92bozpl8A7pk6Fo9TmLH072xuaGXJ9GMyesp/tqPF0q6NO1ss6+XnP70Op+wLS7mdjpTKoIgxDOrnYdpv38r4/Pyzj8DjcnDzcx9y4xmHp4SHFk4dyyCfh0DIJNYNAIl8Qigc4YHLxvNfz6zjhXXbOW1UHXNPOTRj5N+IWh8NrUENrSiWaMxesZ18NL2KRAw79vpp8Yf5dMdeyt2OlFrwh2dOYMriN1M+c/Swan529le4IqlfzcKpY5n3x/cTJZhxXrrmBAAGVLoIhYVzFrzGQ5ePZ1OsiqbWV8ZNZ42iv9fN9j3+lPANwIv/9jWuezTaY2bF7Imct/CNxL7jzdfmP70upUzz2tNHpuQTFk0by6BKD2FDQuiT9/HAZeNTktTaWKzvoTF7paDko+mVwyEIwtTfRsXu4ZkTUsTQqs98fbOfZn8oUfs+sNKTmD+bzNABXjZsb2b+0+tYHAu/bG5oxekUavuVc/8lx+L1ONm8qyXh1aeHbz6u35uo8knvIRNffBUvv4wncdOrjmYtW80f5hxHKEtYbHuTv8d5EaX3omKv5IV8JPSScwPp4r7wlY8tV9Qmt/B9+ZoTWPbGZyllmo+t3sSVJ43gpic/YHNDKzOXRccJ3jT5MAIhk1iwVeV1J/riQGr4ptzt4Nbn1yd66v/65Q0pdif3xIm3JjDGZM1zZAuLWV1EtF+NEkfFXuk1JItg+krQ+mY/Xk90OEiV180n9XsTSVeINR1zCl8bOTglObtw6lgGVrqYfeIhiaZfbpdwzEGDuCgpZHLfJdbJ3i/VVOB0CP9x5uHU9S/HHwxx5UkjWLelKeMOIG6H1+1M/L9VnsOqzj3etjkZXeikJKNir/QakkUwvhL0gcvG43QIIsJNT77PC+u2c8HYoUyd+KVEuGboAC93nH8kgmSETmb/fjXLLxvP25/tTAwVaQ1EePKdzSnvi7dLThfnDdubY1U3Tn7w4BrWbGrk6GHVzJs8isOG9MMY+Pmz+1b3Ji9O6miwenqS+upTR6ZcRHShk5KMJmiVXkW2xT7pSeJZ/3Ig0yYdxBcNrQTDEZwOYUhVOSfd/mrGPp+48jgArnwgdfDIXn8IhwiNrUE+2rKHkw8fnOihn7xA6/pvHE65Wzj+llcS+xw6wMsjsybiFAgbMCazH05XFy7pQidFE7RKnyFbbiDdG3a7HPiDEX77t0/43qSDuG5FtLwy2Ts/elg1c08ZQZXXzac79qa0U56z/O3ECt3TRtVx1ckjuOsvH6Ukex9+ayMzjjuIvf4gZS5PYt9xr3tI//J2xbireQ5d6KS0h+2evYjMAa4D9gM+AH5ojPm/9j6jnr2SD0KhCF/sbk3E3o8eVs1NZ32FKx94m1pfWUbP9jvOP5KIMQlv/uBBFXyyo4VD63yJRVRxhg7w8tDMCbgcwtY9bTy39p9cNOFAnA5Rr1uxjfY8e1u7XorIFOAu4OfA0cDrwHMiMtzO4ypKZ3C5HLgckhDpNZsaWf7mRu6/5Fh+MeXIjJ7t1zz6Lm3BCFMWv8n8p9fRFjQ8tnpTSsljnM0NrexsDhCJGAb3L+PSrx3C0AEVOoJPKRh2tzj+N2CpMeZeY8yHxpjvA1uAK2w+rqJ0ivR2vY+s3szNz32IQ8RSwCti05g2N7RyxfLVnDt2WNYZuVVeN06HsF9VBXX92g/ZKIrd2Cb2IuIBxgIvpL30AqAjbZSiwGpw9+X/cjAmNtYwmfQZsfHFUPEyz+R93HPRGB5c+RkOh46MUIoDOxO0gwAnsC1t+zbg6+lvFpGZwEyA4cM1yqPkh4zEbWKAiWHRtLHMWravuia5lQHsE/81mxoTM3K/VFPBhu3N/OqlDVx96kgtfVSKhnxU46RngMViG8aYxcBiiCZo82CXogDZq1gG+cp5/IpJtIXCOEVoDYZTavMXXDSGu2MDwuub/QypKqe6ws0R+/dnzPDRmoRVigo7xX4HEAaGpG2vI9PbV5Siw+EQ6vqXE4kY9rT58ZU7WH7ZeOqb/OzcG+CZd7/gP84cxU8mf4Uyd1KFTWWhLVeUTGwTe2NMQERWA6cCjya9dCrwmF3HVZRc43AI1RXltLWFGFABLkc5tf3KOGL//tR4PZSX63IVpfix+1/pL4BlIvIW8BowG9gfWGjzcRUl55SXuyjHRf+uDXVSlKLAVrE3xjwsIjXAj4kuqnof+KYxZqOdx1UURVFSsf3+0xizAFhg93EURVGU7GgRsKIoSh9AxV5RFKUPoGKvKIrSByjKfvYiUg8UcxJ3ENF1BKWA2moPpWJrqdgJamsu+JIxptbqhaIU+2JHRFZlayNabKit9lAqtpaKnaC22o2GcRRFUfoAKvaKoih9ABX77rG40AZ0AbXVHkrF1lKxE9RWW9GYvaIoSh9APXtFUZQ+gIq9oihKH0DFXlEUpQ+gYt8DROReEflYRFpFpF5EnhCRwwttVzoiMlBEfiUi/4jZuklE7ol1JC06RGSmiLwsIo0iYkTkwELbFEdE5ojIpyLSJiKrReRfCm2TFSLyNRF5UkS+iJ3D6YW2yQoRuVFE/i4ie2K/oadE5IhC22WFiFwpImtjtu4RkTdE5MxC29VZVOx7xipgOnA4cDrRkYsvioi7kEZZsD9wAHA98FVgKvA14MFCGtUOFUQH099UYDtSEJEpwF3Az4GjgdeB50SkGIcm+4i2FP8B0FpgW9rjRKJdcScBJwMhor+hgYU0KgubgRuAMcA44CXgjyIyuqBWdRKtxskhsT/6u8Bhxpj1Hb2/kIjIN4GngWpjzJ5C22OFiIwD/g4cZIz5rMDmICIrgbXGmMuTtm0AVhhjbiycZe0jIs3AVcaYpYW2pSNExAfsBr5ljHmq0PZ0hIjsAm40xiwqtC0doZ59jhCRSmAG8DnwWWGt6RT9AT/QUmhDSgER8QBjid5xJPMCUa9UyQ39iOpSQ6ENaQ8RcYrId4jeQb1eaHs6g4p9D4nFcJuBZuAM4BRjjL/AZrWLiFQD84F7jTGhQttTIgwCnMC2tO3bgCH5N6fXchfwDvBGoQ2xQkS+Gvu9+4mOVz3HGPNegc3qFCr2aYjIf8USWu09Tkz6yHKi8dsTgI+AR0Wkokhtjd+BPAV8QTSGnxe6Y2uRkh73FIttSjcQkV8AxwPnGmPChbYnC+uBo4AJwD3AfcWaUE7H9rGEJcj/Ar/v4D2fx//HGLObaIxxg4i8SfT281xgmW0W7qNLtsbioc/Gnk42xrTZZZgFXbK1CNkBhMn04uvI9PaVLiIidwLfAU4yxnxSaHuyYYwJAP8v9nSViBwDXA1cWjirOoeKfRrGmB10v0+1xB5lubMoO12xVUT6Ac8Rte8bxphmO21Lp4fnteAYYwIisho4FXg06aVTgccKY1XvQETuIir0Jxpj/lFoe7qIgzz93nuKin03EZEvE/XgXwTqgaHAj4jG8p4uoGkZxIT+BaJJ2W8BlbFwDsCumLdSNIjIEKIe9KGxTaNieYbPjTG7CmcZvwCWichbwGvAbKJlrQsLaJMlsbu4L8eeOoDhInIU0b930dxBicivgWlE/102xP72AM35dkg6QkRuBp4BNhFNJF9ItHS0NGrtjTH66MYDGEbUU94OBIj+A1hOtOyy4Pal2Xoi0biy1ePEQttnYe9NWWydXgS2zSFabeUHVgNfK7RNXfybLy20bWl2Zvt3eVOhbbOwdSnRCXr+2O/+ReD0QtvV2YfW2SuKovQBtBpHURSlD6BiryiK0gdQsVcURekDqNgriqL0AVTsFUVR+gAq9oqiKH0AFXulqIk1neruZ5eLyHoReV9Efpc8Z0BEviUiP4n9/02xIR/vxAa83CMijthrt4vIyVn2f6CIvN9N2ypE5JnY8T6ILdhJfv2HInKxiPw6Ztc6iQ6eeSf2OK892xQlHRV7pTezHDiM6MAWL3BZ0mvXEx2aEedOY8xRwKjY+0+Ibf8V0ZXRdnC7MeYwoo30jhORMwBExAVcAjxgjLkyZtc3gY+NMUfFHitstk3pZajYKyWBRLkt5qW/F5sahYg4RGRBzDt+WkSeFZHzAIwxz5oYwFtEW1ogIocCfhPt15OOBygn1k/dGLMRqElaxp/NvnIRWRKzbY2InBTbXiEij0h0nN3DIrJSRMYZY1qMMS/HjhEA3o7bR3Ri09umg/bTnbVNUUDFXikdvk20teyRwNeB20Rkv9j2A4l645cBE9M/GAvfTAOej206jqi4JnO1iLwDbAE+Msa8k/Ta27HPtMeVAMaYrwLfJdr6tpxoe4UGY8xoojMExlrYVw38K/CXJPtWd3C8rtimKCr2SslwPPCgMSZsjNkGvAocE9v+qDEmYozZCrxs8dkFwF+NMf8Xe74f0eZ1ycTDOHVEG8V9J+m17UQbnnVk3zIAE+3cuJFoI7fjgYdi298H1iZ/KBayeRD4pdnX2tfKvmx0xjZFUbFXSgbp4vboiyI/BWqBf0va3Eo0VJOBMSZI9A7ga0mby4FWERmflCA9Kxf2AYuBDcaY/+2MfRaUU9wDxZUiQcVeKRX+CkyR6OzPWqJi/BbwN+DcWOx+MNFujwCIyGXA6cB3jTGRpH19yL72vymIiBCdKftx0uZDgfeNMSuTEqRPWth3UWwfhwLDiU41+htwQWx7PPkbP9Z/AVXAD9P2ldU+Cw4FulURpPQtVOyVUuEPREMg7wIvAdfHwjaPAZuJCt4iYCXRyWEQ7TM/GHgj5o3/JLb9r8DRMWGPE4/Zv090zsMCSMT7vwys6sC+BYBTRN4DHibajtkf214rImuBG2LfYbeIDAX+g2j1z9sx++LVQs+RemdhSRdsUxRtcayUPiLiM8Y0i0gNUW//uNiFoL3P3AU8ZYx5sYP3nQOMMcbM66ZtTsBtjGkTkUOIJmEPNR0MjBGRPxC9oG2wyzalb6GTqpTewNOxihYPML8joY/xc2B8J97nAu7ogW0VwMsxL1yAKzoS+hg/IpqozSr2ObBN6UOoZ68oitIH0Ji9oihKH0DFXlEUpQ+gYq8oitIHULFXFEXpA6jYK4qi9AH+PwWugvEgVj4XAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "figure = sns.scatterplot(x='log2(B)-log2(T)',y='-log p',data=tstats)\n",
    "figure2 = sns.scatterplot(x='log2(B)-log2(T)',y='-log p',data=tstats.loc[altered.index])\n",
    "plt.rc('axes',labelsize=16)\n",
    "plt.xticks(fontsize=14)\n",
    "plt.yticks(fontsize=14)\n",
    "#plt.title(\"B cells vs T cells\", fontsize=20)\n",
    "\n",
    "plt.savefig(\"data/figures_selected_FP/BvsT.png\", dpi=300)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This shows the differentially expressed proteins. Recall that this dataset is restricted to those expressed in at least three of each cell type.\n",
    "\n",
    "The next step is to identify those proteins by name. Here we get the protein names from the fasta file, which remain the same regardless of which software the data was quantified using."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "protein_names = load_data.load_fasta()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "protein_names_by_group={}\n",
    "for ID in data_raw.index:\n",
    "    if ';' in ID:\n",
    "        ids = ID.split(';')\n",
    "        names=''\n",
    "        for i in ids:\n",
    "            if i in protein_names.index:\n",
    "                names=names + protein_names[i] + ';'\n",
    "        names = names[:-1]\n",
    "        protein_names_by_group[ID] = names\n",
    "    else:\n",
    "        if ID in protein_names.index:\n",
    "            protein_names_by_group[ID]=protein_names[ID]\n",
    "protein_names_by_group = pd.Series(list(protein_names_by_group.values()), index=protein_names_by_group.keys())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here we add the names to the <code>tstats</code> dataframe, which contains the fold changes for each protein. Then we can list the proteins that are altered. Optionally, use <code> head()</code> to display the first few or <code>to_csv</code> to write out to a file. We now split up the proteins heightened in B cells from those heightened in T cells."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "tstats['Name']=protein_names_by_group\n",
    "tstats['GeneID'] = tstats.Name.apply(lambda n: str(n).split('_HUMAN')[0])\n",
    "\n",
    "altered_proteins = tstats.loc[altered.index].sort_values('log2(B)-log2(T)',\n",
    "                        ascending=False)[['pvalue','log2(B)-log2(T)','Name','GeneID']]\n",
    "higher_in_B_cells = altered_proteins.loc[altered_proteins['log2(B)-log2(T)'] >0]\n",
    "higher_in_T_cells = altered_proteins.loc[altered_proteins['log2(B)-log2(T)'] <0]\n",
    "\n",
    "#To write out to a file, use:\n",
    "#altered_proteins.to_csv('data/filename.tsv', sep='\\t')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next we do the same for those only identified in one cell type."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "data_in_one_type['Name'] = protein_names_by_group\n",
    "data_in_one_type['GeneID']=data_in_one_type.Name.apply(lambda n: str(n).split('_HUMAN')[0])\n",
    "\n",
    "B_cell_only_geneset = data_in_one_type[(data_in_one_type.Type=='B_')]\n",
    "T_cell_only_geneset= data_in_one_type[(data_in_one_type.Type=='T_')]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Proteins Higher in B Cells:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>pvalue</th>\n",
       "      <th>log2(B)-log2(T)</th>\n",
       "      <th>Name</th>\n",
       "      <th>GeneID</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Protein ID</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>P09960</th>\n",
       "      <td>7.198480e-07</td>\n",
       "      <td>1.972804</td>\n",
       "      <td>LKHA4_HUMAN Leukotriene A-4 hydrolase</td>\n",
       "      <td>LKHA4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P01834</th>\n",
       "      <td>1.092879e-03</td>\n",
       "      <td>1.913174</td>\n",
       "      <td>IGKC_HUMAN Immunoglobulin kappa constant</td>\n",
       "      <td>IGKC</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P29350</th>\n",
       "      <td>5.002525e-07</td>\n",
       "      <td>1.744085</td>\n",
       "      <td>PTN6_HUMAN Tyrosine-protein phosphatase non-re...</td>\n",
       "      <td>PTN6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P0DOY2</th>\n",
       "      <td>5.558317e-04</td>\n",
       "      <td>1.684692</td>\n",
       "      <td>IGLC2_HUMAN Immunoglobulin lambda constant 2</td>\n",
       "      <td>IGLC2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P09525</th>\n",
       "      <td>1.857150e-04</td>\n",
       "      <td>1.525379</td>\n",
       "      <td>ANXA4_HUMAN Annexin A4</td>\n",
       "      <td>ANXA4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P09874</th>\n",
       "      <td>2.778706e-06</td>\n",
       "      <td>1.377227</td>\n",
       "      <td>PARP1_HUMAN Poly [ADP-ribose] polymerase 1</td>\n",
       "      <td>PARP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P27707</th>\n",
       "      <td>2.038183e-04</td>\n",
       "      <td>1.325106</td>\n",
       "      <td>DCK_HUMAN Deoxycytidine kinase</td>\n",
       "      <td>DCK</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q96KP4</th>\n",
       "      <td>2.237074e-05</td>\n",
       "      <td>1.321275</td>\n",
       "      <td>CNDP2_HUMAN Cytosolic non-specific dipeptidase</td>\n",
       "      <td>CNDP2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P58546</th>\n",
       "      <td>4.286378e-03</td>\n",
       "      <td>1.164319</td>\n",
       "      <td>MTPN_HUMAN Myotrophin</td>\n",
       "      <td>MTPN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P22392</th>\n",
       "      <td>2.731861e-03</td>\n",
       "      <td>0.981652</td>\n",
       "      <td>NDKB_HUMAN Nucleoside diphosphate kinase B</td>\n",
       "      <td>NDKB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P59665</th>\n",
       "      <td>3.874038e-04</td>\n",
       "      <td>0.957458</td>\n",
       "      <td>DEF1_HUMAN Neutrophil defensin 1</td>\n",
       "      <td>DEF1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P50395</th>\n",
       "      <td>8.590047e-05</td>\n",
       "      <td>0.941751</td>\n",
       "      <td>GDIB_HUMAN Rab GDP dissociation inhibitor beta</td>\n",
       "      <td>GDIB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q16666</th>\n",
       "      <td>3.008010e-03</td>\n",
       "      <td>0.835062</td>\n",
       "      <td>IF16_HUMAN Gamma-interferon-inducible protein 16</td>\n",
       "      <td>IF16</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UJU6</th>\n",
       "      <td>4.328426e-03</td>\n",
       "      <td>0.596015</td>\n",
       "      <td>DBNL_HUMAN Drebrin-like protein</td>\n",
       "      <td>DBNL</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P41240</th>\n",
       "      <td>6.935452e-03</td>\n",
       "      <td>0.579079</td>\n",
       "      <td>CSK_HUMAN Tyrosine-protein kinase CSK</td>\n",
       "      <td>CSK</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P31150</th>\n",
       "      <td>7.056533e-03</td>\n",
       "      <td>0.568566</td>\n",
       "      <td>GDIA_HUMAN Rab GDP dissociation inhibitor alpha</td>\n",
       "      <td>GDIA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P27695</th>\n",
       "      <td>9.997328e-04</td>\n",
       "      <td>0.562257</td>\n",
       "      <td>APEX1_HUMAN DNA-(apurinic or apyrimidinic site...</td>\n",
       "      <td>APEX1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q06830</th>\n",
       "      <td>9.829618e-04</td>\n",
       "      <td>0.545459</td>\n",
       "      <td>PRDX1_HUMAN Peroxiredoxin-1</td>\n",
       "      <td>PRDX1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P51149</th>\n",
       "      <td>4.379513e-03</td>\n",
       "      <td>0.494201</td>\n",
       "      <td>RAB7A_HUMAN Ras-related protein Rab-7a</td>\n",
       "      <td>RAB7A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P61106</th>\n",
       "      <td>9.252077e-03</td>\n",
       "      <td>0.399791</td>\n",
       "      <td>RAB14_HUMAN Ras-related protein Rab-14</td>\n",
       "      <td>RAB14</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P51159</th>\n",
       "      <td>5.291783e-03</td>\n",
       "      <td>0.376507</td>\n",
       "      <td>RB27A_HUMAN Ras-related protein Rab-27A</td>\n",
       "      <td>RB27A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O14773</th>\n",
       "      <td>9.529256e-03</td>\n",
       "      <td>0.373186</td>\n",
       "      <td>TPP1_HUMAN Tripeptidyl-peptidase 1</td>\n",
       "      <td>TPP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P27797</th>\n",
       "      <td>4.034493e-03</td>\n",
       "      <td>0.354386</td>\n",
       "      <td>CALR_HUMAN Calreticulin</td>\n",
       "      <td>CALR</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                  pvalue  log2(B)-log2(T)  \\\n",
       "Protein ID                                  \n",
       "P09960      7.198480e-07         1.972804   \n",
       "P01834      1.092879e-03         1.913174   \n",
       "P29350      5.002525e-07         1.744085   \n",
       "P0DOY2      5.558317e-04         1.684692   \n",
       "P09525      1.857150e-04         1.525379   \n",
       "P09874      2.778706e-06         1.377227   \n",
       "P27707      2.038183e-04         1.325106   \n",
       "Q96KP4      2.237074e-05         1.321275   \n",
       "P58546      4.286378e-03         1.164319   \n",
       "P22392      2.731861e-03         0.981652   \n",
       "P59665      3.874038e-04         0.957458   \n",
       "P50395      8.590047e-05         0.941751   \n",
       "Q16666      3.008010e-03         0.835062   \n",
       "Q9UJU6      4.328426e-03         0.596015   \n",
       "P41240      6.935452e-03         0.579079   \n",
       "P31150      7.056533e-03         0.568566   \n",
       "P27695      9.997328e-04         0.562257   \n",
       "Q06830      9.829618e-04         0.545459   \n",
       "P51149      4.379513e-03         0.494201   \n",
       "P61106      9.252077e-03         0.399791   \n",
       "P51159      5.291783e-03         0.376507   \n",
       "O14773      9.529256e-03         0.373186   \n",
       "P27797      4.034493e-03         0.354386   \n",
       "\n",
       "                                                         Name GeneID  \n",
       "Protein ID                                                            \n",
       "P09960                 LKHA4_HUMAN Leukotriene A-4 hydrolase   LKHA4  \n",
       "P01834              IGKC_HUMAN Immunoglobulin kappa constant    IGKC  \n",
       "P29350      PTN6_HUMAN Tyrosine-protein phosphatase non-re...   PTN6  \n",
       "P0DOY2          IGLC2_HUMAN Immunoglobulin lambda constant 2   IGLC2  \n",
       "P09525                                ANXA4_HUMAN Annexin A4   ANXA4  \n",
       "P09874            PARP1_HUMAN Poly [ADP-ribose] polymerase 1   PARP1  \n",
       "P27707                        DCK_HUMAN Deoxycytidine kinase     DCK  \n",
       "Q96KP4        CNDP2_HUMAN Cytosolic non-specific dipeptidase   CNDP2  \n",
       "P58546                                 MTPN_HUMAN Myotrophin    MTPN  \n",
       "P22392            NDKB_HUMAN Nucleoside diphosphate kinase B    NDKB  \n",
       "P59665                      DEF1_HUMAN Neutrophil defensin 1    DEF1  \n",
       "P50395        GDIB_HUMAN Rab GDP dissociation inhibitor beta    GDIB  \n",
       "Q16666      IF16_HUMAN Gamma-interferon-inducible protein 16    IF16  \n",
       "Q9UJU6                       DBNL_HUMAN Drebrin-like protein    DBNL  \n",
       "P41240                 CSK_HUMAN Tyrosine-protein kinase CSK     CSK  \n",
       "P31150       GDIA_HUMAN Rab GDP dissociation inhibitor alpha    GDIA  \n",
       "P27695      APEX1_HUMAN DNA-(apurinic or apyrimidinic site...  APEX1  \n",
       "Q06830                           PRDX1_HUMAN Peroxiredoxin-1   PRDX1  \n",
       "P51149                RAB7A_HUMAN Ras-related protein Rab-7a   RAB7A  \n",
       "P61106                RAB14_HUMAN Ras-related protein Rab-14   RAB14  \n",
       "P51159               RB27A_HUMAN Ras-related protein Rab-27A   RB27A  \n",
       "O14773                    TPP1_HUMAN Tripeptidyl-peptidase 1    TPP1  \n",
       "P27797                               CALR_HUMAN Calreticulin    CALR  "
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "higher_in_B_cells"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Proteins Only in B Cells:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Name</th>\n",
       "      <th>GeneID</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Protein ID</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>O15400</th>\n",
       "      <td>TX7_HUMAN Syntaxin-7</td>\n",
       "      <td>TX7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O60262</th>\n",
       "      <td>GBG7_HUMAN Guanine nucleotide-binding protein ...</td>\n",
       "      <td>GBG7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O75223</th>\n",
       "      <td>GGCT_HUMAN Gamma-glutamylcyclotransferase</td>\n",
       "      <td>GGCT</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O75608</th>\n",
       "      <td>LYPA1_HUMAN Acyl-protein thioesterase 1</td>\n",
       "      <td>LYPA1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O75874</th>\n",
       "      <td>IDHC_HUMAN Isocitrate dehydrogenase [NADP] cyt...</td>\n",
       "      <td>IDHC</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P01880</th>\n",
       "      <td>IGHD_HUMAN Immunoglobulin heavy constant delta</td>\n",
       "      <td>IGHD</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q8TDP1</th>\n",
       "      <td>RNH2C_HUMAN Ribonuclease H2 subunit C</td>\n",
       "      <td>RNH2C</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UBR2</th>\n",
       "      <td>CATZ_HUMAN Cathepsin Z</td>\n",
       "      <td>CATZ</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O94874</th>\n",
       "      <td>UFL1_HUMAN E3 UFM1-protein ligase 1</td>\n",
       "      <td>UFL1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UKD2</th>\n",
       "      <td>MRT4_HUMAN mRNA turnover protein 4 homolog</td>\n",
       "      <td>MRT4</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>69 rows × 2 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                         Name GeneID\n",
       "Protein ID                                                          \n",
       "O15400                                  TX7_HUMAN Syntaxin-7     TX7\n",
       "O60262      GBG7_HUMAN Guanine nucleotide-binding protein ...   GBG7\n",
       "O75223             GGCT_HUMAN Gamma-glutamylcyclotransferase    GGCT\n",
       "O75608               LYPA1_HUMAN Acyl-protein thioesterase 1   LYPA1\n",
       "O75874      IDHC_HUMAN Isocitrate dehydrogenase [NADP] cyt...   IDHC\n",
       "...                                                       ...    ...\n",
       "P01880        IGHD_HUMAN Immunoglobulin heavy constant delta    IGHD\n",
       "Q8TDP1                 RNH2C_HUMAN Ribonuclease H2 subunit C   RNH2C\n",
       "Q9UBR2                                CATZ_HUMAN Cathepsin Z    CATZ\n",
       "O94874                   UFL1_HUMAN E3 UFM1-protein ligase 1    UFL1\n",
       "Q9UKD2            MRT4_HUMAN mRNA turnover protein 4 homolog    MRT4\n",
       "\n",
       "[69 rows x 2 columns]"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "B_cell_only_geneset[['Name', 'GeneID']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Proteins higher or only in B cells by GeneID\n",
    "<code>B_cell_set</code> is defined here as both the above lists of names, those higher or only in B cells."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "B_cell_set = list(higher_in_B_cells.GeneID)+list(B_cell_only_geneset.GeneID)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Proteins Higher in T Cells:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>pvalue</th>\n",
       "      <th>log2(B)-log2(T)</th>\n",
       "      <th>Name</th>\n",
       "      <th>GeneID</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Protein ID</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>P26447</th>\n",
       "      <td>0.000013</td>\n",
       "      <td>-3.571146</td>\n",
       "      <td>10A4_HUMAN Protein S100-A4</td>\n",
       "      <td>10A4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q15149</th>\n",
       "      <td>0.001746</td>\n",
       "      <td>-1.865700</td>\n",
       "      <td>PLEC_HUMAN Plectin</td>\n",
       "      <td>PLEC</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P12814</th>\n",
       "      <td>0.000009</td>\n",
       "      <td>-1.710082</td>\n",
       "      <td>ACTN1_HUMAN Alpha-actinin-1</td>\n",
       "      <td>ACTN1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P31949</th>\n",
       "      <td>0.000654</td>\n",
       "      <td>-1.699055</td>\n",
       "      <td>10AB_HUMAN Protein S100-A11</td>\n",
       "      <td>10AB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O15533</th>\n",
       "      <td>0.003014</td>\n",
       "      <td>-1.597947</td>\n",
       "      <td>TPSN_HUMAN Tapasin</td>\n",
       "      <td>TPSN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9Y2Q3</th>\n",
       "      <td>0.003051</td>\n",
       "      <td>-1.525526</td>\n",
       "      <td>GSTK1_HUMAN Glutathione S-transferase kappa 1</td>\n",
       "      <td>GSTK1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q99733</th>\n",
       "      <td>0.000316</td>\n",
       "      <td>-1.454706</td>\n",
       "      <td>NP1L4_HUMAN Nucleosome assembly protein 1-like 4</td>\n",
       "      <td>NP1L4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P08670</th>\n",
       "      <td>0.000004</td>\n",
       "      <td>-1.420240</td>\n",
       "      <td>VIME_HUMAN Vimentin</td>\n",
       "      <td>VIME</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P07195</th>\n",
       "      <td>0.000127</td>\n",
       "      <td>-1.408794</td>\n",
       "      <td>LDHB_HUMAN L-lactate dehydrogenase B chain</td>\n",
       "      <td>LDHB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O60234</th>\n",
       "      <td>0.000007</td>\n",
       "      <td>-1.345577</td>\n",
       "      <td>GMFG_HUMAN Glia maturation factor gamma</td>\n",
       "      <td>GMFG</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P19367</th>\n",
       "      <td>0.000103</td>\n",
       "      <td>-1.323451</td>\n",
       "      <td>HXK1_HUMAN Hexokinase-1</td>\n",
       "      <td>HXK1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q09666</th>\n",
       "      <td>0.001541</td>\n",
       "      <td>-1.280424</td>\n",
       "      <td>AHNK_HUMAN Neuroblast differentiation-associat...</td>\n",
       "      <td>AHNK</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O75083</th>\n",
       "      <td>0.000045</td>\n",
       "      <td>-1.218378</td>\n",
       "      <td>WDR1_HUMAN WD repeat-containing protein 1</td>\n",
       "      <td>WDR1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P48426</th>\n",
       "      <td>0.002104</td>\n",
       "      <td>-1.214703</td>\n",
       "      <td>PI42A_HUMAN Phosphatidylinositol 5-phosphate 4...</td>\n",
       "      <td>PI42A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9BSJ8</th>\n",
       "      <td>0.000016</td>\n",
       "      <td>-1.207312</td>\n",
       "      <td>ESYT1_HUMAN Extended synaptotagmin-1</td>\n",
       "      <td>ESYT1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P09382</th>\n",
       "      <td>0.009102</td>\n",
       "      <td>-1.126815</td>\n",
       "      <td>LEG1_HUMAN Galectin-1</td>\n",
       "      <td>LEG1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q15181</th>\n",
       "      <td>0.008095</td>\n",
       "      <td>-1.124519</td>\n",
       "      <td>IPYR_HUMAN Inorganic pyrophosphatase</td>\n",
       "      <td>IPYR</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q92522</th>\n",
       "      <td>0.000387</td>\n",
       "      <td>-1.078519</td>\n",
       "      <td>H1X_HUMAN Histone H1x</td>\n",
       "      <td>H1X</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P21333</th>\n",
       "      <td>0.000253</td>\n",
       "      <td>-1.027372</td>\n",
       "      <td>FLNA_HUMAN Filamin-A</td>\n",
       "      <td>FLNA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P62891</th>\n",
       "      <td>0.000297</td>\n",
       "      <td>-1.026259</td>\n",
       "      <td>RL39_HUMAN 60S ribosomal protein L39</td>\n",
       "      <td>RL39</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UI08</th>\n",
       "      <td>0.000089</td>\n",
       "      <td>-1.010118</td>\n",
       "      <td>EVL_HUMAN Ena/VASP-like protein</td>\n",
       "      <td>EVL</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P41091</th>\n",
       "      <td>0.005024</td>\n",
       "      <td>-0.967700</td>\n",
       "      <td>IF2G_HUMAN Eukaryotic translation initiation f...</td>\n",
       "      <td>IF2G</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P50238</th>\n",
       "      <td>0.002378</td>\n",
       "      <td>-0.956795</td>\n",
       "      <td>CRIP1_HUMAN Cysteine-rich protein 1</td>\n",
       "      <td>CRIP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P21291</th>\n",
       "      <td>0.008865</td>\n",
       "      <td>-0.934984</td>\n",
       "      <td>CSRP1_HUMAN Cysteine and glycine-rich protein 1</td>\n",
       "      <td>CSRP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P07237</th>\n",
       "      <td>0.006627</td>\n",
       "      <td>-0.901216</td>\n",
       "      <td>PDIA1_HUMAN Protein disulfide-isomerase</td>\n",
       "      <td>PDIA1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P84095</th>\n",
       "      <td>0.008719</td>\n",
       "      <td>-0.892711</td>\n",
       "      <td>RHOG_HUMAN Rho-related GTP-binding protein RhoG</td>\n",
       "      <td>RHOG</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P61224</th>\n",
       "      <td>0.005614</td>\n",
       "      <td>-0.868925</td>\n",
       "      <td>RAP1B_HUMAN Ras-related protein Rap-1b</td>\n",
       "      <td>RAP1B</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P35579</th>\n",
       "      <td>0.001202</td>\n",
       "      <td>-0.844896</td>\n",
       "      <td>MYH9_HUMAN Myosin-9</td>\n",
       "      <td>MYH9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P07355</th>\n",
       "      <td>0.004340</td>\n",
       "      <td>-0.841736</td>\n",
       "      <td>ANXA2_HUMAN Annexin A2</td>\n",
       "      <td>ANXA2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P0DP23</th>\n",
       "      <td>0.009853</td>\n",
       "      <td>-0.830721</td>\n",
       "      <td>CALM1_HUMAN Calmodulin-1</td>\n",
       "      <td>CALM1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P16070</th>\n",
       "      <td>0.000031</td>\n",
       "      <td>-0.821592</td>\n",
       "      <td>CD44_HUMAN CD44 antigen</td>\n",
       "      <td>CD44</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P30626</th>\n",
       "      <td>0.008634</td>\n",
       "      <td>-0.803462</td>\n",
       "      <td>RCN_HUMAN Sorcin</td>\n",
       "      <td>RCN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P60660</th>\n",
       "      <td>0.000241</td>\n",
       "      <td>-0.802568</td>\n",
       "      <td>MYL6_HUMAN Myosin light polypeptide 6</td>\n",
       "      <td>MYL6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q14005</th>\n",
       "      <td>0.002071</td>\n",
       "      <td>-0.790203</td>\n",
       "      <td>IL16_HUMAN Pro-interleukin-16</td>\n",
       "      <td>IL16</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q14847</th>\n",
       "      <td>0.000742</td>\n",
       "      <td>-0.780487</td>\n",
       "      <td>LASP1_HUMAN LIM and SH3 domain protein 1</td>\n",
       "      <td>LASP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P50995</th>\n",
       "      <td>0.000638</td>\n",
       "      <td>-0.759591</td>\n",
       "      <td>ANX11_HUMAN Annexin A11</td>\n",
       "      <td>ANX11</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O15127</th>\n",
       "      <td>0.008961</td>\n",
       "      <td>-0.739596</td>\n",
       "      <td>CAM2_HUMAN Secretory carrier-associated membra...</td>\n",
       "      <td>CAM2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9ULV4</th>\n",
       "      <td>0.004104</td>\n",
       "      <td>-0.680691</td>\n",
       "      <td>COR1C_HUMAN Coronin-1C</td>\n",
       "      <td>COR1C</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UQ80</th>\n",
       "      <td>0.002146</td>\n",
       "      <td>-0.680431</td>\n",
       "      <td>PA2G4_HUMAN Proliferation-associated protein 2G4</td>\n",
       "      <td>PA2G4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O43707</th>\n",
       "      <td>0.002915</td>\n",
       "      <td>-0.678984</td>\n",
       "      <td>ACTN4_HUMAN Alpha-actinin-4</td>\n",
       "      <td>ACTN4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P51608</th>\n",
       "      <td>0.001443</td>\n",
       "      <td>-0.678261</td>\n",
       "      <td>MECP2_HUMAN Methyl-CpG-binding protein 2</td>\n",
       "      <td>MECP2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UII2</th>\n",
       "      <td>0.008367</td>\n",
       "      <td>-0.657541</td>\n",
       "      <td>ATIF1_HUMAN ATPase inhibitor, mitochondrial</td>\n",
       "      <td>ATIF1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q15019</th>\n",
       "      <td>0.004856</td>\n",
       "      <td>-0.645156</td>\n",
       "      <td>EPT2_HUMAN Septin-2</td>\n",
       "      <td>EPT2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P52565</th>\n",
       "      <td>0.003804</td>\n",
       "      <td>-0.624032</td>\n",
       "      <td>GDIR1_HUMAN Rho GDP-dissociation inhibitor 1</td>\n",
       "      <td>GDIR1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P48735</th>\n",
       "      <td>0.004538</td>\n",
       "      <td>-0.620374</td>\n",
       "      <td>IDHP_HUMAN Isocitrate dehydrogenase [NADP], mi...</td>\n",
       "      <td>IDHP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P07741</th>\n",
       "      <td>0.001737</td>\n",
       "      <td>-0.601529</td>\n",
       "      <td>APT_HUMAN Adenine phosphoribosyltransferase</td>\n",
       "      <td>APT</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UHD8</th>\n",
       "      <td>0.004440</td>\n",
       "      <td>-0.600049</td>\n",
       "      <td>EPT9_HUMAN Septin-9</td>\n",
       "      <td>EPT9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P31948</th>\n",
       "      <td>0.002290</td>\n",
       "      <td>-0.543619</td>\n",
       "      <td>TIP1_HUMAN Stress-induced-phosphoprotein 1</td>\n",
       "      <td>TIP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q99439</th>\n",
       "      <td>0.002756</td>\n",
       "      <td>-0.540945</td>\n",
       "      <td>CNN2_HUMAN Calponin-2</td>\n",
       "      <td>CNN2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P62280</th>\n",
       "      <td>0.001128</td>\n",
       "      <td>-0.505126</td>\n",
       "      <td>RS11_HUMAN 40S ribosomal protein S11</td>\n",
       "      <td>RS11</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P78417</th>\n",
       "      <td>0.004740</td>\n",
       "      <td>-0.498805</td>\n",
       "      <td>GSTO1_HUMAN Glutathione S-transferase omega-1</td>\n",
       "      <td>GSTO1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O15144</th>\n",
       "      <td>0.004936</td>\n",
       "      <td>-0.473793</td>\n",
       "      <td>ARPC2_HUMAN Actin-related protein 2/3 complex ...</td>\n",
       "      <td>ARPC2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P26038</th>\n",
       "      <td>0.003892</td>\n",
       "      <td>-0.465095</td>\n",
       "      <td>MOES_HUMAN Moesin</td>\n",
       "      <td>MOES</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P59998</th>\n",
       "      <td>0.008911</td>\n",
       "      <td>-0.421006</td>\n",
       "      <td>ARPC4_HUMAN Actin-related protein 2/3 complex ...</td>\n",
       "      <td>ARPC4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P26583</th>\n",
       "      <td>0.008623</td>\n",
       "      <td>-0.419552</td>\n",
       "      <td>HMGB2_HUMAN High mobility group protein B2</td>\n",
       "      <td>HMGB2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P61158</th>\n",
       "      <td>0.008721</td>\n",
       "      <td>-0.413927</td>\n",
       "      <td>ARP3_HUMAN Actin-related protein 3</td>\n",
       "      <td>ARP3</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "              pvalue  log2(B)-log2(T)  \\\n",
       "Protein ID                              \n",
       "P26447      0.000013        -3.571146   \n",
       "Q15149      0.001746        -1.865700   \n",
       "P12814      0.000009        -1.710082   \n",
       "P31949      0.000654        -1.699055   \n",
       "O15533      0.003014        -1.597947   \n",
       "Q9Y2Q3      0.003051        -1.525526   \n",
       "Q99733      0.000316        -1.454706   \n",
       "P08670      0.000004        -1.420240   \n",
       "P07195      0.000127        -1.408794   \n",
       "O60234      0.000007        -1.345577   \n",
       "P19367      0.000103        -1.323451   \n",
       "Q09666      0.001541        -1.280424   \n",
       "O75083      0.000045        -1.218378   \n",
       "P48426      0.002104        -1.214703   \n",
       "Q9BSJ8      0.000016        -1.207312   \n",
       "P09382      0.009102        -1.126815   \n",
       "Q15181      0.008095        -1.124519   \n",
       "Q92522      0.000387        -1.078519   \n",
       "P21333      0.000253        -1.027372   \n",
       "P62891      0.000297        -1.026259   \n",
       "Q9UI08      0.000089        -1.010118   \n",
       "P41091      0.005024        -0.967700   \n",
       "P50238      0.002378        -0.956795   \n",
       "P21291      0.008865        -0.934984   \n",
       "P07237      0.006627        -0.901216   \n",
       "P84095      0.008719        -0.892711   \n",
       "P61224      0.005614        -0.868925   \n",
       "P35579      0.001202        -0.844896   \n",
       "P07355      0.004340        -0.841736   \n",
       "P0DP23      0.009853        -0.830721   \n",
       "P16070      0.000031        -0.821592   \n",
       "P30626      0.008634        -0.803462   \n",
       "P60660      0.000241        -0.802568   \n",
       "Q14005      0.002071        -0.790203   \n",
       "Q14847      0.000742        -0.780487   \n",
       "P50995      0.000638        -0.759591   \n",
       "O15127      0.008961        -0.739596   \n",
       "Q9ULV4      0.004104        -0.680691   \n",
       "Q9UQ80      0.002146        -0.680431   \n",
       "O43707      0.002915        -0.678984   \n",
       "P51608      0.001443        -0.678261   \n",
       "Q9UII2      0.008367        -0.657541   \n",
       "Q15019      0.004856        -0.645156   \n",
       "P52565      0.003804        -0.624032   \n",
       "P48735      0.004538        -0.620374   \n",
       "P07741      0.001737        -0.601529   \n",
       "Q9UHD8      0.004440        -0.600049   \n",
       "P31948      0.002290        -0.543619   \n",
       "Q99439      0.002756        -0.540945   \n",
       "P62280      0.001128        -0.505126   \n",
       "P78417      0.004740        -0.498805   \n",
       "O15144      0.004936        -0.473793   \n",
       "P26038      0.003892        -0.465095   \n",
       "P59998      0.008911        -0.421006   \n",
       "P26583      0.008623        -0.419552   \n",
       "P61158      0.008721        -0.413927   \n",
       "\n",
       "                                                         Name GeneID  \n",
       "Protein ID                                                            \n",
       "P26447                            10A4_HUMAN Protein S100-A4    10A4  \n",
       "Q15149                                    PLEC_HUMAN Plectin    PLEC  \n",
       "P12814                           ACTN1_HUMAN Alpha-actinin-1   ACTN1  \n",
       "P31949                           10AB_HUMAN Protein S100-A11    10AB  \n",
       "O15533                                    TPSN_HUMAN Tapasin    TPSN  \n",
       "Q9Y2Q3         GSTK1_HUMAN Glutathione S-transferase kappa 1   GSTK1  \n",
       "Q99733      NP1L4_HUMAN Nucleosome assembly protein 1-like 4   NP1L4  \n",
       "P08670                                   VIME_HUMAN Vimentin    VIME  \n",
       "P07195            LDHB_HUMAN L-lactate dehydrogenase B chain    LDHB  \n",
       "O60234               GMFG_HUMAN Glia maturation factor gamma    GMFG  \n",
       "P19367                               HXK1_HUMAN Hexokinase-1    HXK1  \n",
       "Q09666      AHNK_HUMAN Neuroblast differentiation-associat...   AHNK  \n",
       "O75083             WDR1_HUMAN WD repeat-containing protein 1    WDR1  \n",
       "P48426      PI42A_HUMAN Phosphatidylinositol 5-phosphate 4...  PI42A  \n",
       "Q9BSJ8                  ESYT1_HUMAN Extended synaptotagmin-1   ESYT1  \n",
       "P09382                                 LEG1_HUMAN Galectin-1    LEG1  \n",
       "Q15181                  IPYR_HUMAN Inorganic pyrophosphatase    IPYR  \n",
       "Q92522                                 H1X_HUMAN Histone H1x     H1X  \n",
       "P21333                                  FLNA_HUMAN Filamin-A    FLNA  \n",
       "P62891                  RL39_HUMAN 60S ribosomal protein L39    RL39  \n",
       "Q9UI08                       EVL_HUMAN Ena/VASP-like protein     EVL  \n",
       "P41091      IF2G_HUMAN Eukaryotic translation initiation f...   IF2G  \n",
       "P50238                   CRIP1_HUMAN Cysteine-rich protein 1   CRIP1  \n",
       "P21291       CSRP1_HUMAN Cysteine and glycine-rich protein 1   CSRP1  \n",
       "P07237               PDIA1_HUMAN Protein disulfide-isomerase   PDIA1  \n",
       "P84095       RHOG_HUMAN Rho-related GTP-binding protein RhoG    RHOG  \n",
       "P61224                RAP1B_HUMAN Ras-related protein Rap-1b   RAP1B  \n",
       "P35579                                   MYH9_HUMAN Myosin-9    MYH9  \n",
       "P07355                                ANXA2_HUMAN Annexin A2   ANXA2  \n",
       "P0DP23                              CALM1_HUMAN Calmodulin-1   CALM1  \n",
       "P16070                               CD44_HUMAN CD44 antigen    CD44  \n",
       "P30626                                      RCN_HUMAN Sorcin     RCN  \n",
       "P60660                 MYL6_HUMAN Myosin light polypeptide 6    MYL6  \n",
       "Q14005                         IL16_HUMAN Pro-interleukin-16    IL16  \n",
       "Q14847              LASP1_HUMAN LIM and SH3 domain protein 1   LASP1  \n",
       "P50995                               ANX11_HUMAN Annexin A11   ANX11  \n",
       "O15127      CAM2_HUMAN Secretory carrier-associated membra...   CAM2  \n",
       "Q9ULV4                                COR1C_HUMAN Coronin-1C   COR1C  \n",
       "Q9UQ80      PA2G4_HUMAN Proliferation-associated protein 2G4   PA2G4  \n",
       "O43707                           ACTN4_HUMAN Alpha-actinin-4   ACTN4  \n",
       "P51608              MECP2_HUMAN Methyl-CpG-binding protein 2   MECP2  \n",
       "Q9UII2           ATIF1_HUMAN ATPase inhibitor, mitochondrial   ATIF1  \n",
       "Q15019                                   EPT2_HUMAN Septin-2    EPT2  \n",
       "P52565          GDIR1_HUMAN Rho GDP-dissociation inhibitor 1   GDIR1  \n",
       "P48735      IDHP_HUMAN Isocitrate dehydrogenase [NADP], mi...   IDHP  \n",
       "P07741           APT_HUMAN Adenine phosphoribosyltransferase     APT  \n",
       "Q9UHD8                                   EPT9_HUMAN Septin-9    EPT9  \n",
       "P31948            TIP1_HUMAN Stress-induced-phosphoprotein 1    TIP1  \n",
       "Q99439                                 CNN2_HUMAN Calponin-2    CNN2  \n",
       "P62280                  RS11_HUMAN 40S ribosomal protein S11    RS11  \n",
       "P78417         GSTO1_HUMAN Glutathione S-transferase omega-1   GSTO1  \n",
       "O15144      ARPC2_HUMAN Actin-related protein 2/3 complex ...  ARPC2  \n",
       "P26038                                     MOES_HUMAN Moesin    MOES  \n",
       "P59998      ARPC4_HUMAN Actin-related protein 2/3 complex ...  ARPC4  \n",
       "P26583            HMGB2_HUMAN High mobility group protein B2   HMGB2  \n",
       "P61158                    ARP3_HUMAN Actin-related protein 3    ARP3  "
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "higher_in_T_cells.sort_values('log2(B)-log2(T)',ascending=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Proteins Only in T Cells:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Name</th>\n",
       "      <th>GeneID</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Protein ID</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>O00584</th>\n",
       "      <td>RNT2_HUMAN Ribonuclease T2</td>\n",
       "      <td>RNT2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O15117</th>\n",
       "      <td>FYB1_HUMAN FYN-binding protein 1</td>\n",
       "      <td>FYB1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O43665</th>\n",
       "      <td>RGS10_HUMAN Regulator of G-protein signaling 10</td>\n",
       "      <td>RGS10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O60496</th>\n",
       "      <td>DOK2_HUMAN Docking protein 2</td>\n",
       "      <td>DOK2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O60711</th>\n",
       "      <td>LPXN_HUMAN Leupaxin</td>\n",
       "      <td>LPXN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O95571</th>\n",
       "      <td>ETHE1_HUMAN Persulfide dioxygenase ETHE1, mito...</td>\n",
       "      <td>ETHE1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P06127</th>\n",
       "      <td>CD5_HUMAN T-cell surface glycoprotein CD5</td>\n",
       "      <td>CD5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P06729</th>\n",
       "      <td>CD2_HUMAN T-cell surface antigen CD2</td>\n",
       "      <td>CD2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P07305</th>\n",
       "      <td>H10_HUMAN Histone H1.0</td>\n",
       "      <td>H10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P07766</th>\n",
       "      <td>CD3E_HUMAN T-cell surface glycoprotein CD3 eps...</td>\n",
       "      <td>CD3E</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P12544</th>\n",
       "      <td>GRAA_HUMAN Granzyme A</td>\n",
       "      <td>GRAA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P30203</th>\n",
       "      <td>CD6_HUMAN T-cell differentiation antigen CD6</td>\n",
       "      <td>CD6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P32455</th>\n",
       "      <td>GBP1_HUMAN Guanylate-binding protein 1</td>\n",
       "      <td>GBP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P42229</th>\n",
       "      <td>TA5A_HUMAN Signal transducer and activator of ...</td>\n",
       "      <td>TA5A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P46976</th>\n",
       "      <td>GLYG_HUMAN Glycogenin-1</td>\n",
       "      <td>GLYG</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P49863</th>\n",
       "      <td>GRAK_HUMAN Granzyme K</td>\n",
       "      <td>GRAK</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P53634</th>\n",
       "      <td>CATC_HUMAN Dipeptidyl peptidase 1</td>\n",
       "      <td>CATC</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P54577</th>\n",
       "      <td>YYC_HUMAN Tyrosine--tRNA ligase, cytoplasmic</td>\n",
       "      <td>YYC</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q03519</th>\n",
       "      <td>TAP2_HUMAN Antigen peptide transporter 2</td>\n",
       "      <td>TAP2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q13045</th>\n",
       "      <td>FLII_HUMAN Protein flightless-1 homolog</td>\n",
       "      <td>FLII</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q13094</th>\n",
       "      <td>LCP2_HUMAN Lymphocyte cytosolic protein 2</td>\n",
       "      <td>LCP2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q15027</th>\n",
       "      <td>ACAP1_HUMAN Arf-GAP with coiled-coil, ANK repe...</td>\n",
       "      <td>ACAP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q8WWP7</th>\n",
       "      <td>GIMA1_HUMAN GTPase IMAP family member 1</td>\n",
       "      <td>GIMA1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q96B97</th>\n",
       "      <td>H3K1_HUMAN SH3 domain-containing kinase-bindin...</td>\n",
       "      <td>H3K1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q96C23</th>\n",
       "      <td>GALM_HUMAN Aldose 1-epimerase</td>\n",
       "      <td>GALM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9BUH6</th>\n",
       "      <td>PAXX_HUMAN Protein PAXX</td>\n",
       "      <td>PAXX</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9H9B4</th>\n",
       "      <td>FXN1_HUMAN Sideroflexin-1</td>\n",
       "      <td>FXN1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9NYL9</th>\n",
       "      <td>TMOD3_HUMAN Tropomodulin-3</td>\n",
       "      <td>TMOD3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9NZ08</th>\n",
       "      <td>ERAP1_HUMAN Endoplasmic reticulum aminopeptida...</td>\n",
       "      <td>ERAP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9Y3Z3</th>\n",
       "      <td>AMH1_HUMAN Deoxynucleoside triphosphate tripho...</td>\n",
       "      <td>AMH1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O94804</th>\n",
       "      <td>TK10_HUMAN Serine/threonine-protein kinase 10</td>\n",
       "      <td>TK10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P05023</th>\n",
       "      <td>AT1A1_HUMAN Sodium/potassium-transporting ATPa...</td>\n",
       "      <td>AT1A1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P06239</th>\n",
       "      <td>LCK_HUMAN Tyrosine-protein kinase Lck</td>\n",
       "      <td>LCK</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P23743</th>\n",
       "      <td>DGKA_HUMAN Diacylglycerol kinase alpha</td>\n",
       "      <td>DGKA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P41236</th>\n",
       "      <td>IPP2_HUMAN Protein phosphatase inhibitor 2</td>\n",
       "      <td>IPP2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P51812</th>\n",
       "      <td>KS6A3_HUMAN Ribosomal protein S6 kinase alpha-3</td>\n",
       "      <td>KS6A3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q08AF3</th>\n",
       "      <td>LFN5_HUMAN Schlafen family member 5</td>\n",
       "      <td>LFN5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9NUV9</th>\n",
       "      <td>GIMA4_HUMAN GTPase IMAP family member 4</td>\n",
       "      <td>GIMA4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P16150</th>\n",
       "      <td>LEUK_HUMAN Leukosialin</td>\n",
       "      <td>LEUK</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q969E4</th>\n",
       "      <td>TCAL3_HUMAN Transcription elongation factor A ...</td>\n",
       "      <td>TCAL3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q8NHV1</th>\n",
       "      <td>GIMA7_HUMAN GTPase IMAP family member 7</td>\n",
       "      <td>GIMA7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O60831</th>\n",
       "      <td>PRAF2_HUMAN PRA1 family protein 2</td>\n",
       "      <td>PRAF2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O60880</th>\n",
       "      <td>H21A_HUMAN SH2 domain-containing protein 1A</td>\n",
       "      <td>H21A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O95831</th>\n",
       "      <td>AIFM1_HUMAN Apoptosis-inducing factor 1, mitoc...</td>\n",
       "      <td>AIFM1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P01732</th>\n",
       "      <td>CD8A_HUMAN T-cell surface glycoprotein CD8 alp...</td>\n",
       "      <td>CD8A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P02795</th>\n",
       "      <td>MT2_HUMAN Metallothionein-2</td>\n",
       "      <td>MT2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P04234</th>\n",
       "      <td>CD3D_HUMAN T-cell surface glycoprotein CD3 del...</td>\n",
       "      <td>CD3D</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P13164</th>\n",
       "      <td>IFM1_HUMAN Interferon-induced transmembrane pr...</td>\n",
       "      <td>IFM1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P37235</th>\n",
       "      <td>HPCL1_HUMAN Hippocalcin-like protein 1</td>\n",
       "      <td>HPCL1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q15438</th>\n",
       "      <td>CYH1_HUMAN Cytohesin-1</td>\n",
       "      <td>CYH1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q86WV1</th>\n",
       "      <td>KAP1_HUMAN Src kinase-associated phosphoprotei...</td>\n",
       "      <td>KAP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9BU61</th>\n",
       "      <td>NDUF3_HUMAN NADH dehydrogenase [ubiquinone] 1 ...</td>\n",
       "      <td>NDUF3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UHL4</th>\n",
       "      <td>DPP2_HUMAN Dipeptidyl peptidase 2</td>\n",
       "      <td>DPP2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O95834</th>\n",
       "      <td>EMAL2_HUMAN Echinoderm microtubule-associated ...</td>\n",
       "      <td>EMAL2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P09564</th>\n",
       "      <td>CD7_HUMAN T-cell antigen CD7</td>\n",
       "      <td>CD7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O15260</th>\n",
       "      <td>URF4_HUMAN Surfeit locus protein 4</td>\n",
       "      <td>URF4</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                         Name GeneID\n",
       "Protein ID                                                          \n",
       "O00584                            RNT2_HUMAN Ribonuclease T2    RNT2\n",
       "O15117                      FYB1_HUMAN FYN-binding protein 1    FYB1\n",
       "O43665       RGS10_HUMAN Regulator of G-protein signaling 10   RGS10\n",
       "O60496                          DOK2_HUMAN Docking protein 2    DOK2\n",
       "O60711                                   LPXN_HUMAN Leupaxin    LPXN\n",
       "O95571      ETHE1_HUMAN Persulfide dioxygenase ETHE1, mito...  ETHE1\n",
       "P06127             CD5_HUMAN T-cell surface glycoprotein CD5     CD5\n",
       "P06729                  CD2_HUMAN T-cell surface antigen CD2     CD2\n",
       "P07305                                H10_HUMAN Histone H1.0     H10\n",
       "P07766      CD3E_HUMAN T-cell surface glycoprotein CD3 eps...   CD3E\n",
       "P12544                                 GRAA_HUMAN Granzyme A    GRAA\n",
       "P30203          CD6_HUMAN T-cell differentiation antigen CD6     CD6\n",
       "P32455                GBP1_HUMAN Guanylate-binding protein 1    GBP1\n",
       "P42229      TA5A_HUMAN Signal transducer and activator of ...   TA5A\n",
       "P46976                               GLYG_HUMAN Glycogenin-1    GLYG\n",
       "P49863                                 GRAK_HUMAN Granzyme K    GRAK\n",
       "P53634                     CATC_HUMAN Dipeptidyl peptidase 1    CATC\n",
       "P54577          YYC_HUMAN Tyrosine--tRNA ligase, cytoplasmic     YYC\n",
       "Q03519              TAP2_HUMAN Antigen peptide transporter 2    TAP2\n",
       "Q13045               FLII_HUMAN Protein flightless-1 homolog    FLII\n",
       "Q13094             LCP2_HUMAN Lymphocyte cytosolic protein 2    LCP2\n",
       "Q15027      ACAP1_HUMAN Arf-GAP with coiled-coil, ANK repe...  ACAP1\n",
       "Q8WWP7               GIMA1_HUMAN GTPase IMAP family member 1   GIMA1\n",
       "Q96B97      H3K1_HUMAN SH3 domain-containing kinase-bindin...   H3K1\n",
       "Q96C23                         GALM_HUMAN Aldose 1-epimerase    GALM\n",
       "Q9BUH6                               PAXX_HUMAN Protein PAXX    PAXX\n",
       "Q9H9B4                             FXN1_HUMAN Sideroflexin-1    FXN1\n",
       "Q9NYL9                            TMOD3_HUMAN Tropomodulin-3   TMOD3\n",
       "Q9NZ08      ERAP1_HUMAN Endoplasmic reticulum aminopeptida...  ERAP1\n",
       "Q9Y3Z3      AMH1_HUMAN Deoxynucleoside triphosphate tripho...   AMH1\n",
       "O94804         TK10_HUMAN Serine/threonine-protein kinase 10    TK10\n",
       "P05023      AT1A1_HUMAN Sodium/potassium-transporting ATPa...  AT1A1\n",
       "P06239                 LCK_HUMAN Tyrosine-protein kinase Lck     LCK\n",
       "P23743                DGKA_HUMAN Diacylglycerol kinase alpha    DGKA\n",
       "P41236            IPP2_HUMAN Protein phosphatase inhibitor 2    IPP2\n",
       "P51812       KS6A3_HUMAN Ribosomal protein S6 kinase alpha-3   KS6A3\n",
       "Q08AF3                   LFN5_HUMAN Schlafen family member 5    LFN5\n",
       "Q9NUV9               GIMA4_HUMAN GTPase IMAP family member 4   GIMA4\n",
       "P16150                                LEUK_HUMAN Leukosialin    LEUK\n",
       "Q969E4      TCAL3_HUMAN Transcription elongation factor A ...  TCAL3\n",
       "Q8NHV1               GIMA7_HUMAN GTPase IMAP family member 7   GIMA7\n",
       "O60831                     PRAF2_HUMAN PRA1 family protein 2   PRAF2\n",
       "O60880           H21A_HUMAN SH2 domain-containing protein 1A    H21A\n",
       "O95831      AIFM1_HUMAN Apoptosis-inducing factor 1, mitoc...  AIFM1\n",
       "P01732      CD8A_HUMAN T-cell surface glycoprotein CD8 alp...   CD8A\n",
       "P02795                           MT2_HUMAN Metallothionein-2     MT2\n",
       "P04234      CD3D_HUMAN T-cell surface glycoprotein CD3 del...   CD3D\n",
       "P13164      IFM1_HUMAN Interferon-induced transmembrane pr...   IFM1\n",
       "P37235                HPCL1_HUMAN Hippocalcin-like protein 1   HPCL1\n",
       "Q15438                                CYH1_HUMAN Cytohesin-1    CYH1\n",
       "Q86WV1      KAP1_HUMAN Src kinase-associated phosphoprotei...   KAP1\n",
       "Q9BU61      NDUF3_HUMAN NADH dehydrogenase [ubiquinone] 1 ...  NDUF3\n",
       "Q9UHL4                     DPP2_HUMAN Dipeptidyl peptidase 2    DPP2\n",
       "O95834      EMAL2_HUMAN Echinoderm microtubule-associated ...  EMAL2\n",
       "P09564                          CD7_HUMAN T-cell antigen CD7     CD7\n",
       "O15260                    URF4_HUMAN Surfeit locus protein 4    URF4"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "T_cell_only_geneset[['Name', 'GeneID']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Proteins higher or only in T cells by GeneID\n",
    "<code>T_cell_set</code> is defined here as both the above lists of names, those higher or only in T cells."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "T_cell_set = list(higher_in_T_cells.GeneID)+list(T_cell_only_geneset.GeneID)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.8.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}