{
 "cells": [
  {
   "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 tr th {\n",
       "        text-align: left;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th>Category</th>\n",
       "      <th colspan=\"4\" halign=\"left\">Count</th>\n",
       "      <th colspan=\"4\" halign=\"left\">GC content</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Enzyme</th>\n",
       "      <th>ImProm-II 42C</th>\n",
       "      <th>ImProm-II 55C</th>\n",
       "      <th>SuperScriptIV</th>\n",
       "      <th>TGIRT</th>\n",
       "      <th>ImProm-II 42C</th>\n",
       "      <th>ImProm-II 55C</th>\n",
       "      <th>SuperScriptIV</th>\n",
       "      <th>TGIRT</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>Low abundance</th>\n",
       "      <td>288</td>\n",
       "      <td>239</td>\n",
       "      <td>230</td>\n",
       "      <td>211</td>\n",
       "      <td>0.500000</td>\n",
       "      <td>0.500000</td>\n",
       "      <td>0.500000</td>\n",
       "      <td>0.500000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Verrucomicrobiae</th>\n",
       "      <td>1310</td>\n",
       "      <td>620</td>\n",
       "      <td>1391</td>\n",
       "      <td>1029</td>\n",
       "      <td>0.534486</td>\n",
       "      <td>0.533333</td>\n",
       "      <td>0.532315</td>\n",
       "      <td>0.532096</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Unclassified</th>\n",
       "      <td>1211</td>\n",
       "      <td>1210</td>\n",
       "      <td>1292</td>\n",
       "      <td>1216</td>\n",
       "      <td>0.492766</td>\n",
       "      <td>0.507756</td>\n",
       "      <td>0.505267</td>\n",
       "      <td>0.503933</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Alphaproteobacteria</th>\n",
       "      <td>12943</td>\n",
       "      <td>9302</td>\n",
       "      <td>12148</td>\n",
       "      <td>12385</td>\n",
       "      <td>0.554835</td>\n",
       "      <td>0.553426</td>\n",
       "      <td>0.562606</td>\n",
       "      <td>0.563610</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Deltaproteobacteria</th>\n",
       "      <td>4330</td>\n",
       "      <td>2154</td>\n",
       "      <td>3141</td>\n",
       "      <td>3001</td>\n",
       "      <td>0.563930</td>\n",
       "      <td>0.566020</td>\n",
       "      <td>0.565218</td>\n",
       "      <td>0.565658</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Gammaproteobacteria</th>\n",
       "      <td>5967</td>\n",
       "      <td>7610</td>\n",
       "      <td>2793</td>\n",
       "      <td>2576</td>\n",
       "      <td>0.543056</td>\n",
       "      <td>0.536829</td>\n",
       "      <td>0.544228</td>\n",
       "      <td>0.545745</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Planctomycetacia</th>\n",
       "      <td>357</td>\n",
       "      <td>710</td>\n",
       "      <td>247</td>\n",
       "      <td>283</td>\n",
       "      <td>0.598206</td>\n",
       "      <td>0.624758</td>\n",
       "      <td>0.607147</td>\n",
       "      <td>0.619704</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Phycisphaerae</th>\n",
       "      <td>101</td>\n",
       "      <td>128</td>\n",
       "      <td>18</td>\n",
       "      <td>17</td>\n",
       "      <td>0.576535</td>\n",
       "      <td>0.573142</td>\n",
       "      <td>0.572757</td>\n",
       "      <td>0.577386</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Saccharimonadia</th>\n",
       "      <td>130</td>\n",
       "      <td>57</td>\n",
       "      <td>49</td>\n",
       "      <td>30</td>\n",
       "      <td>0.503040</td>\n",
       "      <td>0.501756</td>\n",
       "      <td>0.496003</td>\n",
       "      <td>0.496590</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Longimicrobia</th>\n",
       "      <td>328</td>\n",
       "      <td>217</td>\n",
       "      <td>213</td>\n",
       "      <td>197</td>\n",
       "      <td>0.615191</td>\n",
       "      <td>0.616098</td>\n",
       "      <td>0.615348</td>\n",
       "      <td>0.614556</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Gemmatimonadetes</th>\n",
       "      <td>532</td>\n",
       "      <td>591</td>\n",
       "      <td>518</td>\n",
       "      <td>463</td>\n",
       "      <td>0.626024</td>\n",
       "      <td>0.629010</td>\n",
       "      <td>0.626045</td>\n",
       "      <td>0.628619</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>BD2-11</th>\n",
       "      <td>26</td>\n",
       "      <td>33</td>\n",
       "      <td>48</td>\n",
       "      <td>42</td>\n",
       "      <td>0.609861</td>\n",
       "      <td>0.610253</td>\n",
       "      <td>0.608890</td>\n",
       "      <td>0.610935</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>S0134</th>\n",
       "      <td>116</td>\n",
       "      <td>75</td>\n",
       "      <td>73</td>\n",
       "      <td>114</td>\n",
       "      <td>0.589186</td>\n",
       "      <td>0.603075</td>\n",
       "      <td>0.599525</td>\n",
       "      <td>0.598523</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Nitrospira</th>\n",
       "      <td>32</td>\n",
       "      <td>34</td>\n",
       "      <td>74</td>\n",
       "      <td>70</td>\n",
       "      <td>0.600736</td>\n",
       "      <td>0.599553</td>\n",
       "      <td>0.599522</td>\n",
       "      <td>0.599818</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Bacilli</th>\n",
       "      <td>1356</td>\n",
       "      <td>656</td>\n",
       "      <td>663</td>\n",
       "      <td>520</td>\n",
       "      <td>0.522874</td>\n",
       "      <td>0.533793</td>\n",
       "      <td>0.531269</td>\n",
       "      <td>0.537584</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Fibrobacteria</th>\n",
       "      <td>45</td>\n",
       "      <td>35</td>\n",
       "      <td>46</td>\n",
       "      <td>36</td>\n",
       "      <td>0.554449</td>\n",
       "      <td>0.554841</td>\n",
       "      <td>0.553832</td>\n",
       "      <td>0.553822</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Entotheonellia</th>\n",
       "      <td>125</td>\n",
       "      <td>122</td>\n",
       "      <td>221</td>\n",
       "      <td>164</td>\n",
       "      <td>0.576794</td>\n",
       "      <td>0.576664</td>\n",
       "      <td>0.576262</td>\n",
       "      <td>0.575045</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Deinococci</th>\n",
       "      <td>89</td>\n",
       "      <td>28</td>\n",
       "      <td>44</td>\n",
       "      <td>29</td>\n",
       "      <td>0.565357</td>\n",
       "      <td>0.574220</td>\n",
       "      <td>0.561380</td>\n",
       "      <td>0.564073</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Oxyphotobacteria</th>\n",
       "      <td>2798</td>\n",
       "      <td>962</td>\n",
       "      <td>1711</td>\n",
       "      <td>877</td>\n",
       "      <td>0.546075</td>\n",
       "      <td>0.544738</td>\n",
       "      <td>0.546060</td>\n",
       "      <td>0.545485</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Chloroflexia</th>\n",
       "      <td>441</td>\n",
       "      <td>825</td>\n",
       "      <td>666</td>\n",
       "      <td>789</td>\n",
       "      <td>0.639146</td>\n",
       "      <td>0.644771</td>\n",
       "      <td>0.639664</td>\n",
       "      <td>0.640914</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Anaerolineae</th>\n",
       "      <td>191</td>\n",
       "      <td>175</td>\n",
       "      <td>113</td>\n",
       "      <td>89</td>\n",
       "      <td>0.552375</td>\n",
       "      <td>0.559423</td>\n",
       "      <td>0.561238</td>\n",
       "      <td>0.564171</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Gitt-GS-136</th>\n",
       "      <td>270</td>\n",
       "      <td>171</td>\n",
       "      <td>239</td>\n",
       "      <td>274</td>\n",
       "      <td>0.560948</td>\n",
       "      <td>0.561890</td>\n",
       "      <td>0.561070</td>\n",
       "      <td>0.560207</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>JG30-KF-CM66</th>\n",
       "      <td>76</td>\n",
       "      <td>145</td>\n",
       "      <td>175</td>\n",
       "      <td>175</td>\n",
       "      <td>0.596772</td>\n",
       "      <td>0.597473</td>\n",
       "      <td>0.597070</td>\n",
       "      <td>0.597318</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>KD4-96</th>\n",
       "      <td>298</td>\n",
       "      <td>264</td>\n",
       "      <td>370</td>\n",
       "      <td>406</td>\n",
       "      <td>0.569651</td>\n",
       "      <td>0.569340</td>\n",
       "      <td>0.569242</td>\n",
       "      <td>0.568869</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>TK10</th>\n",
       "      <td>275</td>\n",
       "      <td>583</td>\n",
       "      <td>522</td>\n",
       "      <td>611</td>\n",
       "      <td>0.614815</td>\n",
       "      <td>0.628671</td>\n",
       "      <td>0.618417</td>\n",
       "      <td>0.616980</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Ktedonobacteria</th>\n",
       "      <td>26</td>\n",
       "      <td>41</td>\n",
       "      <td>61</td>\n",
       "      <td>52</td>\n",
       "      <td>0.578646</td>\n",
       "      <td>0.578946</td>\n",
       "      <td>0.580507</td>\n",
       "      <td>0.580245</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Dehalococcoidia</th>\n",
       "      <td>70</td>\n",
       "      <td>98</td>\n",
       "      <td>155</td>\n",
       "      <td>97</td>\n",
       "      <td>0.578148</td>\n",
       "      <td>0.577860</td>\n",
       "      <td>0.577543</td>\n",
       "      <td>0.578580</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Bacteroidia</th>\n",
       "      <td>2825</td>\n",
       "      <td>1817</td>\n",
       "      <td>2307</td>\n",
       "      <td>3005</td>\n",
       "      <td>0.514352</td>\n",
       "      <td>0.514383</td>\n",
       "      <td>0.513738</td>\n",
       "      <td>0.508924</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Fimbriimonadia</th>\n",
       "      <td>69</td>\n",
       "      <td>29</td>\n",
       "      <td>32</td>\n",
       "      <td>43</td>\n",
       "      <td>0.542453</td>\n",
       "      <td>0.543326</td>\n",
       "      <td>0.536859</td>\n",
       "      <td>0.541226</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>uncultured</th>\n",
       "      <td>41</td>\n",
       "      <td>66</td>\n",
       "      <td>31</td>\n",
       "      <td>28</td>\n",
       "      <td>0.608757</td>\n",
       "      <td>0.612651</td>\n",
       "      <td>0.611116</td>\n",
       "      <td>0.616270</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Actinobacteria</th>\n",
       "      <td>5864</td>\n",
       "      <td>11263</td>\n",
       "      <td>10105</td>\n",
       "      <td>11353</td>\n",
       "      <td>0.583321</td>\n",
       "      <td>0.583805</td>\n",
       "      <td>0.584725</td>\n",
       "      <td>0.586320</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Acidimicrobiia</th>\n",
       "      <td>1203</td>\n",
       "      <td>1617</td>\n",
       "      <td>1518</td>\n",
       "      <td>1858</td>\n",
       "      <td>0.585221</td>\n",
       "      <td>0.586662</td>\n",
       "      <td>0.588191</td>\n",
       "      <td>0.587488</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Thermoleophilia</th>\n",
       "      <td>2745</td>\n",
       "      <td>4643</td>\n",
       "      <td>3733</td>\n",
       "      <td>3562</td>\n",
       "      <td>0.589096</td>\n",
       "      <td>0.588953</td>\n",
       "      <td>0.589492</td>\n",
       "      <td>0.589915</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0319-7L14</th>\n",
       "      <td>527</td>\n",
       "      <td>442</td>\n",
       "      <td>332</td>\n",
       "      <td>392</td>\n",
       "      <td>0.580750</td>\n",
       "      <td>0.582987</td>\n",
       "      <td>0.581233</td>\n",
       "      <td>0.582942</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Nitriliruptoria</th>\n",
       "      <td>253</td>\n",
       "      <td>380</td>\n",
       "      <td>410</td>\n",
       "      <td>408</td>\n",
       "      <td>0.576569</td>\n",
       "      <td>0.578654</td>\n",
       "      <td>0.576860</td>\n",
       "      <td>0.577517</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Rubrobacteria</th>\n",
       "      <td>1351</td>\n",
       "      <td>1505</td>\n",
       "      <td>2380</td>\n",
       "      <td>1974</td>\n",
       "      <td>0.597608</td>\n",
       "      <td>0.600372</td>\n",
       "      <td>0.601017</td>\n",
       "      <td>0.601583</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>MB-A2-108</th>\n",
       "      <td>61</td>\n",
       "      <td>101</td>\n",
       "      <td>67</td>\n",
       "      <td>99</td>\n",
       "      <td>0.608917</td>\n",
       "      <td>0.608527</td>\n",
       "      <td>0.608236</td>\n",
       "      <td>0.605956</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Thermoanaerobaculia</th>\n",
       "      <td>117</td>\n",
       "      <td>127</td>\n",
       "      <td>122</td>\n",
       "      <td>146</td>\n",
       "      <td>0.582242</td>\n",
       "      <td>0.582820</td>\n",
       "      <td>0.583003</td>\n",
       "      <td>0.582521</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subgroup</th>\n",
       "      <td>385</td>\n",
       "      <td>284</td>\n",
       "      <td>681</td>\n",
       "      <td>501</td>\n",
       "      <td>0.580126</td>\n",
       "      <td>0.580621</td>\n",
       "      <td>0.579717</td>\n",
       "      <td>0.580671</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Holophagae</th>\n",
       "      <td>235</td>\n",
       "      <td>139</td>\n",
       "      <td>308</td>\n",
       "      <td>191</td>\n",
       "      <td>0.567214</td>\n",
       "      <td>0.566678</td>\n",
       "      <td>0.568034</td>\n",
       "      <td>0.568013</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Blastocatellia</th>\n",
       "      <td>336</td>\n",
       "      <td>243</td>\n",
       "      <td>309</td>\n",
       "      <td>328</td>\n",
       "      <td>0.543584</td>\n",
       "      <td>0.546420</td>\n",
       "      <td>0.548164</td>\n",
       "      <td>0.549827</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Acidobacteriia</th>\n",
       "      <td>257</td>\n",
       "      <td>259</td>\n",
       "      <td>444</td>\n",
       "      <td>359</td>\n",
       "      <td>0.555791</td>\n",
       "      <td>0.561708</td>\n",
       "      <td>0.561813</td>\n",
       "      <td>0.559804</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "Category                    Count                                     \\\n",
       "Enzyme              ImProm-II 42C ImProm-II 55C SuperScriptIV  TGIRT   \n",
       "Low abundance                 288           239           230    211   \n",
       "Verrucomicrobiae             1310           620          1391   1029   \n",
       "Unclassified                 1211          1210          1292   1216   \n",
       "Alphaproteobacteria         12943          9302         12148  12385   \n",
       "Deltaproteobacteria          4330          2154          3141   3001   \n",
       "Gammaproteobacteria          5967          7610          2793   2576   \n",
       "Planctomycetacia              357           710           247    283   \n",
       "Phycisphaerae                 101           128            18     17   \n",
       "Saccharimonadia               130            57            49     30   \n",
       "Longimicrobia                 328           217           213    197   \n",
       "Gemmatimonadetes              532           591           518    463   \n",
       "BD2-11                         26            33            48     42   \n",
       "S0134                         116            75            73    114   \n",
       "Nitrospira                     32            34            74     70   \n",
       "Bacilli                      1356           656           663    520   \n",
       "Fibrobacteria                  45            35            46     36   \n",
       "Entotheonellia                125           122           221    164   \n",
       "Deinococci                     89            28            44     29   \n",
       "Oxyphotobacteria             2798           962          1711    877   \n",
       "Chloroflexia                  441           825           666    789   \n",
       "Anaerolineae                  191           175           113     89   \n",
       "Gitt-GS-136                   270           171           239    274   \n",
       "JG30-KF-CM66                   76           145           175    175   \n",
       "KD4-96                        298           264           370    406   \n",
       "TK10                          275           583           522    611   \n",
       "Ktedonobacteria                26            41            61     52   \n",
       "Dehalococcoidia                70            98           155     97   \n",
       "Bacteroidia                  2825          1817          2307   3005   \n",
       "Fimbriimonadia                 69            29            32     43   \n",
       "uncultured                     41            66            31     28   \n",
       "Actinobacteria               5864         11263         10105  11353   \n",
       "Acidimicrobiia               1203          1617          1518   1858   \n",
       "Thermoleophilia              2745          4643          3733   3562   \n",
       "0319-7L14                     527           442           332    392   \n",
       "Nitriliruptoria               253           380           410    408   \n",
       "Rubrobacteria                1351          1505          2380   1974   \n",
       "MB-A2-108                      61           101            67     99   \n",
       "Thermoanaerobaculia           117           127           122    146   \n",
       "Subgroup                      385           284           681    501   \n",
       "Holophagae                    235           139           308    191   \n",
       "Blastocatellia                336           243           309    328   \n",
       "Acidobacteriia                257           259           444    359   \n",
       "\n",
       "Category               GC content                                        \n",
       "Enzyme              ImProm-II 42C ImProm-II 55C SuperScriptIV     TGIRT  \n",
       "Low abundance            0.500000      0.500000      0.500000  0.500000  \n",
       "Verrucomicrobiae         0.534486      0.533333      0.532315  0.532096  \n",
       "Unclassified             0.492766      0.507756      0.505267  0.503933  \n",
       "Alphaproteobacteria      0.554835      0.553426      0.562606  0.563610  \n",
       "Deltaproteobacteria      0.563930      0.566020      0.565218  0.565658  \n",
       "Gammaproteobacteria      0.543056      0.536829      0.544228  0.545745  \n",
       "Planctomycetacia         0.598206      0.624758      0.607147  0.619704  \n",
       "Phycisphaerae            0.576535      0.573142      0.572757  0.577386  \n",
       "Saccharimonadia          0.503040      0.501756      0.496003  0.496590  \n",
       "Longimicrobia            0.615191      0.616098      0.615348  0.614556  \n",
       "Gemmatimonadetes         0.626024      0.629010      0.626045  0.628619  \n",
       "BD2-11                   0.609861      0.610253      0.608890  0.610935  \n",
       "S0134                    0.589186      0.603075      0.599525  0.598523  \n",
       "Nitrospira               0.600736      0.599553      0.599522  0.599818  \n",
       "Bacilli                  0.522874      0.533793      0.531269  0.537584  \n",
       "Fibrobacteria            0.554449      0.554841      0.553832  0.553822  \n",
       "Entotheonellia           0.576794      0.576664      0.576262  0.575045  \n",
       "Deinococci               0.565357      0.574220      0.561380  0.564073  \n",
       "Oxyphotobacteria         0.546075      0.544738      0.546060  0.545485  \n",
       "Chloroflexia             0.639146      0.644771      0.639664  0.640914  \n",
       "Anaerolineae             0.552375      0.559423      0.561238  0.564171  \n",
       "Gitt-GS-136              0.560948      0.561890      0.561070  0.560207  \n",
       "JG30-KF-CM66             0.596772      0.597473      0.597070  0.597318  \n",
       "KD4-96                   0.569651      0.569340      0.569242  0.568869  \n",
       "TK10                     0.614815      0.628671      0.618417  0.616980  \n",
       "Ktedonobacteria          0.578646      0.578946      0.580507  0.580245  \n",
       "Dehalococcoidia          0.578148      0.577860      0.577543  0.578580  \n",
       "Bacteroidia              0.514352      0.514383      0.513738  0.508924  \n",
       "Fimbriimonadia           0.542453      0.543326      0.536859  0.541226  \n",
       "uncultured               0.608757      0.612651      0.611116  0.616270  \n",
       "Actinobacteria           0.583321      0.583805      0.584725  0.586320  \n",
       "Acidimicrobiia           0.585221      0.586662      0.588191  0.587488  \n",
       "Thermoleophilia          0.589096      0.588953      0.589492  0.589915  \n",
       "0319-7L14                0.580750      0.582987      0.581233  0.582942  \n",
       "Nitriliruptoria          0.576569      0.578654      0.576860  0.577517  \n",
       "Rubrobacteria            0.597608      0.600372      0.601017  0.601583  \n",
       "MB-A2-108                0.608917      0.608527      0.608236  0.605956  \n",
       "Thermoanaerobaculia      0.582242      0.582820      0.583003  0.582521  \n",
       "Subgroup                 0.580126      0.580621      0.579717  0.580671  \n",
       "Holophagae               0.567214      0.566678      0.568034  0.568013  \n",
       "Blastocatellia           0.543584      0.546420      0.548164  0.549827  \n",
       "Acidobacteriia           0.555791      0.561708      0.561813  0.559804  "
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "\n",
    "df1 = pd.read_csv(\"GC_class_forFigure_ordered.csv\", header=[0,1], index_col=0)\n",
    "df1.columns.names = [\"Category\", \"Enzyme\"]\n",
    "df1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "|Enzyme        |TGIRT  |SuperScriptIV|Promega42  |Promega55|\n",
    "|:-------------|:-----:|:-----------:|:---------:|:-------:|\n",
    "|TGIRT         | x     | TvS         | Tv42      | Tv55    |  \n",
    "|SuperScriptIV | TvS   | x           | Sv42      | Sv55    |\n",
    "|Promega42     | Tv42  | Sv42        | x         | P42v55  |\n",
    "|Promega55     | Tv55  | Sv55        | P42v55    | x       |"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "def relativeEnrichment(a, b):\n",
    "    return (a-b)/(a+b)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "idx = pd.IndexSlice\n",
    "\n",
    "TvS = relativeEnrichment(df1.loc[:, idx[\"Count\", \"TGIRT\"]], df1.loc[:, idx[\"Count\", \"SuperScriptIV\"]])\n",
    "df1[\"Comparison\", \"TvS\"] = TvS\n",
    "\n",
    "Tv42 = relativeEnrichment(df1.loc[:, idx[\"Count\", \"TGIRT\"]], df1.loc[:, idx[\"Count\", \"ImProm-II 42C\"]])\n",
    "df1[\"Comparison\", \"Tv42\"] = Tv42\n",
    "\n",
    "Tv55 = relativeEnrichment(df1.loc[:, idx[\"Count\", \"TGIRT\"]], df1.loc[:, idx[\"Count\", \"ImProm-II 55C\"]])\n",
    "df1[\"Comparison\", \"Tv55\"] = Tv55\n",
    "\n",
    "Sv42 = relativeEnrichment(df1.loc[:, idx[\"Count\", \"SuperScriptIV\"]], df1.loc[:, idx[\"Count\", \"ImProm-II 42C\"]])\n",
    "df1[\"Comparison\", \"Sv42\"] = Sv42\n",
    "\n",
    "Sv55 = relativeEnrichment(df1.loc[:, idx[\"Count\", \"SuperScriptIV\"]], df1.loc[:, idx[\"Count\", \"ImProm-II 55C\"]])\n",
    "df1[\"Comparison\", \"Sv55\"] = Sv55\n",
    "\n",
    "P42v55 = relativeEnrichment(df1.loc[:, idx[\"Count\", \"ImProm-II 42C\"]], df1.loc[:, idx[\"Count\", \"ImProm-II 55C\"]])\n",
    "df1[\"Comparison\", \"P42v55\"] = P42v55"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "from matplotlib import cm\n",
    "import numpy as np\n",
    "\n",
    "default_dpi = plt.rcParamsDefault[\"figure.dpi\"]\n",
    "plt.rcParams[\"figure.figsize\"] = [6,8]\n",
    "\n",
    "def barplotFigure(y, GC_content, name):\n",
    "    # Bar position\n",
    "    x_pos = np.arange(len(y))\n",
    "    # \"Stretching\" color values betwenn 0 and 1 to get the best colourfulness\n",
    "    col_norm =(GC_content-min(GC_content))/(max(GC_content-min(GC_content)))\n",
    "    # Getting color pallette\n",
    "    colors = cm.viridis(col_norm)\n",
    "    \n",
    "    # Legend construction\n",
    "    plot = plt.scatter(GC_content, GC_content, c=GC_content, cmap=\"viridis\")\n",
    "    plt.clf()\n",
    "    plt.colorbar(plot)\n",
    "\n",
    "    # Barplot \n",
    "    plt.barh(x_pos, y, color=colors, align=\"center\")\n",
    "    plt.yticks(x_pos, df1.index.values)\n",
    "    plt.title(name)\n",
    "    \n",
    "    axes = plt.gca()\n",
    "    axes.set_xlim([-0.8,0.8])\n",
    "    \n",
    "    plt.savefig(name+\".pdf\", bbox_inches=\"tight\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdUAAAHiCAYAAAC+xlbiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsnXeYnGXZvs8rCU1CEYgISAyGXgNZejF0LEhQisgnRJGInygW4KeCEECqiCgREPgwNAGphiK9S0tCSEIHIQiCkFCEQEBIrt8fzzPZN5OZnS2zu7O793kc7zHvPPV+Z5K952nXLdsEQRAEQdBx+nW3AUEQBEHQWwinGgRBEAR1IpxqEARBENSJcKpBEARBUCfCqQZBEARBnQinGgRBEAR1IpxqEARBENSJcKpBwyFpuqQd2ll3nKT/Spol6U1Jt0pas942ttKWIZIsaUDBtl/VKDurcP2ykF98rtLVv5B/tKR/S7pC0iKF9CUlnS7pn7nOc/n9ch19nnbUny5pdrbjNUl/kjSwPW11FEmjJN1XZlvNf3OS7ij7Tj8l6VJJr0j6j6S/S9q0UF6S/pif9w9lba0g6f8kvSrpXUlPSTpG0uL1fNagawmnGvRGTrE9EPgM8DowrlKh9jqHTmZp2wPzdVxZ3imFvIG25wBIWg3YBPgsMBn4n5y+MHA7sA6wC7AksAXwRi7fHeyav5uNgI2BI8sLZEfUcH+bJO0LlP+bGQhMAIYDywAXADcUfiyUHPVKwABJ2+e2lgEeABYDNre9BLAjsDQwtDOfI+hcGu4fbhAUySOKv0v6raS3JT0vaYuc/pKk1yXtX6mu7feBPwPr5rbGSLpS0sWS3gFGSVokj9xeydfppZGepBGSXpZ0eO7nVUkjJX1R0jN5JPyLLvswqtMfUH4tXQD7AYOB3W0/YXuu7ddtH2f7xo52mkfPZ0r6Wx59/l3Sp/Nn+FYeeW1Yqa7tfwF/o/m7uUvS8ZL+DrwPfE7SipLG58/5OUkHFvoek0flF+dR3jRJq0v6ef6uXpK0U0efsdDfUsDRwOFlz/G87dNsv2p7ju1zgIWBNXKR/qS/s8VXgJ8A7wL/Y3t6busl24fYnlovu4OuJ5xq0BPYFJgKLEtykpeRRjmrkkZlYytNI+a0fUmjtxK7AVeSRgSXAEcAmwHDgA1II7ji6OnTwKKkkcZRwLm5z+HA1sBRkj5Xp+cEeDE78j9VmKL93+xgJkn6WinR9lOkz+fFbNdFOWsH4Cbbs1rbuaTrJf2sDfbuRfq8lgM+JI2+HsnvrwROq9LPysAXmf+7+SYwGlgiP8ulwMvAisAewAmlkV5mV9KzfjK3czPpb9pKwLHAH9vwHLU4ATgL+HdLhSQNIznV53LSLfn9y/n9rfl1B+Bq23PraGPQCNiOK66GuoDpwA75fhTwbCFvPcDA8oW0N4Bh+X4c8AHwNukP4HhgaM4bA9xT1tc/gC8W3u8MTM/3I4DZQP/8fonc96aF8pOAkVWeY0guP6Bg26+qlB0INJGmF5cnOaSbC/kbkX5UDCA5o3eBLWt8jrcCJ9Xxe6n0POcW8n8APFn2Xb1d9r3Oyt/Ni8CZwGI57y7g2ELZlYE5wBKFtBOBcYXv8tZC3q657fLvaukqzzIKuK/Sv7kKZZuAR/NnP99nUFZuSWAa8PNWfJbPAgd19/+1uOp/NeKaUhCU81rhfjaA7fK04kj1VNsLrNVlXip7vyLpD3yJF3NaiTec1y5LfVewZyCApOKIcO0q/VfEaTQ5sdS+pIOBVyUtafsd248Uit8o6RLgq8DfW2j2DWCFttjRDso/i5a+F0g/QG6r0lbxu1kReNP2u4W0F0kOrlrfMyt8VwMlrUeaagZ40fY6VfpfgLy2eyZwiO2PJVUrtxhwHfCg7RNb0XRXfDdBNxDTv0Ffozws0yukDT4lBue0tjc8/yaif7bXwFJz+bXyX/GUXy2vxG3Azj1oN2nxu3kFWEbSEoW0wcC/2tyofW/he2m1Q80sSXLkl0v6N2lTEsDLkrYGyGvw12bbvtvKdm8Ddm/EDVlBx4gvNOjrXAocKWlQXsM8Cri4q42QtKmkNST1k7Qs8HvgLtv/yfl7SBqY83cireuOr9HsRaTR31WS1iy1LekXkr7YqQ/UQWy/BNwPnChpUUnrAweQ1sG7kv+QRs3D8lX63IYDD0laiDRVPxvYz61fIz2N5LAvkPRZAEkrSTotP2vQQwmnGvR1fkWadp1KWg97JKd1NZ8DbiKtlT5G2vSzTyH/ENJI6G3g18CBtu9qqUHbH5I2xDxFWl99B3iYtInooUp18k7eRtjRDOn5h5BGrdcAR9u+tcUadcaJf5cuYEbOes32f0lHlL4M7AS8reYzxFvXaPfNXPcjknN+l3T86T80b3IKeiCyI0h5EARBENSDGKkGQRAEQZ0IpxoEQRAEdSKcahAEQRDUiXCqQRAEQVAnwqkGQRAEQZ0IRaU+zHLLLechQ4Z0txlBEPQAJk2aNNP2oPbW33nbxf3Gm3NqF2wjk6Z+eLPtXerecDsJp9qHGTJkCBMnTqxdMAiCPo+kF2uXqs4bb87h4ZsH18ucefRf4dk2xwbuTMKpBkEQBJ2Ogbn0/qA8saYaBEEQBHUiRqpBEARBF2Dm9IHwsTFSDYIgCII6ESPVIAiCoNNJa6q9X2s+nGoQBEHQJcRGpSAIgiAIWk2MVIMgCIJOx5g5fSDUaIxUgyAIgqBOxEg1CIIg6BJio1IQBEEQ1AEDc/qAU43p3yAIgiCoEzFSDYIgCLqEvjD9GyPVIAiCIKgTMVINgiAIOh1DnzhSE041CIIg6BJ6v55STP8GQRAEQd2IkWoQBEHQ6RjHkZqehKTdJVnSmjXKHStphwrpIyRdX6PuGEmH1sHWkZLWbke9r0j6Wb4/SNJ+HbUlCIqsecxvu9uEIOjR9KaR6j7AfcDXgTHVCtk+qqsMaoGRwPXAE62tIGmA7fHAeADbZ3eSbUEQBPXHMKf3D1R7x0hV0kBgS+AAklMtpR8uaZqkKZJOymnjJO2R73eR9JSk+4CvFuotI+laSVMlPShp/UJ3G0i6Q9Kzkg4s9S/pdkmP5P52K7S1X25niqSLJG0BfAX4taRHJQ3N102SJkm6tzTazraeJulO4GRJoySNzXnzRs2SDpQ0IfdxlaRPdMbnHARBELRMbxmpjgRusv2MpDclbQQsn9M3tf2+pGWKFSQtCpwLbAc8B1xeyD4GmGx7pKTtgAuBYTlvfWAzYHFgsqQbgNeB3W2/I2k54EFJ44G1gSOALW3PlLSM7Tdz3vW2r8y23A4cZPtZSZsCZ2a7AFYHdrA9R9KoKs9/te1zc1u/Iv24OKNSQUmjgdEAgwcPbuEjDYIgqB8pSHnvp7c41X2A0/P9Zfl9P+BPtt8HsP1mWZ01gRdsPwsg6WKyswG2Ar6W690haVlJS+W8v9qeDczOI8hNgBuAEyRtQ/p3sxLJqW8HXGl7ZhUbSqPsLYArJJWSFykUucL2nBrPv252pksDA4GbqxW0fQ5wDkBTU1MfmIwJgqAxEHNQ7WI9nB7vVCUtS3Je60oy0J/0o+iq/NoS1fIrffMuey2m7wsMAobb/kjSdGDR3E4tG/oBb9seViX/vRr1AcYBI21PyaPZEa2oEwRBENSZ3rCmugdwoe3P2h5ie2XgBeBN4Nul9cXy6V/gKWAVSUPz+30KefeQHCWSRgAzbb+T83aTtGh25iOACcBSwOvZoW4LfDaXvR3YK5ct2vAusARAbvcFSXvmMpK0QRs/gyWAVyUtVLI7CIKgkTAw1/W/Go3e4FT3Aa4pS7sKWJG0U3aipEeB+Y7C2P6ANN17Q96o9GIhewzQJGkqcBKwfyHvYdJ074PAcbZfAS7J5SeSnNpTuY/HgeOBuyVNAU7LbVwGHCZpcnbq+wIH5DKPA7vRNn4JPATcWuo7CNrDU0f/uLtNCIIejdwHtBiDyjQ1NXnixIndbUYQBD0ASZNsN7W3/rrrL+y/3DConiYBsM7gVzpkV73p8WuqlZA0B5hGWtOcAxxs+35JQ0i7btdtR5u/sH1Cne0cAfzX9v01yo0CmmwfLGkMMMv2qZKOBe6xfVs97Qr6NsMO7joBiEfHxsi4r5CClPf+jUq9Yfq3ErNtD7O9AfBz4MQ6tPmLOrRRzgjSzt92YfuocKhBEASNQ291qkWWBN4qT5Q0JAstPJKvLXL6CpLuycIMj0naOgtHLJbTLsnlfpLzH5P0o0K784k95LRdJT2U11Bvk7R8HjUfBPw4t7u1pEFZvGFCvrZs6cE0v5DFUbnOY5LOUeF8ThAEQSMw16r71Wj0yulfsgMkHWtZgWYhhSKvAzva/kDSasClQBPwDeBm28dL6g98wva9kg4uHXuRNBz4FrApaYr5IUl3A/+lTOwh93UfsJltS/oOcLjtn0o6mzyVm9v9M/Bb2/dJGkw6b7pWK595rO1jczsXAV8Grmv1JxYEQRB0mN7qVGcXHODmwIWSytdRFwLGShpGWnddPadPAM7Px1Outf1ohfa3Aq6x/V7u42pga9KyQSWxh88Al0taAViYdOSnEjsAaxcGmUtKWqKVz7ytpMOBTwDLkHYRL+BUQ1EpCILuINZUewm2HwCWI4kzFPkx8BqwAWmEunAufw+wDfAv4CJVjgRT7V9GNbGHM0gjyfWA75JG0JXoB2ye14OH2V7J9rtVH67UaZJcPBPYI/dxbrU+bJ9ju8l206BB9d+JFwRBUAkj5tCv7lej0XgW1Rklcfr+wBtlWUsBr9qeC3wzl0HSZ0lCDucC/wdslMt/lEevkMQhRkr6hKTFgd2Be6ku9rAUyUnD/Gde54lAZG4BDi7YXk1lqZySA52ZZQ/3aGW9IAiCoI701unf0poqpNHj/lmQvljmTOCqrGR0J81ygCNIwgwfAbOA0kj1HGCqpEds7ytpHEkIAuA825MBJJXEHuYAk4FRJDGJKyT9iyQasUqudx1wpVJUmx8APwT+kEUnBpCc90G1Htb225LOJR0jmk6awg6CIGgoGnFjUb0J8Yc+TIg/BEHQWjoq/rDW+ot43HUr1tMkADYbMj3EH4IgCIK+RV/ZqNQnnaqkWbYHFt6PIisWtVBnDIXjL23sr911g6CrWeWM33S3Ca3mhR/8tLtNCFqNmONev42n929UCoIgCIKuIpxqGZI+K+n2rIp0exZhKC8zTNKDucw1kj6Z0++SdLqk+7Oy0SaFamvn/Ocl/bDQ1rWSJkl6PJ8hLaUfIOmZXOdcSWNz+gLqTDl9cUnnZ1WlyXnzUxAEQUNgYC796n41Go1nUddQkhx8NO8SPraQN5YUn3V9Uki331eofyHw/3KZacDRhbzFbW8B/C9wfiF9TWBnYBPg6MLxnG/bHk46K/tDSctKWpEUzm0zYMdct0RJnWlDUgi5w3P6EcAdtjcGtgV+nY/7zIek0ZImSpo4Y8aMFj+kIAiCoG30yTVVCopL0Lymmt9uDnw1318EnFKsKGkpYGnbd+ekC4ArCkUuhSQiIWlJSUvn9Btsfwh8KOl1YHngZZIj3T2XWRlYDfg0cHdJkUnSFTQrPlVTZ9oJ+IqkUtzYRYHBwJNF+22fQzoeRFNTU2z9DoKgy4iNSgFUVkhqS/nS+w8LaXOAAUqh33YgqSi9L+kukjNs6V/eGcBptsfn+mNyuoCv2X66jfYGQRAEdaKvTv+2xP3A1/P9vqTp1nnY/g/wlqStc9I3gbsLRfYGkLQV8J9cvhpLAW9lh7omaboXkqjE5yV9UtIA4GtldSqpM90M/KAUnUbShjWfNAiCoIuw0+7fel+NRoxUF+SHJEH9w4AZpGg05ewPnC3pE8DzZWXeknQ/KeTct2v0dRNwUFZQepqktoTtf0k6AXgIeAV4Aig55zFUVmc6DjidpPokkrLSl1v5zEEwjzimEnQWc/vA9G8oKtWRPH17qO0OyxRJGmh7Vh6pXgOcb/uajrZbJBSVgiBoLR1VVFp9vcV8xvhVahdsI7t87slQVOppSDqCFGd1DjAX+K7th6oUHy1pRB2EHsZI2oG0xnoLcG0H2wuCVjH0shM6vY9/fP0Xnd5H0FgkRaXGm66tN+FUa5DjsX4Z2Mj2h5KWI4eJK8f2iKye1JH++tueY/vQ2qWDIAiCRqL3/2zoOCsAM/NxGGzPtP2KpOnZwSKpKU/9lthA0h2SnpV0YC7TT9KZWeThekk3Stoj502XdJSk+4A9a4hLNOX75SRNz/ejJP1V0k2SnpZUPDcbBEHQAMRGpSBxC3CUpGeA24DLC2dUq7E+aSfv4sBkSTcAWwBDgPWAT5HOjxbFIT6wvRVA3rj0A9t3SzqWJC7xoxp9bgKsC7wPTJB0Qz3WdoMgCOpBSVGpt9P7n7CD2J4FDAdGk3YDX57FIlrir7Zn255JitW6CbAVcIXtubb/ndOLXA5VxSW2aYWpt9p+w/Zs4Orc3wKEolIQBEHnEU61FeQ1zrtsHw0cTDo3+jHNn9+i5VUqvK+1l/y9Gvm0o88FsH2O7SbbTYMGDWpFl0EQBPVhjlX3qzVI2iUvjT0n6WdVyuwl6Ym8RPfnnDZM0gM5baqkvWv1FU61BpLWkLRaIWkY8CLpHOjwnPa1smq7SVpU0rLACGACSUTia3ltdfmcvgA1xCWKfe5RVnVHSctIWgwYCfy9tc8YBEHQW5HUH/gD8AVgbWAfSWuXlVkN+Dmwpe11aF5uex/YL6ftApxekJ6tSKyp1mYgcEb+ID8GniNNBa8F/J+kX5BEGoo8DNxA0t49Lm9sugrYHngMeCbXqaa2VE1c4lTgL5K+CdxRVuc+klbxqsCfYz01CIJGwqi7jtRsAjxn+3kASZcBu5FEdUocCPzB9lsAtl/Pr8+UCuS/468Dg4C3q3UWTrUGtieRNhmVcy/NIvfF8mOqtDNX0qFZ0GFZkuOdlvOGlJV9lGbJwmL6U6RNUCWOLNy/3lKQ9SBoLXGGNOgs5nbPbt2VgJcK718GNi0rszqApL8D/YExtm8qFlAK5bkw8I+WOgun2rVcn0e8C5NGsP/uboOCIAh6OMtJKs7MnZOjcZWotPBavudkAClC2AhSJLB7Ja1r+22AHBXsImB/23NbMqbVTjWPrm7Pbz9NUheaQTom8orttatU7VNImmV7YKU82yPq0P444HrbV0o6jxSxZhwwrqNtB32H1a88rmreM3v8sgstCfoKnaioNLOGTOHLpLCaJT5D0lQvL/Og7Y+AFyQ9TXKyEyQtSVrOO9L2g7WMafUT5uMaw3Ic0rOB3+b7YSTpvnaRtW17DEo0xAYv29+x/UTtkkEQBH2WCcBqklaRtDApCtn4sjLXAttCEtYhTQc/n8tfA1xo+wpaQb2cQ39J5+Ztx7fkHahIGppVfiZJulcpvBmSxkk6TdKdwMmSxki6INedLumrkk6RNC3XXyjX217S5Jx+vqRFcvpRkiZIekzSOdK88Gd3STpZ0sOSnintqJU0JNvzSL7mrZlKOiy3NVXSMYXyT0o6E3gEWFnSPtmOxySdXPwwJP0mt3u7pEE57cDc7hRJV+VNSEhaXkk1aUq+tsj9PVZo71BVkD/U/ApLZymdP328ZHcQBEGjYOp/nKY1R2psf0w6CnkzSXTnL7Yfl3SspK/kYjcDb0h6gqQhcJjtN4C9SDoBoyQ9mq9hLfVXL6e6Gmnn1DqkXVGlIybnkJSBhgOHAmcW6qwO7GC7FGdqKPAl0q6si4E7ba8HzAa+JGlR0hTn3jl9APC9XHes7Y1trwssxvwhzwbY3oS0Rbok3/c6sKPtjUjxT38PIGmn/CybkEbgwyWVhBfWIP1a2RD4CDgZ2C6X21jSyFxuceCR3PbdhT6vzjZuQPpiD8jpvwfuzukbAY+38Dm3xBF5CmR9UizW9SsVUog/BEHQx7B9o+3VbQ+1fXxOO8r2+Hxv2z+xvbbt9WxfltMvtr1QaZY2X4+21Fe9nOoLhY4mAUMkDSTtmr1C0qPAH0k6uiWusD2n8P5veT57Gmn3VWnn1TTSuu0auZ/SFuei0tC2kh6SNI3k6NYptHt10a58vxBwbi5/BensEsBO+ZpMGpGuSXKyAC8W5tM3Bu6yPSP/CrqkYMtcsjoS6cdBSdlo3Tw6nkYKfl6ycTvgLJgnMtFSUPOW2EvSI9n2dQrPNB8h/hAEQXcxl351vxqNeq1nfli4n0MaLfYD3s7rrpUoVxAqCdbPlfSRmwO9zs12Vhzn5xHsmUCT7ZfyNGlRbahk2xyan/fHwGvABtnOD0rNASfa/mNZH0PK7G1LpN3Sc4wDRtqeoiRzOKKFOkXlJFhQPWk+JK1CmgnY2PZbeTNTi3WCIAi6EpuGFMCvN532hLbfIe2i2hPmbfDZoANNPkUaAa+a35eUhkrOY2YeHZcrDVViKeDVvDX6m6SRMaR59W/ndpC0kqRPVaj/EGmKdTkltY59aFY96lew4RskUQaAJYBX8/rwvoW2bidPY0vqr7TT7DXgU5KWzevGxensSixJcvr/UVJr+kLNTyAIgiCoO52983Zf4CxJR5KmXC8DprSnIdsfSPoWaTp5AGlH19k5xum5pGni6Tm9FmcCV2WHfyd5FGr7FklrAQ/kvU6zgP8hjXKLtrwq6ee5roAbbf81Z78HrCNpEkkxqaQV+UuSM34x27pETj8EOEfSAbmf79l+QCk6zUPAC6QfFC19NlMkTSatxz5PSBQGNYhjM0HXI+a2aZKvZ6LmWdagr9HU1OSJE0PNMAiC2kiaVOM8aIsMWXcJH3l1ixtn28WBa9zXIbvqTY86I9qZSDJwse1v5vcDgFeBh2x/Oa+D/hr4F2nU/SRJaPn9Ku1tSNrstIvtm3PaysCFJPGMuSTlj99VqX8+adr39byruZS+DGkj1BDSyHyvvI66FGlj1GDS93qq7T+1+wMJ6s7PprZmZaL7OWn9K7vbhKAXYmJNta/xHmmH7mL5/Y4kB1rk8ryleh3gvzRP7VZiH9J66j6FtI+Bn9pei6Tt+32VRUsoMI4UFaGcnwG3216NtB5bCmP0feCJfDRnBPAbpYPLQRAEDcEc+tX9ajQaz6Lu5W+ks7KQnOGllQrlUeziwFtV8kXarDQK2CnvUMb2q7Yfyffvkka7K1Vqw/Y9wJsVsnYjHSciv5bOxxpYIvc9MNf9uMpzBkEQBJ1AONX5uQz4enaC67NgSLe985nbfwHLANdVaWdL0pnafwB3AV8sL5CP6WxYoY9aLG/7VUhOGijtTh5LCkf3Cmkj1CG1hJ+DIAi6CiPmuv5XoxFOtYDtqaS1yn2AGysUuTyfu/00yXEdVqWpfUgOmvxanAImH9m5CvhRPnpUD3YGHgVWJKk8jc3Hc+YjFJWCIAg6j3CqCzKeFAy84tQvJEkr0ih1m3y2tKQJeWw+t/o14ChJ04EzgC9IWgIgn1O9CrjE9tU5beVCGwfVsO81pTBEpXBEr+f0b5GkEG37OdJRnDUr2B6KSkEQdAt9YU01dv8uyPnAf2xPkzSihXJbAf/IUovz9olL2hmYYnvnQtoFwEhJFwP/Bzxp+7RSvu2Xim3UYDywP3BSfi2dj/0nsD0pDuDyJFnH51vZZhAEQadiui1IeZfS+5+wjdh+udoxF/KaqqSppPXQSkEp9yGFCipyFUldaUuSgtN2hZHpAuutAJIuBR4A1pD0chaHgORMd5T0LGmH8kk5/Thgi6wtfDvw/2zPbM0zB0EQBPUhxB/6MCH+EARBa+mo+MPK6y7lQ67YrJ4mAXDY2reE+ENnIGkOafPQANJ64jdtv12jznSSEH+7RnRZvH+W7VPbU7/QzkjgmbYGHM+xANe2fVLNwkEA7HrvD2qWuW7rM7rAkiDonfSm6d/ZWZhhXdIZze/Xo9F8JrWzGUmVUG3VkDTA9vhwqEEQ9ARKa6r1vhqNxrOoPjxAFlWQNELS9aUMSWOz5GCJwyQ9nK9Vc5lxkk6TdCdwsqRlJF0raaqkB8sCgG8g6Q5Jz0o6MNcfKOl2SY9ImiZpt0L/++V2pki6SNIWwFeAX+c11qH5uknSJKUYrGtWsWuUpLE5b1elmLKTJd2WNysFQRAEXUivmf4tkY+0bE/aZdsa3rG9iaT9gNNpDrO2OrCD7TmSzgAm2x4paTuSfm9pt+76JMnBxYHJkm4gHXPZ3fY7kpYDHpQ0njQaPQLY0vZMScvYfjPnXW/7yvwMtwMH2X5W0qakqDrbVbBrVOE57gM2s21J3wEOB37a6g8uCIKgk5nTB6LU9CanulhWOxoCTAJubWW9Swuvvy2kX5GPy0A6PvM1ANt3KMU5XSrn/dX2bGB2HkFuAtwAnCBpG5Jw/krA8iTHeGVpDdf2AjKEWRhiC1KIu1LyIlXsKvIZ4PJ8dnVh0rryAkgaDYwGGDx4cKUiQRAEdcdWQ07X1pve9ISzs9rRZ0lOpbSm+jHzP+eiZfVc5f69wn2ln1cuey2m7wsMAoZnm17L/apC+XL6AW/n9eHStVYVu4qcAYy1vR7wXRZ8zmRciD8EQRB0Gr3JqQJg+z/AD4FDs3rRi8DakhbJo8vty6rsXXh9oEqz95AcJVkQYmZBXnA3SYtKWpYUHWYCsBQpZNtHkrYlOXpI50f3ymVLYdwA3iUHLc/tvqAUQB0lNmjFoy9Fc1Sd/VtRPgiCoEuZ4351vxqN3jT9Ow/bkyVNAb5u+yJJfwGmAs8Ck8uKLyLpIdIPjH2ozBjgT1n04X3md1oPk6Z7BwPH2X5F0iXAdZImkvR4n8p2PS7peODufARoMimSzWXAuZJ+SIpusy9wlqQjSbFbLwOm1HjsMaQp438BDwKr1CgfBEEQ1JkQf+jDhPhDEAStpaPiDyuus7RHX/75epoEwDHrjQ/xhyAIgqCvoYacrq03vdapSjJwmu2f5veHAgNtj8mRYN63fWE+lnKL7VeqtHMscI/t2yTdBRxqe6KkG4Fv1FJtaqWtw4AVbVcKN9dSvSZgP9s/7KgNQfcz5MLG0PGYvt/PutuEIOix9FqnCnwIfFXSieUyhLbPLrwdBTxGCu49H5L62z6qUuO2KwUeF2lKva3BwYcBTVSO4VqRrKg0EYj52yAIGp6kqNT7z6n25rH4x8A5wI/LMySNkXSopD1IzuyabZm0AAAgAElEQVSSrGa0mKTpko6SdB+wZ1Yx2qNCG9MlLSdpiKQnJZ0JPAKsLGlWodweksbl+3GSzs4qSc9I+rKkhYFjaY6As3c1Bads9zmSbgEuVEEtStImku7Pikr3S1qjzp9nEARBUIPePFIF+AMwVdIplTJtXynpYPKULkAWXPjA9lb5/S6t6GcN4Fu2/7fQRjWGAJ8HhgJ3AqsCR5GE/Q/O9VtScBoObGV7tuaP9/oUsI3tjyXtAJxAFqwIgiBoBBoxqHi96dVONcsEXkg6tzq7DVUvb2NXL9p+sJVl/5Knh5+V9DywZoUyLSk4jc8KTuUsBVwgaTXSTMtClToPRaUgCLoDo5j+7SWcDhxA0uZtLdVUi1pbvnhOqSUFp0rvoWUFp2q2HQfcmaP07Fqh39RIKCoFQRB0Gr3eqWZ93b+QHGsl5qkZ1ZHXJK0lqR+we1nenpL6SRoKfA54uoINLSk4VaOoqDSqY+YHQRDUn7n0q/vVaPTq6d8CvwEOrpI3Djhb0mxg8zr19zPgeuAl0s7igYW8p4G7SQL7B9n+IAvx/ywHBDiRlhWcqnEKafr3J8AddXqOoAuJoyxB0PMJRaUuJO8CnhfirbsJRaUgCFpLRxWVll97GX/jzzvV0yQATt/w8lBUamTybt/fAf2B82yfJOn/SEdvBDwDjLI9K4d2O50UU/XrRWcp6WTgS/ntcbYX2PyUxfaL4ebWzO1cWxSaKJRfFrgS2BgYV9otXNbmeOBzeW01CNrM0MuPb1e9f+x9RJ0tCYKeRzjVAkoBzv8A7Ai8DEzITurHpTVNSaeRppJPAv5JWr88tKydLwEbkY7BLEIS0P+b7VHFcrbvzGVKEWueA25pwcQPgF8C6+ar3P6vArPK04MgCBqB2P3b99gEeM7287b/S4oOs1vBoQpYjLwT1/Z021NJgciLrA3cbftj2++RIszUOu+6B/A32+9XK2D7Pdv3kZzrfCgFN/8J8KtWPGcQBEGXko7U9Kv71Wg0nkXdy0qkzUUlXs5pSPoT8G/SFO0ZNdqZAnxB0ickLQdsC6xco87XgUvbY3TmONKGrKpOOQiCIOhcwqnOT9Xzoba/BawIPElzYPOK2L6FpON7P8lRPkCSTazcqbQCsB5wc7uMToL8q9q+phVlR0uaKGnijBkz2tNdEARBu5iD6n41GuFU5+dl5h9RfoaC0L7tOSS1pZryf7aPtz3M9o4kZ/2spE2zvu+jkr5SKL4XcI3tj9pp9+bAcEnTgfuA1fNGp0p2hfhDEARBJxEbleZnArCapFVIQgpfB74haVXbz+U11V1JOrtVyRuelrb9RhbDX58UXu5jmjV8i+wD/Ly9Rts+Czgr9z2EdGxnRHvbC4IgqDd9JUpNONUCWYz+YNI0bH/gfNJ0772SliSNOKcA3wOQtDFwDfBJYFdJx9heh6S7e28W1n8H+J/sUBcgO8GVSYIQ5dwgqTR6fcD2nnk0uiSwsKSRwE62n+joswdBEHQuasiNRfUmnGoZOVB4eVzTLauUnUCaIi5P/4C0A7g1/U0nb4YqSx9RpfyQVrQXZ1SDdhPnTYOg/YRTDYIgCLqEuQ24sajeNLRTzQpCt+e3nwbmADNIMUlfsd2q0WBnUEnxqA5tjqOCjKGkFYHf294jC+wfavvLebPT2rZPqpcNQecw/G89Z/Q36QvtU1QKgqDBnartN2hWHBoDzLJ9amkzTnvblTSg2hpnI2L7FZI4RHn6eGB811sUBEHQNmyY0wc2KvXkVeP+ks6V9LikWyQtBiBpqKSbJE2SdK+kNXP6OEmn5YgwJ0saI+mCXHe6pK9KOkXStFx/oVxve0mTc/r5khYpN0TSTpIekPSIpCuyulHVurm/kyU9nK9VC81tI+l+Sc9L2iOXHyLpsQr9jpI0Nt/vKumh3N9tkpav78cdBEHQMUJRqbFZDfhD3m37Ns1nR88BfmB7OEmT98xCndWBHWz/NL8fShK93w24mBTkez1gNvAlSYuSQsPtndMHkHf+lsiKSUfmdjcCJgI/aUXdd2xvAowlifKXWAHYCvgySV+4tdwHbGZ7Q5K84uFtqBsEQRDUgYae/q3BC7YfzfeTgCF5hLgFcEU+zgJJ0L7EFVnAocTfbH8kaRrpCM1NOX0aad12jdzPMzn9AuD7zO8ENyPt9P177nNhkoJSrbqXFl6LkWqutT0XeKKNo83PAJdndaaFgRcqFZI0GhgNMHjw4DY0HwRB0H6S9m/vn/7tyU71w8L9HJLQfT/gbduVBBYA3qvUhu25kj5yc3DZuaTPpjX/AgTcanuf+RKTdGBLuMp98bna8i/wDOA02+PzZqYxFTu1zyGN5mlqaopgukEQBHWkJ0//LkCOJvOCpD0hRZWRtEEHmnyKNAIurXl+kwVFGh4EtiyVySL6q7ei7t6F1wc6YGOJpUgqUAD716G9IAiCujIX1f1qNHrySLUa+wJnSTqSpGx0GUkFqc3Y/kDSt0jTyQNIMoZnl5WZIWkUcGlhE9ORtp+pUXcRSQ+RftjMN8ptJ2NyX/8iOfpV6tBmUCfimEoQdB+SdgF+R1rmO6/SMURJe5H+jhqYYvsbOX1/0r4ZgF/ZvqDFvppnPIOuIksNNtme2Z12NDU1eeLEuh2zDYKgFyNpku2m9tZfZq1B3vlPu9fTJAAu2/zcFu3KWuzPADuSgqZMAPYpyrtKWg34C7Cd7bckfcr265KWIW0+bSI520nAcNtvVeuvV41UJe0OXA2sZbuq6L2kY4F7bN9Wlj6CDggrSLoR+Ibtt9v1AM3tjCGfyS1LL4pANAH72f5hR/oKehZrXzum0/t4YmTn9xH0TbrpCMwmwHO2nweQdBnpxEdRM/1A0mmStwBsv57TdybtmXkz170V2IUWYl/3KqdKmka9jxRdZky1QraPqtVQe4QVbH+xleWG5CnhNlEUgchKTjHMDIIgaJmVgJcK718GNi0rszqApL+TpojH2L6pSt0FtNqL9JqNSvk4zZbAASSnWko/PIsvTJF0Uk4bVxBW2EXSU5LuA75aqFcUVhgn6SxJd2ZRhs9nMYcnlaQFS3Wm53OrSNpP0tTc70WFdooCFMtIujaXe1ApTFyJDSTdIelZSQfm+vNEICSNkHR9vt8kC0ZMzq9r1P0DDoIg6AhOR2rqfQHLSZpYuEaX9VxpN1P5uucAkvbBCNLg7DxJS7ey7gIN9RZGAjflDUJvStoIWD6nb2r7/Tw/Po8s0HAusB3wHCkAeTU+mct9BbiO5MC/A0yQNKxwZhZJ6wBHAFvanlnWb0mAYo6kM4DJtkdK2g64kOZ4q+uTzsAuDkyWdEMLtj0FbJND1+0AnEArAqkHQRD0AmbWWOt9mRRes8RngFcqlHnQ9kekEyRPk5zsyyRHW6x7V0vG9JqRKunXxWX5/rL8fgfgT7bfByjNixdYkyTQ8Gw+o3pxC+1fl8tMA16zPS2LNDxOEooosh1wZWkjUlm/RQGKrYCLcpk7gGUlLZXz/mp7dm7jTtK6QDWWIu38fYwkJLFOtYKSRpd+0c2YMaOFJoMgCOqH6bYjNROA1SStImlh0kxm+dLetcC2ME8lb3XgeVJs7Z0kfVLSJ4GdclpVesVIVSmazXbAupJMmhM3cBU1huqtyC9REmWYy/wCDSWhiPlMaqHd98rKVbOnvH5Ldh5HkljcXSnYwF3VCob4QxAE3UV3KCrlGbyDSc6wP3C+7cfzhtWJef9MyXk+QRITOiwHdEHScSTHDHBshcHZfPSWkeoewIW2P2t7iO2VSTJ9bwLflvQJgPLpX9K06SqShub39TgvCilc3V7Z2Vfqt8Q9pHO1pZ3HM7OABcBukhbNbYyg+UutRFH4YVSHLA+CIOhl2L7R9uq2h9o+PqcdlR0qTvzE9tq217N9WaHu+bZXzdefavXVW5zqPsA1ZWlXASuShvkTJT1KEtifh+0PSDq4N+SNSi/WwxjbjwPHA3dLmgKcVqXoGKBJ0lSSeH5RCelh4AaSkMNxeedvNU4BTizsXAuCIGgoDJ21UamhCPGHPkyIPwRB0Fo6Kv6w9Jqf8ufP26ueJgEwfus/dMiuetMr1lSDIAiCxqcRR5b1plc41TYoKY0CbilNpUo6jxTZ5YlqdVpoaxxwve0r22V0FZvaUO8g4H3bF3ak/6BjrHH1sd1tQt15+qs1tVGCoM30ldBvvWlNtaSk1BKjSOusANj+Tnscap0ZRcGm1iBpgO2zw6EGQRA0Fj3eqbZWSSkrKDUBl0h6VNJiku7KGrpImiXp+Fz+QeUA4ZI+K+n2rHp0u6RiZO8dJN0r6RlJX87lh+S0R/K1RRttGi7pbkmTJN2sFHScbOsJku4GDpE0RtKhOe9ASRNyu1eVdjsHQRA0En0h9FuPd6oUlJSANyVtJOkLNCspbQCckqdpJwL72h5me3ZZO4uTFDU2IB11OTCnjyUd11kfuAT4faHOEODzwJeAs7NC0+vAjrY3IsVK/T1Aa2wCPiYFG9/D9nDgfNIu4hJL2/687d+U2X617Y1zu0+SfmAEQRAEXUxvWFPdBzg935eUlPrRspJSJf4LXJ/vJ5HCBAFsTrMm8EWk4ysl/pJVlZ6V9DxZoQkYK2kY6RDx6rlsLXUngDWAdYFbJUE6HvNqIb+ajOK6kn4FLA0MpAXFj6yLORpg8ODB1YoFQRDUF8dGpYang0pK5Xzk5vNFc6j+2bjKfen9j4HXgA1Izv2DkrmtsEnA47Y3r5L/XpX0ccBI21PyxqcR1ToIRaUgCILOo6dP/7ZVSeldYIk29nE/zWu1+5I2RJXYU1K/rMj0OeBpkrrRq3kE+02axRhuaYVNTwODJG2eyyyUxflrsQTwqqSFso1BEAQNRV8Rf+jRI1XSVG95EPGrgLVoVlL6L3Aj8AvSiO5sSbNJ07qt4YfA+ZIOA2YA3yrkPQ3cTYqGc5DtDySdCVwlaU+SEP57ALZvylPCtWzaA/h9FtYfQJrafryGjb8EHiIpQk2j7T8cgnYSx0+CoPU0ohOsN6Go1IcJRaUgCFpLRxWVllxjeW98dv0n0u7Y7rehqNQXyWudTbYPljQSeKazzshKGgPMsn1qZ7QfdJzDp+zZ3SZU5ZQNruhuE4JeSIg/BJ3JSGDttlSQFD+AgiAIGpxwqpks2vBY4f2hWWDhLkknS3o4izxsnfP7Szo1izlMlfSDnD5dKcgtkpok3VXWzxbAV4BfZ8GHoWUiFMtJmp7vR0m6QtJ1pI1OSDosCz1MlXRMod0jJD0t6TbS0ZwgCIKGwlbdr0YjRj+tY4DtTSR9ETiadOZ0NLAKsGEOglstZup82L5f0ngKusH5TGo1NgfWt/2mpJ2A1YBNSMdvxkvahrQZ6uvAhqTv9BHSWdsgCIKGoREVkOpNONXWcXV+nURSUYLkWM+2/TG0WmCiPdxaaHunfE3O7weSnOwSwDUlYYnstCsS4g9BEASdRzjVZj5m/unwRQv3H+bXoihENTGHYjuLVsiv1Xd5naLgg4ATbf+xWEDSj6rYsgAh/hAEQXfgPqKoFGuqzbwGfErSspIWAb5co/wtwEGlDUSF6d/pwPB8/7UqdctFKIp19mihz5tJAhIDc58rSfoUSat49yzIvwSwaw3bgyAIgk4gnGrG9kfAsSQRheuBqnFZM+cB/wSmSpoCfCOnHwP8TtK9pJFtJS4DDpM0OasxnQp8T9L9wHIt2HgL8GfgAUnTgCuBJWw/QtIFfpQkfnFvrecNgiDoavrCRqUQf+jDhPhDEAStpaPiDwNXX8Hrjd2/niYB8ODOJzeU+EOMVIMgCIKgTvSojUqS5pC0bUuMJE2X7mf7h21oZzpJ3WhmjXJfAda2fZKkg4D3bV/Ydsvri6RZtgdKWhH4ve2W1mGDoE1s+L+/7W4TGorJZ/64u03oNTTidG296VFOFZidg3kXmU4K9F1XJA2wPZ4kzI/ts+vdR0ex/Qotb2wKgiAIupAeP/0raYSk6/P9GEkXSLolKxt9VdIpWfXophwarcRhWSXpYUmr5vrjJJ0m6U7g5KxoNLbQ9qH5/i5Jv5V0j6QnJW0s6WpJz+Zg4SXbfiLpsXz9KKcNyXXOlfR4tnWxnHdgVkuaIukqNYeJW0XSAznvuEL781Sg8v29kh7J1xad+bkHQRC0hb4S+q2nOdXFsrTfo5KuqVJmKPAlYDfgYuBO2+sBs3N6iXdsbwKMJYVXK7E6sIPtn9aw5b+2twHOBv4KfB9YFxiVj+UMJ4WJ2xTYDDhQ0oa57mrAH2yvA7xN89Gbq21vbHsD4EnggJz+O+As2xsD/65iz+vAjrY3AvYGfl+pkKTRkiZKmjhjxowajxgEQRC0hZ7mVGfbHpav3auU+Vs+HjONFCD8ppw+jWY1JIBLC6/F2KpX2K52FKZISbVoGvC47Vdtfwg8D6wMbEVSOXrP9iySKtPWuc4Lth/N90WVpnXzaHMaKdh4KUD5lgV7L6piz0LAubnuFVQR7Ld9ju0m202DBg1qxWMGQRDUAScBiHpfjUZPW1NtDR8C2J4r6SM3nxmay/zP6yr3RQWjmv3kdj8spJf6aWleolh+DrBYvh8HjLQ9RSlU3IgqNlbixyQBiw1IP5Y+qFE+CIKgS+kL2r89baRaT/YuvD7QCe3fA4yU9AlJiwO7U1uUYQng1bz2W4zm+3eSYD5l6UWWAl61PRf4JmmUHgRBEHQhvXGk2loWkfQQ6YfFPvVu3PYjksYBD+ek82xPljSkhWq/JCk6vUiaVi5JGR4C/FnSISTFpEqcCVwlaU/gTlo/4g6C+YgjJEFnYPrGkZpQVOrDhKJSEAStpaOKSp9YbQWv+bsDahdsI5O/dHxDKSr15ZEqMJ+ghEjrmwfbvr8d7YwEnrH9RJ3smic8USEvxB+CTmPnpjHdbcJ83DxxTHebENSFxjwCU2/6vFOlICghaWfgRODz7WhnJEmIv9VONQtMfFwpryg8UY0QfwiCoCfRFyZG+/JGpUosCbwFIGmgpNuzkMI0SbuVCknaT9LULNJwURZa+Arw63yGdmi+bpI0KR+TWTPXLReYWEbStbm9ByWtn8sVhSdC/CEIgqAHECPVLChBCg6+ArBdTv8A2N32O5KWAx6UNJ50/vMIYEvbMyUtY/vNnHe97SsBJN0OHGT7WUmbkjYSldouCUzMkXQGMNn2SEnbARcC5VKMJfGHCyV9v8pzlMQfPpC0Gulca8OsMwRBEPSFjUrhVOef/t0cuFDSuqQ11hMkbUM6e7oSsDzJMV5ZEuO3/WZ5g0pBxLcArpDm/SNapFCkKDCxFVlRyfYdWY1pqbImt6RZdeki4OQKz7EQMFbSMNLa8OqVHlbSaGA0wODBgysVCYIgCNpJONUCth/Io9JBwBfz63DbHylFtlmU5GxrrQz0A96uIP5fonjcpdJPt0rt10X8wfY5wDmQdv/WaDMIgqAuJAWk3j9SjTXVAnndsz/wBklM4fXsULcFPpuL3Q7sJWnZXGeZnP4u+Vyp7XeAF/KZUZTYoEq395AFHSSNAGbm+kVC/CEIgh5PCOr3DeaJ9AOXA/vnqdlLgCZJE0mO7CkA248DxwN3S5oCnJbbuYwU+WaypKG5zgG5zOMkgf9KjMn9TAVOAvavUOYQ4PuSJpCcZyXOBPaX9CBp6jfEH4IgCLqYEH/ow4T4QxAEraWj4g+LrbqiVzl1dD1NAuDJ3Y9pKPGHGKkGQRAEQZ1o6I1KBbWjhYCPgQuA0/O6YbU6I4BDbX+5Df2Mo3AcphGQ1ATsZ/uHFfKmA035SM/9tuNMai9mlUtO7NL+Xtj3513aX9B36AsblRraqTL/cZdPAX8mrSke3a1WdQG2JwI152bDoQZB0BMw6hNOtcdM/9p+nXS+8uC8m7a/pF9nlaGpkr5bKD5Q0pWSnpJ0ifJhUUlH5fKPSTqnlF5E0vZ5s9E0SedLWiSnbyzp/qyi9LCkJSQtKulPuezkvEuYbNupOX2qpB+0o40Rkq7P98tKuiXn/5HCMRxJs/JrVQWoIAiCoGvoMU4VwPbzJJs/BRwA/Mf2xsDGwIGSVslFNwR+RFI/+hxJPAFgrO2Nba9LCgw+3xSxpEVJgcL3tr0eaST/PUkLk3YGH2J7A2AHYDbw/WzXeqTwcRfkNkYDqwAb2l4fuKQdbRQ5GrjP9oYkPeBKqg0lBaiNgG2B31T50TBa0kRJE2fMmFHxcw6CIOgM3AlXo9GjnGqm5Ch2AvbLR2EeApYFVst5D9t+Oa+9PgoMyenbSnpI0jSSMtI6ZW2vAbxg+5n8/gJgm5z+qu0JkM6hZiH8rUgKR9h+ihQHdXWSwzy7JJafVZfa2kaRbYCLc5kbyPrEFT6XE/LRnNtoVoCaD9vn2G6y3TRo0KAKzQRBEATtpdHXVOdD0udIEnyvk5zID2zfXFZmBPBhIWkOMCCP/s4kbfB5SdIYkkLSfNWrdU3lH0VtKd/WNsqp9aNsXyorQAVBEHQ/oajUWEgaBJxNmsI1cDNpanahnL+6pMVbaKLkYGZmbd5KIdOeAoZIWjW//yZwd05fUdLGua8lJA1gfjWk1UnTsk8DtwAH5TIl1aW2tlGkWOYLwCcr2F5NASoIgiDoIhp9pFqKIFM6UnMRzQpG55GmdR/Ja4czSDFNK2L7bUnnko7oTAcmVCjzgaRvkYTwB+QyZ9v+r6S9gTMkLUZaC92BNPI9O08nfwyMsv2hpPNIU7hTJX0EnGt7bBvbKJp2DHCppEdITv6fFR7xEuA6JQWoR8kKUEHPJ464BL2GRlwErTOhqNSHCUWlIAhaS0cVlRYdupJXPul79TQJgOf2+mVDKSo19EhV0izbA/P9F0lxRbcHNgKesf1Ee9uro43TyUIMHWznF7ZPaEe984DT2vpZBN3P5/7c5q+7S3j+G7/obhOCoMfSI9ZUJW0PnAHsYvufpGnetbvXqrrT5r9kkvrb/k441CAIegIp/Ft9r9YgaRdJT0t6TtLPKuSPkjRDObiKpO8U8k6R9LikJyX9vtJRxSIN71QlbQ2cC3zJ9j8kbQF8Bfh1fvih+bpJ0iRJ9yqFcEPSKpIeUBJ8OK7QppSEIx7LQgl75/QRku5SZeGIiqIQmcOUxBweLm1ykrRrPr4zWdJtkpbP6QPVLPYwVdLXJJ1Ec7ScS3K5/8ntPSrpj5L65/RZko6V9BCweba3KeedpXQG9XFJx3TqFxMEQdADyH87/wB8gTQY20dSpUHZ5baH5eu8XHcLks7B+sC6JE2Ez7fUX6M71UWAvwIj8xlObN9PEkA4LD/8P0hBt39gezhwKGnzD6Tp4rOyQMS/C+1+FRhGCui9A8lBr5DzFhCOUBVRiEJ779jeBBgLnJ7T7gM2y4INlwGH5/RfkkQr1svCEHfY/hlZktH2vpLWAvYGtswyjXNojqO6OPCY7U1t31f2eR2R1xbWBz4vaf3aH3EQBEHnY9KRmnpfrWAT4Dnbz9v+L+nvcWsV50w6ObIwyR8tBLzWUoVGd6ofAfeT1JMqko/HbEHasfso8Eeg5CC3BC7N9xcVqm0FXGp7ju3XSDtqN855lYQjqolClLi08Lp5vv8McHPe1XsYzUITO5B+NQFgu5KQw/bAcGBCfqbtSQ4ekoO9qsrHsVfeITw597fArzGFolIQBN2BAav+V21WAl4qvH85p5XztTx7eKWklQFsPwDcCbyar5ttP9lSZ43uVOcCewEbS6q25tgPeLswbB9me61CflsFFxYQjqhRvryP0v0ZpDO16wHfpfmcbDURiHL7Lig8zxq2x+S8D3IQ9fkrJInGQ4Ht8wj4BiqIP4SiUhAEvYzlSgOFfJUHba3097v8b/B1wJD8t/M20sCJvJy3FmmQtBKwnaRtaIFGd6rYfp+k0buvpNKI9V1giZz/DvCCpD1h3nrpBrnc34Gv5/t9m1vlHmBvJeH7QaRR58MtmFFNFKLE3oXXB/L9UsC/8v3+hbK3AAeX3kgqCTl8pCxkAdwO7KEUmQdJy0iqJeawJPAe8J+8fvuFGuWDIAi6lE7aqDSzNFDI1zll3b4MrFx4/xnglfnt8hu2SwOqc0kzhQC7Aw/anmV7FvA3YLOWnrHhnSrM087dBThSKfrKZaTNQZMlDSU5zAMkTQEep3m+/BDg+5ImkJxciWuAqcAU4A7gcNvFNdfy/j8ASqIQ00gj6LMLRRbJG4cOAX6c08bk8vcCxeM2vwI+mTdJTSGJ30NaF54q6ZK8m/dI4BYlLd9baZ7SrmbjFNK07+PA+aQfFEEQBH2dCcBqeePqwqSB1vhigcKeGkgbYUtTvP8k7U8ZkAc9ny/kVSTEH/owIf4QBEFrUQfFHxb53Epe6Vffr6dJALyw7xE17VLSOTgd6A+cb/t4SccCE22Pl3QiyZl+DLwJfM/2U3nn8Jmk2UwDN9n+SUt9NbT4QxA0Ar97aofuNqFLOWTN27rbhKBX0n1Bym3fCNxYlnZU4f7nwAJ6oHn/ynfL01uiR0z/9laUgo+XDhv/W9K/Cu/fL5T7oqRnJQ2WtI1SIPKPJe1R1t7+udyzkvZfsMcgCIKgM4mRajdi+w3SeVmUQtHNsn1qfj8rv5bUpHay/U9J/YBRpJ2+81CKhHM00ESappgkaXyVIztBEARdTx9YbYyRagOjMjUpANvTbU8lbZYqsjNwq+03syO9lbS5KwiCIOgiYqTauJTUpEaU1KRq0KoDzvkM12iAwYMH18HMIAiCVhBByoNupqaaVBmtOeAc4g9BEASdSDjVxqU1alJFah5wDoIg6FbcCVeDEU61gamiJlWNm4GdJH0yqzTtlNOCIAgaBHXC1VjEmmqDY/tNSbsA90iaSRp9XgN8EthV0jG218nljiOphwAcm5Wogg4S5zaDIGgt4VQbhIJgfun9wML9S8AqhezPVGnjfJJEYRAEQePRgNO19Sacap2RNKvkELM01u9Iodu+DRwIzCDFRJ0GHJl1fov1zwC+VXSqZQCsetsAACAASURBVPl7A0eQ5LZusH14IW8vkuawgSm2v1Hfpwv6AkPOOrVV5aZ/79DahYKgjxFrqp1EQbRhF9v/zMm/zaHcVgMuB+7IUXJKdZqApVto8/+zd97hdlXV+n4/irRAFIhKCUSaNEOAJEiTgOi9+lMIAhcxAhEUQRABUVG8VFGwgFeKGLhSFBSp0gSkd0hISAIIiJQLghCkGenh+/0xx05WdvY+JWefknPG+zznOWvPNdtaBzL2nHOMbywH/ISS3m094AMxDpLWpMhsbR73DuyO50qSJJlv0lEpmR8aiTbUY/t8Shq4L0SbhSkG89uN6gerAY/YrmUXvw7YMa6/ApxSU1Cy/XxXnyNJkqRl9F6S8h4ljWrrqYk2jO2AaMNkYO243h+4zPazbdR/FFhb0jBJiwBjmRNGsxawlqTbJd0Vzk1JkiRJD5Jnqq2nKtrwjXbqCkDSisDOwJi2Ktt+SdK+lK3jd2Oc1eL2IsCa0cfKwK2S1rf98lwDpqJSkiS9xEDINJor1dbTGdGGDSkJbzcE1gAelfQEsKSkRyUtXMlaczSA7cttb2J7U+Bh4K/R19PAH22/bfvxuLdm/YCpqJQkSdJ95Eq1G7D9mqTPUFaLz9n+3/o6knakCDR8M85IP1i5N9P2GvFxRF2799t+PgQevkYx4ACXArsCZ0lanrId/Firny1JkmS+GQAr1TSq3UQD0QaAgyR9kRJScz+wTcXpqKP8j6QN4vpo24/EdU1R6UFgFvCtSC2XJJ0iQ2WSbqMPOha1mjSqLaYN0YY/UmJIO9VHg3u7Nik3cHD8JEmSJL1AvzGqNdGFiNc8EVgHeBl4FTjC9i2StgeOoZx7vgMcaPu2aL8H8P3o7ge2z24wxhjgENufic8/AEYB21FWiisAr1f6uLBBH5+KOSxFcVS6wvYhkaT8CGBN249G3YOAE4BRtidJeg9wMsUZ6V3gMNsXRd0UfugFPnVLe75oCx5/+tj/9PYUkn6Kcvt3wULS4sCVFMN3WZStD4wEbgGup4StWNJw4A+UEJVlKQZtJMUo3SvpslrMZ5OxDgM2Bz5t+01JAONsT2qjzfoUo/j/bD8UYTF7V6pMBz4P/CA+7wRUFZcOA563vZakhYBlo9+q8MNLkt7f7stKkiRJWk5/8/4dB9xZM6gAtu+3fVZcz4xtUigrxdr1fwB/tv1iGNI/A03jPCV9E/g08Fnbrzer14BvA8fW4ldtv2P71Mr9S4HtY4zVgFcosoY19gR+FG3ftV07q03hhyRJ+jbdoabUB1e+/c2orkcRVGiKpB0kPURZ0e4ZxSsBT1WqPR1ljdgc2Af4lO2ZdffOrYTALNeg7frAvW1M71XgqVjR7kqJR63NuyZfeIykyZIukPSBKEvhhyRJkj5AfzOqcyHpEkn3S7q4Vmb7EttrU9SIjqlVbdC82XegR6P+JxvcGxfaviO64Hn7e8oW8FhKircai1BEHW63vRFwJ/DTyr2a8MOuwBkVIzwXkvaWNEnSpBkzOut4nCRJMr90g0RhH/Qm7m9G9QFgo9oH2zsA44mzxyq2bwFWj5jOp5kj9wfFeD0Tq9raynNk3HuOsvV7oqSt25qMpP0q7VeM+W3czjNcDuwG/J/tVyvl/wReY46hvaDyrB0SfojnTvGHJEl6h9z+XeA4D9hc0naVsiVrF5LWUHgUSdoIeA/FWNViPN8XogqfBK6JVW1t5TnbASliQz8H/FbSXOIMVWyfUmn/DEUw/3uS1oo5LCTp4Lo2rwPfAY6tKzfF4I6Joo8zx4npUmDr6DOFH5IkSXqJfuH9G160b9p+PZSMTpD0c8qq8l/M8abdEdhd0tuU0Jddwli9KOkYYGLUO9r2i22NaXuipC8Bl7W3Yq20mSbpQOB3kpakfM+6skG93zfp4jvAb+LZZgBfivIUfkiSpO/TB1eWrUbuBwrHoTB0uu3RvT2XBYmRI0d60qSmEUBJkiSzkXSv7ZHt12zMYqsO9QqHtj6u+8mvfatL82o1C/xKVdI+wAFkUu4kSZK+zYK/hmuXBd6o2j4NOG1+2kZIyonAR4GXgLeAH1PCa3a3fUCoKL1l+45oM5aSKPzBzvRp+5LY8j0dGE7xIH4Z+M/60BxJawNnUhyRDrP90yhfnCJisRjlb3eh7SPinijb3DtTtoB/afsX8/NekoHNGsef2NtTmIdHv3NQb08h6Sq1JOX9nAXeqM4vYYQuBc6uSfpJWhXYzvYlQG1fdAwwk5K7FEqoyxXMrXTUbp9R5RvAc7Y/Evc+TMm/Ws+LlNX32LryNyki/DMlLQrcJulPtu+ieDkPBda2/W6qKiVJkvQ8/c37tzNsQ1mBzl7l2n7S9kmSxki6QtIwitDDQREWsxXFQP4kPq/e0T7j4wrA3yv3Hrb9Zv3EbD9veyJ1BteF2qp20fipbajsS3GwerfWR6feRpIkSTcjt/6nrzGQjWq76ku2n6BsLZ8YYTE3A5dRvGtH2P5bJ/v8NfAdSXdK+kFo9naKWuJy4HmKtOLdcWt1YJcQdvjT/PSdJEmSdI2BbFTnQtIpkqZKmth+7fnr0/Z9wGqUeNVlgYmS1ulMn7Zn2R5BEagYHZKGUM5Z3wgvuNMpBrzRnFJRKUmS3iHFH/o19epL+1EEFTosMyRpaEUxaZ+O9Bmi/hfb/hrwW+DTDZSX2sX2y8BNzBH+fxq4KK4voThDNWqXikpJkiTdxEA2qjcAi0vat1K2ZIN6/wKWbvTZ9lMVxaTT2utT0uah2IRKbtR1gScbKC81RNKQmqavpCWAbYGH4vallDNdgK2AR9p5/iRJkqTFDFjv38ipOpai4fttikLRvymqRVUuBy5USXD+dYrg/emSDgB2qp6rdqDP1YFfhpfwQhQ1pYuoQ9IHKd7HywDvhgrTuhRHp7MlLRzt/2D7imh2HCVLzkEUb+Uvd+H1JAOYDF9Juou+6FjUagasUQWw/SwlI0wjboo6jzDvVuq689On7XOAczowr39QzkzrmQZs2KTNy8D/a6/vJEmSpPvo90ZV0g7AxcA6teTgvTCH8cBI2/vH2etrYWCTpM+x7mHNxR8ePDZXsUkXGADiDwPhTHVX4Daar0hbQoj6t4vt09KgJkmS9E/6tVGVNAjYHNiLMKoh7HCTpAslPSTp3Eo6uMMlTVRJbD6hUr66pKsl3Svp1pARRNJZkk6QdCNwvKRlJV0qaZqkuyTN44Er6UhJh8T1TZKOl3SPpEckbRnlC0v6ScxlmqSv1p5H0vWSJkuaHue8tX6/GP3cJ+lXce6aJEnSN+iOcJo+eEbbr40qRebv6jgXfVElhyqUc8ma889qFMMLcLLtUbbXB5YAPhPlE4Cv294YOAQ4tTLGWsC2tr8JHAVMsT0c+B4dOD8FFonsOgcCR0TZXsArtkcBo4CvSPoQ8Aawg+2NKPlTf6bCOsAuwOYRwzoLGNfBd5QkSdIzDACj2t/PVHcFfh7Xv4/PVwL32H4aINSJhlG2iLcOr90lKeIMD8QqdDPggli4QhFaqHGB7VlxvQUlZyu2b5C0nKTB7czx4vh9b8wDSpL04ZJ2is+DgTUpsag/lPQx4F1gJeADlFjYjSliElC+EDSUKZS0N7A3wCqrrNLO1JIkSZLO0G+NqqTlKHGb60sysDDle81VFGH6GrOARVQywJxKcSh6StKRwOKU1fzLsQJsxL+rwza43953qdpcZjHn7yHKyviaumcaTxGS2Nj225KeiDmKIuL/3XbGwvYEysqbkSNH9sHveUmS9FcGQkhNf97+3Qk4x/aqtofZHgo8TllNNmLx+P1CnMXuBGD7VeBxSTtDyUSjkhS9EbcQ264qKeNeiPad5RpgX5VMNEhaS9JSlBXr82FQtwZWjfrXAzspMtPE2e6qjTpOkiRJuo9+u1KlbPUeV1d2ESWbS70QPrZflnQ6MB14AqhqAI+jiDZ8n5IZ5vfA1AZjHgmcKWka8Bqwx3zO/QzKVvDkcJaaQTkfPhe4XNIk4D5CTcn2gzG3ayUtRMlusx/w5HyOnyRJ0noGwEpV9gB4yqQhI0eO9KRJk9qvmCTJgEfSvZGwY75YbOhQr/yN1sc5P/atb3ZpXq2mP2//JkmSJEmP0p+3f9sk9HV/TglZeZOy5XspsJ3tzzSofxNwiO0uL+0iHvU0yjbtrhQP4vXbbtWwn5HA7rYP6Oqckq7x7j/W6u0ptIyFPpi5GJLW01eTireaAblSjXPKS4CbbK9ue11KXOkHWjhGW+IL44Cfhkfx6/M7hu1JaVCTJEn6DgPSqFKEE96OdG3A7ATitwKDGqktVZG0ayga3S/p+Er5TElHS7ob2FTSxyVNibq/lrSYpC8D/wUcLuncun6bKSntIOm68DxeIdSXPqiiDnVF1Bkt6Y4Y7w5JH+6OF5ckSTLfWK3/6WMMVKO6PkVsoRHN1JYAUEkifjwlBnYEMEol3RvAUsD9tjehpG47C9jF9kcoW+372j4DuAz4lu161aOGSkq2LwH+QfHoPR04IjLZVHkI+JjtDYHDgR929GUkSZL0CL2kqCTpPyU9LOlRSYc2uD9e0oyQeb0vFj+1e6tIulbSXyQ9KGlYW2MNVKPaFvfYftr2u5SwlWF190dRto1n2H6HEubysbg3izn5UT8MPB4SiQBnV+o145PA7qHydDewHEVJCUou1+8Cb9r+XYO2gymqT/cDJwLrNRpA0t6SJkmaNGPGjHamkyRJsmATR3GnAJ+iLJZ2ldQofef5tkfEzxmV8nOAn9heBxhNE7W6GgPVqD5AkfVrxDxqS3X329pveKMiWTg/+xI1JaXaH/ZDtq+NeytRpAk/ELGo9RwD3BgOT59ljpjFXNieYHuk7ZFDhgyZjykmSZLMHzVnpVb+dIDRwKO2H7P9FkVnYPt22pT5FuO7iO0/A9ieafu1ttoMVKN6A7CYpK/UCiSNArbqQNu7ga0kLR/fgHYFbm5Q7yFgmKQ14vNuTepVaaikpJJW7kzgC8BfgIMbtB0M/D2ux3fgOZIkSQYCKwFPVT4/HWX17Bi+LBdKGhplawEvS7o4/FV+0o4T6sAMqbFtleTlP4/99TeYE1LTXttnJX0XuJGysrzK9h8b1HtD0pcoW7KLUBSaTquvV0czJaVvArfavjW2hidKurKu7Y+BsyUdTPnSkPQgGYaSJB2ge0Jqlg+VuRoTQuO8Rkc02S8Hfmf7TUn7UI7rtqHYyC0pvjb/B5xPWbT8b7PJpKLSACYVlZIk6ShdVVRafOWhHrpfo022rvHo9w5uc16SNgWOtP0f8fm7ALZ/1KT+wsCLtgdL+ihwnO0xcW834KO292s2Xq+vVCWtTDlEXpeyHX0FxTP2rRb1P56SeWb/TrT5nu02vWfDA+yK+RFtaNDPZrbPm4+2d9jerCvjJ0k9ww86sbensMAz7cTWy/Et8PSe+MNEYE2VnNR/Bz5POUqbjaQVbD8bH7ejHLPV2r5P0hDbMyir1zZXIr16phpbnBcDl9pek7J/PQg4tjfnRRGC6CmGUfcHbo/ann4a1CRJFih6IaQmojT2p/is/AX4g+0HQlNgu6h2gKQHJE0FDiD8UsLx9BDgeknTKVvJp7c1Xm87Km1D8Zg9E2Y/wEHAnpK+J+nXAJI+EkILS0r6q6QhUb5QxB0tL+ksSadJujXEEapSgytKujra/rhWqAYiDpKOA5aIWKVzo+zgqHO/pAMr/S4i6ezK4faSUf9wFQGH+yVNiC8PSFpDRcRhqqTJklanZNLZMsY7SM0FIMZIulHSeZRMOkiaGb8HSbo++pwuqUOebUmSJAMB21fZXisU9I6NssNtXxbX37W9nu0NbG9t+6FK2z/bHm77I7bHt7eL2ttGdT3qRBgi/+j/UbaB1wiHojOBr4Yr82+JnKXAtsBU2y/E52EUD97/B5ymkngcikjDLsBHgF0kDVUTEQfbhwKvR0jLOEkbA18CNgE+ShFk2DD6/TDlUHw48CrwtSg/2fao2BpeAqgZ+HOBU2xvAGwGPAscSnFCGmH7RJoIQET70cBhIatY5Q1gB9sbUdSiflYz5EmSJH2GXhJ/6El626iKxq+lVj4e+A1ws+3b496vgd3jek+Kwa3xB9vv2v4r8BiwdpRfb/sV228AD1KSe7cl4lBlC+AS2/+2PZOyXb1l3HuqMq/fMicB+taS7o7tgm2A9SQtDawU6kjYfqNJvFNbAhD32H68yfv6oUoe1+so7uINdYyV4g9JkiTdRm8b1QeAuby2JC0DDKUkEl8TmAmsWLtv+yngOUnbUFaPf6o0rzfQtc+NBB06upJrq94848Xq+FRgp5AnPJ0ixNCZ8ZoJQPy7SZtxwBBg4xDpf44Uf0iSpI/RS+IPPUpvG9XrgSUl7Q6zHXB+RtHMXRT4H8rqcTlJO1XanUFZGf6homAEsHOcs65O0e19uI2x2xJxeFshwADcAoyN89ylgB0owvsAq4S7NtH+NuYYsxckDQJ2gtnb2k8rdIJVxPWXBP4FLF2ZV0MBiDaeA4rww/O235a0NWUlniRJkvQwvWpUXYJkd6AYw78Cj1DOB79H0a89NbRz9wKOk/T+aHoZxUv4zLouH6YYxj8B+8R2b7Oxn6Vo6d4ITAUmV0QcJgDTJJ1rezLFyN9DMcRn2J4S9f4C7BHbrssCv7T9MmV1Op0iJjGxMuxuFC+zacAdwAeBacA74bx0EOULw4MUAYj7gV/RfujTucBIlQDocRQ1pyRJkqSHWSDFH1SSc59oe8tK2VmUuNELe21iCxgp/pAkSUdRV8UfVhrqVb/aevGHR45oW/yhp+l18YfOoiIruC9zPICTJEmSpE+wwBlV28dRYjvry8d3pV9JsyhbtosC71C0H3/ukgKuWZsVgV/Y3qlZnd4gAprXjXeVJJ1itZ+f0NtT6Lc8dmDrV2oLDH3UsajVLHBGtRt5PTxnibPb8ygOQEc0a2D7GcIRqS8RAc2X9fY8kiRJ5mIAGNXe9v7tk9h+Htgb2F+FZipHw8KZqJY5/uKOKjdF+X+GCtJUSddH2bKSLo1x7pI0PMoHSToz+pkmacc2+hgv6eSeel9JkiRJIVeqTbD9mEoy8PdTEtq+YnuUpMWA2yVdy7zfu0ZQUgS9CTws6SRKXOzxlKToLwHXRljN7RQv4Y/ZflzSstHHUcAU22MjFvec6Pe/Yw4fAZD0PhW5xkZ9NEXS3pQvDKyyyirz93KSJEnmhwGwUk2j2jY1wYZPAsMrsbKDKcIU9Uk0r7f9CoCkmnLTcoRyU5TXlJtmAbfUFJJsvxh9bAHsGGU3SFpO0mCKJOPnawPZfknSZ5v00ZTIMzgBivdvR19EkiRJ0j5pVJsgaTWK4XueOSpH19TVGVbXrDPKTW1JNNbjJvWb9ZEkSdKnEAPDUSnPVBsQ26qnUYTxzfypHNVoptx0Z5R/KPqsbd3eQoQLSRoDvBBqTNdS0hfV5vi+NvpIkiRJeoFcqc5hiRCxr4XU/AaoxRacQcmAMzmyv8wAxnakU9vPqmSav5HyZe2qmnJTnG9eHGe3zwOfAI4EzgzVpdeAPaKrHwCnhGPULOAo2xc36SNJ5psBHfaRdC8DYKW6QCoqJa0hFZWSJOkoXVVUWmKFof7QXq3/wvaXY/uWolK/2P6VNEslyXft59B26o+RtFkH+h0rad3K55tCIrHHUUnCvlP9PCRdJem9vTGnpH8y7JzjGHZO6oYkyfzQX7Z/Zws3dJAxlJRyd7RTbywlWfqD8zmvbsf2p3t7DkmSJB1iAGyM9ouVajMkPSHpqBBHmC5p7fDY3Qc4KFa1W0paVdL1IapwvaRVYiW7HfCTqLd6dLuzpHskPSJpyxinmTiEovz+GH+XKB8Tq80LJT0k6dw4q0XSxpJulnSvpGskrdCBZ1w+ri+Ndg/EWWuSJEnSg/SXlWrNyajGj2yfH9cv2N5I0teAQ2x/WdJpwEzbPwWQdDlwju2zJe1J0fMdK+kyKplvwu4tYnu0pE9TJAy3paSmayQOsRFFuGEDYHlgoqRbYl4bAusBz1CEIDaXdDdwErC97RlhhI8F9uzge9jT9ouSloixLrL9z868yCRJkm5jAKxU+4tRbWv79+L4fS/wuSZ1Nq3c+w3w4yb16vsbFtfNxCG2AH4XidSfk3QzMAp4FbjH9tMA8YVgGPAysD7w5zDgCwPPtjGXeg6QtENcD405zGVUU1EpSZLeYiDEqfYXo9oWNUGGmhhDR2jrT9+ov2biEG2ddzYTinjA9qYdnGd1rDGUVfOmtl+TdBOweH29VFRKkiTpPvr1mWob/AtYuvL5DuZIAI4DbmtSrxnNxCFuAXaJM9chFHnCe9ro52FgiKRNo59FJa3XwWcaDLwUBnVt4KMdbJckSdIzuBt++hj9xaguURdS0148wOXADjVHJeAA4EshuLAb8I2o93vgW5KmVByVGnEGxUN4cogz/Iqy8rwEmAZMBW4Avm37H806sf0WJZXc8ZKmAvcB7Yb+BFcDi8QzHAPc1cF2SZIkSYtI8YcBTIo/JEnSUbos/vDBoV5tj9aLPzz4474l/jAQzlSTJEmSPkA6Ki3gSJoFTK8UjaWEtuxu+wBJR1IJrZmP/odRQm7W7+I8hwGb2T5vPtreYbujW8RJ0i7rHH5iy/v8y9EHtbzPJOmL9GujSuNQmyeADu95Slo4QmK6k2HAF4AOG9XavNKgJkmywDAAVqr9xVGpw4Sa0RWVog0k3SDpr5K+Uqlzo6TziJWupINDGel+SQdW2i8i6exQUrpQ0pJR//BQWLpf0oSKYtIakq6TNFVF6Wl14Dhgy3CcOqgNhaZG85oZvwepqEHV1KO27943mSRJktTT31eqVaWlx23v0KDOcEr4yVLAFElXRvloYH3bj0vaGPgSsAkllvTuEHJ4CfgwsJft2yX9Gvga8FNKLtajAST9BvgMxev4XOA425dIWpzyxeZQitrTZ6L+3jRWaJprXnXP8Qawg+1XQ7bwLkmXuc4TLcUfkiTpLQbCmWp/X6m+bntE/DQyqAB/tP267RcoOU9HR/k9FcO1BXCJ7X/bnklRVdoy7j1l+/a4/m3UBdha0t2SpgPbAOtJWhpYyfYlALbfsP1agzl9Etg9vhDcDSxHUUeqn1cVAT+MkJrrgJWAD9RXsj3B9kjbI4cMGdLklSRJkiTzQ39fqXaE+u9Otc//rpSpM+1jBXoqMNL2U+EQtXg7/VRpptA0pm5eVcYBQ4CNbb8t6QkaKColSZL0GrlSHRBsL2lxSctRUsJNbFDnFmCspCVDKWkH4Na4t0pNAQnYlaLGVDNmL0gaRBF0wParwNOSxgJIWizOYOuVm5opNLXFYOD5MKhbA6t28PmTJEm6n+5QU+qDRjpXqkU28EpgFeAY289IWqtawfZkSWcxR2LwDNtTIhTmL8Aekn4F/BX4ZUgFnk5xJnqCuQ31bsCvJB0NvA3sTFFdeidUlM4C/ofiETw5HJxmUMKB2uJc4HJJkyhKTA917jUkSSHDX5Jk/klFpQFMKiolSdJRuqqotOQHhnqNca1XVJp+YioqDVgqYhSiZKbZ3/Yd89HPPsBrts+JFfQVti+MzDSH2J4k6SrgC7Zfbt0TJAOBrT95fG9PoUe48drv9PYUkn5IGtWeZbYYhaT/AH4EbNXZTmyf1oE6baWdS5Ik6XkGwMZoOir1HstQ4lzbFG6QtHsIQEyNeFckHSnpkLY6l/RExKsmSZL0CeTW//Q1cqXas9TEKBYHVqDEr0IT4QZgXeAwYHPbL0hatldmnSRJknSINKo9S3X7d1PgHEnrM0e44WPAu8wRbtgGuDCEKbD9YlcnkIpKSZL0Gn1wZdlqcvu3l7B9JyVjzhDmFm4YATzHHLGIlv5nmIpKSZIk3Uca1V5C0trAwsA/aS7ccD3wXyFMQW7/JkmyQJPiD0mLqQr8C9jD9ixJDYUbbD8g6Vjg5gjHmQKM74V5J0mSdI0+6ljUatKo9iC2F25S/gKwaZN7ZwNn15UdWbkeX7keU7ke1pW5JgOXjN9MkvknjWqSJEnSM+RKNZFk4ATb34zPhwCDbB9Zp2w0HrjW9jPdOJcVgV/Y3qm7xkgGFrvcuc88Zedv2q62SJIkTUhHpfZ5E/hcIyEF26fZPic+jgdWbNSBpIbbvp3F9jONDKqk/HKUJEmfZyCIP6RRbZ93gAnAPKk7aspGknYCRgLnSrpP0hKhaHS4pNuAnSWNkHRXqCNdIul90ccBkh6M8t9X+v2NpBsk/VXSV6J8mKT743q8pAskXQ5c25YqU5IkSdIz5AqnY5wCTJP040Y3Q8x+f0LMHqBkbOMN21vE52mUxOM3R9q3I4ADgUOBD9l+U9J7K90OBz4KLAVMkXRlg6E3BYbbfjFWq/OoMrkuDVGKPyRJ0mv0wZVlq8mVageI5OLnAAd0sun5AJIGA++1fXOUnw18LK6nUVa4X6Ssimv80fbr4Rl8IzC6Qf9/rqgs1VSZpgHXMUeVqf5ZUvwhSZJeIbd/kyo/B/airBw7yr87UOf/UVbCGwP3Vs5H6/9zafSfT7X/ZqpMSZIkSQ+RRrWDxIrwDxTD2oh/AUs3afsK8JKkLaNoN4qgw0LAUNs3At8G3gsMijrbS1o81JTGABPbmWIzVaYkSZLepzvUlDq4UpX0n5IelvSopEMb3B8vaUb4xNwn6ct195eR9HdJJ7c3Vp6pdo6fAfs3uXcWcJqk12ks5LBH3F8SeAz4EkWm8LexPSzgRNsvx3nsPcCVwCrAMbafkTSsjbk1VGVKkrbI8JmkvxPRF6cAnwCeBiaGv8mDdVXPt93s3/djgJub3JuLNKrtYHtQ5fo5YMnK5yMr1xcBF1WaDqvr5z6K41E9WzQZ+hHbe9f18QSwflyfRTHktXtNVZmSJEn6BL1zBjoaeNT2YwARZbE9UG9UGyJpY4p/ytWUKI82SaPaQ0g6DPgCMIuS3u2rwPPA74FlgcnAbrbfomzdUJ6vCgAAIABJREFUfkLSccDnbV8YfawKXExZ4S4KnGT7tLpxLgNWs71+jzxY0u/Y+E+H9fYU5uLeTx3b21NIWoDoNsei5WOHrsYE2xMqn1cCnqp8fhrYpEE/O0b6zUeAg2w/FUd0P6Mc2X28I5PJM9UeIHKnfgbYyPZwYFvKH/l4ypbvmsBLzDmvPQr4NHBeXVfPApuFI9ImwKGhslQb53PAzO58liRJkj7GC7WIhviZUHdfDdrUm/fLgWHx7/N1zNFb/xpwle2n6CBpVHuGFSh/+Ddh9lbts0QS8qhzNjA27j9hexplRTsb22/V+gAWo/L3kzQIOBj4QTc+R5IkyfzTO45KTwNDK59XBuaSk7X9z8q/radTojGgHKntL+kJ4KfA7rGD2JQ0qj3DtcBQSY9IOlXSVsBywMu2a7GpT1O2KdpE0tCIRX0KOL6iNXwMZZvitdZPP0mSZIFlIrCmpA9Jeg/weeCyagVJK1Q+bgf8BcD2ONurRNavQ4BzbM/jPVwlz1R7ANsz47B7S2BriijEzxpV7UBfTwHDY9v3UkkXUlbCa9g+qB0P4VRUSpKk15B73lPJ9juheHcNxR/l15Gr+mhgku3LgAMkbUcR4HmRLuStTqPaQ9ieBdwE3CRpOiXE5r2SFonV6jxbEu3094ykByiGegiwcWxRLAK8X9JN1fyqlXYTKFrGjBw5sg/qkSRJ0i/pRFxpy4e2rwKuqis7vHL9XeC77fRxFpWIi2bk9m8PIOnDktasFI0AnqTID9ayzuwB/LGdflaWtERcvw/YHHjY9i9trxhbFFtQwnHGtPYpkiRJkvbIlWrPMAg4KQTz3wEepWzBLgP8XtIPgCnA/wJIGgVcArwP+Kyko2yvB6wD/CxyvAr4qe3pPf40SZIk80Ff1OptNWlUewDb9wKbNbj1Ag2E8m1PpGwH15f/mZK9pq2xniAEIpJkfsi40CSZf9KothBJs4DplFXkLGB/23eE89BfKNKBi1N0gk+xfXa0Gwd8J7qZCexre2qD/tcGzgQ2Ag6z/dPKvV9TYmGfT+GHgcVh0z7X0v6OHX5xS/tLktnkSjXpJK+HMAOS/gP4EbBV3Pub7Q3j3mrAxZIWsn0m8Diwle2XJH2K4kjUSPHjRUr6ubEN7p0FnExJUZckSdLnGAjbv+mo1H0sQ1FJmofQoDyYyM9q+w7btbp30WDrN+o9H1vDbze4dwvF6CZJkiS9RK5UW8sSku6jbPGuQFFMasZkYO0G5XsBf+qGuSVJkvQuA2Clmka1tVS3fzcFzpHU7HxzHj3KyIO6F80z13SZFH9IkiTpPnL7t5uwfSewPEWYoREbElJYAJKGA2cA29v+Z5TtV0mau2KTfjo7rwk14ekhQ5pNLUmSpMW4nKm2+qevkSvVbiI8dRcG/kklB2vcG0YRZz4pPq9CSem2m+1HavVsn0JJrpskSZIsAKRRbS21M1Uo27t72J4lCWB1SVOYE1JzUnj+AhxOEdg/Neq+Y3ueZLiSPghMojhBvSvpQGBd269K+h0whpJb8GngCNv/210PmiRJ0mn64Mqy1aRRbSG2F25S/gSwRBvtvgx8uQP9/4PmnsG7dmyWSX8j40qTBYFuTFLep8gz1SRJkiRpES1dqUr6AHAi8FFKjOZbwI9tX9LKceZjXgcCE2y/Fp+vAr5g++Ventd4YKTt/duoMwZ4y/YdPTWvpP9yx5OrtVtns1Uf64GZJAOSXkj91tO0bKWqchh4KXCL7dVsb0xJBttwu7KHOZCKs5DtT/e2Qe0EY2isG5wkSZL0MVq5/bsNZUV1Wq3A9pO2T5K0sKSfSJooaZqkr0JZhUm6WdIfJD0i6ThJ4yTdI2m6pNWj3lmSfinpRkmPSdpK0q8l/UXSWbXxos4kSQ9IOirKDgBWBG6UdGOUPSFpeUnDJD0k6QxJ90s6V9K2km6X9FdJo6P+spIujbnfFeEvSDoy5nFTzOuAylwulXRvzGXvSvmX4llvpqRuq5UPkXRRvKOJkjYPL+F9gIMirGbLRvWi/VaV8JspkpZu4d82SZKky2RITedYj6IS1Ii9gFdsj5K0GHC7pGvj3gaUlGYvAo8BZ9geLekbwNcpq0woadC2AbYDLqcYpC8DEyWNsH0fRWT+RUkLA9dLGm77F5IOBra2/UKDua0B7EwRRJgIfIEivrAd8D2Kzu5RwBTbYyVtQ9HXHRHt1wa2BpYGHpb0S9tvA3vGXJaIOV4EvCf62hh4hZJPdUr08z/AibZvixCba2yvI+k0YGZNPF/SefX14v0dAuxn+3ZJg4A3mvwtkiRJep5eTFLek3Sb96+kUyjG6S1KQu7hkmoJuQcDa8a9ibafjTZ/A2rGdjrFWNW43LYlTQeeq+URlfQAMAy4D/ivWBUuQpEJXBeY1s5UH6/r6/rKOMOizhbAjgC2b5C0nKTBce9K228Cb0p6HvgA8DRwgKQdos7QeN4PAjfZnhHjnQ+sFXW2BdaNkBqAZZqsNpvVux04QdK5wMW2n270sEpFpSRJkm6jlUb1AcLwANjeT9LylLjK/wO+bvuaaoNwwnmzUvRu5fO7dfN7s0Gd2fUkfYiyWhsV2V7OosSEtkdHxp9HUpA537mq7WfFXMZQjN+mtl+TdFNlLs2+qy0U9V+vFlaMZ5v1gOMkXQl8GrhL0ra2H5pn0vYEShYcRo4cOQC+NyZJ0lfQu709g+6nlWeqNwCLS9q3UlZzDroG2FfSogCS1pK0VAvHhiKI8G/gFRUv5E9V7v2Lsj07v9wCjIPZXwResP1qG/UHAy+FQV2b4g0NcDcwJla6i1K2nWtcC8z2ApZU216un3vDepJWtz3d9vGULzKNxPqTJEmSbqRlK9XYMh0LnCjp28AMipH7DnABZSt1ssrSawaNc4J2ZfypKopFD1DOZm+v3J4A/EnSs7a3bthB2xwJnClpGvAasEc79a8G9on6D1PSuWH7WUlHAncCz1LOoGuCEQcAp0SbRSiGfB/K+fGFkrannDE3q3egiiD/LOBBMtNN0oAMl0l6lQGwNyYPgLihpDEjR470pEmTensaSZIsAEi6t5F8akcZtOxQb/DxA9uv2EnuuPCQLs2r1fQbmUJJM20P6sb+zwBOsP1gB+uPBHa3fUC7ldvv6wmKSMQLdeX7AK/ZPqerYyRJjZMfaisNcM+x/9o39PYUkqTT9Buj2t2EPm9n6k+inG12CEkL257VyTFOa79WkiRJH8CkotKCjqRVJV0fog3XR1xnTUziF5LuCNGGnaJ8IUmnhmDDFZKuqty7KVafSJop6fgQd7hO0uiKAMR2UWeMpCviepCkM1UELaZJ2rHSz9GS7gY2lfTxEG6YriIqsVjlcb6lIopxj6Q1ov2Rkg6J66+oiEFMVRGHmCvdXJIkSdL99GujCpwMnGN7OHAu8IvKvRUo8aefAY6Lss9RHKo+QhGW2LRJv0tR4k03pnjn/gD4BLADcHSD+v9NEb/4SMzlhko/99vehLKqPQvYxfZHKLsIVU/qV22Pjmf6eYMxLrY9yvYGlOTnezWZe5IkSa8wEBSV+rtR3RQ4L65/QzGiNS61/W6ckX4gyrYALojyf1AUjxrxFsXDF4pIxc2holQVjKiyLZVk47ZfistZwEVx/WGKEEUtSfnZwMcqffyu8ruRsV9f0q0hWjGOonA1D5L2VpFynDRjxowmj5ckSZLMD/3dqNZT/V5TFW1Q3e/2eNtz3KZnC0bYrhesqPbf6DvVG5Vz1PbGdpPrGmcB+8cq9yiaCF/YnmB7pO2RQ4YMaWfIJEmSFuJu+Olj9HejegclUw6U1dtt7dS/DdgxzlY/QMkQ0wrqBRve16DOQ8Cw2nkpsBtwc+X+LpXfdzZovzTwbIhKjOvyjJMkSVpILUl5f9/+7U/ev0tKqurdnkARSvi1pG9RBCe+1E4fFwEfB+4HHqEoIL3Sgrn9gCLYcD9ly/co4OJqBdtvSPoScIGkRSji/lXv3sXCoWkhYNcGY/x3zPdJyjZ0ZqlJkiTpYVL8oQ5Jg2zPlLQccA+weZyv9jtS/CFJko7SVfGHpd+7skeM+UYrpwTAbX/8doo/9HGukPReSpq2Y/qrQU2SJElaT58yqpIOo+QznUVxAPqq7btb0O9NwCEhyNAmtsc06eMO25t1dS5dJbSDZ9r+qaSjgVtsX9fL00r6EZ9a6eu9PYVe509/P6m3p9Av6YtnoK2mzxhVSZtSYkY3sv1mpI17Ty9Pa7bSUV8wqPXYPry355AkSdJhBoBR7UvevytQUqrVwlNesP2MpMNDKeh+SRMiyw2S1gg1o6mSJktaPcq/HYpEUyUdV+l/51AjekTSllF3WMR2To6fzaJ8jKQbJZ1HcfpB0szKvZsl/SH6Ok7SuOh7emUenVVzGhT1Jkc/29cmLukwSQ9Luo4Sz0qlr1r7hu8pSZIk6Tn6klG9FhgahupUSVtF+cmhFLQ+sARlNQtFIemUUBDajBJO8ilKSrlNovzHlf4XCUWiA4Ejoux54BO2N6KEqlQVl0YDh9let8FcNwC+QVFe2g1YK/o+g5KeDTqv5vQGsEPMZWvgZypsTAkL2pCi+DSqyftr9p6SJEn6BAMhpKbPGFXbM4GNgb0p4S/nSxoPbC3p7lAK2gZYT9LSwEq2L4m2b9h+jaJcdGZcY/vFyhC1EJZ7maN6tChwevR9AVA1oPfYfrzJdCfafjZW1X+jfCGAuRWVOqvmJOCHKnlSrwNWintbApfYfi0So1/WZE7zvKdGlVJRKUmSpPvoM2eqAKEudBNwUxiHrwLDKWnPngonncVprj7UTLkI5igozWLOcx8EPEdZeS5EWS3W+HcbU62qMb1b+dxMUYm6eTVScxoHDAE2tv22Srq3xRu0nQdJiwOnMu97mncS9gRK0nZGjhzZB7/nJUnSLzHwbv//J6fPrFQlfVjSmpWiEcDDcf2CpEHATgCxYnta0thou1hkZbkW2LOWoUXSsu0MOxh4NuQFdwMWbtkDdV7NaTDwfBjUrYFVo/wWYAdJS8QK/bMN2tYM6FzvKUmSpE8xAGQK+9JKdRBwUsSIvgM8StkKfpmyrfoERWWoxm7AryKs5G1gZ9tXSxoBTJL0FnAV8L02xjwVuEjSzhTx/LZWp52ls2pO5wKXS5oE3EeRLcT2ZEnnR9mTwK31DW2/LOl0Gr+nJOkUGU6SJPNPKioNYFJRKUmSjtJlRaXBK3vjzQ5o5ZQAuPnq76SiUiuRNIuyQluEkkd0D+D9wBXhCduVvkcCu9vu9H8JkmbaHtSV8ZOkIww757j2K3WCJ3Y/tKX9JclAos+cqXaB122PCAP6FrBPqzq2PWl+DGp3oSK0nyRJsmBit/6nj9EfjGqVW4Fa6rSFJZ0u6QFJ14ajz+qSJtcqS1pT0r1xPSoEGaaGkMPSIfRwRdzfStJ98TOlcv8WSZdIelDSaZIWqvR/bPR3l0oqOSR9NkJfpqiIV9TKl5L06xBwmFITf5A0XtIFki4Hrm1HJOKLMff7JP1KUisdr5IkSZJ26DdGNVZxnyIUkIA1KeIQ61GcnXa0/TfglXBmguI8dJak9wDnA98I0YhtgdfrhjgE2M/2CErsaO3+aOCbFCGI1SkCDQBLAXdFf7cAX4ny24CP2t4Q+D3w7Sg/DLjB9iiK+MNPJC0V9zYF9rC9Dc1FItahCFhsHnOcReZVTZKkDzEQxB/6w3biEpLui+tbgf8FVgQet10rrwo+nAF8SdLBFCM0miL996ztiTA7ZIc6pb/bgRMknQtcbPvpuH+P7cei/u8oIg8XUrair6iM/4m4XpkibLECRdu4JjDxSWA7SYfE58WBVeL6zxUhi5pIxMcocbE1kYiPU8QzJsa8lqAoRs2FpL0pXtWsssoq9beTJEm6hz4aAtNq+oNRfT1WZrMJo1IVWJhFMTJQEpEfAdwA3Gv7n5JWop0/t+3jJF0JfBq4S9K2tVv1VeP3257jWl0VnDgJOMH2ZZLGAEfWpk1ZTddic2vPsglzh/o0E4kQcLbt77bzHCn+kCRJ0k30m+3fjmL7DeAa4JfAmVH8ELCipFEAcV461xcOSavbnm77eGASsHbcGi3pQ3GWugsdE3n4e1zvUSm/Bvi6NDthwIZttG8kEnE9sJOk90f7ZSWt2qSPJEmSHkWA7Jb/9DUGnFENzqWsKK8FsP0WxSCeJGkq8Gfmlfk7UCUDzFTKeeqfovxOiij+/ZSt3EvaGftI4AJJtwIvVMqPoWgRT5N0f3xuNveRIRIxjjkiEQ8C36c4M02LZ1ihnbkkSZIkLWRAij/EueVg2//dxX7GUJKfL5AZYVL8IUmSjtJV8YdlllnZI0ft38opAXDjDd/tU+IPA26lKukSYHfgf3p7LkmSJAOJ3tr+lfSfKjmpH5U0j7pJhC7OqIRNfjnKR0i6M0Izp0napb2xet1RqZEiku3XWq1IJGkYsJntHVrVp+2bKFl15ovKnM5rp2qz9kcDt9i+bn7nkHSe7W9r/bftvsQftzi5t6eQJC0j4vVPoURgPE2JkLgsjsyqnG+7/n/u1yiqen+VtCJwr6RrbL/cbLy+sFLtNkWkOoYBX+imvueXYXRhTrYPT4OaJMkCQXdkqOnYQnU08Kjtx8J/5vfA9u20KVO2H7H917h+hhKmOKStNn3BqFapKiIB0ExBSNIwSX9RnWpS3Fsj1IqmRrvVKc5EW8bS/iBJi0s6M/qcEp60tW2ASyVdLulxSftLOjjq3BVetZ1VZlpY0k9U1JKmSfpqNK2f0zBJt8acJ0varDLGt2OuUyUdF2VnSdoprg+P/u+XNKHmRZwkSTLAWQl4qvL56SirZ8f49/lCSUPrb0oaTdEW+Ftbg/UZo6p5FZFqNFQQinvzqCZF+blRvgGwGfAscChwa6yKTwT2A7D9EWBX4GyVZN8A61NWkKOBY4HXQgHpTspWQGeVmfYCXgm1pFHAVyR9qMGcngc+Ec+6C/CLeDefAsYCm0S/P27wCk+2PSpW/EsAC6TzVJIk/ZVu0P0tZ6rLS5pU+dm7buBGC4z6Ne7lwDDbw4HrgLPn6qCI9fwG+FLk325Kr5+p0lgRqUozBSFooJqkksh7JduXwOy41Hp1JCjKRydFnYckPQmsFfdutP0v4F+SXqG8cCgGf3hcd0aZ6ZPA8NqqkhJruiZlu7vKosDJYaxnVeazLXCm7dei3xeZl60lfRtYElgWeKAy79koFZWSJOkluklW8IV2vH+fBqorz5WBZ6oVbP+z8vF04PjaB0nLAFcC37d9V3uT6QtGdR5FpDqaKQhBY9Wkjm57tlWv2u+7lc/vMueddUaZScDXbV8zV2EJyalyEPAcsAFlF+GNSvum/znGCvtUYKTtpyQdybxxtkAqKiVJMuCYCKwZu4N/Bz5PnS+LpBVsPxsft6M4zRK7j5cA59i+oCOD9Znt3zZopiDUkFgdPi1pLICkxSQtCfwLWLpS9RZCcF7SWhSd3YfpIJ1UZroG2FfSorXxVMTy6+c0mLLSfRfYDahlmbkW2DOeA0nL1k2nZkBfkDQI2IkkSZK+Ri+kfrP9DrA/5d/hvwB/sP2ApKMlbRfVDgjfnKnAAcD4KP8v4GPAeM0Jt2lrEdgnVqrtcS5wuYqC0H2EglA77Ab8SiXk5G1gZ2Aa8E68tLMoK7vTJE0H3gHG236zk/4951Ky0sxWZlKJYzopnKZep2zdnkHx9J0c58EzKGekjeZ0kaSdgRsJzV/bV8cfcpKkt4CrgO/VJmH7ZUmnU7ann6B8M0u6iQw5SZIFC9tXUf7drJYdXrn+LjCPbrrt3wK/7cxYA1JRqVWoRcpMvUUqKiVJ0lHUVUWlQSt5kxFfa+WUALju9u/3KUWlDq1UVRJpnwh8FHiJ4mDz45oz0IJKnGm+ZfuO+Wh7CcWz+JQWzGMs8EiDYOT22m0HrGv7uK7OIVkw2XC/E1ve55RTDmp5n0kyUGj3TDW2Ky+lKPesZntjykHvyt09uVagumwzdYyhGMZOE8pMv6QobnSVscC6nWkgaRHbl6VBTZJkgaEXzlR7mo44Km1DWc2dViuw/aTtk5qJFUgaI+lmSX+Q9Iik4ySNCzGE6SpiDDXxgl9KulHSY5K2kvRrFVGHs2rjRZ1JcZB8VKX8CUnHR7/3SFqj0u8Jkm4EjlcRbLg0AnvvkjRcRSJwH+CgOHzeUtIQSRepiChMlLR59DdP+8r72UDSDZL+KukrUb+hYEXc2z36mSrpN/HOtgN+EvNYPX6ulnRvvN+1mzzXeEknx73PSrpbRaTiuthdSJIk6Tv0jqJSj9KR7d/1gMlN7tXECt6QtCbwO6C2t70BsA7wIvAYcIbt0ZK+AXwdODDqvY9iuLejxFVuDnyZos84IuJQD7P9ooqG4/WShtueFu1fjX53B37OHNGDtYBtbc+SdBIwxfZYSdtQ3KNHSDoNmGn7pwCSzgNOtH2bpFUo3mLrAEfVtwdqHmDDKdviSwFTVBKZP08RrHhV0vKUpOaXUVajhwGb235B0rLxXJcBV9i+MOZxPbBP6E1uQnFg2qbBc42v/C1uAz5q2ypi0N8Gvtnk75YkSZJ0A532/pV0CkU44S2KZ2sjsQKAibW4H0l/IzxkKR6qW1fqXR6GYDrwnO3p0eYBisfsfcB/qYgWLELJEbouxXMWiiGv/a4eMF1ge1Zcb0GoLdm+QdJykgY3eLxtgXU1xwN4GRUxibba/9H268DrsYIcTQkUbiRYsQ1woe0Xoq95RBxUQmI2o+RcrRUv1uS5qqwMnK+i/PEeSm7XeVCKPyRJ0kv0xaTiraYjRvUB5sj/YXu/WH1NorlYAXRMQIG68vo2i6gE7B4CjLL9UmwLV4UN3OT635XrjshUEc+waRjJOY3VMM7Gdb+r5c0EK9oUcajM4eU2BDH+3aT8JOAE25epOGAd2ahSij8kSZJ0Hx05U70BWFzSvpWyJeN3M7GCVrIMxZC8EueEn6q7v0vl951N+qgKPYyhyFq9yrziC9dSgoSJuiPaaQ+wvYo4/3IUx6eJNBesuJ6y6l4u+qqJOMyeR/T7uEqsKips0OzlVBhMUQsB2KMD9ZMkSXqWAeCo1O5KNbZmxwInqmjLzqAYue9QzlrnEStoJbanSppCWTE/BtxeV2UxSXdTviDs2qSbI4EzJU2jeOvWjM7lwIXhSPR1ipLGKVFvEYox3aeN9gD3ULZ7VwGOsf2MpIaCFaHicSxws0oe2SkU5Y7fA6dLOoCihjQO+KWk71P0gH8PTG3nVR1J2TL+O3AX8KF26idJkvQcpuw/9nMWaPGH2FYdWTujTDpHij8kSdJR1FXxh6VW8kfX+2r7FTvJnyceseCJPyRJkiRJVxBOR6X5JbY2p1O2Lt+h5Kb7eVt56CJu9Arb68dZ5oqh19gU28NaMNfv2f5hF9rPtD2oi3N4L/AF26fOR9urou3LXZlDktRY/YQTensK3crfDj64t6eQ9GO6K0vN65F4ez3gE8CnKWnSOsqIaNNlwtGnref8Xhv3eor3Ap0Sxaw9l+1Pp0FNkmSBYAA4KnV76jfbz1PiIvcPQ7CwpJ+oKBZNkzTXJrtK/rqjgV1CYWgXSaMl3RFqQXdI+nDUHS/pjyrqQw9LOiLKh6moMp1KcaYaKmlXFXWj+yUdH/WOI5Kkh3MRkr6oos50n6RfqQhO0Kh9Zc4/U1FPul7SkCj7SjzjVBWVplratg9IuiTKp6ooKh0HrB5j/iTqfavyjo5q47meiBAnVFSf7lVRntq7xX/KJEmSrpFGtTXYfizGej+wF/CK7VHAKOArEYtaq/sWcDhwfqx2z6d4z37M9oZxr7pdO5riLTsC2FlS7cD6wxTlpA0p6d+Op4gvjABGSRpr+1DmrKrHSVqHEpqzecSJzgLGSVqxUfsYZylgsu2NgJuZsyK/2PYo2xtQcvjtFeW/AG6O8o0oXs2HAn+LeXxL0ieBNePZRgAbqwhJzPVctp+se9V7hjbzSEp+wOXq/xaS9laRfJw0Y8aM+ttJkiRJF+hJR6WagMIngeGSaom0B1MMyCNttB0MnK0ihWjKWW2NP9v+J4CkiynqR5cCT9q+K+qMAm6yPSPqnUtJPHtp3TgfBzamSCQCLEGRHGyr/bvA+dH+t8DFcb2+pB9QtnYHUSQPoRjm3QFCGekVSe+rm8cn42dKfB4U7+j/6p6rngMk7RDXQ6PNP6sVUvwhSZJeYYCE1PSIUZW0GmXV9zzFuH7d9jV1dYa10cUxwI22d4h6N1XuNVI0gvYVlRpOFTg7EtZW5za2Sf1G1MY/CxgbcbbjKcIQHUXAj2z/qm4ew2gSC6wiSrEtRRHqNUk3MbfyVJIkSdLNdPv2b5wxngac7BIUew2wr6RF4/5akpaqa1avdFRVCxpfV/cTKllklqCkUKsXhwC4G9hK0vJxRrorZasW4O3aXCiKRztJen/MbVlJq7bTfiGKYAPAFyjC9sT8n42+x1Xmcj2wb/S/sKRlGjzvNcCeKjrASFqpNqc2GAy8FAZ1bYrIf5IkSZ9Bdst/+hrdtVJdQtJ9zAmp+Q1Q89M/gyKUP1llj3UGxRhWuRE4NPr4EfBjyvbvwRTZxCq3Rf9rAOfZnlS/6rX9rKTvRr8CrrL9x7g9AZgmaXKcq34fuFbFY/htYD/bd7XR/t/AepLuBV5hjmzif1OM8ZOU8KKa0fwGMEHSXpTV+76275R0u6T7gT/Fueo6wJ2xDT0T+GLUb8bVwD4qqk8PU1SVkqTTZMhJksw/C7qi0niKotL+7dVN5iUVlZIk6SjqoqLS4CVX9KZr7tV+xU5yzbQf9G9FpXCUuRhYx/ZDqog6tNGm3TqtJsbczPZ589l+DHCI7c+0V7edfjokdNGg3YrAL2zv1G7lJOkEax9xYvuVWsRDRx3UY2MlvU3fDIFpNd1xprorZUv2893Q91zYPqutVaqktr40DKOcgfY2nRa6kLSI7WfSoCZJkvQtWmq11yh1AAAUo0lEQVRUw7Fmc0pM5jxGVU3EGoKFJZ0ewgXXhuNRWyIKZ0k6TdKtkh6R9JnKGBdIupxyNioVsYn7VcQbameexwFbhuDCQWoiStFGeyhJzC+R9GDMZaFo88uIBX1AIdwQ5aNUxCumqghMDGZeoYulJP065jFFJYNOo+caFmewNVGIW1UEKCarCEokSZL0HcyAEH9o9fbvWOBq249IelHSRsCLdXVGA+tTUqhNlHQl8AIlpnJX21+R9AdKYvTfUkQUTgdQifvci5KQG8pqcytgdeBGSWtE+abAcNsvStqRshrcAFg+xryFIrgwe/tWRYHoFdujJC0G3C7pWopAQ6P2tWdZl+KMdDXwOeBC4LAYe2HgeknDKQIW5wO72J4YXr+vUcQsZp8LS/ohcIPtPVU0ge+RdF2D5xpWeafPA5+w/YZKLO/vKAIQSZIkSQ/S6u3fXSm5P4nfjfKb/tn2P22/Tjl73SLKH7d9X1zfSzGYUEQUbpU0nRKasl6lrz/Yftf2Xym5VteujFEz5lsAv7M9y/ZzlFCYUQ3m9Ulg9/A4vhtYjmLo22p/j+3HQsThd5Vn+S9JkyniDetRDO+HKQndJ0JJRm77nSbzqHk+30SJNV2lwXNVWZSSj3U6cEGM1xClolKSJL3Fu93w08do2UpVRRJvG4oRNLAwZcFfn3mlmVjDm5WyWRQ1I2hbRKHVwg+NRCnaOu+cZ3wVycVDgFG2X5J0FsUwqkH9ZvPY0fbDdfPYhOZJ4A8CnqOsphcC3mg64VRUSpKkl+iLcaWtppUr1Z0omrSr2h5meyjwOLByXb2OiDVUaSaiAEXrdyHp/7d37tFSVucZ/z0C3qIYLEhARdCISqQRuSimKgoammYZiVixaqA1UnRFa2xszIrxmq6I8UrQZZBSvETi8oJaNQgaUWIAUREQL4iECEoUvKU0UVHe/rHf8XwOczjDOd+ZmcO8v7W+Nd/s2Zdnvplz9rz72/vZ2gfYm7Q+s5gnSfcs2ykZURwBPE1pw4VSphSNlQcYJKmX30s9iTRBqyOp8/tAUlfg7z3vy0B3SQO9/p2VJlKV0nG2lBaoSurXxPWBZPywxrfWO430gyYIgiCoMHneUz2ZNPknyz1surVak2YNRTRmogCpE30C6AqM83uKxeWnk+5FLiJFiv9hZn+S9A7wiaRFpGj4ekqbUjRWfn9grr/nvqTOd7qZbZS0kGSUvwL/0WBmH/skp1/4D4q/kmwFi40uLgeuIxlSCFgJNLVs50bgHkknen2NRbRBEATVow4i1YqaPyhHswYfVn3QzO5uaV31Spg/BEFQLmqp+cMO3eywnmNyVJSY8fIVW7f5QxAEQRBsggEbt/5ItSL7qRZoyqxhC+sa05woNbu+M5N2iaQfNKOu2WrYv7VFSBogaYKfbyfp0cza1cmSGp3RW6KuIZIezENXEARBPrTCGtUaHE6OSLVGMLNngMJYbD+gg2+UDg37tQZBEAQ1TEUj1VrHI8/x7na0TNLhnt5O0lVKjkqLJZ1domxjLkpXKDkuLZZ0laedqOTQtKhgJFGILpW2eLsdOMgj1X2yEbGkYyXNdeeku9SwPdxwSS9L+h3JhCIIgqC2iEi1LmlvZoN8ferFpBm6Y4FeQD8z+0TSriXKlXJRWg2MAPY3M3OHJEguSl83szcyaQCY2duSvsvn3Z7wx87AhcAwM/s/ST8EzpN0JXAzaZ3wcjYT2So5R40F6NGjR2PZgiAIgmZQj5FqYz9tCun3+mPW1WkYcFPBAakRV6NSLkp/JhkxTJb0bZItIaRlNlMlncGWrSk91Ot9ypfgjAb2IjlJ/cHMXvWN4G9vrAIzm2RmA8xsQJcuXbag6SAIghYSkepWyTtAp6K0XUlGFdDg7PQpDddns25IjbkoeVQ7CBhK2mDge8DRZjbOHZL+AXheafu3chDJqvBz9o9evva+XUEQBHVG3UWqZrae5NA0FMCHcoeTTCkaYyYwzh2QKDH8W9JFye937uJ7pZ5LMuZH0j5mNt/MLiJtJrBnmfLnAV+TbxwgaUdJvUluTb3cWQpKey4HQRBUj8KSmryPGqMeI1WA7wA3SLran19qZq+VcGMqMBnoTXI52kC6fzmx8KL7Em/iokRyf7pfUsH7t7Aj88+VdpMR8BjJrenIpkSb2Vo30JimtJMOwIW+K9BY4CFJ60g/ECq24XsQBEHTGFgNOuDnTEUdlYLaIhyVgiAolxY7Km3X1Q7rXmzf3nJmrLw2HJUKSJoN/Cy7M4ykc4HeZnZW1YSViS9z+Y6ZnZNDXStJFo7ritLHAX8xs1tb2kYQlMOhp15TbQkVZd7t51VbQv1QB0FctYd/p5Em8GS3WxsFnN9UQTebl+/MUkhr53ubVoQiw4YmaY4+M7tpi4UFQRAEVaHaE5XuBr5ZuD/ou9V0B34n6XxJC9w04dLC65JeknQj8Bywp6T1ki6TNB8YLGmlr+csWP/N9vOdJP13xsDhBE8/2dNekDS+IMzrHS/pWbcMHOQmDCskHed5PrMD3Ez9xfqGSlro+aZk7o0CnK9kPPF0ZjLSZxaKks7wa7JI0j2SdvT0Lv58gR9fa4XPKgiCoPnUyUSlqnaqZvYOaW/S4Z40imRccAywLzCINGO2v6QjPM9+pH1b+5nZH4EvAC+Y2SFmtrkZvD8BPjCzvmb2t8BvJXUHxpNMEw4CBko63vN/AZhtZv1Je57+1HWNAC4rp/5MPS+Y2SGkqHYqcJKZ9SWNFJyZqePPZjaINAnquhJt3GtmA83sq8BLwOmefj1wrZkNBE4gTawKgiCoLepgnWq1I1VoGALGH6cBx/qxkBSR7k/qZAH+aGbzMuU/Je3b2hTDgBsKT8zsPWAgqeNc68YOvyJtQg7wMTDDz5cAT5jZBj/vWWb9xfr2I5k0LPPnt2Tag/TeC4+DS7RxoKQ5kpaQNmz/SqbtiW4I8QDQUdLOJcojaaySneIza9euLZUlCIIgaCbVvqcKcB9wjaSDgR3M7DlJp5AmMP0ym9GHh4s34P6w6D7lJzT8WNg+W5xNDRIaXUMDbLCGqdEbcVMI34S81HVrzCAiq29z7VFUvlRdU4HjfQnPGGCIp28DDDazvzZRP2Y2CZgEafZvU/mDIAhyowYjy7ypeqTqZgyzgSk0RGqPAP+iBrP43ZWM5sthJdDfz0/IpM8kORrhdXYC5gNHSuqs5Nl7MvBE895JyfqLeRnoWbhfCpxW1N5Jmce5JcrvTDKu6ECKVBtru1yHpiAIgiBHqt6pOtOArwK/BjCzmcAdwFwf6ryb1KGUw6XA9ZLmkIZeC/wU6OQTkhYBR5nZGuBHwOMkA4bnzOz+Zr6HTeovzmBmHwL/DNzl72sjkJ3du51PaPo3GowisvyE9ENgFqmDLnAOMMAnSL0IjGvmewiCIGglWuF+ag1GvmH+UMeE+UMQBOXSYvOHDrvZYZ1PzFMSADP+dGOYPwRBULsM7zy22WVnrJuUo5IgaHvUyvBvk0haX+H2PluD2gp1T5U0sjXqDoIgqFnqYPi3zXSqQRAEQVDrtOlOVdJekh7zCTqPSeohqZ27HknSFyVtLBhH+BrPLxfV0dPTn/PjsMzLHSVNl/SipJskbeNl1mfKj1TaP7UQgU6Q9HvXMNLTJWmi1/MQsFum/EXugvSCpElS2irH3ZvGu7vSMkmHe3o7SVdlnJvO9vT+kp5QcoB6RFK3VrjkQRAEzadKkaqk4ZJekbRc0gUlXh8jaa2k5/34bua10ZJe9WN0U2216U6V5Dx0qzsY/QqY4GtClwF9gL8DngUOV7ID3MPMlhfV8TZwjJkdTFrKMiHz2iDg34G+wD7At8vQ1M3b/SZwhaeNIBk/9AXOALId90R3SToQ2MHLFWjvDkvnAhd72ligF9Cv8L59ic0vgJHuADUF+M9S4sL8IQiCesKXS95A2ue6D3CypD4lst5pZgf5MdnL7kr633sIqT+4uJHlkp/R1jvVwaSlNwC3kTozgDkkp6IjgJ95+kBgQYk6OgA3+xKXu0gXvcDTZrbCO+ppmfo3x31mttHMXgS6etoRwDQz+9TM3qTBwhDgKEnzvf2jaXBJArjXH5+lwcVpGHCTO0BhZu+SOuwDgVnuqnQhsEcpcWY2ycwGmNmALl26lPF2giAI8qAVfH/L8/4dBCz3/+Ufk5ZufqtM0V8HZpnZu+6SN4sGW92SbG2zfwtXeA5prWZ34CLSrjdDgCdLlPk+8BZpnew2wIcl6it+nk3fvijPR5nzrIPSJp++0ublN5K2fFsl6ZKi+gp1fUrDZ9WYM9RSMytlbRgEQVB9DKx1NinvLCm7NnCSO8cV2B1YlXm+mhR5FnOC3ypcBnzfzFY1Unb3zYlp65Hq72nwDT4FKBjqzycNsW50w4XngX8ldbbF7AKs8S3kTgPaZV4bJKmX30s9KVP/W5IO8PQRZeh8Ehjl90O70WAMUehA17l7VDkzgmcC4+RWiT488QrQRdJgT+sg6SubqSMIgmBrYV1h9M2P4nVdpexhiwOT/wF6+i21R0m+7OWW/RxtKVLdUdLqzPNrSE5CUySdD6wluRVhZh9JWgUUjPfnkCwIl5So90bgHkknkpyVst7Cc0n3RfuSOsbpnn4B8CDpF8wLwE5NaJ9OGtpdQvoV9ITrfF/SzZ6+ktLD08VMBnoDiyVtAG42s4k+KWqCpF1In+t1wNIy6guCzxFrTYNWozpbta0G9sw83wN4M5vBd0wrcDNp97JC2SFFZWdvrrFwVKpjwlEpCIJyabGjUvsuNrjj8U1n3EIeeW/yZnX5qN4yYCjwBil4+SczW5rJ081ta5E0AvihmR3qI4HPAgd71ueA/j6XpSRtKVINgiAI2jJVCOLM7BNJ3yNt1NIOmGJmSyVdBjxjZg8A50g6jrTL2bvAGC/7rqTLaRhFvGxzHSpEpxoEQRBUAjPY2CoTlcpo2h4GHi5Kuyhz/iPS5iqlyk4hLVMsi7Y+USkIgiAIaoaIVIMgCILKUAdzeCJSDYIgCIKciEg1CIIgqAhWpXuqlSQ61SAIgqAC1OZWbXkTw79BEARBkBMRqQZBEAStj1EtR6WKEpFqEARBEORERKpBEARBZWidXWpqiohUgyAIgiAnIlINgiAIWh0DrA7uqUanGgRBELQ+ZjH8GwRBEARB+USkGgRBEFSEehj+jUg1CIIgCHIiItUgCIKgMtTBPVVZHXgxBqWR9L/AK1WW0RlYV2UNUBs6akED1IaOWtAAtaGjFjQA7GdmOze3sKQZpPeSN+vMbHgr1NssolOtYyQ9Y2YD6l1DreioBQ21oqMWNNSKjlrQUEs6ap24pxoEQRAEORGdahAEQRDkRHSq9c2kagugNjRAbeioBQ1QGzpqQQPUho5a0AC1o6OmiXuqQRAEQZATEakGQRAEQU5Ep1pHSNpV0ixJr/pjp0byXSlpqaSXJE2QpCpo6CFppmt4UVLPvDRsiQ7P21HSG5ImVlqDpIMkzfXPY7Gkk3Jqe7ikVyQtl3RBide3k3Snvz4/7+u/BTrO889/saTHJO1VaQ2ZfCMlmaRWmQFbjg5J/+jXY6mkOyqtwf8uH5e00D+Tb+Stoc1jZnHUyQFcCVzg5xcA40vkOQx4Cmjnx1xgSCU1+GuzgWP8fCdgx0pfi0ze64E7gIlV+Dx6A/v6eXdgDfDFFrbbDngN2BvYFlgE9CnKcxZwk5+PAu5she9jOTqOKnz2wJl56yhHg+fbGXgSmAcMqNK12BdYCHTy57tVQcMk4Ew/7wOszPtatPUjItX64lvALX5+C3B8iTwGbE/6o9oO6AC8VUkNkvoA7c1sFoCZrTezv+SooSwdrqU/0BWYmXP7ZWkws2Vm9qqfvwm8DXRpYbuDgOVmtsLMPgZ+7Voa03Y3MDTPEYtydZjZ45nPfh6wR6U1OJeTfgR9mHP7W6LjDOAGM3sPwMzeroIGAzr6+S7AmzlraPNEp1pfdDWzNQD+uFtxBjObCzxOiojWAI+Y2UuV1ECKzt6XdK8PM/1cUrscNZSlQ9I2wNXA+Tm3XbaGIj2DSD92Xmthu7sDqzLPV3tayTxm9gnwAfA3LWy3OTqynA78ptIaJPUD9jSzB3Nue4t0kP4uekt6StI8SXm7CJWj4RLgVEmrgYeBs3PW0OYJ79+tDEmPAl8q8dKPyyz/ZeAAGiKCWZKOMLMnK6WB9L08HOgHvA7cCYwB/qtcDTnpOAt42MxWNTdIy0FDoZ5uwG3AaLMWG6iWejPFywDKydNSym5D0qnAAODISmrwH1bXkr5/rUk516I9aQh4COnvc46kA83s/QpqOBmYamZXSxoM3OYatn5T3zKJTnUrw8yGNfaapLckdTOzNf5PutTw0Qhgnpmt9zK/AQ4l3U+qlIbVwEIzW+Fl7nMNW9Sp5qBjMHC4pLNI93W3lbTezBqdzNIKGpDUEXgIuNDM5pXb9mZYDeyZeb4Hmw7jFfKsltSeNNT3bg5tb6kOJA0j/Qg50sw+qrCGnYEDgdn+w+pLwAOSjjOzZyqoo5BnnpltAP4g6RVSJ7ugghpOB4ZDGtWStD3Jzzfvoeg2Swz/1hcPAKP9fDRwf4k8rwNHSmovqQMpMshz+LccDQuATpIK9w6PBl7MUUNZOszsFDPrYWY9gR8At25Jh5qHBknbAtO97btyancBsK+kXl7/KNfSmLaRwG/NLO9ItUkdPvT6S+C4VriH2KQGM/vAzDqbWU//HsxzLXl2qE3qcO4jTdxCUmfScPCKCmt4HRjqGg4gzb9Ym6OGtk+1Z0rFUbmDdE/sMeBVf9zV0wcAk/28Hemf2EukjuyaSmvw58cAi4ElwFRg22royOQfQ/6zf8v5PE4FNgDPZ46Dcmj7G8Ay0v3ZH3vaZaQOA9I/y7uA5cDTwN6t9J1sSsejpIlyhff+QKU1FOWdTSvM/i3zWgi4xv8ulwCjqqChD2l1wCL/PI5tjWvRlo9wVAqCIAiCnIjh3yAIgiDIiehUgyAIgiAnolMNgiAIgpyITjUIgiAIciI61SAIgiDIiehUgyAIgiAnolMNgiAIgpyITjUIgiAIcuL/ATvXTP4ngyqjAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x576 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "name = \"SuperScriptIV : TGIRT            \" \n",
    "barplotFigure(y=df1[\"Comparison\", \"TvS\"], GC_content=df1[\"GC content\", \"TGIRT\"], name=name)\n",
    "\n",
    "name = \"ImProm-II 42°C : TGIRT              \" \n",
    "barplotFigure(y=df1[\"Comparison\", \"Tv42\"], GC_content=df1[\"GC content\", \"ImProm-II 42C\"], name=name)\n",
    "\n",
    "name = \"ImProm-II 55°C : TGIRT               \" \n",
    "barplotFigure(y=df1[\"Comparison\", \"Tv55\"], GC_content=df1[\"GC content\", \"ImProm-II 55C\"], name=name)\n",
    "\n",
    "name = \"ImProm-II 42°C : SuperScriptIV  \" \n",
    "barplotFigure(y=df1[\"Comparison\", \"Sv42\"], GC_content=df1[\"GC content\", \"ImProm-II 42C\"], name=name)\n",
    "\n",
    "name = \"ImProm-II 55°C : SuperScriptIV  \" \n",
    "barplotFigure(y=df1[\"Comparison\", \"Sv55\"], GC_content=df1[\"GC content\", \"ImProm-II 55C\"], name=name)\n",
    "\n",
    "name = \"ImProm-II 55°C : ImProm-II 42°C\" \n",
    "barplotFigure(y=df1[\"Comparison\", \"P42v55\"], GC_content=df1[\"GC content\", \"ImProm-II 42C\"], name=name)"
   ]
  }
 ],
 "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.6.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}