5_characterize_iModulons/4_motif_search.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from pymodulon.io import *\n",
    "from pymodulon.motif import *"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Load Data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "ica_data = load_json_model(os.path.join('..','data','processed_data','bsu.json.gz'))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "fasta_file = os.path.join('..','data','external','genome.fasta')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Search for motifs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Finding motifs for 16 sequences\n",
      "Found 5 motifs across 35 sites\n",
      "Finding motifs for 82 sequences\n",
      "Found 3 motifs across 188 sites\n",
      "Finding motifs for 13 sequences\n",
      "Found 5 motifs across 32 sites\n",
      "Finding motifs for 149 sequences\n",
      "Found 2 motifs across 161 sites\n",
      "Finding motifs for 3 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 4 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 5 sequences\n",
      "Found 2 motifs across 9 sites\n",
      "Finding motifs for 5 sequences\n",
      "Found 4 motifs across 18 sites\n",
      "Finding motifs for 82 sequences\n",
      "Found 5 motifs across 142 sites\n",
      "Finding motifs for 164 sequences\n",
      "Found 3 motifs across 239 sites\n",
      "Finding motifs for 16 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 3 sequences\n",
      "Found 3 motifs across 6 sites\n",
      "Finding motifs for 48 sequences\n",
      "Found 4 motifs across 139 sites\n",
      "Less than two sequences found for iModulon: PyrR\n",
      "Finding motifs for 17 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 37 sequences\n",
      "Found 5 motifs across 105 sites\n",
      "Finding motifs for 2 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Less than two sequences found for iModulon: SG_1\n",
      "Finding motifs for 14 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 108 sequences\n",
      "Found 5 motifs across 301 sites\n",
      "Finding motifs for 2 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 6 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 34 sequences\n",
      "Found 5 motifs across 115 sites\n",
      "Finding motifs for 6 sequences\n",
      "Found 4 motifs across 17 sites\n",
      "Finding motifs for 5 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 50 sequences\n",
      "Found 5 motifs across 138 sites\n",
      "Finding motifs for 41 sequences\n",
      "Found 5 motifs across 88 sites\n",
      "Finding motifs for 20 sequences\n",
      "Found 2 motifs across 23 sites\n",
      "Finding motifs for 9 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Less than two sequences found for iModulon: MtlR\n",
      "Finding motifs for 17 sequences\n",
      "Found 5 motifs across 49 sites\n",
      "Finding motifs for 52 sequences\n",
      "Found 3 motifs across 68 sites\n",
      "Finding motifs for 18 sequences\n",
      "Found 5 motifs across 35 sites\n",
      "Less than two sequences found for iModulon: SwrA\n",
      "Finding motifs for 7 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 21 sequences\n",
      "Found 5 motifs across 40 sites\n",
      "Finding motifs for 98 sequences\n",
      "Found 5 motifs across 244 sites\n",
      "Finding motifs for 5 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 10 sequences\n",
      "Found 3 motifs across 15 sites\n",
      "Finding motifs for 14 sequences\n",
      "Found 1 motif across 9 sites\n",
      "Finding motifs for 101 sequences\n",
      "Found 3 motifs across 187 sites\n",
      "Finding motifs for 54 sequences\n",
      "Found 5 motifs across 111 sites\n",
      "Finding motifs for 69 sequences\n",
      "Found 3 motifs across 189 sites\n",
      "Less than two sequences found for iModulon: FruR\n",
      "Less than two sequences found for iModulon: AcoR\n",
      "Finding motifs for 8 sequences\n",
      "Found 1 motif across 7 sites\n",
      "Finding motifs for 19 sequences\n",
      "Found 5 motifs across 33 sites\n",
      "Finding motifs for 39 sequences\n",
      "Found 5 motifs across 71 sites\n",
      "Finding motifs for 24 sequences\n",
      "Found 4 motifs across 47 sites\n",
      "Finding motifs for 53 sequences\n",
      "Found 4 motifs across 152 sites\n",
      "Finding motifs for 4 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 4 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 10 sequences\n",
      "Found 3 motifs across 14 sites\n",
      "Finding motifs for 4 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 3 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 35 sequences\n",
      "Found 2 motifs across 49 sites\n",
      "Finding motifs for 33 sequences\n",
      "Found 2 motifs across 36 sites\n",
      "Finding motifs for 74 sequences\n",
      "Found 4 motifs across 195 sites\n",
      "Finding motifs for 80 sequences\n",
      "Found 5 motifs across 146 sites\n",
      "Finding motifs for 8 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 12 sequences\n",
      "Found 2 motifs across 9 sites\n",
      "Finding motifs for 23 sequences\n",
      "Found 2 motifs across 24 sites\n",
      "Finding motifs for 5 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 2 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 6 sequences\n",
      "Found 1 motif across 5 sites\n",
      "Finding motifs for 15 sequences\n",
      "Found 5 motifs across 33 sites\n",
      "Finding motifs for 17 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 2 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 54 sequences\n",
      "Found 2 motifs across 81 sites\n",
      "Finding motifs for 4 sequences\n",
      "No motif found with E-value < 1.0e-03\n",
      "Finding motifs for 71 sequences\n",
      "Found 2 motifs across 112 sites\n",
      "Finding motifs for 11 sequences\n",
      "Found 5 motifs across 37 sites\n"
     ]
    }
   ],
   "source": [
    "for k in ica_data.imodulon_names:\n",
    "    find_motifs(ica_data, k, fasta_file, force=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'S-box': <MotifInfo with 5 motifs across 35 sites>,\n",
       " 1: <MotifInfo with 3 motifs across 188 sites>,\n",
       " 'WalR': <MotifInfo with 5 motifs across 32 sites>,\n",
       " 3: <MotifInfo with 2 motifs across 161 sites>,\n",
       " 6: <MotifInfo with 2 motifs across 9 sites>,\n",
       " 7: <MotifInfo with 4 motifs across 18 sites>,\n",
       " 'CsoR/Fnr': <MotifInfo with 5 motifs across 142 sites>,\n",
       " 9: <MotifInfo with 3 motifs across 239 sites>,\n",
       " 'EAR_riboswitch': <MotifInfo with 3 motifs across 6 sites>,\n",
       " 'SigB-1': <MotifInfo with 4 motifs across 139 sites>,\n",
       " 'CodY': <MotifInfo with 5 motifs across 105 sites>,\n",
       " 'RhgR': <MotifInfo with 5 motifs across 301 sites>,\n",
       " 'stringent_response': <MotifInfo with 5 motifs across 115 sites>,\n",
       " 'SG_3': <MotifInfo with 4 motifs across 17 sites>,\n",
       " 'SigB-2': <MotifInfo with 5 motifs across 138 sites>,\n",
       " 26: <MotifInfo with 5 motifs across 88 sites>,\n",
       " 'TnrA-1': <MotifInfo with 2 motifs across 23 sites>,\n",
       " 'ImmR': <MotifInfo with 5 motifs across 49 sites>,\n",
       " 'SigD': <MotifInfo with 3 motifs across 68 sites>,\n",
       " 'TnrA-2': <MotifInfo with 5 motifs across 35 sites>,\n",
       " 'ComA': <MotifInfo with 5 motifs across 40 sites>,\n",
       " 'SigK': <MotifInfo with 5 motifs across 244 sites>,\n",
       " 'Zur': <MotifInfo with 3 motifs across 15 sites>,\n",
       " 'Thi-box': <MotifInfo with 1 motif across 9 sites>,\n",
       " 'NusA': <MotifInfo with 3 motifs across 187 sites>,\n",
       " 'CsoR': <MotifInfo with 5 motifs across 111 sites>,\n",
       " 'SigB-3': <MotifInfo with 3 motifs across 189 sites>,\n",
       " 'LnrK': <MotifInfo with 1 motif across 7 sites>,\n",
       " 'DegU': <MotifInfo with 5 motifs across 33 sites>,\n",
       " 'Xpf': <MotifInfo with 5 motifs across 71 sites>,\n",
       " 'Fur': <MotifInfo with 4 motifs across 47 sites>,\n",
       " 'CcpA-1': <MotifInfo with 4 motifs across 152 sites>,\n",
       " 52: <MotifInfo with 3 motifs across 14 sites>,\n",
       " 'ResD': <MotifInfo with 2 motifs across 49 sites>,\n",
       " 'SigW': <MotifInfo with 2 motifs across 36 sites>,\n",
       " 'SigE': <MotifInfo with 4 motifs across 195 sites>,\n",
       " 'Rok': <MotifInfo with 5 motifs across 146 sites>,\n",
       " 60: <MotifInfo with 2 motifs across 9 sites>,\n",
       " 61: <MotifInfo with 2 motifs across 24 sites>,\n",
       " 'FapR': <MotifInfo with 1 motif across 5 sites>,\n",
       " 'AbrB': <MotifInfo with 5 motifs across 33 sites>,\n",
       " 68: <MotifInfo with 2 motifs across 81 sites>,\n",
       " 70: <MotifInfo with 2 motifs across 112 sites>,\n",
       " 71: <MotifInfo with 5 motifs across 37 sites>}"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ica_data.motif_info"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Compare motifs to known databases using TOMTOM"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [],
   "source": [
    "for info in ica_data.motif_info.values():\n",
    "    compare_motifs(info)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1\n",
      "3\n",
      "CsoR/Fnr\n",
      "9\n",
      "SigB-1\n",
      "CodY\n",
      "RhgR\n",
      "stringent_response\n",
      "SigB-2\n",
      "SigK\n",
      "NusA\n",
      "SigB-3\n",
      "Fur\n",
      "CcpA-1\n",
      "SigW\n",
      "SigE\n",
      "Rok\n"
     ]
    }
   ],
   "source": [
    "for k in ica_data.motif_info.keys():\n",
    "    if not ica_data.motif_info[k].matches.empty:\n",
    "        print(k)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "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>motif</th>\n",
       "      <th>target</th>\n",
       "      <th>target_id</th>\n",
       "      <th>pvalue</th>\n",
       "      <th>Evalue</th>\n",
       "      <th>qvalue</th>\n",
       "      <th>overlap</th>\n",
       "      <th>optimal_offset</th>\n",
       "      <th>orientation</th>\n",
       "      <th>database</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>MEME-1</td>\n",
       "      <td>PhoP</td>\n",
       "      <td>MX000057</td>\n",
       "      <td>9.800910e-07</td>\n",
       "      <td>0.000579</td>\n",
       "      <td>0.000571</td>\n",
       "      <td>29</td>\n",
       "      <td>-3</td>\n",
       "      <td>+</td>\n",
       "      <td>prodoric</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>MEME-1</td>\n",
       "      <td>OxyR_P.aeruginosa</td>\n",
       "      <td>EXPREG_00001560</td>\n",
       "      <td>1.040420e-06</td>\n",
       "      <td>0.000615</td>\n",
       "      <td>0.000571</td>\n",
       "      <td>21</td>\n",
       "      <td>19</td>\n",
       "      <td>+</td>\n",
       "      <td>collectf</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>MEME-1</td>\n",
       "      <td>lrp</td>\n",
       "      <td>lrp</td>\n",
       "      <td>1.831820e-06</td>\n",
       "      <td>0.001083</td>\n",
       "      <td>0.000670</td>\n",
       "      <td>23</td>\n",
       "      <td>2</td>\n",
       "      <td>+</td>\n",
       "      <td>dpinteract</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>MEME-1</td>\n",
       "      <td>rpoS17</td>\n",
       "      <td>rpoS17</td>\n",
       "      <td>3.302720e-06</td>\n",
       "      <td>0.001952</td>\n",
       "      <td>0.000749</td>\n",
       "      <td>23</td>\n",
       "      <td>6</td>\n",
       "      <td>-</td>\n",
       "      <td>dpinteract</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>MEME-1</td>\n",
       "      <td>rpoD19</td>\n",
       "      <td>rpoD19</td>\n",
       "      <td>4.501930e-06</td>\n",
       "      <td>0.002661</td>\n",
       "      <td>0.000823</td>\n",
       "      <td>25</td>\n",
       "      <td>-7</td>\n",
       "      <td>+</td>\n",
       "      <td>dpinteract</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    motif             target        target_id        pvalue    Evalue  \\\n",
       "0  MEME-1               PhoP         MX000057  9.800910e-07  0.000579   \n",
       "1  MEME-1  OxyR_P.aeruginosa  EXPREG_00001560  1.040420e-06  0.000615   \n",
       "2  MEME-1                lrp              lrp  1.831820e-06  0.001083   \n",
       "3  MEME-1             rpoS17           rpoS17  3.302720e-06  0.001952   \n",
       "4  MEME-1             rpoD19           rpoD19  4.501930e-06  0.002661   \n",
       "\n",
       "     qvalue  overlap  optimal_offset orientation    database  \n",
       "0  0.000571       29              -3           +    prodoric  \n",
       "1  0.000571       21              19           +    collectf  \n",
       "2  0.000670       23               2           +  dpinteract  \n",
       "3  0.000749       23               6           -  dpinteract  \n",
       "4  0.000823       25              -7           +  dpinteract  "
      ]
     },
     "execution_count": 33,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ica_data.motif_info[1].matches"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "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>regulator</th>\n",
       "      <th>pvalue</th>\n",
       "      <th>qvalue</th>\n",
       "      <th>precision</th>\n",
       "      <th>recall</th>\n",
       "      <th>f1score</th>\n",
       "      <th>TP</th>\n",
       "      <th>regulon_size</th>\n",
       "      <th>imodulon_size</th>\n",
       "      <th>n_regs</th>\n",
       "      <th>single_gene</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>SigB-1</th>\n",
       "      <td>SigB</td>\n",
       "      <td>4.367180e-66</td>\n",
       "      <td>6.987488e-65</td>\n",
       "      <td>0.962264</td>\n",
       "      <td>0.230769</td>\n",
       "      <td>0.372263</td>\n",
       "      <td>51.0</td>\n",
       "      <td>221.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CodY</th>\n",
       "      <td>CodY</td>\n",
       "      <td>2.915811e-73</td>\n",
       "      <td>6.997947e-72</td>\n",
       "      <td>0.821918</td>\n",
       "      <td>0.320856</td>\n",
       "      <td>0.461538</td>\n",
       "      <td>60.0</td>\n",
       "      <td>187.0</td>\n",
       "      <td>73.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SigW</th>\n",
       "      <td>SigW</td>\n",
       "      <td>8.102748e-84</td>\n",
       "      <td>7.292473e-83</td>\n",
       "      <td>0.872727</td>\n",
       "      <td>0.607595</td>\n",
       "      <td>0.716418</td>\n",
       "      <td>48.0</td>\n",
       "      <td>79.0</td>\n",
       "      <td>55.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CcpA-1</th>\n",
       "      <td>CcpA</td>\n",
       "      <td>7.893839e-56</td>\n",
       "      <td>2.683905e-54</td>\n",
       "      <td>0.719512</td>\n",
       "      <td>0.225191</td>\n",
       "      <td>0.343023</td>\n",
       "      <td>59.0</td>\n",
       "      <td>262.0</td>\n",
       "      <td>82.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Fur</th>\n",
       "      <td>Fur</td>\n",
       "      <td>1.213043e-74</td>\n",
       "      <td>1.334347e-73</td>\n",
       "      <td>0.950000</td>\n",
       "      <td>0.622951</td>\n",
       "      <td>0.752475</td>\n",
       "      <td>38.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       regulator        pvalue        qvalue  precision    recall   f1score  \\\n",
       "SigB-1      SigB  4.367180e-66  6.987488e-65   0.962264  0.230769  0.372263   \n",
       "CodY        CodY  2.915811e-73  6.997947e-72   0.821918  0.320856  0.461538   \n",
       "SigW        SigW  8.102748e-84  7.292473e-83   0.872727  0.607595  0.716418   \n",
       "CcpA-1      CcpA  7.893839e-56  2.683905e-54   0.719512  0.225191  0.343023   \n",
       "Fur          Fur  1.213043e-74  1.334347e-73   0.950000  0.622951  0.752475   \n",
       "\n",
       "          TP  regulon_size  imodulon_size  n_regs  single_gene  \n",
       "SigB-1  51.0         221.0           53.0     1.0          NaN  \n",
       "CodY    60.0         187.0           73.0     1.0          NaN  \n",
       "SigW    48.0          79.0           55.0     1.0          NaN  \n",
       "CcpA-1  59.0         262.0           82.0     1.0          NaN  \n",
       "Fur     38.0          61.0           40.0     1.0          NaN  "
      ]
     },
     "execution_count": 58,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ica_data.imodulon_table.loc[['SigB-1','CodY','SigW','CcpA-1','Fur']]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "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>rel_position</th>\n",
       "      <th>pvalue</th>\n",
       "      <th>site_seq</th>\n",
       "      <th>genes</th>\n",
       "      <th>locus_tags</th>\n",
       "      <th>start_pos</th>\n",
       "      <th>strand</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>operon</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>Op1002</th>\n",
       "      <td>163</td>\n",
       "      <td>3.03e-06</td>\n",
       "      <td>CCTATTTTAGAAAAGTGAAATAAAATA</td>\n",
       "      <td>steT</td>\n",
       "      <td>BSU_12860</td>\n",
       "      <td>1352591</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op1009</th>\n",
       "      <td>455</td>\n",
       "      <td>6.59e-08</td>\n",
       "      <td>AATTTGTTAGAATATTCATAATTTAGT</td>\n",
       "      <td>dppA,dppB,dppC,dppD,dppE</td>\n",
       "      <td>BSU_12920,BSU_12930,BSU_12940,BSU_12950,BSU_12960</td>\n",
       "      <td>1359901</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op1010</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>ldcA,aeeB,eepC,ykfD</td>\n",
       "      <td>BSU_12970,BSU_12980,BSU_12990,BSU_13000</td>\n",
       "      <td>1365388</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op1028</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>ispA</td>\n",
       "      <td>BSU_13190</td>\n",
       "      <td>1386883</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op1094</th>\n",
       "      <td>126</td>\n",
       "      <td>5.89e-08</td>\n",
       "      <td>ACGAAGATTGAATTTTCAGTAAATATA</td>\n",
       "      <td>ykwB</td>\n",
       "      <td>BSU_13940</td>\n",
       "      <td>1463393</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op1196</th>\n",
       "      <td>446</td>\n",
       "      <td>6.77e-07</td>\n",
       "      <td>AAGTTTTATGAAAATTAAAACAATTAT</td>\n",
       "      <td>ylmA</td>\n",
       "      <td>BSU_15340</td>\n",
       "      <td>1606060</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op1296</th>\n",
       "      <td>427</td>\n",
       "      <td>3.29e-08</td>\n",
       "      <td>TAAAAGTTTGAATATTATGAAAATATA</td>\n",
       "      <td>acpK,pksF,pksG,pksH</td>\n",
       "      <td>BSU_17130,BSU_17140,BSU_17150,BSU_17160</td>\n",
       "      <td>1787969</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op1505</th>\n",
       "      <td>406</td>\n",
       "      <td>1.50e-06</td>\n",
       "      <td>AATTATTCAGAAAATTATTTATAATAA</td>\n",
       "      <td>yocS</td>\n",
       "      <td>BSU_19350</td>\n",
       "      <td>2105990</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op1800</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>yppF</td>\n",
       "      <td>BSU_22260</td>\n",
       "      <td>2338082</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op181</th>\n",
       "      <td>451</td>\n",
       "      <td>3.05e-07</td>\n",
       "      <td>AAAAAGACTGAATATTTAAACAATTAT</td>\n",
       "      <td>ilvE</td>\n",
       "      <td>BSU_02390</td>\n",
       "      <td>258516</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op226</th>\n",
       "      <td>191</td>\n",
       "      <td>1.13e-07</td>\n",
       "      <td>ACATATTGCTACTATTCAGAATAATTA</td>\n",
       "      <td>amhX</td>\n",
       "      <td>BSU_03010</td>\n",
       "      <td>325089</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op227</th>\n",
       "      <td>441</td>\n",
       "      <td>1.13e-07</td>\n",
       "      <td>ACATATTGCTACTATTCAGAATAATTA</td>\n",
       "      <td>ycgA</td>\n",
       "      <td>BSU_03020</td>\n",
       "      <td>324839</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op2347</th>\n",
       "      <td>454</td>\n",
       "      <td>9.40e-07</td>\n",
       "      <td>TAATTGACAGAATATTTTAACATAATT</td>\n",
       "      <td>braB</td>\n",
       "      <td>BSU_29600</td>\n",
       "      <td>3027797</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op243</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>putB,putC,putP</td>\n",
       "      <td>BSU_03200,BSU_03210,BSU_03220</td>\n",
       "      <td>344051</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op2486</th>\n",
       "      <td>155</td>\n",
       "      <td>1.29e-06</td>\n",
       "      <td>AAAAATAGATAATTTTTAAAAAATTAT</td>\n",
       "      <td>nupN,nupO,nupP,nupQ</td>\n",
       "      <td>BSU_31540,BSU_31550,BSU_31560,BSU_31570</td>\n",
       "      <td>3239430</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op2519</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>spsC</td>\n",
       "      <td>BSU_32070</td>\n",
       "      <td>3299140</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op2520</th>\n",
       "      <td>308</td>\n",
       "      <td>4.66e-10</td>\n",
       "      <td>CCTTTTTCTTAATTTTCTGATAATATA</td>\n",
       "      <td>yuiB</td>\n",
       "      <td>BSU_32080</td>\n",
       "      <td>3299564</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op2521</th>\n",
       "      <td>111</td>\n",
       "      <td>4.66e-10</td>\n",
       "      <td>CCTTTTTCTTAATTTTCTGATAATATA</td>\n",
       "      <td>yuiA</td>\n",
       "      <td>BSU_32090</td>\n",
       "      <td>3299761</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op2553</th>\n",
       "      <td>442</td>\n",
       "      <td>2.47e-06</td>\n",
       "      <td>ACTAATTTACAATTTTAATGAAAAATA</td>\n",
       "      <td>pucG,pucF</td>\n",
       "      <td>BSU_32520,BSU_32530</td>\n",
       "      <td>3343571</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op2556</th>\n",
       "      <td>185</td>\n",
       "      <td>9.40e-07</td>\n",
       "      <td>CATTTTTGCTATATTTCTAAAATAAAA</td>\n",
       "      <td>frlP</td>\n",
       "      <td>BSU_32550</td>\n",
       "      <td>3346016</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op2558</th>\n",
       "      <td>184</td>\n",
       "      <td>6.59e-08</td>\n",
       "      <td>TTATTTTCTGAAATTTCTCAATAATTT</td>\n",
       "      <td>frlD,frlM,frlN,frlO,frlB</td>\n",
       "      <td>BSU_32570,BSU_32580,BSU_32590,BSU_32600,BSU_32610</td>\n",
       "      <td>3351996</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op2935</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>rocC,rocB</td>\n",
       "      <td>BSU_37760,BSU_37770</td>\n",
       "      <td>3878792</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op2952</th>\n",
       "      <td>427</td>\n",
       "      <td>2.89e-09</td>\n",
       "      <td>GATTTTTTTGAATTTTCAAGAAATATA</td>\n",
       "      <td>vpr</td>\n",
       "      <td>BSU_38090</td>\n",
       "      <td>3907344</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op3042</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>hutP</td>\n",
       "      <td>BSU_39340</td>\n",
       "      <td>4040992</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op3065</th>\n",
       "      <td>161</td>\n",
       "      <td>1.89e-07</td>\n",
       "      <td>ATATTGACAGAATTATCAGAGGATTAT</td>\n",
       "      <td>yxbD,yxbC</td>\n",
       "      <td>BSU_39870,BSU_39880</td>\n",
       "      <td>4096807</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op3066</th>\n",
       "      <td>52</td>\n",
       "      <td>1.89e-07</td>\n",
       "      <td>ATATTGACAGAATTATCAGAGGATTAT</td>\n",
       "      <td>yxbB,yxbA,yxnB,asnH</td>\n",
       "      <td>BSU_39890,BSU_39900,BSU_39910,BSU_39920</td>\n",
       "      <td>4096916</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op4</th>\n",
       "      <td>306</td>\n",
       "      <td>5.32e-06</td>\n",
       "      <td>CCTTTTCTAGAAATTTCTGGGTTGATT</td>\n",
       "      <td>guaB</td>\n",
       "      <td>BSU_00090</td>\n",
       "      <td>15415</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op459</th>\n",
       "      <td>362</td>\n",
       "      <td>1.77e-08</td>\n",
       "      <td>GATTTTTCTGAATTATCAGATGATATA</td>\n",
       "      <td>gmuF,gmuG</td>\n",
       "      <td>BSU_05870,BSU_05880</td>\n",
       "      <td>631308</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op600</th>\n",
       "      <td>170</td>\n",
       "      <td>4.80e-07</td>\n",
       "      <td>CTTAATTTTGATTTTTCTGAAAATAGG</td>\n",
       "      <td>yfmB</td>\n",
       "      <td>BSU_07530</td>\n",
       "      <td>827111</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op732</th>\n",
       "      <td>440</td>\n",
       "      <td>7.37e-08</td>\n",
       "      <td>ATGAATTCAGAAAATTCAAATAAAATA</td>\n",
       "      <td>bcaP</td>\n",
       "      <td>BSU_09460</td>\n",
       "      <td>1022850</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op750</th>\n",
       "      <td>178</td>\n",
       "      <td>1.19e-06</td>\n",
       "      <td>ATGTTTTTTCACAATTATGAATATAAA</td>\n",
       "      <td>nhaC</td>\n",
       "      <td>BSU_09680</td>\n",
       "      <td>1044146</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op818</th>\n",
       "      <td>157</td>\n",
       "      <td>1.89e-07</td>\n",
       "      <td>AAATTTTTATAATTTTTAGACAACTAA</td>\n",
       "      <td>yhzC</td>\n",
       "      <td>BSU_10410</td>\n",
       "      <td>1116716</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op822</th>\n",
       "      <td>234</td>\n",
       "      <td>2.52e-07</td>\n",
       "      <td>GTATAGTGTGAATTGTCTGAAATTAAT</td>\n",
       "      <td>yhjB,yhjC</td>\n",
       "      <td>BSU_10450,BSU_10460</td>\n",
       "      <td>1120728</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op887</th>\n",
       "      <td>400</td>\n",
       "      <td>1.02e-06</td>\n",
       "      <td>AAATTTTTCGATAATTCAATATTATTT</td>\n",
       "      <td>appD,appF,appAn,appAc,appB,appC</td>\n",
       "      <td>BSU_11360,BSU_11370,BSU_11381,BSU_11382,BSU_11...</td>\n",
       "      <td>1210977</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op890</th>\n",
       "      <td>431</td>\n",
       "      <td>9.17e-09</td>\n",
       "      <td>CATTTGTCTGATATTTCTGAGGATTTA</td>\n",
       "      <td>oppA,oppB,oppC,oppD,oppF</td>\n",
       "      <td>BSU_11430,BSU_11440,BSU_11450,BSU_11460,BSU_11470</td>\n",
       "      <td>1219349</td>\n",
       "      <td>+</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op925</th>\n",
       "      <td>113</td>\n",
       "      <td>3.71e-08</td>\n",
       "      <td>CTATTGTTAGAAAATTAAGACATATTA</td>\n",
       "      <td>yjcL</td>\n",
       "      <td>BSU_11900</td>\n",
       "      <td>1262516</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Op966</th>\n",
       "      <td>168</td>\n",
       "      <td>2.47e-06</td>\n",
       "      <td>AAATTTTTTGAATTTTTGCACGTTAGT</td>\n",
       "      <td>yjnA</td>\n",
       "      <td>BSU_12400</td>\n",
       "      <td>1313515</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       rel_position    pvalue                     site_seq  \\\n",
       "operon                                                       \n",
       "Op1002          163  3.03e-06  CCTATTTTAGAAAAGTGAAATAAAATA   \n",
       "Op1009          455  6.59e-08  AATTTGTTAGAATATTCATAATTTAGT   \n",
       "Op1010          NaN       NaN                          NaN   \n",
       "Op1028          NaN       NaN                          NaN   \n",
       "Op1094          126  5.89e-08  ACGAAGATTGAATTTTCAGTAAATATA   \n",
       "Op1196          446  6.77e-07  AAGTTTTATGAAAATTAAAACAATTAT   \n",
       "Op1296          427  3.29e-08  TAAAAGTTTGAATATTATGAAAATATA   \n",
       "Op1505          406  1.50e-06  AATTATTCAGAAAATTATTTATAATAA   \n",
       "Op1800          NaN       NaN                          NaN   \n",
       "Op181           451  3.05e-07  AAAAAGACTGAATATTTAAACAATTAT   \n",
       "Op226           191  1.13e-07  ACATATTGCTACTATTCAGAATAATTA   \n",
       "Op227           441  1.13e-07  ACATATTGCTACTATTCAGAATAATTA   \n",
       "Op2347          454  9.40e-07  TAATTGACAGAATATTTTAACATAATT   \n",
       "Op243           NaN       NaN                          NaN   \n",
       "Op2486          155  1.29e-06  AAAAATAGATAATTTTTAAAAAATTAT   \n",
       "Op2519          NaN       NaN                          NaN   \n",
       "Op2520          308  4.66e-10  CCTTTTTCTTAATTTTCTGATAATATA   \n",
       "Op2521          111  4.66e-10  CCTTTTTCTTAATTTTCTGATAATATA   \n",
       "Op2553          442  2.47e-06  ACTAATTTACAATTTTAATGAAAAATA   \n",
       "Op2556          185  9.40e-07  CATTTTTGCTATATTTCTAAAATAAAA   \n",
       "Op2558          184  6.59e-08  TTATTTTCTGAAATTTCTCAATAATTT   \n",
       "Op2935          NaN       NaN                          NaN   \n",
       "Op2952          427  2.89e-09  GATTTTTTTGAATTTTCAAGAAATATA   \n",
       "Op3042          NaN       NaN                          NaN   \n",
       "Op3065          161  1.89e-07  ATATTGACAGAATTATCAGAGGATTAT   \n",
       "Op3066           52  1.89e-07  ATATTGACAGAATTATCAGAGGATTAT   \n",
       "Op4             306  5.32e-06  CCTTTTCTAGAAATTTCTGGGTTGATT   \n",
       "Op459           362  1.77e-08  GATTTTTCTGAATTATCAGATGATATA   \n",
       "Op600           170  4.80e-07  CTTAATTTTGATTTTTCTGAAAATAGG   \n",
       "Op732           440  7.37e-08  ATGAATTCAGAAAATTCAAATAAAATA   \n",
       "Op750           178  1.19e-06  ATGTTTTTTCACAATTATGAATATAAA   \n",
       "Op818           157  1.89e-07  AAATTTTTATAATTTTTAGACAACTAA   \n",
       "Op822           234  2.52e-07  GTATAGTGTGAATTGTCTGAAATTAAT   \n",
       "Op887           400  1.02e-06  AAATTTTTCGATAATTCAATATTATTT   \n",
       "Op890           431  9.17e-09  CATTTGTCTGATATTTCTGAGGATTTA   \n",
       "Op925           113  3.71e-08  CTATTGTTAGAAAATTAAGACATATTA   \n",
       "Op966           168  2.47e-06  AAATTTTTTGAATTTTTGCACGTTAGT   \n",
       "\n",
       "                                  genes  \\\n",
       "operon                                    \n",
       "Op1002                             steT   \n",
       "Op1009         dppA,dppB,dppC,dppD,dppE   \n",
       "Op1010              ldcA,aeeB,eepC,ykfD   \n",
       "Op1028                             ispA   \n",
       "Op1094                             ykwB   \n",
       "Op1196                             ylmA   \n",
       "Op1296              acpK,pksF,pksG,pksH   \n",
       "Op1505                             yocS   \n",
       "Op1800                             yppF   \n",
       "Op181                              ilvE   \n",
       "Op226                              amhX   \n",
       "Op227                              ycgA   \n",
       "Op2347                             braB   \n",
       "Op243                    putB,putC,putP   \n",
       "Op2486              nupN,nupO,nupP,nupQ   \n",
       "Op2519                             spsC   \n",
       "Op2520                             yuiB   \n",
       "Op2521                             yuiA   \n",
       "Op2553                        pucG,pucF   \n",
       "Op2556                             frlP   \n",
       "Op2558         frlD,frlM,frlN,frlO,frlB   \n",
       "Op2935                        rocC,rocB   \n",
       "Op2952                              vpr   \n",
       "Op3042                             hutP   \n",
       "Op3065                        yxbD,yxbC   \n",
       "Op3066              yxbB,yxbA,yxnB,asnH   \n",
       "Op4                                guaB   \n",
       "Op459                         gmuF,gmuG   \n",
       "Op600                              yfmB   \n",
       "Op732                              bcaP   \n",
       "Op750                              nhaC   \n",
       "Op818                              yhzC   \n",
       "Op822                         yhjB,yhjC   \n",
       "Op887   appD,appF,appAn,appAc,appB,appC   \n",
       "Op890          oppA,oppB,oppC,oppD,oppF   \n",
       "Op925                              yjcL   \n",
       "Op966                              yjnA   \n",
       "\n",
       "                                               locus_tags  start_pos strand  \n",
       "operon                                                                       \n",
       "Op1002                                          BSU_12860    1352591      -  \n",
       "Op1009  BSU_12920,BSU_12930,BSU_12940,BSU_12950,BSU_12960    1359901      +  \n",
       "Op1010            BSU_12970,BSU_12980,BSU_12990,BSU_13000    1365388      +  \n",
       "Op1028                                          BSU_13190    1386883      -  \n",
       "Op1094                                          BSU_13940    1463393      -  \n",
       "Op1196                                          BSU_15340    1606060      +  \n",
       "Op1296            BSU_17130,BSU_17140,BSU_17150,BSU_17160    1787969      +  \n",
       "Op1505                                          BSU_19350    2105990      +  \n",
       "Op1800                                          BSU_22260    2338082      +  \n",
       "Op181                                           BSU_02390     258516      +  \n",
       "Op226                                           BSU_03010     325089      -  \n",
       "Op227                                           BSU_03020     324839      +  \n",
       "Op2347                                          BSU_29600    3027797      +  \n",
       "Op243                       BSU_03200,BSU_03210,BSU_03220     344051      +  \n",
       "Op2486            BSU_31540,BSU_31550,BSU_31560,BSU_31570    3239430      +  \n",
       "Op2519                                          BSU_32070    3299140      -  \n",
       "Op2520                                          BSU_32080    3299564      -  \n",
       "Op2521                                          BSU_32090    3299761      -  \n",
       "Op2553                                BSU_32520,BSU_32530    3343571      -  \n",
       "Op2556                                          BSU_32550    3346016      -  \n",
       "Op2558  BSU_32570,BSU_32580,BSU_32590,BSU_32600,BSU_32610    3351996      -  \n",
       "Op2935                                BSU_37760,BSU_37770    3878792      -  \n",
       "Op2952                                          BSU_38090    3907344      +  \n",
       "Op3042                                          BSU_39340    4040992      +  \n",
       "Op3065                                BSU_39870,BSU_39880    4096807      -  \n",
       "Op3066            BSU_39890,BSU_39900,BSU_39910,BSU_39920    4096916      +  \n",
       "Op4                                             BSU_00090      15415      +  \n",
       "Op459                                 BSU_05870,BSU_05880     631308      +  \n",
       "Op600                                           BSU_07530     827111      -  \n",
       "Op732                                           BSU_09460    1022850      +  \n",
       "Op750                                           BSU_09680    1044146      -  \n",
       "Op818                                           BSU_10410    1116716      -  \n",
       "Op822                                 BSU_10450,BSU_10460    1120728      -  \n",
       "Op887   BSU_11360,BSU_11370,BSU_11381,BSU_11382,BSU_11...    1210977      +  \n",
       "Op890   BSU_11430,BSU_11440,BSU_11450,BSU_11460,BSU_11470    1219349      +  \n",
       "Op925                                           BSU_11900    1262516      -  \n",
       "Op966                                           BSU_12400    1313515      -  "
      ]
     },
     "execution_count": 63,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ica_data.motif_info['CodY'].sites.loc['MEME-1']"
   ]
  },
  {
   "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.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}