{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 27,
   "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\">Average</th>\n",
       "      <th colspan=\"4\" halign=\"left\">StdDev</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",
       "      <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>68.925165</td>\n",
       "      <td>52.142589</td>\n",
       "      <td>35.392582</td>\n",
       "      <td>32.702259</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>84.882271</td>\n",
       "      <td>43.312816</td>\n",
       "      <td>112.839266</td>\n",
       "      <td>49.726251</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>40.947527</td>\n",
       "      <td>21.505813</td>\n",
       "      <td>33.328666</td>\n",
       "      <td>40.133527</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>481.747652</td>\n",
       "      <td>232.488279</td>\n",
       "      <td>303.496787</td>\n",
       "      <td>382.188435</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>440.766945</td>\n",
       "      <td>215.612384</td>\n",
       "      <td>78.349856</td>\n",
       "      <td>48.864097</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>556.201672</td>\n",
       "      <td>552.653599</td>\n",
       "      <td>80.198504</td>\n",
       "      <td>162.080844</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>20.329781</td>\n",
       "      <td>50.906778</td>\n",
       "      <td>16.577093</td>\n",
       "      <td>16.531788</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>5.630275</td>\n",
       "      <td>11.058933</td>\n",
       "      <td>3.049590</td>\n",
       "      <td>1.341641</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>19.798990</td>\n",
       "      <td>4.393177</td>\n",
       "      <td>16.452963</td>\n",
       "      <td>9.110434</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>36.956731</td>\n",
       "      <td>12.198361</td>\n",
       "      <td>25.510782</td>\n",
       "      <td>19.165072</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>111.576431</td>\n",
       "      <td>77.212046</td>\n",
       "      <td>105.991981</td>\n",
       "      <td>106.924272</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>5.761944</td>\n",
       "      <td>7.987490</td>\n",
       "      <td>8.820431</td>\n",
       "      <td>6.730527</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>32.645061</td>\n",
       "      <td>13.619838</td>\n",
       "      <td>18.796276</td>\n",
       "      <td>30.654527</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>4.037326</td>\n",
       "      <td>3.768289</td>\n",
       "      <td>9.808160</td>\n",
       "      <td>8.689074</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>133.784155</td>\n",
       "      <td>60.870354</td>\n",
       "      <td>74.751589</td>\n",
       "      <td>65.076878</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>3.000000</td>\n",
       "      <td>3.082207</td>\n",
       "      <td>2.489980</td>\n",
       "      <td>2.949576</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>23.430749</td>\n",
       "      <td>22.074873</td>\n",
       "      <td>34.629467</td>\n",
       "      <td>23.488295</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>19.917329</td>\n",
       "      <td>2.966479</td>\n",
       "      <td>4.147288</td>\n",
       "      <td>3.962323</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>523.532520</td>\n",
       "      <td>243.335776</td>\n",
       "      <td>212.430459</td>\n",
       "      <td>139.263419</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>39.232639</td>\n",
       "      <td>12.389512</td>\n",
       "      <td>33.714982</td>\n",
       "      <td>33.729809</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>7.563068</td>\n",
       "      <td>8.093207</td>\n",
       "      <td>12.700394</td>\n",
       "      <td>5.805170</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>32.893768</td>\n",
       "      <td>25.674890</td>\n",
       "      <td>16.513631</td>\n",
       "      <td>14.872794</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>11.278298</td>\n",
       "      <td>14.352700</td>\n",
       "      <td>12.845233</td>\n",
       "      <td>7.245688</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>21.173096</td>\n",
       "      <td>12.194261</td>\n",
       "      <td>35.686132</td>\n",
       "      <td>26.470739</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>25.951879</td>\n",
       "      <td>45.970643</td>\n",
       "      <td>38.448667</td>\n",
       "      <td>43.888495</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>3.563706</td>\n",
       "      <td>4.438468</td>\n",
       "      <td>4.086563</td>\n",
       "      <td>5.319774</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>5.873670</td>\n",
       "      <td>11.104053</td>\n",
       "      <td>14.525839</td>\n",
       "      <td>9.476286</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>141.663333</td>\n",
       "      <td>188.279845</td>\n",
       "      <td>143.224649</td>\n",
       "      <td>194.803747</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>7.596052</td>\n",
       "      <td>3.834058</td>\n",
       "      <td>4.098780</td>\n",
       "      <td>3.646917</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>4.438468</td>\n",
       "      <td>3.420526</td>\n",
       "      <td>1.923538</td>\n",
       "      <td>2.701851</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>Acidobacteriia</th>\n",
       "      <td>257</td>\n",
       "      <td>259</td>\n",
       "      <td>444</td>\n",
       "      <td>359</td>\n",
       "      <td>29.364945</td>\n",
       "      <td>6.942622</td>\n",
       "      <td>18.088670</td>\n",
       "      <td>28.455228</td>\n",
       "      <td>0.555791</td>\n",
       "      <td>0.561708</td>\n",
       "      <td>0.561813</td>\n",
       "      <td>0.559804</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>109.917242</td>\n",
       "      <td>47.300106</td>\n",
       "      <td>46.301188</td>\n",
       "      <td>55.410288</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>254.825430</td>\n",
       "      <td>135.417133</td>\n",
       "      <td>190.525064</td>\n",
       "      <td>174.316666</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>85.207394</td>\n",
       "      <td>81.041964</td>\n",
       "      <td>50.257338</td>\n",
       "      <td>56.389715</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>17.840964</td>\n",
       "      <td>14.747881</td>\n",
       "      <td>12.409674</td>\n",
       "      <td>4.560702</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>67.310475</td>\n",
       "      <td>71.333723</td>\n",
       "      <td>193.530359</td>\n",
       "      <td>158.676400</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>14.754660</td>\n",
       "      <td>18.686894</td>\n",
       "      <td>3.435113</td>\n",
       "      <td>4.147288</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>11.436783</td>\n",
       "      <td>10.285913</td>\n",
       "      <td>11.567195</td>\n",
       "      <td>16.269604</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>51.200586</td>\n",
       "      <td>37.036469</td>\n",
       "      <td>44.549972</td>\n",
       "      <td>42.710654</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>25.416530</td>\n",
       "      <td>14.307341</td>\n",
       "      <td>12.095454</td>\n",
       "      <td>36.833409</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>31.546791</td>\n",
       "      <td>14.876155</td>\n",
       "      <td>33.899853</td>\n",
       "      <td>36.929663</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>Actinobacteria</th>\n",
       "      <td>5864</td>\n",
       "      <td>11263</td>\n",
       "      <td>10105</td>\n",
       "      <td>11353</td>\n",
       "      <td>445.351771</td>\n",
       "      <td>581.639321</td>\n",
       "      <td>223.156895</td>\n",
       "      <td>274.296919</td>\n",
       "      <td>0.583321</td>\n",
       "      <td>0.583805</td>\n",
       "      <td>0.584725</td>\n",
       "      <td>0.586320</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "Category                  Average                                     \\\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",
       "Acidobacteriia                257           259           444    359   \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",
       "Actinobacteria               5864         11263         10105  11353   \n",
       "\n",
       "Category                   StdDev                                          \\\n",
       "Enzyme              ImProm-II 42C ImProm-II 55C SuperScriptIV       TGIRT   \n",
       "Low abundance           68.925165     52.142589     35.392582   32.702259   \n",
       "Verrucomicrobiae        84.882271     43.312816    112.839266   49.726251   \n",
       "Unclassified            40.947527     21.505813     33.328666   40.133527   \n",
       "Alphaproteobacteria    481.747652    232.488279    303.496787  382.188435   \n",
       "Deltaproteobacteria    440.766945    215.612384     78.349856   48.864097   \n",
       "Gammaproteobacteria    556.201672    552.653599     80.198504  162.080844   \n",
       "Planctomycetacia        20.329781     50.906778     16.577093   16.531788   \n",
       "Phycisphaerae            5.630275     11.058933      3.049590    1.341641   \n",
       "Saccharimonadia         19.798990      4.393177     16.452963    9.110434   \n",
       "Longimicrobia           36.956731     12.198361     25.510782   19.165072   \n",
       "Gemmatimonadetes       111.576431     77.212046    105.991981  106.924272   \n",
       "BD2-11                   5.761944      7.987490      8.820431    6.730527   \n",
       "S0134                   32.645061     13.619838     18.796276   30.654527   \n",
       "Nitrospira               4.037326      3.768289      9.808160    8.689074   \n",
       "Bacilli                133.784155     60.870354     74.751589   65.076878   \n",
       "Fibrobacteria            3.000000      3.082207      2.489980    2.949576   \n",
       "Entotheonellia          23.430749     22.074873     34.629467   23.488295   \n",
       "Deinococci              19.917329      2.966479      4.147288    3.962323   \n",
       "Oxyphotobacteria       523.532520    243.335776    212.430459  139.263419   \n",
       "Chloroflexia            39.232639     12.389512     33.714982   33.729809   \n",
       "Anaerolineae             7.563068      8.093207     12.700394    5.805170   \n",
       "Gitt-GS-136             32.893768     25.674890     16.513631   14.872794   \n",
       "JG30-KF-CM66            11.278298     14.352700     12.845233    7.245688   \n",
       "KD4-96                  21.173096     12.194261     35.686132   26.470739   \n",
       "TK10                    25.951879     45.970643     38.448667   43.888495   \n",
       "Ktedonobacteria          3.563706      4.438468      4.086563    5.319774   \n",
       "Dehalococcoidia          5.873670     11.104053     14.525839    9.476286   \n",
       "Bacteroidia            141.663333    188.279845    143.224649  194.803747   \n",
       "Fimbriimonadia           7.596052      3.834058      4.098780    3.646917   \n",
       "uncultured               4.438468      3.420526      1.923538    2.701851   \n",
       "Acidobacteriia          29.364945      6.942622     18.088670   28.455228   \n",
       "Acidimicrobiia         109.917242     47.300106     46.301188   55.410288   \n",
       "Thermoleophilia        254.825430    135.417133    190.525064  174.316666   \n",
       "0319-7L14               85.207394     81.041964     50.257338   56.389715   \n",
       "Nitriliruptoria         17.840964     14.747881     12.409674    4.560702   \n",
       "Rubrobacteria           67.310475     71.333723    193.530359  158.676400   \n",
       "MB-A2-108               14.754660     18.686894      3.435113    4.147288   \n",
       "Thermoanaerobaculia     11.436783     10.285913     11.567195   16.269604   \n",
       "Subgroup                51.200586     37.036469     44.549972   42.710654   \n",
       "Holophagae              25.416530     14.307341     12.095454   36.833409   \n",
       "Blastocatellia          31.546791     14.876155     33.899853   36.929663   \n",
       "Actinobacteria         445.351771    581.639321    223.156895  274.296919   \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",
       "Acidobacteriia           0.555791      0.561708      0.561813  0.559804  \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",
       "Actinobacteria           0.583321      0.583805      0.584725  0.586320  "
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "import math\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": 41,
   "metadata": {},
   "outputs": [],
   "source": [
    "def relativeEnrichment(x, y, stdDevX, stdDevY):\n",
    "    valueVect = []\n",
    "    errorVect = []\n",
    "    \n",
    "    for a, b, stdDevA, stdDevB in zip(x, y, stdDevX, stdDevY):\n",
    "        value = (a-b)/(a+b)\n",
    "        \n",
    "        ab_1 = 2/(math.pow((a+b), 2))\n",
    "        ab_2 = math.sqrt(math.pow(b,2)*math.pow(stdDevA,2) + math.pow(a,2)*math.pow(stdDevB,2))\n",
    "        error = ab_1*ab_2\n",
    "        \n",
    "        valueVect.append(value)\n",
    "        errorVect.append(math.fabs(error))\n",
    "    \n",
    "    return valueVect, errorVect"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [],
   "source": [
    "idx = pd.IndexSlice\n",
    "\n",
    "TvS, TvSerror = relativeEnrichment(df1.loc[:, idx[\"Average\", \"TGIRT\"]],\n",
    "                                   df1.loc[:, idx[\"Average\", \"SuperScriptIV\"]],\n",
    "                                   df1.loc[:, idx[\"StdDev\", \"TGIRT\"]],\n",
    "                                   df1.loc[:, idx[\"StdDev\", \"SuperScriptIV\"]])\n",
    "df1[\"Comparison\", \"TvS\"] = TvS\n",
    "df1[\"Comparison\", \"TvSerror\"] = TvSerror\n",
    "\n",
    "Tv42, Tv42error = relativeEnrichment(df1.loc[:, idx[\"Average\", \"TGIRT\"]],\n",
    "                                     df1.loc[:, idx[\"Average\", \"ImProm-II 42C\"]],\n",
    "                                     df1.loc[:, idx[\"StdDev\", \"TGIRT\"]],\n",
    "                                     df1.loc[:, idx[\"StdDev\", \"ImProm-II 42C\"]])\n",
    "df1[\"Comparison\", \"Tv42\"] = Tv42\n",
    "df1[\"Comparison\", \"Tv42error\"] = Tv42error\n",
    "\n",
    "Tv55, Tv55error = relativeEnrichment(df1.loc[:, idx[\"Average\", \"TGIRT\"]],\n",
    "                                     df1.loc[:, idx[\"Average\", \"ImProm-II 55C\"]],\n",
    "                                     df1.loc[:, idx[\"StdDev\", \"TGIRT\"]],\n",
    "                                     df1.loc[:, idx[\"StdDev\", \"ImProm-II 55C\"]])\n",
    "df1[\"Comparison\", \"Tv55\"] = Tv55\n",
    "df1[\"Comparison\", \"Tv55error\"] = Tv55error\n",
    "\n",
    "Sv42, Sv42error = relativeEnrichment(df1.loc[:, idx[\"Average\", \"SuperScriptIV\"]],\n",
    "                                     df1.loc[:, idx[\"Average\", \"ImProm-II 42C\"]],\n",
    "                                     df1.loc[:, idx[\"StdDev\", \"SuperScriptIV\"]],\n",
    "                                     df1.loc[:, idx[\"StdDev\", \"ImProm-II 42C\"]])\n",
    "df1[\"Comparison\", \"Sv42\"] = Sv42\n",
    "df1[\"Comparison\", \"Sv42error\"] = Sv42error\n",
    "\n",
    "Sv55, Sv55error = relativeEnrichment(df1.loc[:, idx[\"Average\", \"SuperScriptIV\"]], \n",
    "                                     df1.loc[:, idx[\"Average\", \"ImProm-II 55C\"]],\n",
    "                                     df1.loc[:, idx[\"StdDev\", \"SuperScriptIV\"]], \n",
    "                                     df1.loc[:, idx[\"StdDev\", \"ImProm-II 55C\"]])\n",
    "df1[\"Comparison\", \"Sv55\"] = Sv55\n",
    "df1[\"Comparison\", \"Sv55error\"] = Sv55error\n",
    "\n",
    "P42v55, P42v55error = relativeEnrichment(df1.loc[:, idx[\"Average\", \"ImProm-II 42C\"]],\n",
    "                                         df1.loc[:, idx[\"Average\", \"ImProm-II 55C\"]],\n",
    "                                         df1.loc[:, idx[\"StdDev\", \"ImProm-II 42C\"]],\n",
    "                                         df1.loc[:, idx[\"StdDev\", \"ImProm-II 55C\"]])\n",
    "df1[\"Comparison\", \"P42v55\"] = P42v55\n",
    "df1[\"Comparison\", \"P42v55error\"] = P42v55error"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "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, y_stdDev, 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, xerr=y_stdDev, color=colors, align=\"center\")\n",
    "    #plt.errorbar(y, x_pos, 0, y_stdDev)\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": 44,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdUAAAHiCAYAAAC+xlbiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsnXeYXVW5xn8voRM6iIDE0HsPINVQBESRIE1AIIrk4gXlisBVQQhNELmI1CiIoQlID72HIjUhJKGDELqQ0AMBIXnvH2udZOfknDlnZs7MnJn5fs+zn7P32qvtfZL5zlrrW+8n2wRBEARB0H5m6+oOBEEQBEFPIYxqEARBEDSIMKpBEARB0CDCqAZBEARBgwijGgRBEAQNIoxqEARBEDSIMKpBEARB0CDCqAZNh6QJkrZpY9nhkv4jabKk9yTdIWmVRvexzr70l2RJsxf6dkKNvJMLx28L94vPVTr6FO4fI+nfkq6UNFchfQFJp0t6NZd5MV8v1t7naUP5CZKm5H68Lelvkvq2pa72ImmwpAfK+lbz35yku8u+069IukzSm5I+lPRPSRsV8kvSn/Pznl1W15KS/irpLUkfS3pW0rGS5mvkswadSxjVoCdyiu2+wNeAd4DhlTK11Th0MAvZ7puP48vunVK419f2VABJKwIbAl8HxgA/zOlzAncBqwPbAwsAmwDv5vxdwY75u1kP2AA4qjxDNkRN97dJ0t5A+b+ZvsBjwPrAIsCFwE2FHwslQ700MLukrXNdiwAPAfMAG9ueH/gWsBCwfEc+R9CxNN0/3CAokkcU/5T0R0kfSHpJ0iY5/TVJ70jar1JZ258CfwfWyHUNlXSVpEskfQQMljRXHrm9mY/TSyM9SQMlvS7piNzOW5IGSdpB0vN5JPybTnsZ1ekDKH+WDoB9gX7Azraftj3N9ju2j7d9c3sbzaPncyTdkkef/5T01fwO388jr3UrlbX9BnALM76bkZJOlPRP4FNgOUlLSRqR3/OLkg4otD00j8ovyaO88ZJWkvTr/F29Jmnb9j5job0FgWOAI8qe4yXbp9l+y/ZU238B5gRWzln6kP7OFj8BDgU+Bn5oe0Ku6zXbh9ge16h+B51PGNWgO7ARMA5YlGQkLyeNclYgjcrOqjSNmNP2Jo3eSuwEXEUaEVwKHAl8A1gHWJs0giuOnr4KzE0aaRwNnJfbXB/YHDha0nINek6AV7Ih/1uFKdr/zgZmtKRdSom2nyW9n1dyvy7Ot7YBbrU9ud7GJd0o6Vet6O/upPe1GPA5afT1eL6+CjitSjvLADsw83ezDzAEmD8/y2XA68BSwK7A70ojvcyOpGddONdzG+lv2tLAccCfW/EctfgdcC7w75YySVqHZFRfzEm35+vX8/Ud+XMb4Brb0xrYx6AZsB1HHE11ABOAbfL5YOCFwr01AQNLFNLeBdbJ58OBz4APSH8ARwDL53tDgfvK2voXsEPhejtgQj4fCEwB+uTr+XPbGxXyjwYGVXmO/jn/7IW+nVAlb19gAGl6cQmSQbqtcH890o+K2UnG6GNg0xrv8Q7g5AZ+L5We57zC/Z8Bz5R9Vx+Ufa+T83fzCnAOME++NxI4rpB3GWAqMH8h7SRgeOG7vKNwb8dcd/l3tVCVZxkMPFDp31yFvAOAJ/K7n+kdlOVbABgP/LqOd/kCcGBX/1+Lo/FHM64pBUE5bxfOpwDYLk8rjlRPtT3LWl3mtbLrpUh/4Eu8ktNKvOu8dllqu0J/+gJIKo4IV6vSfkWcRpOjSvVLOhh4S9ICtj+y/Xgh+82SLgW+D/yzhWrfBZZsTT/aQPm7aOl7gfQD5M4qdRW/m6WA92x/XEh7hWTgqrU9qcJ31VfSmqSpZoBXbK9epf1ZyGu75wCH2P5SUrV88wA3AA/bPqmOqjvjuwm6gJj+DXob5WGZ3iQ5+JTol9NaX/HMTkSvtrWDperyZ+W/4ul+tXsl7gS260bepMXv5k1gEUnzF9L6AW+0ulL7/sL3UrdBzSxAMuRXSPo3ySkJ4HVJmwPkNfjrct/+q8567wR2bkaHrKB9xBca9HYuA46StHhewzwauKSzOyFpI0krS5pN0qLAGcBI2x/m+7tK6pvvb0ta1x1Ro9qLSaO/qyWtUqpb0m8k7dChD9RObL8GPAicJGluSWsB+5PWwTuTD0mj5nXyUXpv6wOPSJqDNFU/BdjX9a+RnkYy2BdK+jqApKUlnZafNeimhFENejsnkKZdx5HWwx7PaZ3NcsCtpLXSJ0lOP3sW7h9CGgl9APwBOMD2yJYqtP05ySHmWdL66kfAoyQnokcqlcmevM3g0Qzp+fuTRq3XAsfYvqPFEg3GiX+XDmBivvW27f+Qtih9F9gW+EAz9hBvXqPe93LZL0jG+WPS9qcPmeHkFHRDZEeQ8iAIgiBoBDFSDYIgCIIGEUY1CIIgCBpEGNUgCIIgaBBhVIMgCIKgQYRRDYIgCIIGEYpKvZjFFlvM/fv37+puBEHQDRg9evQk24u3tfx2W87nd9+bWjtjKxk97vPbbG/f8IrbSBjVXkz//v0ZNWpU7YxBEPR6JL1SO1d13n1vKo/e1q9R3ZlOnyVfaHVs4I4kjGoQBEHQ4RiYRs8PyhNrqkEQBEHQIGKkGgRBEHQCZmovCB8bI9UgCIIgaBAxUg2CIAg6nLSm2vO15sOoBkEQBJ1COCoFQRAEQVA3MVINgiAIOhxjpvaCUKMxUg2CIAiCBhEj1SAIgqBTCEelIAiCIGgABqb2AqMa079BEARB0CBipBoEQRB0Cr1h+jdGqkEQBEHQIGKkGgRBEHQ4hl6xpSaMahAEQdAp9Hw9pZj+DYIgCIKGESPVIAiCoMMx7hVbanqEUZW0M3ANsKrtZ1vINxi43fab+fp84DTbT7ehzeHAjbavalOnq/SpFeUOBD61fVF72g/ax8rXHNfVXWg4z33/6K7uQhB0W3qEUQX2BB4AfgAMbSHfYOBJ4E0A2z/p6I7VwWAKfaoHSbPbHtZhPQqCIGg0hqk9f6Da/Y2qpL7ApsCWwAiyUZV0BLAPaW38FmAUMAC4VNIUYOOcfpjtUZImA38CvgtMAXay/bakrwMXAIsDE4Ef2X41N7+NpEOAJYBDbd8oqT9wMTBfznOw7Qdb0afVgNOAvsAkYLDttySNBB7MzzpC0vzAZNunSjoAGALMCbwI7GP70wa83gAYOHBg1XuvTXql6r1ljvtRB/QmCIJmpic4Kg0CbrX9PPCepPUkfTunb2R7beCUPE07Ctjb9jq2p5TVMx/wcM5/H3BATj8LuMj2WsClwBmFMv2BbwLfAYZJmht4B/iW7fWAPUr56+kT8CVwJrCr7fVJxvzEQnsL2f6m7f8r6/s1tjfI9T4D7F/tZUkaImmUpFETJ06sli0IgqChpCDljT+ajW4/UiVN/Z6ezy/P17MBfyuN1my/V0c9/wFuzOejgW/l842B7+fzi4FTCmX+YXsa8IKkl4BVgJeBsyStA0wFVsp5t6mjTysDawB3SALoA7xVuH9Flb6vIekEYCHSCPe2ag9p+y/AXwAGDBjQCyZj2s/IkSOr3uuJa6pB0DGIqairO9HhdGujKmlRYCuSUTHJCBm4On+2hi/s6TuTp1L93bjKeen6F8DbwNok4/5Zqbt19EnAU7Y3rnL/kyrpw4FBtsdmx6eBNdoJgiAIOoDuPv27K2lq9uu2+9tehjRSfA/4saR5ASQtkvN/DMzfyjYeJDlAAexNcogqsZuk2SQtDywHPAcsCLyVR7D7kAw9wO119Ok5YHFJG+c8c0havY4+zg+8JWmO3McgCIKmwsA0N/5oNrr1SJU01XtyWdrVwKokp6VRkv4D3Az8hjSiG1ZwCqqHnwMXSDqc7KhUuPcccC/JUelA259JOge4WtJuwD3k0aXtW/OUcK0+7QqcIWlB0vdzOvBUjT7+FngEeAUYT+t/OARtJLafBEFQRO4FWoxBZQYMGOBRo0Z1dTeCIOgGSBpte0Bby6+x1pz+x02LN7JLAKze78129avRdPeRakUkTSWN2ERaHz3Y9oN5u8uNttdoQ52/sf27BvdzIPCf0pabFvINBgbYPljSUGZspTkOuM/2nY3sV9C7WefgP3ZaW0+c9YtOayvoWlKQ8p7vqNTd11SrMSVvm1kb+DVwUgPq/E0D6ihnILBJWwvbPjoMahAEQfPQU41qkQWA98sTJfWXdL+kx/OxSU5fUtJ9kp6Q9KSkzSWdDMyT0y7N+Q7N95+U9D+FeveVNE7SWEkX57QdJT0iaYykOyUtkUfNBwK/yPVuLmlxSVdLeiwfm7b0YJKGS9o1nx+dyzwp6S/Ke3KCIAiahWlWw49mo0dO/5INIDA3sCRp2005JZGGzyStCFxGUjfaC7jN9omS+gDz2r5f0sFZoAFJ65McljYiTTE/Iule0l7XI4FNbU8qePg+AHzDtiX9BDjC9i8lDSNP5eZ6/w780fYDkvqR9puuWuczn2X7uFzPxSRlqBvqfmNB0EBevPbsuvINfPL6mnla2iccBM1GTzWqUwoGcGPgIknl66hzUFmk4TGSt+8cwHW2n6hQ/2bAtbY/yW1cA2xOWja4yvYkmEng4WvAFZKWJEkJvlyl39sAqxUGmQtkOcJ62DLLIM4LLELyGJ7FqEoaQpI0pF+/fnVWHQRB0D56y5pqTzWq07H9kKTFSNq9RSqKNNi+T9IWJOnBiyX9oUIkmGr/MqoJPJxJioYzIjsnDa1SfjZg43IJxVozuVke8RySM9Nr2Zlp7kp5Q1Ep6AxW2PmguvKNDEelXoMRU3vBimOPf0JJq5AEGN4tu1VRpCEL6L9j+zzgr8B6Of8XefQKSRt4kKR5Jc0H7AzcD9wF7J6VnooCDwsCb+Tz/Qp9KBejuB04uND3dep8zJIBnaQUYGDXOssFQRAEDaSnjlRLa6qQRo/72Z5aNuKrKNJA8sg9XNIXwGRg35z+F2CcpMdt760UT/XRfO9822MAJJ0I3Ju39YwhhXYbClwp6Q3gYWDZXO4G4CpJOwE/IwlNnC1pHOm7uY/kzNQitj+QdB5pG9EE0hR2EARBU9GMjkWNJsQfejEh/hAEQb20V/xh1bXm8vAblmpklwD4Rv8JIf4QBEEQ9C7CUakHI2my7b6F68FkxaIWygylsP2lle21uWwQdDbLnlkerrd5eflnv+zqLgR1I6a6x7vx9HxHpSAIgiDoLHrlSLUlsvfvBaQtOBOBH9l+tSzPOsAw0p7QfwE/tv2+pJHAE8CGJCWnH9suOTOtlu/3A063fUau6zpgGZIH75/ylhck7Q/8L/Am8ALwedb+3RE4irTf9V1gb9tvZy/kM4E1Sd/rUNu1d9YHQRPx1hnntCr/wKvr1zcJEYmuxcC0XjCO6/lPWJmS5OAT2Uv4uMK9s0gxWtcCLgXOqFD+IuB/c57xwDGFe/PZ3gT4b5JxLrEKsB3J4B5T2J7zY9vrk9Scfi5pUUlLkcK5fQP4Vi5boqTOtC5wOXBETj8SuNv2BsCWwB+yoZ0JSUMkjZI0auLEiS2+pCAIgqB19NaR6nTFJZixppovNwa+n88vBk4pFsxxTheyfW9OuhC4spDlMpguIrGApIVy+k22Pwc+l/QOKQbr6yRDunPOswywIvBV4N6SIpOkK5mh+FRNnWlb4HuSDsvXc5NGxc8U+x/iD0Ezs+TP/7tV+UfGmmq3IhyVAqiskNSa/KXrzwtpU4HZs7rSNiQVpU/z9PDcVFdsgurqTAJ2sf1cK/sbBEEQNIjeOv3bEg8CP8jne5OmW6dj+0PgfUmb56R9gHsLWfYAkLQZ8GHOX40FgfezQV2FNN0LSVTim5IWljQ7sEtZmUrqTLcBPytFp5G0bs0nDYIg6CTs5P3b6KPZiJHqrPycJKh/ONlRqUKe/YBhkuYFXirL876kB8mOSjXauhU4MCsoPUdSW8L2G5J+BzxCclR6GigZ56FUVmc6HjidpPokkrLSd+t85iCYTmxTCTqKab1g+jcUlRpInr49zHa7ZYok9bU9OY9UrwUusH1te+stEopKQRDUS3sVlVZacx6fOWLZ2hlbyfbLPROKSt0NSUeS4qxOBaYB/2X7kSrZh0ga2AChh6GStiGtsd4OXNfO+oKgLpa//Hcd3sa/fvCbDm8jaC6SolLzTdc2mjCqNcjxWL8LrGf78xxGbs5KeW0PzOpJ7Wmvj+2ptg+rnTsIgiBoJnr+z4b2syQwKW+HwfYk229KmpANLJIG5KnfEmtLulvSC5IOyHlmk3SOpKck3SjpZkm75nsTJB0t6QFgN0nrSHpY0jhJ10paOOcbKWlAPl9M0oR8PljS9ZJulfScpOK+2SAIgiYgHJWCxO3A0ZKeB+4ErijsUa3GWiRP3vmAMZJuAjYB+pMUj75C2j9aFIf4zPZmANlx6We275V0HElc4n9qtLkhsAbwKfCYpJsasbYbBNV487jz21x24LDb21w2lJG6J6GoFABgezKwPjCE5A18RRaLaInrbU+xPYkUq3VDYDPgStvTbP87pxe5AqqKS2xRR1fvsP2u7SnANbm9WQhFpSAIgo4jRqp1YHsqMBIYKWk8aUvNl8z4UTJ3eZEK17V8yT+pcZ82tDkLoagUNIqljv5Jm8uODEelXsnULgpSLml74E9AH+B82ydXyLM7acuigbG298o67+eStkhOBU60fUVLbcVItQaSVpa0YiFpHeAV0j7Q9XPaLmXFdpI0t6RFgYHAYyQRiV3y2uoSOX0WaohLFNvctazotyQtImkeYBDwz3qfMQiCoKciqQ9wNvBtYDVgT0mrleVZEfg1sKnt1Zmx3PYpsG9O2x44vSA9W5EYqdamL3BmfpFfAi+SpoJXBf4q6TckkYYijwI3kbR3j8+OTVcDWwNPAs/nMtXUlqqJS5wK/EPSPsDdZWUeIGkVrwD8PdZTgyBoJoy6akvNhsCLtl8CkHQ5sBNJVKfEAcDZtt8HsP1O/ny+lCH/HX+HFMHsg2qNhVGtge3RJCejcu5nhsh9Mf/QKvVMk3RYFnRYlGR4x+d7/cvyPsEMycJi+rMkJ6gSRxXO32kpyHoQ1EvsIQ06imld4627NPBa4fp1YKOyPCsBSPonaYp4qO1bixkkbUjaTvmvlhoLo9q53JhHvHOSRrD/7uoOBUEQdHMWk1ScmftLKS51ptJCbrk/yeykCGEDSZHA7pe0hu0PAHJUsIuB/WxPa6kzdRvVPLq6K19+lbRoO5G0TeRN26tVKdqrkDTZdt9K92wPbED9w4EbbV8l6XxSxJrhwPD21h30Hla66viq957f9bed2JOgt9CBikqTasgUvk4Kq1niayRN9fI8D9v+AnhZ0nMkI/uYpAVIy3lH2X64VmfqfsK8XWOdHId0GPDHfL4OSbqvTWRt226DEk3h4GX7J7afrp0zCIKg1/IYsKKkZSXNSYpCNqIsz3XAlpCEdUjTwS/l/NcCF9m+kjpolEHrI+k80trjG8BOtqdIWp7kdbU4yYvqANvP5tHWe8C6wOOSPiZFW1kyP8yhpDXFb+f6drT9haStSc46s5Ne1E+zdODRwI7APKTQbf9l21nl6BHSy1oI2N/2/ZL6k4by8+X+H2z7QYAcnWZ3YC7gWtvH5Py3kPaWbgwMkrQJ8BvS1MJNtv+39DIk/V9u833gB7YnZmWlIaSp3xeBfXLItyVIP1KWy8V/SvoVdaPtNXJ9hwF9y9driwL+ks4FNsjv4CrboarUyxg4cGDdeV+f+Er1es66a/p5CC0EjcKoS7bU2P5S0sGk8Jh9SMFJnsrCOqNsj8j3tpX0NGkW9nDb70r6IUknYNGCPsHg7PdSkUaNuFYkeU6tTvKKKm0x+QtJGWh94DDgnEKZlYBtbJfiTC0PfIfklXUJcI/tNYEpwHckzU2a4twjp89OMkAAZ9neIBuheZg55NnstjckuUiXDM07wLdsr0eKf3oGgKRt87NsSBqBry+pJLywMunXyrrAF8Dvga1yvg0kDcr55gMez3XfW2jzmtzHtUlqSvvn9DOAe3P6esBTLbznljgyT4GsRYrFulalTCH+EARBb8P2zbZXsr287RNz2tHZoOLEobZXs72m7ctz+iW25yjN0uajqkGFxo1UXy40NBroL6kvaeR6pTT918lchTJXZlGFErfk0eh40q+JkufVeNK67cq5nZKL84XAQaQYoltKOgKYF1iEZJhuyPmuKfYrn88BnJU39k5lhhfvtvkYk6/7kozsq8Arhfn0DYCRticCSLqU9GvmOtJUeGlz8CWF9teQdAJpxNyX9MsIkmHeF6aLTHxY0vptJbtLGkL6Tpck7ccaV54pxB96Lq0ZVba0pjoy1lSDDqI3yBQ2yqh+XjifShotzgZ8kNddK1GuIFQSrJ8m6Qt7eqDXabmfFecN8gj2HGCA7deUosQU1YZKfZvKjOf9BfA2sHbu52el6oCTbP+5rI3+Zf1tzRxG6TmGA4Nsj83TCANbKFNUToJZ1ZNmQtKypJmADWy/n6fXWywTBEHQmdg0pQB+o+mwJ7T9EcmLajeY7uCzdjuqfJY0Al4hX5eUhkrGY1IeHZcrDVViQeCt7Bq9D2lkDGn0+ONcD5KWlvSVCuUfIU2xLpbVOvZkhurRbIU+7EUSZQCYH3hL0hzA3oW67iJPY0vqkz3N3ga+ImlRSXMx83R2JRYgGf0P8xrtt2u+gSAIgqDhdLTn7d7AuZKOIk25Xg6MbUtFtj+T9CPSdHLJUWlYdlQ6jzRNPCGn1+Ic4Ops8O8hj0Jt3y5pVeChPGU9GfghaZRb7Mtbkn6dywq42fb1+fYnwOqSRpMUk/bI6b8lGeNXcl/nz+mHAH+RtH9u56e2H8qL6I8AL5N+ULT0bsZKGkOa9n6JkCgMahDbZoLOR0xr1SRf90QzZlmD3saAAQM8alSoGQZBUBtJo2vsB22R/mvM76OuqbYa2HYOWPmBdvWr0XSrPaIdiSQDl9jeJ1/PDrwFPGL7u3kd9A+kLT5zkDx497X9aZX61gUeB7a3fVtOWwa4iCSeMY2k/PGnKuUvIE37vlPaWpPTFyE5QvUnjcx3z+uoC5Ico/qRvtdTbf+tzS8kaDi/GlfPykTXc/JaV3V1F4IeiIk11d7GJyQP3Xny9bdIBrTIFdmlenXgP8yY2q3EnqT11D0LaV8Cv7S9Kmkf7kHl0RIKDCdFRSjnV8Bdtlckrcf+KqcfBDydt+YMBP4vb1wOgiBoCqYyW8OPZqP5etS13ELaKwvJGF5WKVMexc5HEneodF8kZ6XBpA3Fc0Nai7X9eD7/mDTaXbpSHbbvIwlklLMTaTsR+bO0P9bA/Lntvrnsl1WeMwiCIOgAYvp3Zi4HjpZ0I0lE4QJg88L9PSRtRtoH+jwz9sKWsylpT+2/surRDszYrwpM36azLrOGjavFErbfgukOUyXv5LNI0ltvkpyg9qgl/BzMTGsUidrCS5+0VdejPvb66zc7tP4gaA9GTOuiIOWdSYxUC9geR1qr3BO4uUKWK/K+26+SPHgPr1LVniQDTf4sTgGTt+xcDfxP3nrUCLYDngCWIqk8nZW358xEKCoFQRB0HDFSnZURJH3hgcCilTJkXeEbgJ9J+gNJralU9liSTOP3JB1J2nKzqKT5bX+c96leDVxq+xqY7sBUGvUOsz2shf69LWnJPEpdkiS5CCmQ+clZNONFSS8Dq5Dithb7HopKVehondvu4qgUBB1FM66BNpowqrNyAfCh7fGSBraQbzPgX1lacLqfuKTtgLG2tyukXUgS4b8E+CvwjO3TSvdtv1asowYjgP2Ak/NnaX/sq8DWpDiAS5BkHV+qs84gCIIOxXRZkPJOpec/YSux/Xq1bS6kNdUnJI0jrYdWElDdkxQqqMjVJHWlTUkKTlvlep6QtEOlhiRdBjwErCzp9SwOAcmYfkvSCyQP5ZNz+vHAJlk7+S7gf21PqueZgyAIgsYQ4g+9mBB/CIKgXtor/rDMGgv6kCu/0cguAXD4areH+ENHIGkqyXlodpK03z62P6hRZgJJiL9NI7os3j/Z9qltKV+oZxDwvFsZcFzS94DVbJ9cM3MQADve/7OaeW7Y/MxO6EkQ9Ex60vTvlCzMsAZpj+ZBjag070ntaAaRQrXVjaTZbY8IgxoEQXegtKba6KPZaL4eNYaHyKIKkgbmfafk67MKEdwBDpf0aD5WyHmGSzpN0j3A7yUtIuk6SeMkPVwWAHxtSXdLekHSAbl8X0l3SXpc0nhJOxXa3zfXM1bSxZI2Ab4H/CGvsS6fj1sljZZ0v6RVqvRrsKSz8r0dJT0iaYykO7OzUhAEQdCJ9Jjp3xI5FNvWJC/bevjI9oaS9iUFPC+FWVsJ2Mb2VElnAmNsD5K0FUm/t+StuxZJcnA+YIykm0jbXHa2/ZGkxYCHJY0gjUaPBDa1PUnSIrbfy/dutH1Vfoa7gANtvyBpI1JUna0q9Gtw4TkeAL6Rt/v8BDgC+GXdLy7oUVQTshj/4QuzpG1yxvc7uDdBkJjaC6LU9CSjOo+kJ0jiDaOBO+osd1nh84+F9CvzdhlI22d2AbB9t1Kc0wXzvettTwGm5BHkhsBNwO8kbUESzl8aWIJkGK8qreHankWGMAtDbEIKcVdKnqtKv4p8Dbgi712dk7SuPAuShgBDAPr161cpSxAEQcOx1ZTTtY2mJxnVKbbXycbuRtKa6hkk/dviNzl3WTlXOf+kcF7p55XLPovpewOLA+vb/iI7RM2d66nlbj0b8EFWbqrEJ1XSzwROsz0i768dWilTiD/0DqoJWdTjqBQEQdvpcT8bbH8I/Bw4LKsXvQKsJmmubHC3LiuyR+HzoSrV3kcylGSDNakgL7iTpLklLUpSYXoMWJAUsu0LSVsCX8957wJ2z3lLYdwAPiYHLc/1vqwUQB0l1q7j0RdkRlSd/erIHwRB0KlM9WwNP5qNnjRSnY7tMZLGAj+wfbGkfwDjgBeAMWXZ55L0COkHxp5UZijwtyz68CkzG61HSdO9/YDjbb8p6VLgBkmjSHq8z+Z+PSXpRODevAVoDCmSzeXAeZJ+TopuszdwrqSjSLFbLwfG1njsoaQp4zeAh4Fla+QPgiAIGkyIP/RiQvwhCIJ6aa/4w1KrL+QhVzQ+ktKHlKIYAAAgAElEQVSxa44I8YcgCIKgt6GmnK5tND3WqEoyyXHnl/n6MKCv7aGSDgQ+tX1R3pZyu+03q9RzHHCf7TuVYqMeZnuUpJuBvWqpNtXZ13WApWxXCjfXUrkBwL62f97ePgRdT/+LmkPHY8K+v+rqLgRBt6XHGlXgc+D7kk4qlyEsC602GHiSFNx7JiT1sX10pcptzyKEr7QHRm0IDr4OMIDKMVwrkhWVRgExfxsEQdOTFJV6/j7VnjwW/5K0deQX5TckDZV0mKRdScbs0qxmNI+kCZKOlvQAsFtWMZolEGbOt5ik/pKekXQO8DiwjKTJhXy7Shqez4dLGpZVkp6X9F1JcwLHMSMCzh7VFJxyv/8i6XbgIhXUoiRtKOnBrKj0oKSVG/w+gyAIghr05JEqwNnAOEmnVLpp+ypJB5OndAGy4MJntjfL19vX0c7KwI9s/3ehjmr0B74JLA/cA6wAHE0S9j84l29JwWl9YDPbUzRzvNdngS1sfylpG+B3ZMGKoLmopnb077dfrVrmq78e0kG9CYLOI4KUd3OyTOBFpH2rU1pR9IpWNvWK7YfrzPuPPD38gqSXgFUq5GlJwWlEVnAqZ0HgQkkrkmZa5qjUeCgqBUHQFRj1iunfHm1UM6eTpmX/1ooy1VSL6s1f3KfUkoJTpWtoWcGpWt+OB+6xvbOk/sDISplCUanrqaZ21CyOSkEQtJ0ePxbP+rr/APavkmW6mlEDeVvSqpJmA3Yuu7ebpNkkLQ8sBzxXoQ8tKThVo6ioNLh93Q+CIGg805it4Uez0RtGqgD/Bxxc5d5wYJikKcDGDWrvVyT94ddInsV9C/eeA+4lCewfaPuzLMT/qxwQ4CRaVnCqximk6d9Dgbsb9BxBJxJbWYKg+xOKSp1I9gKeHuKtqwlFpSAI6qW9ikpLrLaI9/r7to3sEgCnr3tFKCo1M9nb909AH+B82ydL+itp642A54HBtifn0G6nk2Kq/qBoLCX9HvhOvjze9izOT1lsvxhubpVcz3VFoYlC/kWBq4ANgOElb+GyOkcAy9leo80vIejVLH/FiW0q9689jmxwT4Kg+xFGtYBSgPOzgW8BrwOPZSP1i9KapqTTSFPJJwOvktYvDyur5zvAeqRtMHORBPRvsT24mM/2PTlPKWLNi8DtLXTxM+C3wBr5KO//94HJ5elBEATNQG/w/m2+Vd6uZUPgRdsv2f4PKTrMTgWDKmAesieu7Qm2x5ECkRdZDbjX9pe2PyFFmKm133VX4Bbbn1bLYPsT2w+QjOtMKAU3PxQ4oY7nDIIg6FTSlprZGn40GzFSnZmlSc5FJV4HNgKQ9DdgB+Bp4Jc16hkLHJNHtfMCW+ZyLfED4LQ29LnE8SSHrKpGOQiKVBOhePOdV9pW37l3tKM3lam2/SgImpXmM/NdS9X9obZ/BCwFPMOMwOYVsX07Scf3QeAyUvDzL6s2Ki0JrAnc1qZOJ0H+FWxfW0feIZJGSRo1ceLEtjQXBEHQJqaihh/NRoxUZ+Z1YJnC9dcoCO3bnirpCuBwaohJ2D4ROBFA0t9JCkobAX/OWY62PSKf7w5ca/uLNvZ7Y2B9SRNI3+lXJI20PbBCv0L8IQCqjwLb6qg0MhyVgiCMahmPAStKWpYkpPADYC9JK9h+Ma+p7kjS2a1KdnhayPa7WQx/LVJ4uS+ZoeFbZE/g123ttO1zgXNz2/1J23YGtrW+IAiCRtNbotSEUS2QxegPJk3D9gEuIE333i9pAdL08FjgpwCSNgCuBRYGdpR0rO3VSbq792dh/Y+AH2aDOgvZCC5DEoQo5yZJpdHrQ7Z3y6PRBYA5JQ0CtrVda702CIKgi1FTOhY1mjCqZeRA4eVxTTetkvcx0hRxefpnJA/getqbQHKQKk8fWCV//zrqiz2qQZuJ/aZB0HbCqAZBEASdwrQmdCxqNE1tVLOC0F358qvAVGAiKSbpm7brGg12BJUUjxpQ53AqyBhKWgo4w/auWWD/MNvflfQ9YDXbEd6kyVn/lu4z+hv97bY5KgVB0ORG1fa7zFAcGgpMtn1qyRmnrfVKmr3aGmczYvtNkjhEefoIYMSsJYIgCJoLG6b2Akel7rxq3EfSeZKeknS7pHkAJC0v6VZJoyXdL2mVnD5c0mk5IszvJQ2VdGEuO0HS9yWdIml8Lj9HLre1pDE5/QJJc5V3RNK2kh6S9LikK7O6UdWyub3fS3o0HysUqttC0oOSXpK0a87fX9KTFdodLOmsfL6jpEdye3dKWqKxrzsIgqB9hKJSc7MisKftAyT9A9gFuIS0B/NA26V9oecAW+UyKwHb5P2mQ4HlSWpHq5EEGnaxfYSka4HvSLqVFBpua9vPS7qI5Pl7eqkTkhYDjsr1fiLpf4FDJZ1So+xHtjeUtG9O+25OXxLYjCSuP4IkoF8PDwDfsG1JPwGOoLbyU9AgqqkTlXj+vZer3lvp9z9scG+CIOgqurNRfdn2E/l8NNA/jxA3Aa7M21kgCdqXuNL21ML1Lba/kDSetIXm1pw+nrRuu3Ju5/mcfiFwEAWjCnyDZJT/mduck2Sga5W9rPBZjFRzne1pwNOtHG1+DbgiqzPNCVT8Ky5pCDAEoF+/fq2oPgiCoO0k7d+eP/3bnY3q54XzqSSh+9mAD2xXElgA+KRSHbanSfrCM4LLTiO9m3r+BQi4w/aeMyUm6cCWcJXz4nO15l/gmcBptkdkZ6ahFRsNRaUOoZZGbXdyVAqCoO0034R0O8jRZF6WtBukqDKS1m5Hlc+SRsClNc99mFWk4WFg01IeSfNKWqmOsnsUPh9qRx9LLEhSgQLYrwH1BUEQNJRpqOFHs9GdR6rV2Bs4V9JRJGWjy0kqSK3G9meSfkSaTp6dJGM4rCzPREmDgcsKTkxH5XXUlsrOJekR0g+bmUa5bWRobusNkqFftgF1Bg0itqkEQdchaXvgT6RlvvMrbUOUtDvp76iBsbb3yun7kfxmAE6wfWGLbc2Y8Qw6iyw1OMD2pK7sx4ABAzxqVMO22QZB0IORNNr2gLaWX2TVxb3d33ZuZJcAuHzj81rsV9Zifx74FiloymMkJ9enC3lWBP4BbGX7fUlfsf2OpEWAUcAAkrEdDaxv+/1q7fWokaqknYFrgFVtVxW9l3QccJ/tO8vSB9IOYQVJNwN72f6gTQ8wo56h5D25ZelFEYgBwL62f96etoLuxWrXDe3wNp4e1PFtBL2TLtoCsyHwou2XACRdDuzEzDGuDwDOLhlL2+/k9O1IPjPv5bJ3ANszw9F0FnqUUSVNoz5Aii4ztFom20fXqqgtwgq2d6gzX/88JdwqiiIQWckphplBEAQtszTwWuH6dWCjsjwrAUj6J2mKeKjtW6uUnUWrvUiPcVTK22k2BfYnGdVS+hFZfGGspJNz2vCCsML2kp6V9ADw/UK5orDCcEnnSronizJ8M4s5PKMkLVgqMyHvW0XSvpLG5XYvLtRTFKBYRNJ1Od/DSmHiSqwt6W5JL0g6IJefLgIhaaCkG/P5hlkwYkz+XLnhLzgIgqA9OG2pafQBLCZpVOEYUtZyJW+m8nXP2UnaBwNJg7PzJS1UZ9lZKuopDAJuzQ5C70laD1gip29k+9M8Pz4dSXMD55HEIV4Ermih/oVzvu8BN5AM+E+AxyStU9gzi6TVgSOBTW1PKmu3KEBxJjDG9iBJWwEXMSPe6lqkPbDzAWMk3dRC354Ftsih67YBfkcSwwi6GbVEJF6ZNKHF+18/YXDD+hIE3YRJNdZ6XyeF1yzxNeDNCnketv0FaQfJcyQj+zrJ0BbLjmypMz1mpEr6dXF5Pr88X28D/M32pwClefECq5AEGl7Ie1QvaaH+G3Ke8cDbtsdnkYanSEIRRbYCrio5IpW1WxSg2Ay4OOe5G1hU0oL53vW2p+Q67iGtC1RjQZLn75MkIYnVq2WUNKT0i27ixIktVBkEQdA4TJdtqXkMWFHSspLmJM1kli/tXUdS1yup5K0EvESKrb2tpIUlLQxsm9Oq0iNGqkrRbLYC1pBk0py4gaupMVSv436JkijDNGYWaCgJRczUpRbq/aQsX7X+lJdvqZ/HA/fY3lkp2MDIahlD/KG5qSUi0RmOSkHQUXSFolKewTuYZAz7ABfYfio7rI7K/jMl4/k0SUzo8BzQBUnHkwwzwHEVBmcz0VNGqrsCF9n+uu3+tpchyfS9B/xY0rwA5dO/pGnTZSUtn68bsV8UUri63bOxr9RuiftI+2pLnseTsoAFwE6S5s51DGTGl1qJovDD4Hb1PAiCoIdh+2bbK9le3vaJOe3obFBx4lDbq9le0/blhbIX2F4hH3+r1VZPMap7AteWpV0NLEUa5o+S9ARwWDGD7c9IOrg3ZUelVxrRGdtPAScC90oaC5xWJetQYICkccDJzKyE9ChwE0nI4fjs+VuNU4CTCp5rQRAETYWhoxyVmooQf+jFhPhDEAT10l7xh4VW+Yq/ef7ujewSACM2P7td/Wo0PWJNNQiCIGh+mnFk2Wh6yvQvknaWZOWg5C3kOy5vOylPn77vs4WyQyUd1lKeOvs6SNJqbSj3PUm/yucHKsViDYKGscqxf6ydKQjaQCn0W0+f/u1JI9WGqSl1AoOAG5lZJqtFJM1eVHmyPaxGkSAIgqCT6REj1Q5QU2qt0lFfSXdJejy3t1OhrpmUlSRtQhKQ+IOkJyQtn49bJY2WdH9ptF1Bgamo8jR91CzpAEmP5TauLnk7B0EQNBMR+q370Gg1pWNpndLRO8DOtj/KG4cfljQCWI0yZSXb7+V7N9q+KvflLuBA2y9I2gg4J/cLZlZgGlzl+a+xfV6u6wTSj4szW/sSg55LLaWmEq9OeJ2B91xfM1+t/bRB0FvpKUZ1T+D0fF5SU5qNOtWUACRdQtpeA0npaJdc7m5JsygdAVPyCHJD0taX30nagiQGsTTJqLekrERuty+wCUkRqZQ8VyFLUYGpGmtkY7oQ0JcWFD+UdDGHAPTr169GtUEQBA3CvcNRqdsb1Q5SU2qt0tHewOKkOHtfKMVLnZuWlZVKzAZ8YHudKvc/qZJeZDgwyPbYPJodWC1jKCr1TuodWa5y7B8ZecwvOrYzQdCD6Qlrqh2hptRapaMFgXeyQd0S+HrOW01Z6WNgfoBc78uSdst5JGntVr6D+YG3JM1R6ncQBEEz0VvEH3qCUe0INaWhtE7p6NKcfxTJqD2b26imrHQ5cLhSqLblc5n9c56nSAF0W8NvgUeAO0ptB0FbeDZGqUEH0huMaigq9WJCUSkIgnppr6LSAisv4Q2GNX4i7e6t/hiKSr2RvNY5wPbBkgYBz9uue59qK9saCky2fWpH1B+0nyPG7tbVXajKKWtf2dVdCHogJfGHnk5PmP7tjgwibbepG0nxAygIgqDJCaOakdRfKch36fqwLLAwUtLvJT0q6XlJm+f7fSSdmsUexkn6WU6fkPeqImmApJFl7VQSfxgpaUC+v1j2HiaLPVwp6Qbg9px2eBZ6GCfp2EK9R0p6TtKdwMod96aCIAjahq2GH81GjH7qY3bbG0raATgG2Ibk5LQssG4OglstZupM2H6wgvhDS0U2BtbKohHbAiuS9sYKGJH3xn5CUpJal/SdPg6MbsNzBu2kXpGFlyY/VTPPD/5aX11B0F1oRgWkRhNGtT6uyZ+jgf75fBtgmO0vobKwQ4O4o1D3tvkYk6/7kozs/MC1JaGLbLQrEuIPQRAEHUcY1Rl8yczT4XMXzj/Pn1OZ8c6qCTsU65m7wv1abZeXKYo/CDjJ9p+LGST9T5W+zEKIP3Qs9YosNLOjUhB0BO4likqxpjqDt4GvZEnCuYDv1sh/O3BgyYGoMP07AVg/n+9Spex08YcKZXZtoc3bSIIWfXObS0v6CkmsYmdJ80iaH9ixRt+DIAiCDiCMasb2F8BxJBGFG6ktonA+8CowLos27JXTjwX+JOl+0si2EuXiD6cCP5X0ILBYC328Hfg78JCk8cBVwPy2HycFBHiCJHxxf63nDYIg6Gx6g6NSiD/0YkL8IQiCemmv+EPflZb0mmftVztjK3l4u983lfhDjFSDIAiCoEF0K0clSVOB8YWkQaTp0n1t/7wV9UwgqRtNqpHve8Bqtk+WdCDwqe2LWt/zxiJpsu2+kpYCzrDd0jpsELSKdf/7j13dhaZizDmhh9womnG6ttF0K6MKTKkQIm0C0PA5TEmz2x5BEuXH9rBGt9Fesph/GNQgCIImobsZ1VnIodkOs/3drHm7LLAksBJwKPAN4NvAG8CO2SEJkqPQlvl8L9svShpOChm3LvB4dgYq6fUOJevpZpWkMSSP3cWBfYFfA2sCV9g+KvftUODHuY3zbZ8uqT9wC/AAKTj5G8BOtqdIOoC0h3RO4EVgH9ufSlqW5KA0O3Br4dn7k0Qk1sjnFwPz5dsH236wbW81CHouL1x/dqvyD3z6+lblr3dbVW+jFPqtp9Pd1lTnydJ+T0gqD/dWYnngO6TwaZcA99heE5iS00t8ZHtD4Czg9EL6SsA2tn9Zoy//sb0FMAy4HjgIWAMYnLflrA/8CNiIZNgPkLRuLrsicLbt1YEPmLH15hrbG9heG3gG2D+n/wk41/YGwL+r9Ocd4Fu21wP2AM6olEnSEEmjJI2aOHFijUcMgiAIWkN3G6lWmv4t55YcLHw80IcZI7vxzFBDAris8FlcRLrSdrWtMEVKqkXjgadsvwUg6SVgGWAzksrRJzn9GmDzXO5l20/k8kWVpjUknQAsRFJLui2nb8oMw3sx8PsK/ZkDOEvSOqStPCtV6nSIPwS9nRV3OqhV+UfGmmpjcBKA6Ol0N6NaD58D2J4m6QvP2DM0jZmf11XOiwpGNdvJ9X5eSC+109I8RzH/VGCefD4cGGR7bA4VN7BKHyvxC5KAxdqkGYjPauQPgiDoVHqD9m93m/5tJHsUPh/qgPrvAwZJmlfSfMDO1BZlmB94S9IcQDGa7z9JgvmUpRdZEHjL9jRgH9IoPQiCIOhEeuJItV7mkvQI6YfFno2u3Pbj2fHp0Zx0vu0x2aGoGr8lKTq9QppWLkkZHgL8XdIhJMWkSpwDXC1pN+Ae6h9xB8FMxBaSoCMwvWNLTSgq9WJCUSkIgnppr6LSvCsu6VX+tH/tjK1kzHdObCpFpd48UgVmEpQQaX2zTVtRJA0Cnrf9dIP6NV14osK9EH8IOoztBgzt6i7MxG2jhnZ1F4KGoF6xpabXG1UKHsWStgNOAr7ZhnoGkYT46zaqWWDiy0r3isIT1QjxhyAIuhO9YWK0NzsqVWIB4H0ASX0l3SXpcUnjJe1UyiRpX0njJI2VdLGkTYDvAX/Ie2iXz8etkkZLul/SKrnscEmnSboH+L2kRSRdl+t7WNJaOd9gSWfl82UlPSTpMUnHF/rRX9KThfP7c38fz30KgiAIOpEYqWZBCVJw8CWBrXL6Z8DOtj+StBjwsKQRwGrAkcCmtidJWsT2e/nejbavApB0F3Cg7RckbURyJCrVXRKYmCrpTGCM7UGStgIuAsr34pbEHy6SVG2TXUn84TNJK5L23zbNOkMQVOKx54bXzDNw4MiaeULFqHvQGxyVwqjOPP27MXCRpDVIa6y/k7QFae/p0sASJMN4VUmM3/Z75RXmIOKbAFdK0/8RzVXIUhSY2Iws7GD77qzGtGBZlQ0Tf5A0hCSFSL9+/SplCYIgCNpIGNUCth/Ko9LFgR3y5/pZoWkCaTQragsxzAZ80IL6U3G7S6WfbpXqb4j4QygqBc3EBisPrpnntpFDO7wfQcdj946RaqypFsjrnn2Ad0liCu9kg7ol8PWc7S5gd0mL5jKL5PSPyftKbX8EvJz3jKLE2lWavY8s6JCDA0zK5YuE+EMQBN2eaVbDj2YjjGpBpB+4AtgvT81eCgyQNIpkyJ4FsP0UcCJwr6SxwGm5nstJkW/GSFo+l9k/53mKJPBfiaG5nXHAycB+FfIcAhwk6TGS8azEOcB+kh4mTf2G+EMQBEEnE+IPvZgQfwiCoF7aK/4wzwpLedlThzSySwA8s/OxTSX+ECPVIAiCIGgQTe2oVFA7mgP4ErgQOD2vG1YrM5ActLwV7QynsB2mGZA0ANjX9s8r3JtACp4+SdKDtmNPag9m2UtP6tT2Xt77153aXtB76A2OSk1tVJl5u8tXgL+T1hSP6dJedQK2RwE152bDoAZB0B0wCqPaTNh+J++xfEzSUNLU9cmkmKNzAWfb/nPO3lfSVcAapCDgP7RtSUcDO5Lilz4I/JfLFpUlbQ2cSno3jwE/tf25pA1IIgzzkeKhbg18AZxLEln4EjjU9j2S+pD2km5H2gpznu0zW1nHQPKIO3saX0ba4vMohW04BR3gvsD1wMKkkf1Rtq9v+xsPOoqBAwe2Kv9bb79ad94ljzqglb0JgqCRdKs1Vdsvkfr8FWB/4EPbGwAbAAdIWjZnXRf4H5L60XIk8QSAs2xvYHsNkmGdaYpY0tykQOF72F6TZFh/KmlOkmfwIbbXBrYBpgAH5X6tSQofd2GuYwiwLLCu7bWAS9tQR5FjgAdsr0vSA66k2lBSgFoP2BL4PxWUJwrPOETSKEmjJk6cWPE9B0EQdATugKPZ6DYj1QIlQ7EtsJakkqD8gsCKwH+AR22/DpC3yvQHHgC2lHQEMC+wCGmryw2FulcGXrb9fL6+kGT07iLtAX0Mpu9DRdJmwJk57VlJr5AlCIFhJbH8LGO4ZivrKLIF8P2c5yZJ71d5L5UUoP5dzBTiD11PayX1OntNNQiCttOtjKqk5UgSfO+QjMjPbN9WlmcgaWq1xFRg9jz6O4fk4PNankIuHxFWm/CvpqLUmvytraOcWgZwbyorQAVBEHQ9oajUXEhaHBhGmsI1cBtpanaOfH8lSfO1UEXJwEzK64+VQqY9C/SXtEK+3ge4N6cvlddEkTS/pNmZWQ1pJdK07HPA7cCBOU9Jdam1dRQp5vk2ad20nGoKUEEQBEEn0ewj1VIEmdKWmouZoWB0Pmla9/G8djiRFNO0IrY/kHQeaYvOBJITUnmezyT9iCSEX3JUGmb7P5L2AM6UNA9pLXQb0sh3mKTxuX+Ds1PT+aQp3HGSviA5Kp3VyjqKXTsWuEzS4yQjX8lz5VLghqwA9QRZASro/sQWl6DH0AsWnEJRqRcTikpBENRLexWV5l5+aS9z8k8b2SUAXtz9t02lqNTUI9XSdpF8vgNpO8rWwHrA87afbmt9DezjBLIQQzvr+Y3t37Wh3PnAaa19F0HXs9zfW/11dwov7fWbru5CEHRbusWaat47eiawve1XSdO8q3VtrxpOq/+SSepj+ydhUIMg6A6k8G+NPepB0vaSnpP0oqRfVbg/WNLEUnAVST8p3DtF0lOSnpF0RqWtikWa3qhK2hw4D/iO7X9J2gT4HvCH/PDL5+NWSaMl3a8Uwg1Jy0p6SNJjko4v1ClJf5D0pKTxea0TSQMljZR0laRnJV1aeoGStlaKQDNe0gWSikHHD5f0aD5WyPl3lPRILnOnpCVyel9Jf8v1jJO0i6STmREt59Kc74e5vick/VlJUAJJkyUdJ+kRYOPc3wH53rlKe1CfknRsh34xQRAE3YD8t/Ns4NukwdiekioNyq6wvU4+zs9lNyHpHKxFEhPaAPhmS+019fQvSSnpemCg7VLotQcljaCg1SvpLuBA2y9I2ojk/LMVabr4XNsXSTqoUO/3gXVIAb0XI6k03ZfvrQusDrxJimO6aXb+GQ5sbft5SRcBPwVOz2U+sr2hpH1z2ndJ+2K/kZWcfgIcAfwS+C1JtGLN3PeFbV8t6eCCJOOqwB7Aptmb9xyS9+9FJDWmJ20fnfMW39eReU9sH+AuSWvZHtfGdx/USWsVkkq89U79SklFljzqJ7UzBUGTYbpsS82GwItZPAhJl5NCcdYzw2fSzpE5Sdsf5wDebqlAs49UvyDJCe5fLYPS9phNSB67TwB/BpbMtzclyftB8hwusRlwme2ptt8medRukO89avv1LNpfEo6oJAqxRaG+ywqfG+fzrwG3Za/ew0mGGpLH79mlgrYrCTlsDaxPMvZP5Ovl8r2pwNVVXsfu2UN4TG5vll9jCkWlIAi6AgNW44/aLA28Vrh+PaeVs0uePbxK0jIAth8C7gHeysdttp9pqbFmH6lOA3YH7mzBkWc24IPSKK8CrRVcmEU4okb+8jZK52eSHIhGKAlSDC20XWslQMCFtivtpfgsB1GfuUCSaDwM2MD2+0qRd2YRfwhFpcbTWoWkEs3qqBQE3YzF8mxiib/kv3MlKv39Lv/bdwNpoPW5pANJA6et8nLeqqRBEsAdkrawfR9VaPaRKrY/JU2n7i2pNGL9GJg/3/8IeFnSbjB9vXTtnO+fwA/y+d6Fau8D9pDUR0lUYguSUH01qolClNij8PlQPl8QeCOf71fIeztwcOlCUknI4QtlIQuSLOKuSpF5kLSIpFpiDgsAnwAf5vXbb9fIHwRB0Kl0kKPSJNsDCsdfypp9HVimcP010vJeoV9+13ZpQHUeaaYQYGfgYduTbU8GbgG+0dIzNr1RhaSdC2wPHCVpJ+ByknPQGEnLkwzm/pLGkvR8d8pFDwEOkvQYyciVuBYYB4wF7gaOsD2TRm5Z+58BJVGI8aQR9LBClrmy49AhwC9y2tCc/36guN3mBGDh7CQ1liR+D2n0OE7Spdmb9yjgdknjgDuYMaVdrY9jSdO+TwEXkH5QBEEQ9HYeA1bMjqtzkgZaI4oZJBX/vn4PKE3xvgp8U9LsedDzzcK9ioT4Qy8mxB+CIKgXtVP8Ya7llvbSJxxUO2MreXnvI2v2S0nn4HSgD3CB7RMlHQeMykt0J5GM6ZfAe6SQn89mp89zSLOZBm61fWhLbTX7mmoQdDl/enabru5Cp3LIKnd2dReCHknXBSm3fTNwc1na0YXzXwOz+LBk/5X/ak1b3WL6t6ciaVHN2Gz8b0lvFB3OhioAACAASURBVK4/LeTbQdILkvpJ2kLS45K+1Iywd6V8++V8L0jab9YWgyAIgo4kRqpdiO13SftlUQpFN9n2qfl6cv4sqUlta/tVSbMBg0mevtNRioRzDDCANE0xWtKIKlt2giAIOp9esNoYRrWJ0Qw1qR1s/wvA9oR8b1pZ9u2AO7JTF5LuIDl3XUZQlXqEG17/dGzNPD+7aP2aeYIg6PmEUW1eZlGTqkFdG5wlDQGGAPTr168B3QyCIKiDXhKkPIxq81JUkzqkjvz1bHAO8Ycy6hFu6G2OSkEQtJ1wVGpeSmpSG0iqJ4JNzQ3OQRAEXYo74Ggywqg2MVXUpKpxG7CtpIWzStO2OS0IgqBJUAcczUVM/zY5OerM9sB9kiaRRp/XAgsDO0o61vbqOd/xJPWQ/2fvvMPsrKr2fT8UCRCIUlR6pEkTggSQJkPRn/AhBAERIxBAKYIIiFQ/qmhABJRePkoQBIEEQq8JvSQkJCEIoUsEIaEaOsnz+2Pvk7w5c87MmZkzfd3Xda55z353e88ks87ee61nAZxUcloK2kbEbQZBUCthVLsItk8oe9+3cP0a8I3C7WWpgO1LSRKFQRAEXY8uuF1bb8Ko1hlJM0oGMUtj/YWUum1v4BfANFJO1EnA77LOb7H92cBeRaNadn9X4FiS3Natto8o3PsxSXPYwATbP63v0wW9gf7nn15TvVcOOLz5SkHQy4gz1XaiINrwA9ulbNRn5qzyqwDXAvflLDmlNgOBLzfR5+LAn0jJ0tcEvpbHQdIqJJmtTfK9Q9rjuYIgCFpNOCoFraEg2vA/JdGGcmxfS0oD99PcZl6SwTyiUv3MisAU26Xs4vcAO+XrXwDnlhSUbL/V1ucIgiCoG52XpLxDie3f+tMS0YZxwGr5+iBgpO03pKr/UF4AVpPUnxRCMwj4Ur63KoCkh0lbwyfYvqOVzxAEAPznzPOr3mu49paq91qbuD0IujthVOtPS0QbBCBpaWAXoKGpyrbflXQAaet4Vh5nxXx7PmCV3MeywIOS1rL93lwDhqJSEASdRG/INBpGtf6URBvukXSM7T80UXddYGz+uTLwQl6lLiTpBeCbwJO57kjbx9m+GbgZZhvImfn+VFKG+s+BlyU9RzKypRAbIBSVgpbx9UMPqHpvdDgqBUEjwqi2A7Y/krQdabX4pu3/K68jaSeSQMNv8hnp1wv3ZtheOb8dUNbuq7bfygIPvyQZcIAbgd2AyyUtQdoOfqnezxYEQdBqesHX+DCq7UQF0QaAQyX9jBRS8zSwZcHpqFb+ImmdfH2S7Sn5uqSo9Axp9frbnFouCFpEhMoE7UYXdCyqN2FU60wTog03kWJIW9RHhXu7VSk3cFh+BUEQBJ1AjzGqJdGFHK95JrA68B7wAXC87Qck7QCcTDr3/AI4xPZDuf2ewO9yd7+3fUWFMRqAw21vl9//Hlgf2J60UlwK+LjQx/UV+tgmz2FhkqPSLbYPz0nKjwdWsf1CrnsocAawvu2xkr4EnENyRpoFHGv7hlw3hB86gW0eqCWBUPfi9u/+pbOnEPRQFNu/3QtJfYBbSYZvZC5bCxgIPADcS3L4saS1gX+QQlQWIxm0gSSj9KSkkaWYzypjHQtsQkog/ml2MBpse2wTbdYiGcX/sf2spPnInriZScBPgN/n9zsDRcWlY4G3bK8qaR5gsdxvUfjhXUlfbfbDCoIgCOpOTxN/GAw8WjKoALaftn15vp6Rt0khrRRL1/8PuNv2O9mQ3g38oNogkn4DbAv80PbH1epV4AjglFL8qu0vbJ9XuH8jsEMeY0XgfZKsYYm9gT/mtrNsl85qQ/ghCIKuTXuoKXXBlW+PWqkCa5IEFaoiaUeSYfoq8D+5eBngtUK1qbmsEpuQQl3Wsz2j7N5VkkpGdqsKjkJrAX9uYnofAK/lFe0OpHjUvfK8S/KFJ+dt6BeBg2y/SQg/1IWGhoYWt5n4/gutGmvDv+zUfKUgCLodPW2lOheSRkh6WtLwUpntEbZXI6kRnVyqWqF5te9AL+T6369wb3DW9h3QBs/ba0hbwINIKd5KzEcSdXjY9reBR4HTC/dKwg+7AZcUjPBcSNpX0lhJY6dNa6njcRAEQWtpB4nCLuhN3NNWqpOB75be2N4xi9Q3SruRHZdWyjGdU5lbzWhZYHRe1R6fy36ef75J2ma+V9LbtkdVm4ykA0lbs5C2iycD6wETmniGm0kawGNtf1CQLHwb+Ig5hvY6kmoT1Cj8kJ87xB+q0BppvZ7oqBQE7UYv+IvT01aqVwObSNq+ULZQ6ULSyspWStK3Sbq5bzMnxvMrWVTh+8CdeVVbWnnOdkDKsaE/Av4maS5xhiK2zy20f51kLI+RVNqunUfSYWVtPgaOBE4pKzfJ4Dbkoq2Y48R0I7BF7jOEH4IgCDqJHrFSzV60n9r+OCsZnSHpLNKq8r/M8abdCdhD0uek0Jdds7F6R9LJzFnZnWT7nabGtD1G0l7ASElb1DJP2xMlHQL8XdJCpO9tt1aod02VLo4ErszPNo183koIPwRB0B3oBStVuQcoHGeFoYttb9DZc+lODBw40GPHVo0ACoIgmI2kJ20PbG37BVZYzksdVf/jkld/+ds2zavedPuVqqT9gYOJpNxBEARdm+6/hmuWbm9UbV8AXNCatpK+RlJf+g7wLvAZcBopvGYP2wfn8JXPbD+S2wwiJQp/piV92h6Rt3wvBtYmeRC/B/ygPDRH0mrAZcC3SapJp+fyPiQRiwVIv7vrbR+f74m0zb0LaQv4fNt/bc3nEvRuVj71zM6eQiNeOPLQzp5C0FZKScp7ON3eqLaWbIRuBK4oSfpJWgHY3vYIUko2SI5BM0i5SyGFutzC3EpHzfaZq/waeNP2t/K9b5Lyr5bzDmn1Pais/FOSCP8MSfMDD0m63fZjwBBgOWA127NCVSkIgqDj6Wnevy1hS9IKdPYq1/arts+W1CDpFkn9gf1J2WWekrQ5yUD+Kb9fqdY+89ulgH8X7j1n+9Pyidl+y/YYygyuE6VV7fz5VdpQOYDkYDWr1EeLPo0gCIJ2Rq7/q6vRa1eq1KC+ZPsVSRcAMwpbsCNJIviNxPJr6PNS4C5JO5N0iK+w/XxLJi1pXlLi8pVJ0oSP51srAbvm2NppwMEt7TsI2pN/X3huq9s23H5Tm8ZuTQxyELSG3rxSnQtJ50qaIKmRYEK9+rT9FLAiKV51MWCMpNVb0qftmbYHkAQqNsiShpDOWT/JXnAXkwx4pTmFolIQBJ1DaP/2aCaT4lYBsH1gFk6oOcZE0nIkQQZIzlLN9pm3b4cDwyXNAraVtCUF5aUsFNEktt+TNJok/P80SVXphnx7BMnRqVK7UFQKOoVl9juw1W1Hh6NS0E3ozSvV+4A+kg4olC1Uod5/gUUqvbf9WkEx6YLm+pS0SVZsIudGXQN4tYLyUkUkLVnS9JW0ILA18Gy+fSPpTBdgc2BKM88fBEEQ1Jleu1LNOVUHAWdKOoJ0DvkhSbWoyM3A9TnB+a9IgvcXSzoY2Nn2iy3ocyXg/OwlPA9JTekGypD0ddLqdlFgVlZhWoPk6HRFPledB/iH7Vtys6GkLDmHkryVf17ebxDUQoSvBO1FV3Qsqje91qgC2H6DlBGmEqNznSmkuNIia7SmT9vDgGE1zOs/pDPTciYC61Zp8x5zUtkFQRAEnUCPN6rZG3Y4sHopOXgnzGEIMND2QVkB6qNsYIOgy7HGsdXFH545JVaxQRvoBeIPveFMdTfgIaqvSOtCFvVvFtsXhEENgiDomfRooyqpL7AJKe/oT3JZg6TRkq6X9Kykqwrp4I6TNCYnNr+oUL6SpDskPSnpwSwjiKTLJZ0haRRwqqTFJN0oaaKkxySVbxsj6QRJh+fr0ZJOlfSEpCmSNsvl80r6U57LREn7lZ5H0r2SxkmalM95S/3+LPfzlKQL87lrEARB16A9wmm64BltT9/+HQTcYXuKpHeUcqhCOpdcE3gdeJhkeB8CzrF9EoCkK4HtSI5KFwH7235e0obAeczxtF0V2Nr2TElnA+NtD8phMsOAqvlWM/PZ3kDStqSE6FuTvgS8b3t9SQsAD0u6i6RJvGNOXr4E8FgWo1gN2BXYxPbnks4jJVKPFXHQal79W2OxhoaHG4swhLBCUDNd0AjWm55uVHcDzsrX1+T3twJP2J4KIOkpoD/JqG6RvXYXIokzTM6r0I2B6/LCFZLQQonrbM/M15uS41Rt3ydpcUn9mpnj8PzzyTwPSEnS187KSwD9gFVIsah/kPRdYBawDPA1UsLy9UhiEgALAhVlCiXtC+wLsPzyyzcztSAIgqAl9FijKmlx0mpyLUkG5iV9T7qNJExfYiYwn1IGmPNIDkWvSToB6EPaIn8vqxhV4sPisBXuN/fdrDSXmcz5fQj4le07y55pCLAksF5ekb6S5yiS5OHRzYwV4g9Bzazws8ZiDaPDUSloA70hpKYnn6nuDAyzvYLt/raXA14mrSYr0Sf/nJ7PYncGsP0B8LKkXSBlolFKil6JB0jbriiljJue27eUO4EDlDLRIGlVSQuTVqxvZYO6BbBCrn8vsLNyZpp8trtCpY6DIAiC9qPHrlRJW71Dy8puIGVzebG8cpb9uxiYBLwCFDWAB5NEG35HygxzDTChwpgnAJdJmgh8BOzZyrlfQtoKHpedpaaRzoevAm6WNBZ4iqymZPuZPLe7JM1Dym5zIPBqK8cPgiCoP71gpSq7FzxlUJGBAwd67NiapY6DIOjFSHoyJ+xoFQsst5yX/XX9jw9e+u1v2jSvetOTt3+DIAiCoEPpydu/TZL1dc8C1ic5C71CEqXf3vZ2FeqPBg633ealXY5HvYC0TbsbyYN4raZbVexnILCH7YPbOqegbcz6z6qdPYW6Mc/XIxdDUH+6alLxetMrV6r5nHIEMNr2SrbXAI4hhafUa4ymxBcGA6dnj+KPWzuG7bFhUIMgCLoOvdKoAlsAn+d0bcDsBOIPAn0rqS0VkbRbVjR6WtKphfIZkk6S9DiwkaStJI3PdS+VtICknwM/Bo6TdFVZv9WUlHaUdE/2PF4qqy99XUkd6pZcZwNJj+TxHpH0zfb44IIgCFqNVf9XF6O3bv+uRRJbqEQ1tSUAJC0NnEoSW3iX5HE7yPaNwMLA07aPy3GvzwNbZUWnYcABts+StClwi+3rJfUvjF1RScn2CEk7kTx6fwAcb/s/ynKJmWeB79r+QtLWwB8oJEwP2k5DQ0P1m59NrXrrvuGVEg4FQS+kk7Z/Jf0A+AtJr+AS20PL7g8B/gT8OxedY/uSfG95UkTGcqQn2Nb2K9XG6q0r1aZ4wvZU27NIYSv9y+6vT9o2nmb7C1KYy3fzvZnMyY/6TeDlnDoO4IpCvWp8H9gjqzw9DixOUlKClMv1aOBT23+v0LYfSfXpaeBM0heDRkjaV9JYSWOnTZvWzHSCIAi6N/ko7lxgG1Lazt0kVUrfea3tAfl1SaF8GPAn26sDG1BFra5Eb12pTiaLO1SgkdpS2f2m9hs+KUgWtmZfoqKSUmYZkjTh1yTNk41+kZOBUbZ3zKvf0ZUGCEWl1tOUxm1PclQKgvaikxyVNgBesP0SgKRrgB2AZ5prmI3vfLbvBrA9o7k2vXWleh+wgKRflAokrQ9sXkPbx4HNJS2RvwHtBtxfod6zQH9JK+f3u1epV6SikpJSWrnLgJ8C/wQOq9C2H3O2LobU8BxBEAS9gWVIyUhKTM1l5eyUfVmul7RcLlsVeE/S8Oyv8qdmnFB750rVtpWSl58l6SjgE+aE1DTX9g1JRwOjSCvL22w3St1h+xNJe5G2ZOcjKTRdUF6vjGpKSr8BHrT9YN4aHiPp1rK2pwFXSDqM9KUh6EAiDCUIaqB9VqpLZJW5EhflHbkStWiy3wz83fankvYnHddtSbKRm5F8bf4FXEtatPxftcmEolIvJhSVgiColbYqKvVZdjkvd2ClTba28cIxhzU5L0kbASfY/n/5/dEAtv9Ypf68wDu2+0n6DjDUdkO+tzvwHduNs01kOn2lKmlZ0iHyGqTt6FuA39r+rE79DyFlnjmoBW2Osf2HZur0J3nwtli0oUI/G9u+uhVtH7G9cVvGD4Jy1j70zM6eQrdn4pmRzacRnSf+MAZYRdI3SEdkPyEdpc1G0lK238hvtycds5XafkXSkrankVavTa5EOvVMNW9xDgdutL0Kaf+6L3BKZ86LJATRUfSn7BfcHKU9/TCoQRB0K9wOr+aGTFEaB5F8Vv4J/MP25KwpsH2udrCkyZImAAeT/VKy4+nhwL2SJpG2ki9uarzOdlTakuQxexnMfoBDgb0lHSPpUgBJ38pCCwtJel7Skrl8HkkvZKehyyVdIOnBLI5QlBpcWtIdue1ppUJVEHGQNBRYUNJTJXEGSYflOk9LOqTQ73ySrigcbi+U6x+nJODwtKSL8pcHJK2sJOIwQdI4SSuRMulslsc7VNUFIBokjZJ0NSmTDpJm5J99Jd2b+5wkaYc6/56CIAi6LbZvs71qVtA7JZcdZ3tkvj7a9pq217G9he1nC23vtr227W/ZHtLcLmpnb/+uSZkIg+0PJP2LtA18TnYoOhbYz/ZHkv5Gkvk7C9gamGB7erZb/UkevCsBowqetwNIB82fAs9JOpsULlNJxOEoSQeVkpJLWg/YC9iQ9C3lcUn35zbfBPax/XD+AvBL4HRS4PBJuf2VwHakg/CrSPvzI5TEIeYBjiJpCm+X6+9LBQGI/BwbAGvZfrnsc/wE2DF/dksAj0ka6TgwD4KaeOm6c+vaX8P4Rr6LbaapkK5uQy/4i9TZK1VR+WMulQ8BrgTut/1wvncpsEe+3psUalLiH7Zn2X4eeAkoKQ7da/t925+QYpNWoGkRhyKbAiNsf5hjlIaTvMEAXivM62/MSYC+haTH83bBlsCakhYBlrE9ApJ3sO2PKozXlADEExUMaunz+oNSHtd7SO7iFXWMFeIPQRAE7UZnr1QnUyalJ2lRkhzUiyTjNQNYunTf9muS3pS0JWn1OLjQvNxAl95XEnSoVZyhqXqNxssr0PNIzlGvSToB6NPC8RoJQEhqAD6s0mYwsCSwnu3PJb2Sx2w84RB/CIJGrLhLVWfOVjE6HJUqEllq2p97gYUk7QGzHXD+DFwOzE/SavwusLikogLSJaSV4T8KCkYAu+Rz1pWAFYHnmhi7KRGHz5UFGIAHgEH5PHdhYEeS8D7A8tldm9z+IeYYs+mS+pKVm2x/AEyVNCg/6wL5DPa/wCKFeVUUgGjiOSAJP7yVDeoWpJV4EARB0MF0qlHNZ347kozh88AU0vngMST92vOydu4+wFBJX81NR5K8hC8r6/I5kmG8Hdg/b/dWG/sNkpbuKGACMK4g4nARMFHSVbbHkYz8EyRDfInt8bneP4E987brYsD5tt8jeYdNIolJjCkMuzvJy2wi8AjwdWAi8EV2XjqU9IXhGZIAxNPAhTS/o3AVMFApAHowSc0pCIIg6GC6pfiDUnLuM21vVii7nJz5pdMm1s0I8YcgCGpFbRV/WGY5r7Bf/cUfphzftPhDR9PZZ6otRklW8ADmPksNgiAIgk6n2xnVnAdvaIXyIW3pV9JM0pbt/MAXJO3Hsypkgym2WRr4q+1qGW86BaWA5jXKcwYGQS2seNYZnT2FHstLh9R/pdZt6DxFpQ6l2xnVduTjQmzqV4GrSQ5Ax1drYPt1qqeQ6zRyQPPIzp5HEATBXIRR7Z3YfiuLMIzJITHzkFbHDcACwLm2L1RB/1dJY3h7YCGS+MQI20dAUm4iOV8JuNX2kbn8B8AfSNnop9veStJipFjcFYGPgH1tT8yexGcDA0n/NE+0fUOVPobQQr3jIAhq441zzmt124YbW/9dt0eIP/QCwqhWwfZLkuYBvkpKaFtJ5aj8e1fNyk3AwyQv4e/afjkbU4ATgfG2B+VY3GG53//Nc/gWgKSvKMk1VuqjKvnLwr4Ayy+/fOs+nCAIgtYQK9VeT0mw4fvA2oVY2X4klaPyJJr32n4fQFJJuWlxsnJTLi8pN80EHigpJNl+J/exKVkQw/Z9khaX1I8kyfiT0kC235X0wyp9VCXEH4KgbSx10C9b3XZ0bz5T7SWEUa2CpBVJhu8tqqsc9S9r1hLlpqYkGstxlfrV+giCIOhSiN7hqNTZikpdkrytegFJGN+0TuWoRDXlpkdz+Tdyn6Wt2wfI4UJZmnB6VmO6i5S+qDTHrzTRRxAEQdAJxEp1DgtmEftSSM2VQCm24BJSBpxxSulwpgGDaunU9htKmeZHkb6s3VZSbsrnm8Pz2e1bwPeAE4DLsurSR8CeuavfA+dmlaWZJEel4VX6CIJW06vDPoL2pResVLulolJQH0JRKQiCWmmrotKCSy3nb+xT/y9s/zylaykq9YjtX0kzlZJ8l15HNVO/QdLGNfQ7SNIahfejs0Rih6OUhH3n8nlIuk3SlztjTkHPpP+wofQfFrohQdAaesr272zhhhppIKWUe6SZeoNIydKfaeW82h3b23b2HIIgCGqiF2yM9oiVajUkvSLpREnjJE2StFr22N0fODSvajeTtIKkeyVNzD+XzyvZ7YE/5Xor5W53kfSEpCmSNsvjzCvpT5LG5D72y+XK5U/n8XfN5Q15tXm9pGclXZXPapG0nqT7JT0p6U5JS9XwjEvk6xtzu8n5rDUIgiDoQHrKSrXkZFTij7avzdfTbX9b0i+Bw23/XNIFwAzbpwNIuhkYZvsKSXuT9HwHSRpJIfNNtnvz2d5A0rYkCcOtSanpKolDfJsk3LAOsARJoemBPK91gTWB10lCEJtIepykmrSD7WnZCJ8C7F3j57C37XckLZjHusH22y35IIPuR0NDQ137+8+b/0r9XnpHXfstEcpAvZhesFLtKUa1qe3f4fnnk8CPqtTZqHDvSuC0JsYq9tc/X1cTh9gU+HtOpP6mpPuB9YEPgCdsTwXIXwj6A+8BawF3ZwM+L/BGE3Mp52BJO+br5fIc5jKqoagUBEFn0RviVHuKUW2KkiBDSYyhFpr61Vfqr5o4RFPnndWEIibb3qjGeRbHaiCtmjey/ZGk0UCf8nqhqNTzqPfKr+SkNHqPJv39giCoQI8+U22C/wKLFN4/whwJwMHAQ1XqVaOaOMQDwK75zHVJkjzhE0308xywpKSNcj/zS1qzxmfqB7ybDepqwHdqbBcEQdAxuB1eXYyeYlQXLAupaS4e4GZgx5KjEnAwsFcWXNgd+HWudw3wW0njC45KlbiE5CE8LoszXEhaeY4AJgITgPuAI2z/p1ontj8jpZI7VdIE4Cmg2dCfzB3AfPkZTgYeq7FdEARBUCdC/KEXE+IPQRDUSpvFH76+nFfcs/7iD8+c1rXEH3rDmWoQBEHQBQhHpW6OpJnApELRIFJoyx62D1ZKQD47tKYV/fcnJylv4zz7AxvbvroVbR+xXesWcRA0y+rHnVn3Pv950qF17zMIuiI92qhSOdTmFaDmPU9J8+aQmPakP/BToGajWppXGNQgCLoNsVLteeTQk8Ntb5eL1pF0Hymu8zTbF+c6x5NiRAcAa0g6jDkiDJfYPitfzyfpCpKYwxTSKvgjSccBPwQWJHkX72fbklYmpZVbkhRKswswFFg9x6teAfw1lzUACwDn2r6wyrxm2O4rqS9wE/AVUqad35Wy4QRBe/LqsHObrdPwQPP/FEMUIugJ9HSjWlRaetn2jhXqrE0KP1kYGC/p1ly+AbCW7ZclrQfsBWxIiiV9PAs5vAt8E9jH9sOSLgV+CZxOysV6EoCkK4HtSF7HVwFDbY+Q1IfkgX0UBUOfBRoqKTTNNa+y5/gE2NH2B1m28DFJI13miRbiD0EQdBZxptr9qUVo/ybbHwMfSxpFMlrvkRSPSoZrU2CE7Q8BJA0HNgNGAq/ZfjjX+xspPOd0YAtJRwALAYsBk7MgwzK2RwDY/iT3Vz6nagpNn5XNq4iAP0j6LjALWAb4GjBXCE+IPwT1ZoU9Dmy2zug4Uw16CT3dqNZCuWEpvf+wUNbI6jXVPq9AzwMG2n4tO0T1aaafItUUmhrK5lVkMGlLeT3bn0t6hQqKSkEQBJ1GL/ga31PEH9rCDpL6SFqcdIY5pkKdB4BBkhbKSkk7Ag/me8uXFJCA3UhqTCVjNj2fde4MYPsDYKqkQQCSFpC0EI2Vm6opNDVFP+CtbFC3AFao8fmDIAjan/ZQU+qCRjpWqkk28FZgeeBk269LWrVYwfY4SZczR2LwEtvjcyjMP4E9JV0IPA+cnx2VLiaF87zC3IZ6d+BCSScBn5MclSYCX2QVpcuBv5A8gsfllHDTSOFATXEVcLOksSQlpmdb9jEEQSLCX4Kg9YSiUi8mFJWCIKiVtioqLfS15bzy4PorKk06MxSVei0FMQqRwmkOsv1IK/rZH/jI9rC8gr7F9vXZEepw22Ml3Qb81PZ79XuCoDewxfdP7ewpdAij7jqys6cQ9EDCqHYss72RJf0/4I/A5i3txPYFNdRpKu1cEARBx9MLNkbDUanzWJQU54qkvpLulTRO0iRJO5QqSdpD0kRJE3K8K5JOkHR4U51LeiXHqwZBEHQJ5Pq/uhqxUu1YSmIUfYClgC1zeUXhBmAN4FhgE9vTJS3WKbMOgm7E+LEX1lSvoeH2muqF0lPQEsKodizF7d+NgGGS1qK6cMOWwPW2pwPYfqetEwhFpSAIOo0uuLKsN2FUOwnbj+ZV6ZLAtlQWbhB1/mcYikpBT2fdgfvVVC8clYL2IM5UOwlJqwHzAm9TXbjhXuDHWZiC2P4NgqBbE+IPQZ0pCvwL2NP2TEkVhRtsT5Z0CnB/DscZDwzphHkHQRC0jS7qWFRvwqh2ILbnrVI+Hdioyr0rSOngimUnFK6HFK4bCtf92zLXoPcS26JB0HrCqAZBEAQdQ6xUA0kGzrD9m/z+cKCv7RPKlI2GAHfZfr0dg8YBugAAIABJREFU57I08FfbOzdbOQhqYNdH929Udu1GzWqLBEFQhXBUap5PgR9VElKwfYHtYfntEGDpSh1Iqrjt21Jsv17JoEqKL0dBEHR5QvwhAPiCFIJyKEmIYTY5T+oMUiaagcBVkj4mnY/+E7iUlHD8HEnPAheQkpa/COxt+11JBwP753Gesf2T3O9KpHjV5YDTbF+cs+LcYnutvDL+H1LozcKStgduAr4CzA/8zvZN7fB5BF2QhoaGVrV75oMpjftatPkERyGIEASVCaNaG+cCEyWdVulmFrM/iCxmD5AytvGJ7U3z+4mkxOP357RvxwOHAEcB37D9qaQvF7pdG/gOsDAwXtKtFYbeCFjb9jt5tdpIlcllaYhC/CEIgk6jC64s600Y1RrIhmoYcDDwcQuaXgsgqR/wZdv35/IrgOvy9UTSCvdG4MZC25tsfwx8LGkUsAEp3KbI3QWVpWqqTP8pe5YQf+iBtHblGGeqQUfSFbdr602cqdbOWcA+pJVjrXxYQ53/Ia2E1wOeLJyPlv/zq/TPsdj/YOaoMg0A3iRtDQdBEAQdRBjVGskrwn+QDGsl/gssUqXt+8C7kjbLRbuTBB3mAZazPQo4Avgy0DfX2UFSn6ym1ACMaWaK1VSZgiAIOp/2UFOqceUr6QeSnpP0gqSjKtwfImmapKfy6+dl9xeV9G9J5zQ3Vmz/tow/AwdVuXc5cEHBUamcPfP9hYCXgL1IMoV/y9vDAs60/V4+j30CuBVYHjjZ9uvZUakaFVWZgqApYqs36Onk6Itzge8BU4Ex2d/kmbKq19qu9vf9ZOD+KvfmIoxqM9juW7h+k+S9W3p/QuH6BuCGQtP+Zf08RXI8KmfTKkNPsb1vWR+vAGvl68tJhrx0r6oqUxAEQZegc85UNwBesP0SgKRrgB2AcqNaEUnrkfxT7iBFeTRJGNUOQtKxwE+BmSRHov2At4BrgMWAccDutj8jbd1+T9JQ4Ce2r899rAAMJ61w5wfOtn1B2TgjgRVtr9UhDxb0ONa7/djmK3UgT25zSmdPIagDot0clZbIO3QlLsoOmSWWAV4rvJ8KbFihn52yo+cU4FDbr+Ujuj+Tjuy2qmUycabaAeTcqdsB37a9NrA16Zd8KmnLdxXgXeac155ISgd3dVlXbwAbZ0ekDYGjsspSaZwfkeJmgyAIegvTbQ8svC4qu68KbcrN+81A//z3+R7m6K3/ErjN9mvUSBjVjmEp0i/+U5i9VfsGOQl5rnMFMCjff8X2RNKKdja2Pyv1ASxA4fcnqS9wGPD7dnyOIAiC1tM5jkpTSSI6JZYF5pKTtf124W/rxaRoDEhHagflHNenA3vkHcSqxPZvx3AXcJykKaRvQdcCk4H3bH+R60wlbVM0iaTlSA5MKwO/LWgNn0zapvioznMPgkZMOfJvHTZWw6kPd8g4oRLVYxkDrCLpG8C/gZ+QjuJmI2kp22/kt9uTFPGwPbhQZwgw0HYj7+EiYVQ7ANsz8mH3ZsAWJKP650pVa+jrNWDtvO17o6TrSSvhlW0f2oyHcCgqBUHQacgd76lk+4useHcnyR/l0pyr+iRgrO2RwMFZ6vUL4B3akLc6jGoHYXsmMBoYLWkSKcTmy5Lmy6vVRlsSzfT3uqTJJEO9JLBe3qKYD/iqpNHF/KqFdqGoFLSZVU/9WYeNNToclXoGLYgrrfvQ9m3AbWVlxxWujwaObqaPyylEXFQjzlQ7AEnflLRKoWgA8CowCihlndmTJIjfVD/LSlowX38F2AR4zvb5tpfOick3JYXjNNT3KYIgCILmiJVqx9AXODsL5n8BvEDagl0UuEbS74HxwP8BSFofGEHKOPNDSSfaXhNYHfhzzvEq4HTbkzr8aYIgCFpBb9D+DaPaAdh+Eti4wq3ppMDk8vpjSNvB5eV3k7LXNDXWK2SBiCBoDREXGgStJ4xqHZE0E5hEWkXOBA6y/Uh2HvonSTqwD0kn+FzbV+R2g4EjczczgANsT6jQ/2rAZcC3gWNtn164dykpFvatEH7oXRw78Ud17e+UtYfXtb8gmE2sVIMW8nEWZkDS/wP+CGye771oe918b0VguKR5bF8GvAxsnpOWb0NyJKqk+PEOKf3coAr3LgfOAYbV8XmCIAjqRm/Y/g1HpfZjUZJKUiOyBuVhJAOJ7Udsl+o+RoWt31zvrbw1/HmFew+QjG4QBEHQScRKtb4sKOkp0hbvUiTFpGqMA1arUL4PcHs7zC3oQBoaGjpsrJc/nFzX/h5euKHqvRBICNpEL1iphlGtL8Xt342AYZKqnW820qPMeVD3oXrmmjYT4g9BEATtRxjVdsL2o5KWIAkzVGJdshQWgKS1gUuAbWy/ncsOBH6Rq2xbkCRsy7xC/KED6MgVXTgqBd0C944z1TCq7UT21J0XeJtCDtZ8rz9JnPns/H55Ukq33W1PKdWzfS4puW4QBEHQDQijWl9KZ6qQtnf3tD1TEsBKksYzJ6Tm7Oz5C3AcsDhwXq77he1GyXAlfR0YS3KCmiXpEGAN2x9I+jvQQMotOBU43vb/tdeDBkEQtJhYqQYtwfa8VcpfARZsot3PgZ/X0P9/qO4ZvFttswx6GrFdG3QH2jFJeZciQmqCIAiCoE7UdaUq6WvAmcB3SDGanwGn2R5Rz3FaMa9DgItsf5Tf3wb81PZ7nTyvIaT8fAc1UacB+Mz2Ix01r6Dn8sirKzZbZ+MVXuqAmQS9kk5I/dbR1G2lqnQYeCPwgO0Vba9HSgZbcbuygzmEgrOQ7W0726C2gAYq6wYHQRAEXYx6bv9uSVpRXVAqsP2q7bMlzSvpT5LGSJooaT9IqzBJ90v6h6QpkoZKGizpCUmTJK2U610u6XxJoyS9JGlzSZdK+qeky0vj5TpjJU2WdGIuOxhYGhglaVQue0XSEpL6S3pW0iWSnpZ0laStJT0s6XlJG+T6i0m6Mc/9sRz+gqQT8jxG53kdXJjLjZKezHPZt1C+V37W+0mp20rlS0q6IX9GYyRtkr2E9wcOlfSUpM0q1cvtN891npI0XtIidfzdBkEQtBm5/q+uRj23f9ckqQRVYh/gfdvrS1oAeFjSXfneOqSUZu8ALwGX2N5A0q+BX5FWmZDSoG0JbA/cTDJIPwfGSBpg+ymSyPw7kuYF7pW0tu2/SjoM2ML29ApzWxnYhSSIMAb4KUl8YXvgGJLO7onAeNuDJG1J0tcdkNuvBmwBLAI8J+l8258De+e5LJjneAPwpdzXesD7pHyq43M/fwHOtP1QDrG50/bqki4AZpTE8yVdXV4vf36HAwfaflhSX+CTKr+LoJvQHqpMH3zyRrN1Fu3TsnFDZSmoiU5MUt6RtJv3r6RzScbpM1JC7rUllRJy9wNWyffG2H4jt3kRKBnbSSRjVeJm25Y0CXizlEdU0mSgP/AU8OO8KpyPJBO4BjCxmam+XNbXvYVx+uc6mwI7Adi+T9Likvrle7fa/hT4VNJbwNeAqcDBknbMdZbLz/t1YLTtaXm8a4FVc52tgTVySA3AolVWm9XqPQycIekqYLjtqZUeVqGoFARB0G7U06hOJhseANsHKikKjQX+BfzK9p3FBtkJ59NC0azC+1ll8/u0Qp3Z9SR9g7RaWz9ne7mcFBPaHLWM30hSkDnfuYrtZ+a5NJCM30a2P5I0ujCXat/V5sn1Py4WFoxnk/WAoZJuBbYFHpO0te1nG006FJW6De2xAqzNUan+4wYBgGZ19gzan3qeqd4H9JF0QKGs5Bx0J3CApPkBJK0qaeE6jg1JEOFD4H0lL+RtCvf+S9qebS0PAINh9heB6bY/aKJ+P+DdbFBXI3lDAzwONOSV7vykbecSdwGzvYAllbaXy+desZ6klWxPsn0q6YtMJbH+IAiCoB2p20o1b5kOAs6UdAQwjWTkjgSuI22ljlNaek2jck7Qtow/QUmxaDLpbPbhwu2LgNslvWF7i4odNM0JwGWSJgIfAXs2U/8OYP9c/zlSOjdsvyHpBOBR4A3SGXRJMOJg4NzcZj6SId+fdH58vaQdSGfM1eodoiTIPxN4hsh0E1QgwmWCTqUX7I3JvSBuKKjMwIEDPXbs2M6eRhAE3QBJT1aST62Vvost53W2OqT5ii3kkesPb9O86k2PkSmUNMN233bs/xLgDNvP1Fh/ILCH7YObrdx8X6+QRCKml5XvD3xke1hbxwiCEuc821Qa4I7joNXu6+wpBEGL6TFGtb3J+rwtqT+WdLZZE5LmtT2zhWNc0HytIAiCLoAJRaXujqQVJN2bRRvuzXGdJTGJv0p6JIs27JzL55F0XhZsuEXSbYV7o/PqE0kzJJ2axR3ukbRBQQBi+1ynQdIt+bqvpMuUBC0mStqp0M9Jkh4HNpK0VRZumKQkKrFA4XF+qySK8YSklXP7EyQdnq9/oSQGMUFJHGKudHNBEARB+9PTV6rnAMNsXyFpb+CvzHGQWooUf7oaMBK4HvgRyaHqW8BXSUnEL63Q78KkeNMjJY0Afg98jxQXe0Xur8j/ksQvvgUg6SuFfp62fZykPsDzwFa2p0gaBhwAnJXrfpBFMfbIZduVjTHc9sW5/9+TBDfOru1jCoKO4y97jG++EnD9Qg011Qvxie5DV1RAqjc9eqUKbARcna+vJBnREjfanpXPSL+WyzYFrsvl/yEpHlXiM5KHLySRivuzilJRMKLI1hSSjdt+N1/OBG7I198kCVGUkpRfAXy30MffCz83qjDGWpIezKIVg0kKV42QtK+SlOPYadOmVXm8IAiCoDX09JVqOcXvSUXRBpX9bI7PPcdterZghO1Zkip9pqKyM/knhXPU5sZ2lesSlwODcmjREJIQf+NOQvwh6GR+PWzdmuqFo1IPpBf8xenpK9VHSJlyIK3eHmqm/kPATvls9WtUMUytoFyw4SsV6jwL9C+dlwK7A/cX7u9a+PlohfaLAG9kUYnBbZ5xEARBHSklKQ9B/e7DQpKKerdnkIQSLpX0W5LgxF7N9HEDsBXwNDCFpID0fh3m9nuSYMPTpC3fE4HhxQq2P5G0F3BdXu2OAYrevQtkh6Z5gN0qjPG/eb6vkrahI0tNEARBBxPiD2VI6mt7hqTFgSeATfL5ao8jxB+CIKiVtoo/LPLlZT2g4df1nBIAD910RIg/dHFukfRlUpq2k3uqQQ2CIAjqT5cyqpKOJeUznUlyANrP9uN16Hc0cHgWZGgS2w1V+njE9sZtnUtbydrBM2yfLukk4AHb93TytIIexDbL/Kqzp9Dp3P7viEZrD7riGWi96TJGVdJGpNjLb9v+NKeN+1InT2u20lFXMKjl2D6us+cQBEFQM73AqHYl79+lSCnVSuEp022/Lum4rBT0tKSLcpYbJK2c1YwmSBonaaVcfkRWJJogaWih/12yGtEUSZvluv1zbOe4/No4lzdIGiXpapLTD5JmFO7dL+kfua+hkgbnvicV5tFSNae+ud643M8OpYlLOlbSc5LuIcWzUuir1L7i5xQEQRB0HF1mpUoKOzlO0hTgHuBa2/cD59g+CUDSlaTV7M3AVcBQ2yOyGtE8krYhKSZtmHOZLlbof76sSLQtcDxJkOEt4HvZ83YVkrBC6cB7A2At2y9XmOs6wOrAO6Q0c5fkvn9NSs92CC1Xc/oE2NH2B3mV/pikkcC3SWFB65J+X+OAJyvMqdrnFAS9nsenD2++UoGGhkktqh+qTrXRG7Z/u8xK1fYMYD1gX1L4y7VZxGALSY9npaAtgTUlLQIsY3tEbvuJ7Y9IhvKyfI3tdwpDlP5XPckc1aP5gYtz39eRZAZLPFHFoAKMsf1GXlW/SPpCAHMrKrVUzUnAH5TypN4DLJPvbQaMsP1RToxeLoFYotHnVKlSKCoFQRC0H11ppUpWFxoNjM7GYT9gbVLas9eyk04fqqsPVVMugjkKSjOZ89yHAm+SVp7zkFaLJT5sYqpFNaZZhfezqP6ZNqfmNBhYEljP9udK6d76VGjbiLxSP4/Gn1PjSYSiUtAL2XCJH7Wo/u2jw1Gp7hiY1fP/5HSZlaqkb+Yt2BIDgOfy9XRJfYGdAfKKbaqkQbntAkpZWe4C9s7XlG3/VqIf8IbtWSQFo3nr9kAtV3PqB7yVDeoWwAq5/AFgR0kL5hX6Dyu0LRnQuT6nIAiCLoXb4dXF6Eor1b7A2TlG9AvgBdJW8HukbdVXSCpDJXYHLsxhJZ8Du9i+Q9IAYKykz4DbgGOaGPM84AZJu5DE85tanbaUlqo5XQXcLGks8BRJthDb4yRdm8teBR4sb2j7PUkXU/lzCoIWEeEkQdB6QlGpFxOKSkEQ1EqbFZX6Lev1Nj64nlMC4P47jgxFpXoiaSZphTYfKf/pnqRcqLfYXquNfQ8E9rDd4n8JkmbY7tuW8YOgFvoPG9p8pRbwyh5H1bW/IOhNdJkz1Tbwse0B2YB+Buxfr45tj22NQW0vVDmtXBAEQffArv+ri9ETjGqRB4FS6rR5JV0sabKku7Kjz0qSxpUqS1pF0pP5ev0syDAhCzkskoUebsn3N5f0VH6NL9x/QNIISc9IukDSPIX+T8n9PaaUSg5JP8yhL+OVxCtK5QtLujQLOIwviT9IGiLpOkk3A3c1IxLxszz3pyRdKKmejldBEARBM/SYlU9exW0D3JGLVgF2s/0LSf8AdrL9N0nvSxpg+ymS89Dlkr4EXAvsanuMpEWBj8uGOBw40PbD2cO2FH6zASm+9dU89o9IYg4LA4/ZPlbSacAvSCngHgK+Y9uSfg4cAfwGOBa4z/be2VnrCSUFJUgxr2vbfic/ZyWRiNVIuVY3yR7E55G8jofV4eMNOpGGhoYm7//nzX81ef/rR+9bx9kEQevpDeIPPcGoLijpqXz9IPB/wNLAy9lwwtyCD5cAe0k6jGSENiBJ/71hewzMDtlBcyv9PQycIekqYLjtqfn+E7ZfyvX/ThJ5uJ60FX1LYfzv5etlScIWS5G0jUsCE98Htpd0eH7fB1g+X99dELIoiUR8lxQXWxKJ2IoknjEmz2tBkmLUXEjal+RVzfLLL19+OwiCoH3ooiEw9aYnGNWPbQ8oFmSjUhRYmEkyMpASkR8P3Ac8afttScvQzK/b9lBJtwLbklaHW5dulVfNPz/3HNfqouDE2cAZtkdKagBOKE2btJouxeaWnmVD5g71qSYSIeAK20c38xwh/tDNaE4Cr96OSkEQtJ6edqbaLLY/Ae4Ezgcuy8XPAktLWh8gn5fO9YVD0kq2J9k+FRhL2m4F2EDSN/JZ6q7UJvLw73y9Z6H8TuBX0uyEAes20b6SSMS9wM6SvprbLyZphSp9BEEQdCgCZNf91dXodUY1cxVpRXkXgO3PSAbxbEkTgLtpLPN3iFIGmAmk89bbc/mjwFDgadJW7ohmxj4BuE7Sg8D0QvnJJC3iiZKezu+rzX1gFokYzByRiGeA35GcmSbmZ1iqmbkEQRAEdaRXij/kc8t+tv+3jf00kJKfb1eXiXUwIf4QBEGttFX8YdFFl/XA9Q+q55QAGHXf0V1K/KHXrVQljQD2AP7S2XMJgiDoTXTW9q+kHyjlpH5BUiN1kxy6OK0QNvnzXD5A0qM5NHOipF2bG6vTHZVUQREp50KtqyKRpP7AxrZ3rFeftkeTsuq0isKcrm6marX2JwEP2L6n2cpB3djhofp/2+5K3LTpOZ09hSCoGzle/1xSBMZUUoTEyHxkVuRa2+X/uT8iqeo9L2lp4ElJd9p+r9p4XWGl2m6KSGX0B37aTn23lv60YU62jwuDGgRBt6A9MtTUtlDdAHjB9kvZf+YaYIdm2qQp21NsP5+vXyeFKS7ZVJuuYFSLFBWRAKimICSpv6R/qkw1Kd9bOasVTcjtViI5E22Wl/aHSuoj6bLc5/jsSVvaBrhR0s2SXpZ0kKTDcp3HsldtS5WZ5pX0JyW1pImS9stNy+fUX9KDec7jJG1cGOOIPNcJkobmsssl7Zyvj8v9Py3popIXcRAEQS9nGeC1wvupuaycnfLf5+slLVd+U9IGJG2BF5sarNO3f0uosSJSiU+orCAEFVSTgL+RPGSH2h6hlMB7HuAoCk5Fkn4DYPtbklYjec2umvtdC1iX5AH8AnCk7XUlnUnaCjhLLVNm2gd43/b6khYAHpZ0V4U5LQR8z/YnSrll/07y9N0GGARsmLfGK+WJPcf2SbmfK4HtgJtb+GsIClRTMnr6/eerttnk7LqdLgRBD6PdtHqXyNEQJS7K8fglKi0wyidyM/B3259K2h+4AthydgdJrOdK0vHkrKYm0xWMaiVFpCLVFISggmqSUiLvZWyPgNlxqeXqSJCUj87OdZ6V9CpQMqqjbP8X+K+k95ljnCYBa+frligzfR9Yu7SqJMWarkLa7i4yP3COUk7YmYX5bA1cZvuj3O87NGYLSUcACwGLAZOpYFQVikpBEHQS7SRTOL0Z79+pQHHluSzwerGC7bcLby8GTi29yYujW4Hf2X6sucl0BaPaSBGpjGoKQlBZNanWbc+m6hX7nVV4P4s5n1lLlJkE/Mr2nXMVppCcIocCbwLrkFbXJX1hVem31E8fUsL1gbZfk3QCjeNsgVBUagnVlIx6uqNSEPQwxgCrSPoGSXjnJ5T5skhayvYb+e32JKdZ8u7jCGCY7etqGayrnalWopqCUEXy6nCqpEEAkhbI26r/BRYpVH2AZLDJ277LA89RIy1UZroTOEDS/KXxJC1cYU79SCvdWcDuQCnLzF3A3vk5qLD9WzKg05XE/ncmCIKgq9EJqd9sfwEcRPo7/E/gH7YnSzpJ0va52sHZN2cCcDAwJJf/GPguMERzwm2aWgR2iZVqc1wF3Jz3zJ8iKwg1w+7AhUohJ58DuwATgS/yh3Y5aWV3gaRJwBfAkLyf3tK5/YiCMpNSHNPZSk5TH5O2bi8hefqOyw5E00hnpJXmdIOkXYBRZM1f23fkX+RYSZ8BtwHHlCZh+z1JF5O2p18hfTML2okIOQmC7oXt20h/N4tlxxWujwYa6abb/hvJT6dmeqWiUr1QnZSZOotQVAqCoFbUVkWlvst4wwG/rOeUALjn4d91KUWlmlaqSom0zwS+A7xLcrA5reQM1F3JZ5qf2X6kFW1HABuTgorbOo9BwJQKwcjNtdseWMN2pCnppax74Jl173P8uYfWvc8g6C00e6aatytvJCn3rGh7PdJB77LtPbl6oLJsM2U0kAxji8nKTOeTFDfayiBSovOakTSf7ZFhUIMg6DZ0wplqR1OLo9KWpNXcBaUC26/aPruaWIGkBkn3S/qHpCmShkoarCSGMElJjKEkXnC+pFGSXpK0uaRLlUQdLi+Nl+uMzQfJJxbKX5F0au73CUkrF/o9Q9Io4FQlwYYblQJ7H5O0tpJE4P7AofnweTNJS0q6QUlEYYykTXJ/jdoXPp91JN0n6XlJv8j1KwpW5Ht75H4mSLoyf2bbA3/K81gpv+6Q9GT+fFer8lxDJJ2T7/1Q0uNKIhX35N2FIAiCrkPnKCp1KLVs/64JjKty7y0qiBXke+sAqwPvAC8Bl9jeQNKvgV8Bh+R6XyEZ7u1JcZWbAD8n6TOWxBWOtf2OkobjvZLWtj0xt/8g97sHcBZJ9ABSjOfWtmdKOhsYb3uQpC1J7tEDJF0AzLB9OoCkq4EzbT8kaXmSt9jqwInl7YGSB9japG3xhYHxSonM36KyYMUawLHAJranS1osP9dI4Bbb1+d53Avsn/UmNyQ5MJUCkYvPNaTwu3gI+I5tK4lBHwH8psrvLegGVBOfKPL881Nr6muVQQe2cTZBENRCi71/JZ1LEk74jOTZWkmsAGBMKe5H0otkD1mSh+oWhXo3Z0MwCXjT9qTcZjLJY/Yp4MdKogXzkXKErkHynIVkyEs/iwdM19mema83JaktYfs+SYtL6lfh8bYG1tAcD+BFlcQkmmp/k+2PgY/zCnIDUqBwJcGKLYHrbU/PfTUScVAKidmYlHO1VLxAlecqsixwrZLyx5dIuV0boRB/CIKgk+iKScXrTS1GdTLZoADYPjCvvsZSXawAahNQoKy8vM18SgG7hwPr2343bwsXhQ1c5frDwnUtMlXkZ9goG8k5jVUxzsZlP4vl1QQrmhRxKMzhvSYEMT6sUn42cIbtkUoOWCdUqhTiD92HauITRdrDUSkIgtZTy5nqfUAfSQcUyhbKP6uJFdSTRUmG5P18TrhN2f1dCz8frdJHUeihgSRr9QGNxRfuIgUJk+sOaKY9wA5K4vyLkxyfxlBdsOJe0qp78dxXScRh9jxyvy8rxaqixDrVPpwC/UhqIQB71lA/CIKgY+kFjkrNrlTz1uwg4EwlbdlpJCN3JOmstZFYQT2xPUHSeNKK+SXg4bIqC0h6nPQFYbcq3ZwAXCZpIslbt2R0bgauz45EvyIpaZyb681HMqb7N9Ee4AnSdu/ywMm2X5dUUbAiq3icAtyvlEd2PEm54xrgYkkHk9SQBgPnS/odSQ/4GmBCMx/VCaQt438DjwHfaKZ+EARBx2HS/mMPp1uLP+Rt1YGlM8qgZYT4QxAEtaK2ij8svIy/s+Z+zVdsIXePOb77iT8EQRAEQVsQDkel1pK3NieRti6/IOWmO6upPHQ5bvQW22vls8yls15jVWz3r8Ncj7H9hza0n2G7bxvn8GXgp7bPa0Xb23Lb99oyhyAosdIZZ3T2FNqVFw87rLOnEPRg2itLzce2B9heE/gesC0pTVqtDMht2kx29GnqOY9p4l5H8WWgRaKYpeeyvW0Y1CAIugXhqNR2bL+VYyPHKOX5nAcYSvKUXQA41/aFpfpK+etOIiUv3xT4Iynm8ixSvtSPgb1sP5fFD3bM/XwDuNr2iXnVezvJeWojYFBWLjqGFNZyq+0jJQ1lTpL0ybYHS/oZyWHpS8DjwC+z0MJu5e0Lc/4zKfb2XeAntqcpqSvtm/t5Adjd9kfZg/kCYMXc/IA83kp5Hnfb/q2k35LSDi0AjLCjAjfgAAAY+0lEQVR9fJXnup98rizpRlIy3j7AX3L4TBD0eF4/r/ZNnoaRI2uuW0tYU9ACuqARrDcdkk/V9kt5rK8C+wDv214fWB/4RY5FLdX9DDgOuDavdq8lec9+1/a6+V5xu3YDkrfsAGAXSaUD62+SlJPWJaV/O5UkvjAAWF/SINtHMWdVPVjS6qTQnE1ynOhMYLCkpSu1z+MsDIyz/W3gfuasyIfbXt/2OqQcfvvk8r8C9+fyb5O8mo8CXszz+K2k7wOr5GcbAKyXhSTmei7br5Z91HtnbeaBpPyAi5f/LiTtqyT5OHbatGnlt4MgCII20JGOSiUBhe8Da0sqJdLuRzIgU5po2w+4QkkK0aSz2hJ3234bQNJwkvrRjcCrth/LddYHRtuelutdRUo8e2PZOFsB65FW1ZBWxm81034WcG1u/zdgeL5eS9LvSVu7fUmSh5AM8x4AWRnpfUlfKZvH9/NrfH7fN39G/yp7rnIOlrRjvl4ut3m7WCHEH4KeyNK/rP30ZHScqXYOvSSkpkOMqqQVSau+t0jG9Ve27yyr07+JLk4GRtneMdcbXbhXSdEImldUqjhV4IqcsLY4t0FV6leiNP7lwKAcZzuEtN1dKwL+WNwWz/PoT5VY4CxKsTVJEeojSaOZW3kqCIIgaGfafftX0pKkM8RznIJi7wQOkDR/vr+qpIXLmpUrHRXVgoaU1f2eUhaZBUkp1MrFISCdjW4uaQklUf7dSFu1AJ+X5kJSPNpZ0lfz3BaTtEIz7echCTYA/JQkbE+e/xu578GFudxLOkdF0rySFq3wvHcCeyvpACNpmdKcmqAf8G42qKuRRP6DIAi6DLLr/upqtNdKteT8UwqpuRIo+elfQhLKH5c1daeRjGGRUcBRuY8/AqeRtn8PI8kmFnko978yyVFpbPmq1/Ybko7O/Qq4zfZN+fZFwERJ4/K56u+Au7LH8OfAgbYfa6L9h8Cakp4E3meObOL/kozxq6TwopLR/DVwkaR9SKv3A2w/KulhSU8Dt+dz1dWBR/M29AzgZ7l+Ne4A9s+qT8+RVJWCoMVEyEkQtJ7urqg0hOT5elBzdYPGhKJSEAS10lZFpX4LLe2NVtmn+Yot5M6Jv+/ZikrZUWY4sLrtZ4uiDk20abZOvcljbmz76la2bwAOt71dc3Wb6acmoYsK7ZYG/mp752YrB0ELWO34jst88+yJh3bYWEFn0zXjSutNe5yp7kbakv1JO/Q9F7Yvb2qVKqmpLw39SWegnU2LhS4kzWf79TCoQRAEXYu6GtXsWLMJKSazkVGVNETSTZLukPSc/n975x5uVVnn8c+Xi6glhqGGF4LMG2p5QRJLPYoXpnFM01ImU8rRzKkmGy17cjTNmbw82aRoJuSgNpGJomRGoImQclUEBBOVGCVJQcyGKU3lN3/83s1ZbtY+Z3POOnvvw/59nmc9e6213/W+37X3Pudd7+37k7IuSz0ljZW0RNLUNPEISWdLmidpoaS7JG2dzo+XdJOkmZKWSTo+U8adkn6Bj41K0jWSnpS0WFJpzPNK4DBJT0g6P00auiaVtUjSF1J+la4HD2I+SdLSpKVHuuaHaS3oEkmXZe7/YEmPpnuZKw90fjlwatJxqqR3Sbol6Vggj6CTd1+D0hgsaX+mpMfTdminv8wgCIIiMcJRqQOcCEwxs2WS1ko6EFhblmYYsC8eQm2epF8Ca/A1laPM7GxJP8cDo/8EN1EYCyBf93kWHpAbvLV5BLAb8JCkD6bzw4EPmdlaSSfjrcEPA/1TmTNww4UN3bdy16fXzOxgSX2ARyRNxQ0a8q4v3csQfDLSFOCTwETgW6nsnsCDkj6EG1jcAZxqZvPSrN+/4GYWG8aFJf0H8Bsz+7zcE3iupAdy7mtQ5jN9GTjGzF6Xr+WdgBtABEFD8fytN7zjuOXhe99xHA5GQXen6O7fUXjsT9JrXnzTaWb2ipn9FR97/Vg6/3szeyLtP4ZXmOAmCjMlLcaXpuyTyevnZrbezJ7BY63ulSmjVJl/DJhgZm+b2Uv4UpiDc3QdC5yRZhzPAd6LV/RtXT/XzJYnE4cJmXv5tKTHcfOGffCKd088oPs88GDkZvZWBR2lmc/T8bWmA3PuK0tvPB7rYuDOVF4uCkelIAjqxfou2BqMwlqqcku8o/BK0ICeeIO/3JSzklnDG5lzb+NuRtC2iULRxg95phRtjXduVL7ccvEC4GAze1XSeLxiVE76SjpONrOny3R8hMpB4M8HXsJb0z2A1ysKDkeloI4MPPOf33E8PSYqNRWNuK60aIpsqZ6Ce9K+38wGmdmuuBH+LmXpqjFryFLJRAHc67eHpN1wg/qnN7oaZuBjlj3lRhSHA3PJN1zIM6WodD3AMEmD01jqqfgErb545fea3Dz/71La3wE7STo45b+NfCJVno4vS75AVdIB7Xw+4MYPq1Jovc/iDzRBEARBjSlyTHUUPvkny11sHFqtXbOGMiqZKIBXog8DOwLnpjHF8usn4WORC/GW4tfN7I+SXgHekrQQbw3/gHxTikrX7wXMSve8H175TjKz9ZIW4Eb5y0kPDWb2tzTJ6fr0QPFX3Faw3OjiO3hEnkVJxwqgvWU7NwJ3SfpUyq9SizYIgqB+NEFLtabmDyrQrCF1q95nZhM7m1ezEuYPQRBUizpr/rDVADt00OgCFTlTfnfl5m3+EARBEAQbYcD6zb+lWpN4qiXaM2vYxLxGd6SVml3fmTn3bUkXdCCv6WqN39opJA2VdF3a7yPpgcza1XGSKs7ozcmrRdJ9RegKgiAohi5Yo9qA3cnRUm0QzGw+UOqLPQDonQKlQ2u81iAIgqCBiUo1gzwG6RzgSDy4+FlmNjOZOFwFHId3Yow1s+vLrv0hvn51K2CimV2azl8JnIBH65lqZhekCUWX4kuHXjOzw5W8hIHP46YX26fJSycDP8aNKuZLOha4DOgDPAd8zszWSRqJT3BaAzzeJR9QsNnS0tKSez7MGIJCacCWZdFEpboxvcxsWFqfeik+Q/ccYDBwgJm9JWm7nOvyXJRWAicBe5mZJYckcBel48zsD5lzAJjZy5L+iXe6PZFe+wMXA0eb2f9J+gbwNUlXA2PxdcLP0kbLNjlHnQMwcODASsmCIAiCDtCMlWqlR6XS+bvTa9bV6WjgppIDUgVXo0+nCqsXMAB3NVqKGzGMS3aMpXHOR4DxyY7x7py8KnFIyveRVNFugS/r2Qt3pHoGQNJPSBXnRjcZ5g9BDtEiDWpCtFQ3S14B+pWd2w43qoBWZ6e3af182nRDquSilFq1w4AReICBLwFHmdm5ySHp74En5OHfqkG4VeE77B/T9Zv/rzUIgqDBqens30bAzNbhDk0jAFJX7kjclKISU4FzkwMSOd2/uS5K8qg926ZYqV/FjfmRtJuZzTGzS/Ax0F2rlD8b+GgpcICkrSXtgbs1DU7OUpDvuRwEQVA/Sktqit4ajGZsqQKcAdwg6Xvp+DIzey7HjanEOGAP3OXoTXz8ckzpzeRLvJGLEu7+dK+kkvdvyej0mhRNRsCDuFvTEe2JNrPVyUBjgjySDsDFKSrQOcAvJa3BHxBqFvA9CIKgfQysAR3wC6amjkpBYxGOSkEQVEunHZX67GiH7lRu3955pqz4fjgqlUhLWL6bjQwj6avAHmZ2Xt2EVUkyfjjDzL5SQF4rcAvHNWXnzwX+Yma3dbaMIKiGQ06/tt4Sasrsn3yt3hKahyZoxNW7+3cCPoEnG27tNODC9i5MZvNKkVlK53qm2KY1ocywoV06os/MbtpkYUEQBEFdqPdEpYnA8aXxwRStZifgt5IulDRP0iJJl5Xel/SUpBtxg4NdJa2TdLmkOcBwSSvSes6S9d/0tP9uSf8laXHK8+R0flQ696Skq0rCUr5XSXosWQYOS7aEyyWdkNJssANsI/9yfSMkLUjpbsmMjQJcKGlu2kqTkTZYKEo6O30mCyXdJWnrdH77dDwvbR/tgu8qCIKg48REpa7HzF6RNBeffXsv3kq9AzgG2B0Yhk/mmSzpcOB5YE/cReg8AHnM0yfTTNoNRgk5/BvuXrRfStdP0k64U9JBwKvAVEknmtk9wLuA6Wb2DUmTgCuSriHArcDk9vJP5zfoSxOWngFGpMlFtwFfxJ2QAP6cjCfOSOfKQ77dbWZjU/5XAGcB1+Nh675vZr+VNBBv+e9d6YMIgmZj6bQbK77X0lL+p9xKrN8tmCbo/q13SxVau4BJrxOAY9O2AG+R7oVXsgD/Y2azM9e/jcdtbY+jgRtKB2b2Km4rON3MVidjh//Gg5AD/A2YkvYXAw+b2Ztpf1CV+Zfr2xM3aViWjm/NlAd+76XX4Tll7CtppqTFeMD2fTJlj5HbGk4G+kraJud6JJ0jab6k+atXr85LEgRBEHSQeo+pAtwDXCvpQGArM3tc0mfwCUw/yiZM3cPlAbhfLxunfIvWh4Uts5ezsUFCxWYt8Ka1To1eTzKFSEHI8z63SgYRWX1tlUfZ9Xl5jQdOTEt4RgMt6XwPYLiZ/bWd/MNRKWhKhhxTed7j9JioVDuipdr1JDOG6cAttLbUfg18PpknIGlnSTtUmeUKvDsX3Iy+xFTc0YiUZz/cPP8ISf3lnr2jgIc7die5+ZfzO2BQabwU+GxZeadmXmflXL8NblzRG2+pViq7WoemIAiCoEDqXqkmJgAfBn4GYGZTgZ8Cs1JX50S8QqmGy4AfSJqJd72WuALolyYkLQSONLNVwDeBh3ADhsfN7N4O3sNG+ZcnMLPXgc8Bd6b7Wg9kZ/f2SROa/oVWo4gs/4Y/CEzDK+gSXwGGpglSS4FzO3gPQRAEXURzxFMN84cmJswfgiColk6bP/TewQ7t/6kiJQEw5Y83hvlDEASNy8j+uQGOqmLKmpsLVBIE3Y9G6f5tF0nralzehjWoXZD3eEmndEXeQRAEDUsTdP92m0o1CIIgCBqdbt39K+n9+Kzh7YHV+CSgP+AGC7sB2wJrgRYzm5EmL33OzJ7N5DEIuB03aQD4kpk9mvb7JuOHPYEZwHlpSc06MyvNTD4FON7MRsvjqP4ZGAq8D/i6mU1MlorXA0fhcVs3LK2RdAnwD8BWwKPAF8zMkhPUHHzC03uAs8xsZpqlfBVwHL7sZqyZXS/pIOBa4N14OLnRaSJW0OS0tLRsUvpFry1rP1FBZXWGMGbohtSpZSlpJG6S0xMYZ2ZXlr0/GrgGrz8AxpjZuPTemcDF6fwVZnZrW2V195bqGOA2M/sQbtxwXVoTugx3PvoY8BhwWLID3CVboSZeBo4xswPxpSzXZd4bBvwrsB9eSX+yCk0DUrnHA6Uv7iS8Yt4POBs4NHsPZnawme2LV6xZF6VeZjYMj8V6aTp3DjAYOKB032mJzfXAKWZ2EP6g8e954sL8IQiCZiI1RG7A41wPAUZJGpKT9A4z2z9tpQp1O/x/70fw+uDSCsslN9CtW6q461CporsduDrtz8SdigYD38UrsoeBeTl59MbdiPbHl+DskXlvrpktB5A0Aa8sJ7aj6Z5k8r9UHrCcpGVCqvBflPSbTPojJX0d2BrYDo/J+ov03t3p9TFaXZyOBm5KDlCY2VpJ++LxU6clm8aeQG4rNcwfmo9NbdF1aqLS9JioFFSibl69w4BnM//LfwZ8AlhaxbXHAdPMbG26dhpuqzuh0gXdvVItp/SNzcTXau4EXIJHvWnBu3DLOR94CV8n2wN4PSe/8uPs+S3L0ryR2c86KG30a0pewDfiId9ekPTtsvxKeb1N63dVyRlqiZnlWRsGQRDUHwPrmiDl/SVl1wbenBoPJXYGXsgcr8RbnuWcnDzmlwHnm9kLFa7duS0x3b3791FafYM/A/w27c/Bu1jXJ8OFJ4Av4JVtOdsCq1Lr8rN4K6/EMEmDJfXAu4ZL+b8kae90/qQqdM4ATpPUU9IAWo0hShXomuQeVc2M4KnAuSWrxNQ98TSwvaTh6VxvSfu0kUcQBMHmwhozG5rZyrtL8uxhyxsmvwAGpSG1B3Bf9mqvfQfdqaW6taSVmeNrcSehWyRdSOtEJczsDUkvACXj/Zm4BeHinHxvBO6S9CncWSnrLTwLHxfdD68YJ6XzFwH34U8wT+KTg9piEj5JaTH+FPRw0vknSWPT+RXkd0+XMw7vol4k6U18otKYNGHqOknb4t/rf+JdyUGwScRa06DLqE/370pg18zxLsCL2QRm9krmcCw+GbR0bUvZtdPbKiwclZqYcFQKgqBaOu2o1Gt7G973xCIlAfDrV8e1qSv16i0DRuCze+cB/2hmSzJpBpRWS0g6CfiGmR2SegIfAw5MSR8HDiqNsebRnVqqQRAEQXemDo04M3tL0pfwQC09gVvMbImky4H5ZjYZ+IqkE/AoZ2uB0enatZK+Q2sv4uVtVagQlWoQBEFQC8xgfZdMVKqiaLsfuL/s3CWZ/W/iwVXyrr0FX6ZYFd19olIQBEEQNAzRUg2CIAhqQxPM4YmWahAEQRAURLRUgyAIgppgdRpTrSVRqQZBEAQ1oDFDtRVNdP8GQRAEQUFESzUIgiDoeox6OSrVlGipBkEQBEFBREs1CIIgqA1dE6WmoYiWahAEQRAURLRUgyAIgi7HAGuCMdWoVIMgCIKuxyy6f4MgCIIgqJ5oqQZBEAQ1oRm6f6OlGgRBEAQFES3VIAiCoDY0wZiqrAm8GIN8JP0v8HSdZfQH1tRZAzSGjkbQAI2hoxE0QGPoaAQNAHua2TYdvVjSFPxeimaNmY3sgnw7RFSqTYyk+WY2tNk1NIqORtDQKDoaQUOj6GgEDY2ko9GJMdUgCIIgKIioVIMgCIKgIKJSbW5urrcAGkMDNIaORtAAjaGjETRAY+hoBA3QODoamhhTDYIgCIKCiJZqEARBEBREVKpNhKTtJE2T9Ex67Vch3dWSlkh6StJ1klQHDQMlTU0alkoaVJSGTdGR0vaV9AdJY2qtQdL+kmal72ORpFMLKnukpKclPSvpopz3+0i6I70/p+jPfxN0fC19/4skPSjp/bXWkEl3iiST1CUzYKvRIenT6fNYIumntdaQ/i4fkrQgfScfL1pDt8fMYmuSDbgauCjtXwRclZPmUOARoGfaZgEttdSQ3psOHJP23w1sXevPIpP2B8BPgTF1+D72AHZP+zsBq4D3dLLcnsBzwAeALYCFwJCyNOcBN6X904A7uuD3WI2OI0vfPfDFonVUoyGl2waYAcwGhtbps9gdWAD0S8c71EHDzcAX0/4QYEXRn0V336Kl2lx8Arg17d8KnJiTxoAt8T+qPkBv4KVaapA0BOhlZtMAzGydmf2lQA1V6UhaDgJ2BKYWXH5VGsxsmZk9k/ZfBF4Gtu9kucOAZ81suZn9DfhZ0lJJ20RgRJE9FtXqMLOHMt/9bGCXWmtIfAd/CHq94PI3RcfZwA1m9iqAmb1cBw0G9E372wIvFqyh2xOVanOxo5mtAkivO5QnMLNZwEN4i2gV8Gsze6qWGvDW2Z8k3Z26ma6R1LNADVXpkNQD+B5wYcFlV62hTM8w/GHnuU6WuzPwQuZ4ZTqXm8bM3gJeA97byXI7oiPLWcCvaq1B0gHArmZ2X8Flb5IO/O9iD0mPSJotqWgXoWo0fBs4XdJK4H7gywVr6PaE9+9mhqQHgPflvPWtKq//ILA3rS2CaZION7MZtdKA/y4PAw4AngfuAEYDP65WQ0E6zgPuN7MXOtpIK0BDKZ8BwO3AmWadNlDNu5nyZQDVpOksVZch6XRgKHBELTWkB6vv47+/rqSaz6IX3gXcgv99zpS0r5n9qYYaRgHjzex7koYDtycNm7+pb5VEpbqZYWZHV3pP0kuSBpjZqvRPOq/76CRgtpmtS9f8CjgEH0+qlYaVwAIzW56uuSdp2KRKtQAdw4HDJJ2Hj+tuIWmdmVWczNIFGpDUF/glcLGZza627DZYCeyaOd6FjbvxSmlWSuqFd/WtLaDsTdWBpKPxh5AjzOyNGmvYBtgXmJ4erN4HTJZ0gpnNr6GOUprZZvYm8HtJT+OV7LwaajgLGAneqyVpS9zPt+iu6G5LdP82F5OBM9P+mcC9OWmeB46Q1EtSb7xlUGT3bzUa5gH9JJXGDo8ClhaooSodZvYZMxtoZoOAC4DbNqVCLUKDpC2ASansOwsqdx6wu6TBKf/TkpZK2k4BfmNmRbdU29WRul5/BJzQBWOI7Wows9fMrL+ZDUq/g9lJS5EVars6EvfgE7eQ1B/vDl5eYw3PAyOShr3x+RerC9TQ/an3TKnYarfhY2IPAs+k1+3S+aHAuLTfE/8n9hRekV1baw3p+BhgEbAYGA9sUQ8dmfSjKX72bzXfx+nAm8ATmW3/Asr+OLAMH5/9Vjp3OV5hgP+zvBN4FpgLfKCLfpPt6XgAnyhXuvfJtdZQlnY6XTD7t8rPQsC16e9yMXBaHTQMwVcHLEzfx7Fd8Vl05y0clYIgCIKgIKL7NwiCIAgKIirVIAiCICiIqFSDIAiCoCCiUg2CIAiCgohKNQiCIAgKIirVIAiCICiIqFSDIAiCoCCiUg2CIAiCgvh/QsmcGBpWUq4AAAAASUVORK5CYII=\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\"],\n",
    "              y_stdDev=df1[\"Comparison\", \"TvSerror\"],\n",
    "              GC_content=df1[\"GC content\", \"TGIRT\"], name=name)\n",
    "\n",
    "name = \"ImProm-II 42°C : TGIRT              \" \n",
    "barplotFigure(y=df1[\"Comparison\", \"Tv42\"],\n",
    "              y_stdDev=df1[\"Comparison\", \"Tv42error\"],\n",
    "              GC_content=df1[\"GC content\", \"ImProm-II 42C\"], name=name)\n",
    "\n",
    "name = \"ImProm-II 55°C : TGIRT               \" \n",
    "barplotFigure(y=df1[\"Comparison\", \"Tv55\"],\n",
    "              y_stdDev=df1[\"Comparison\", \"Tv55error\"],\n",
    "              GC_content=df1[\"GC content\", \"ImProm-II 55C\"], name=name)\n",
    "\n",
    "name = \"ImProm-II 42°C : SuperScriptIV  \" \n",
    "barplotFigure(y=df1[\"Comparison\", \"Sv42\"],\n",
    "              y_stdDev=df1[\"Comparison\", \"Sv42error\"],\n",
    "              GC_content=df1[\"GC content\", \"ImProm-II 42C\"], name=name)\n",
    "\n",
    "name = \"ImProm-II 55°C : SuperScriptIV  \" \n",
    "barplotFigure(y=df1[\"Comparison\", \"Sv55\"],\n",
    "              y_stdDev=df1[\"Comparison\", \"Sv55error\"],\n",
    "              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\"],\n",
    "              y_stdDev=df1[\"Comparison\", \"P42v55error\"],\n",
    "              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
}