{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Predict expression modulation targets\n", "\n", "Cameo provides algorithms to search for genes or reactions that can be over- or down-regulated in order to achieve a given biological objective." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from cameo import models\n", "from cameo.visualization.plotting.with_plotly import PlotlyPlotter" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Load the E. coli core model." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "model = models.bigg.e_coli_core\n", "plotter = PlotlyPlotter()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Flux Scanning based on Enforced Objective Flux" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "from cameo.strain_design.deterministic.flux_variability_based import FSEOF" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "fseof = FSEOF(model)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "Model: e_coli_core
\n", "Enforced objective: EX_succ_e
\n", "Primary objective: Maximize\n", "1.0*BIOMASS_Ecoli_core_w_GAM - 1.0*BIOMASS_Ecoli_core_w_GAM_reverse_712e5
\n", "
\n", "Reaction fluxes

\n", "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
12345678910
PFK7.7569418.0364998.2425158.3610488.4795828.5981158.7477939.3931029.8307879.910940
PGI5.4856306.1103996.5131736.6518766.7905796.9292827.1611088.8748829.9651269.981645
PGK-16.431132-16.838737-17.173618-17.421989-17.670359-17.918730-18.198431-18.976723-19.547784-19.761991
PGM-15.238892-15.761645-16.212407-16.577535-16.942662-17.307790-17.704415-18.602292-19.293292-19.628048
PPC3.7583205.0123305.8788756.2845166.6901577.0957997.5329968.4723679.61532311.108854
PPS0.0000000.0000000.0000000.0000000.0000000.0000000.0000000.0000000.4103381.527283
ADK10.0000000.0000000.0000000.0000000.0000000.0000000.0000000.0000000.4103381.527283
ENO15.23889215.76164516.21240716.57753516.94266217.30779017.70441518.60229219.29329219.628048
SUCCt31.4745752.9491504.4237255.8983007.3728758.84745010.32202511.79660013.27117514.745750
THD20.0000000.0000000.0000000.0000000.0000000.0000000.1196742.1438482.9168071.535161
TPI7.7569418.0364998.2425158.3610488.4795828.5981158.7477939.3931029.8307879.910940
FBA7.7569418.0364998.2425158.3610488.4795828.5981158.7477939.3931029.8307879.910940
FRD70.0000000.0000000.0000001.4532353.1753024.8973696.5867037.7551419.12784210.852285
GAPD16.43113216.83873717.17361817.42198917.67035917.91873018.19843118.97672319.54778419.761991
ICL0.0000000.0000000.3860571.2313422.0766272.9219123.7353224.0414594.1433333.893465
MALS0.0000000.0000000.3860571.2313422.0766272.9219123.7353224.0414594.1433333.893465
PDH8.8092708.3360088.1798688.4009578.6220458.8431349.0645399.2909849.1077608.219084
EX_succ_e1.4745752.9491504.4237255.8983007.3728758.84745010.32202511.79660013.27117514.745750
\n", "
\n" ], "text/plain": [ "" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "fseof.run(target=model.reactions.EX_succ_e)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Differential flux variability analysis\n", "\n", "Compares flux ranges of a reference model to a set of models that have been parametrized to lie on a grid of evenly spaced points in the *n*-dimensional production envelope (*n* being the number of reaction bounds to be varied)." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "from cameo.flux_analysis.analysis import phenotypic_phase_plane\n", "from cameo.strain_design.deterministic import DifferentialFVA" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Succinate production" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The production envelope looks like this." ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/html": [ " \n", " " ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/vnd.plotly.v1+json": { "config": { "linkText": "Export to plot.ly", "plotlyServerURL": "https://plot.ly", "showLink": false }, "data": [ { "fill": "toself", "fillcolor": "#B3E2CD", "hoverinfo": "none", "marker": { "line": { "color": "#B3E2CD" }, "opacity": 0.3 }, "mode": "lines", "name": "WT", "opacity": 0.3, "type": "scatter", "x": [ 0, 0.04599586878781212, 0.09199173757562425, 0.13798760636343638, 0.1839834751512485, 0.2299793439390606, 0.27597521272687275, 0.32197108151468484, 0.367966950302497, 0.41396281909030913, 0.4599586878781212, 0.5059545566659334, 0.5519504254537455, 0.5979462942415577, 0.6439421630293697, 0.6899380318171818, 0.735933900604994, 0.7819297693928061, 0.8279256381806183, 0.8739215069684304, 0.8739215069684304, 0.8279256381806183, 0.7819297693928061, 0.735933900604994, 0.6899380318171818, 0.6439421630293697, 0.5979462942415577, 0.5519504254537455, 0.5059545566659334, 0.4599586878781212, 0.41396281909030913, 0.367966950302497, 0.32197108151468484, 0.27597521272687275, 0.2299793439390606, 0.1839834751512485, 0.13798760636343638, 0.09199173757562425, 0.04599586878781212, 0, 0, 0 ], "y": [ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.8811810012836991, 1.762362002567392, 2.643543003851088, 3.524724005134787, 4.396864526652468, 5.265895211723508, 6.134925896794555, 7.003956581865596, 7.8729872669366365, 8.742017952007677, 9.611048637078722, 10.474279162179368, 11.32275599919326, 12.17123283620715, 13.017389427255484, 13.859083737108278, 14.700778046961073, 15.542472356813864, 16.38416666666667, 0, 16.38416666666667 ] } ], "layout": { "height": 400, "template": { "data": { "bar": [ { "error_x": { "color": "#2a3f5f" }, "error_y": { "color": "#2a3f5f" }, "marker": { "line": { "color": "#E5ECF6", "width": 0.5 } }, "type": "bar" } ], "barpolar": [ { "marker": { "line": { "color": "#E5ECF6", "width": 0.5 } }, "type": "barpolar" } ], "carpet": [ { "aaxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "baxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "type": "carpet" } ], "choropleth": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "choropleth" } ], "contour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "contour" } ], "contourcarpet": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "contourcarpet" } ], "heatmap": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmap" } ], "heatmapgl": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmapgl" } ], "histogram": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "histogram" } ], "histogram2d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2d" } ], "histogram2dcontour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2dcontour" } ], "mesh3d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "mesh3d" } ], "parcoords": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "parcoords" } ], "pie": [ { "automargin": true, "type": "pie" } ], "scatter": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatter" } ], "scatter3d": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatter3d" } ], "scattercarpet": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattercarpet" } ], "scattergeo": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergeo" } ], "scattergl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergl" } ], "scattermapbox": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattermapbox" } ], "scatterpolar": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolar" } ], "scatterpolargl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolargl" } ], "scatterternary": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterternary" } ], "surface": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "surface" } ], "table": [ { "cells": { "fill": { "color": "#EBF0F8" }, "line": { "color": "white" } }, "header": { "fill": { "color": "#C8D4E3" }, "line": { "color": "white" } }, "type": "table" } ] }, "layout": { "annotationdefaults": { "arrowcolor": "#2a3f5f", "arrowhead": 0, "arrowwidth": 1 }, "autotypenumbers": "strict", "coloraxis": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "colorscale": { "diverging": [ [ 0, "#8e0152" ], [ 0.1, "#c51b7d" ], [ 0.2, "#de77ae" ], [ 0.3, "#f1b6da" ], [ 0.4, "#fde0ef" ], [ 0.5, "#f7f7f7" ], [ 0.6, "#e6f5d0" ], [ 0.7, "#b8e186" ], [ 0.8, "#7fbc41" ], [ 0.9, "#4d9221" ], [ 1, "#276419" ] ], "sequential": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "sequentialminus": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ] }, "colorway": [ "#636efa", "#EF553B", "#00cc96", "#ab63fa", "#FFA15A", "#19d3f3", "#FF6692", "#B6E880", "#FF97FF", "#FECB52" ], "font": { "color": "#2a3f5f" }, "geo": { "bgcolor": "white", "lakecolor": "white", "landcolor": "#E5ECF6", "showlakes": true, "showland": true, "subunitcolor": "white" }, "hoverlabel": { "align": "left" }, "hovermode": "closest", "mapbox": { "style": "light" }, "paper_bgcolor": "white", "plot_bgcolor": "#E5ECF6", "polar": { "angularaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "radialaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "scene": { "xaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "yaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "zaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" } }, "shapedefaults": { "line": { "color": "#2a3f5f" } }, "ternary": { "aaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "baxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "caxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "title": { "x": 0.05 }, "xaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 }, "yaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 } } }, "title": { "text": "Phenotypic Phase Plane (flux)" }, "width": 700, "xaxis": { "autorange": true, "range": [ 0, 0.8739215069684304 ], "title": { "text": "BIOMASS_Ecoli_core_w_GAM [h^-1]" }, "type": "linear" }, "yaxis": { "autorange": true, "range": [ -0.9102314814814817, 17.29439814814815 ], "title": { "text": "DM_succ_e: succ_e --> [mmol gDW^-1 h^-1]" }, "type": "linear" } } }, "image/png": "iVBORw0KGgoAAAANSUhEUgAABSEAAAGQCAYAAABRUqVJAAAgAElEQVR4nOzd+X9U5f2/8e/fE7XUpa1+2jrIIgiiIi6txQ5bZJFVNhcExF0QcUGtFSkUwX0pFVFUJOxbEkIgQFgSlgCBkJCVrK/vD+mZeYdkQiYnM/eZ5Loej+cPTKI5wK2VVzPn/D8RERERERERERERJbD/5/oCiIiIiIiIiIiIqHvHCElEREREREREREQJjRGSiIiIiIiIiIiIEhojJBERERERERERESU0RkgiIiIiIiIiIiJKaIyQRERERERERERElNAYIYmIiIiIiIiIiCihMUISERERERERERFRQmOEJCIiIiIiIiIiooTGCElEREREREREREQJjRGSiIiIiIiIiIiIEhojJBERERERERERESU0RkgiIiIiIiIiIiJKaIyQRERERERERERElNAYIYmIiIiIiIiIiCihMUISERERERERERFRQmOEJCIiIiIiIiIiooTGCElEREREREREREQJjRGSiIiIiIiIiIiIEhojJBERERERERERESU0RkgiIiIiIiIiIiJKaIyQRERERERERERElNAYIYmIiIiIiIiIiCihMUISERERERERERFRQmOEJCIiIiIiIiIiooTGCElEREREREREREQJjRGSiIiIiIiIiIiIEhojJBERERERERERESU0RkgiIiIiIiIiIiJKaIyQRERERERERERElNAYIYmIiIiIiIiIiCihMUISERERERERERFRQmOEJCIiIiIiIiIiooTGCElEREREREREREQJjRGSiIiIiIiIiIiIEhojZCf6YNVapYXCulBS5vpSenzpMxZq0PDZri8jYX20Zp3SQmGdLS5xfSkd6lJZuf5070QNGj5bFZXVkdebmpr0r8/Wa+AjM9Wr7yjdPnSy6urqdfOAdM15dVlCr+nr7zcrLRTWsjXfJfTrEBEREREREVHsGCFN8xctV1ooHHHTnekaNmauPv32F9U3NEQ+L1VGyO9+2qEVn61Pytd6bekapYXCyj5wtEv+foOGz27xe/GHweM1ZvpCbd6Z0+LzUnWE7OhZS7URctq8d3TLwLE6cfJsi9fXbtimtFBYjz+1RB+sWqvnFq9QXX1DUkZISXr57Y/Vq99oHS04k/CvRUREREREREStY4Q0ecPQW8u+1HsrvtXC9z7R3Y8+qbRQWFPnvh35vFQZISc+vUR9H3oiKV/ry+82adzsxTpeWNQlf79Bw2crdP8UvbfiW7274hs9t3iF/njPRKWFwlr15YbI56X6CHmts5ZKI2T2gaNKC4X1rzaG7/QZCzV01JxWrydrhKy5Uqs+D0zTxGfeTPjXIiIiIiIiIqLWMUKavGGoqrom8lpjY5PSZyxUWiisXVl5khghk9Gg4bP1QPq8Fq9dvHRZvYdN0c0D0lVdUysp9UfIa521VBohJz7zpkL3T1FtbV2rjw0aPlszF7zX6vVkjZCStOrLDUoLhXWsi4ZyIiIiIiIiIup4jJCmtoYhSfpx026lhcJa/sn3kqIj5NGCM3r57Y91+9DJ+v2gcRo3e7GKzrceiw4cPqHxTy7W7weN02/7j9GDj83Tz1syW3yON6bl5Rdq4jNv6paBY3X70Ml6/o2VqrlS2+Jzm5qatOrLH3X3o09G7q/33OIVKiuvjHzOqGmvtni7r+dCSZnuH/2s+j44TY2NTa2u9ZEJz+uuv82K/LhXv9Gav2i51m/cqftHP6te/Uar38PTtfrrn1r8dbHGMu/n/ofB43XLwLF66LH5Wv7J920OVba2RkhJevHNfystFFbuoRMtft2OFpzRlGffbvfXbfE/PtOQ8NP6w+DxuqHPSPUeNkWvvLNaZZcrW3zepbJyvbBkpQb8daZu7D9GA/46U9PmvaOcvONt/tza+32NVUfPWlu/rv/9abseemy+/nTvRN3QZ6Ruu3uCps17R/nHT7f4e3nn9PCxU3r+jZX6832TfJ/TWJWVV+q63iP00lurWn3s94PGtTqH3m0Crh4hM3bkKC0UVsaOlm+7v1BSprRQWB+sWitJOl5YpJvuTNd9I59RXV19i8995pUP9Zs+o5Rz8FiL10tKy3X9HSP0xgefd+jnRERERERERERdFyOkKdYw5A0jH61ZJyk67vzxnonq//B0PfvaMj02c5HSQmH9dfyCFn/tnn2H9dv+YzRszFx9tGad/v3Fj3pkwvNKC4W1aXt25PPSZyzUTXemq1ffUfrL+AWat3C57h/9rNJC4VajySvvrFZaKKzHZi7Sh6u/03OLV6hX31EaOmpOZHjbsmu/7h/9rG4bMkFfrcuIqLlSq39/8aPSQmFt253b4u974uRZpYXC+ufH/4281qvfaN08IF2/7T9Gk+a8qade+kC3D52stFBYa3/cFvm8tsaynZl5+m3/Mbrjgala+N4n+mDVWo2Z3vydfoWnz7X7exFrhHz9/U+VFgpr/6HjkV+3WwY+phv7j7nmr9vEp5do+nPv6o0PPo/cl/CWgY8pfcbCyOfUNzTogfR5unXweL3xwef6+Kuf9OKbq/THeybqtaVr4v59jVVHz1pbv65L//WNxs56PfJruui9T9Xv4ekK3T9FpWUVkc/zzukfBo/vsnMaqw0Ze1p8B6dt7YZtum3IBD342LzIOTx6onkw7ewIKUmfrd2otFBYr7//aeS1nzfvbTHiXt2jk17Uw+Oeu+bPh4iIiIiIiIi6NkZIU6xhyHvde+iKN+4sXf616uqjDxF58c1VSguFdaqoWFLzdywOfvRJPTLh+RYPG2loaNSQ8NP62+MvRF7zxrQNGXsir9U3NGjYmLnq++C0yGvHCouUFgrrqZc+aHGNP2XsbfUE4Fhvxy4tq1CvvqM06/n3W7y++B+f6YY+I1V8sTTyWq9+o3X/6GcjPyepeRC67e4JuvvRJyOvXT2WNTY2aeAjM3XHA1NbDGOStDfniEpKy1tdl62tEbK+oUFDwk/r94PG6cr/vpOyo79usVr1ZfMge/LMeUlSTt5xpYXC+mpdRovPq6q+ooNHCiTF9/saq46etY6+HfvQ0ZNKC4X15XebIq8l4pzG6rV31+i3/ceooaGxzY/3fegJTXn27Vav+xkhJWnynLd0Xe8R2ptzRBcvXdYf75mosbNeV1NT6+/ylaSly7/W9XeMaPXdk0RERERERESU2BghTd4AdOjoSRWcOqes3PzIa5PmRB9oEeuekOs37lJaKKydmc3fDZaXXxj5jsErtXUtvPLOavXqOyoylsS6t+Hzb6zU9XeMiIw7H67+LnKNV9f/4el6ZMLzkR+3d0/IqXPf1k13pquislpS8+AUun+KHn9qSYvP896OfXXzFjb/unhj4tVjWe6hE0oLhfXh6u9a/bUdadDw2bpnxNMqPH1OxwqLtGXXfoWnvKy0UFj//uLHyOd19NdNkioqq/Xuim/08LjndNvdE3T9HSMibw/O3H+kxXXPfvEfkftOXl08v6+x6uhZa2uEbGxs0lfrMjR6+mu6fehk9eo7KvLz+Me//xP5vESc01jNeuF93TvimZgfT9QIWXa5Ur2HTdGdf5mhsbMWqfewKbpUFnvgXvfzDqWFwjp/oTTm5xARERERERFR15fUEbJXv9Gdkqy8EehqMxa8q8qq6HesxRp3tu7e32JA+f6XnW3+/axrPWDljQ8+V1ooHPnOP+8a2xrIvFHKq70RctP2bKWFwvp87a+SpIwd+5QWCre6B2CsEfKfH/83cr9BqfVY5o09G7dmtfn1r9Wg4bNb/VrddGe63l/5nxaDWEd/3crKKzXwkVlKC4U1bd5Srf76J23I2KOly79WWiisvTmHI3/t408tiXy9Rye9qNff/1R5+YWRj8fz+xqrjp61tkbIqXPfjryl+qM167R+487IvSTfW/Ft5PMScU5jNXbWIoWnvBzz44kaISVpR+ZBXdd7hK7rPUK7sw+1e51bdu2POeITERERERERUeJK6gjpvc00HmmhcNKuzxuG/vPDVv3w625t3b2/zSdgxxp3tu3ObTGgfPdT8xC34rP1MX9+3sNhYo1pS/75RYdHyFFPvNrhEdL7zkfvrbbT5r2j24dObvF2XCn2CPneim+VFoo+aTjWCPnrtmvfT7CtBg2frTv/MkPrN+7Show92p19SFXVV1p9Xkd/3ZZ/8r3SQmF9/8vOFp/nvY3djpCNjU3atD1b81//l+4f/Wxk4Ppm/RZJ8f2+xqqjZy3Wd5g+/8bKFmPs5YqqDo+Qfs5prMbOWqSRU1+J+fH4R8h9LT6vvRFyz77Dke9qXb9xZ6uP2zbvzGkxnhMRERERERFRckr6CJmMv6azxbpP39V1dNzZ978RNdZDMmwdHdO870C035knNd/Xr+9DT7R4O/bkOW+pzwOx74u4+B+fKS0U1r4DR9Wr3+gWD/jwijVCzlzwnnr1HRW5rqvHMu/nbu9RGU+xHkxzdR39dXv2tWVKC4Vb3bOwrRHy6s4Wl6jvg9MibzeO5/c1Vh09a1f/uq79cZvSQmHtyDzY4vP8jJBd8fOZ/eI/NGzM3Jgf7+gIuTMz739j4q4Wnxfz7djllerzwDSNnv6aHh73nP4weLyKzl2MeR3e4Grve0pEREREREREiS+pI+Tzb6xMyl/T2bp6hKxvaFDfB6fpz/dNajV6NDU1aevu/ZEfd3RMO1pwRt49C23eW2rtPRgXLF6h3/QZFbnv49V5T8Pu9/B0pYXCOnHybKvPaWuEPHritHr1HdViVLp6LKtvaFDvYVPU98FpKrtc2eKv7+yDadqqo79uby37stV3wF2uqNKkOW+2GCEPHD6hnIPHWv39HpnwvAb/70E88fy+xqqzI+T2vQeUFgpr1Zc/Rj6nsbEp8p2pnRkhu+Lns/C9T/S7u8bG/HhHR8hTRcVKC4W1YPGKyGv1DQ16f+V/2hwhJ815U7cOHq9zxZd0vLBIN/Yfo0cnvRjzOzd5MA0RERERERGRm3gwjamrR0ip+R50vfqO0u8HjdNzi1foozXr9NrSNbrrb7M04K8zI5/X0TFNij7deMz0hfpw9Xea//q/9Js+o3TPiKdbvE17Q8YepYXCenTSi1q2pvnzyspbDoJ/e/wFpYXCGj6x7Scge/flnPLs23p72Veat3C5fnfXWN06eHzkidJS2/cu3JCxR9ffMUJ9HpimRe99qg9WrdVjMxcpLRRW4elz7f4ad/UIebywSL36jtLtQyfrhSUrNX/Rcv35vkm6eUB6ixHS+065sbNe17I132nVlxv0xPylSguF9a/P1kf+/h39fY1VZ0fI2to63fW3WerVb7SeeukDvfLOat0z4unIz6MzI2RX/Hy87yjdf+h4mx/v6AgpRc/kE/OX6rnFKzTwkZn6/aBxrUbINd/8rLRQWN/9tCPy2r+/+LHVr4PtkQnP6y/jF1zz50NEREREREREXRsjpCkRI6QkHTxSoGnz3tGf7p2oXn1Hqf/D0zV17tst7pcYzwjZ1NSklZ//oMGPPqnf9BmlP983SfMXLVdpWUWrv37p8q91xwNT1avfaA0aPrvVCLn665+UFgrrq3UZbf5ce/UbreETX9BjMxfp1sHjdfOAdE146g0d/9+9IL3aGiGl5rfXjpr2qn5311jdMnCsHnpsvlZ+/oNqzc+nrbp6hJSaH8hy38hn9Nv+YzTwkVl6d8U3+u9P21uMkCWl5Xrjg89174jmz/v9oHF6eNxz+u9P21t9jY78vsaqsyOkJBWePq/R01/TzQPSdfvQyXr2tWXKP37a1wjp9+dzuaJK198xQov/8VmbH49nhDxbXKIJT72h3w8apz/eM1HzFy2PPMHbGyHzj5/Wjf3H6In5S1v8tU1NTRox9WXd0Geksg8cbfGx4oulSguFteSfX1zz50NEREREREREXVvgRsjqmtpWwwIlrjmvLtMtA8e2+dAXKfY9IYmubtKcN9X3oSdUV99w7U920Aer1uq63iNaDehERERERERElPgCN0J6D9igxFdRWa2bB6Rr7sKPYn4OIyR1NO8BN//+4sdrf3KSu1RWrt8PGqdJc950fSlEREREREREPbKkj5DPv7GyXXMXfsQImaRWfbkh8nTsWDFCUjw9MX+pbhsyodVbwF33zCsf6rf9x+gY3wVJRERERERE5KSkj5CPTnpR1/UeoVHTXm3T3ye/xAiZpO4Z8bSG/P2pdj+HEZLi6VJZuf5070QNn/hCi/txuuyT/z3AZtma7679yURERERERESUkJI+QhadL1GvfqNjfpy3YxMREREREREREXWvnNwT8v7Rz8b8WGVVTZtPOyYiIiIiIiIiIqLULHAPpiEiIiIiIiIiIqLuFSMkERERERERERERJTTnI2RDQ6N+ytirsvJK15dCRERERERERERECcj5CFlVXaO0UFjZB466vhQiIiIiIiIiIiJKQIyQRERERERERERElNAYIYmIiIiIiIiIiCihBWCEvKKbB6QrJ++460shIiIiIiIiIiKiBOR8hCQiIiIiIiIiIqLuHSMkERERERERERERJTTnI+SPm3a7vgQiIiIiIiIiIiJKYM5GyCu1dXpu8QrdNmSCq0sgIiIiIiIiIiKiJORkhDxacEb3jnhG4Skv61zxJReXQEREREREREREREkq6SPkF//dpN/dNVbvrvhGjY1Nyf7yRERERERERERElOSSOkJ+/NVP+t1dY5W5/0gyvywRERERERERERE5LKkj5MEjBer38HTNWPCuKiqrk/mliYiIiIiIiIiIyFFJfzv25YoqTZ37tvo/PF1ZufnJ/vJERERERERERESU5Jw9HfvTb3/RLQPH6v2V/3F1CURERERERERERJSEnI2QknTk+CkNCT/t8hKIiIiIiIiIiIgowTkdISWp5kqt60sgIiIiIiIiIiKiBOZ8hCQiIiIiIiIiIqLundMRsvD0ef33p+365JufWyEiIiIiIiIiIqLukbsH0/xno66/Y4R69Rut24dOVu9hU1ogIiIiIiIiIiKi7pGzEfKOB6Zq/uv/Ul1dvatLICIiIiIiIiIioiTkbIS8eUC6Mnbsc/XliYiIiIiIiIiIKEk5GyGfmL9US5d/7erLExERERERERERUZJK6gi5cWtWxKff/qJbBo7VD7/ubvG6h4iIiIiIiIiIiLpHSR0hr79jRIcRERERERERERFR98jZ27GJiIiIiIiIiIioZ8QISURERERERERERAmNEZKIiIiIiIiIiIgSGiMkERERERERERERJTRGSCIiIiIiIiIiIkpojJBERERERERERESU0BghiYiIiIiIiIiIKKEFboSsq6vXJ9/87PoyiIiIiIiIiIiIqIsK3Ah5uaJKaaGw68sgIiIiIiIiIiKiLirpI+T6jbva9c36LYyQRERERERERERE3aikj5B/vGei0kJhXX/HiJgYIYmIiIiIiIiIiLpPSR8hcw+dUK9+o2N+nLdjExERERERERERda+SPkI2Njbpxv5jYn6cEZKIiIiIiIiIiKh75eTBNBcvXY75saamJl0oKUvi1RAREREREREREVEiC9zTsYmIiIiIiIiIiKh7xQhJRERERERERERECc35CFlzpVZDR83R4WOnXF8KERERERERERERJSDnI2RVdY3SQmFlHzjq+lI63ZmL1YAqa+pVVlnn/DoQDNVXGnSpotb5dSAYamobVHL5ivPrQDDU1jfqQhnnAc3qG5p0vrTG+XUgGBqbpLMl7q8DwcCfNWFxHoKB/MUI2QW5/ocAwcAICYsREhYjJCxGSFiMkLAYIWHxZ01YnIdgIH8xQnZBrv8hQDAwQsJihITFCAmLERIWIyQsRkhY/FkTFuchGMhfzkfIuvoGfbRmnYovlrq+lE7n+h8CBAMjJCxGSFiMkLAYIWExQsJihITFnzVhcR6CgfzlfITsDrn+hwDBwAgJixESFiMkLEZIWIyQsBghYfFnTVich2AgfzkfIevq6l1fgu9c/0OAYGCEhMUICYsREhYjJCxGSFiMkLD4syYszkMwkL+cjpCFp8/pgfR5Li+hS3L9DwGCgRESFiMkLEZIWIyQsBghYTFCwuLPmrA4D8FA/nI2Qq7dsE2/HzROi9771NUldFmu/yFAMDBCwmKEhMUICYsREhYjJCxGSFj8WRMW5yEYyF9JHyGrqq/o6Zf/qduHTta23bnJ/vIJyfU/BAgGRkhYjJCwGCFhMULCYoSExQgJiz9rwuI8BAP5K6kjZP7x0xo0fLbGznpdJaXlyfzSCc31PwQIBkZIWIyQsBghYTFCwmKEhMUICYs/a8LiPAQD+SupI+SKz9brd3eNVV5+YTK/bELLOdCo7Vk1yjpYo9yjNco/Wa1TxfyHZE/ECAmLERIWIyQsRkhYjJCwGCFhSYxOiOI8BAP5K+lvx/aGyFVf/pjsL52QsnIa9fPWmlY27qhhnOxhGCFhMULCYoSExQgJixESFiMkLInRCVGch2Agfzl5MM2Bwyc08JFZGv/k4pR/W3ZWTqOyCkq0+3CpduRe1ubMCv2yvbrNYZJxsntjhITFCAmLERIWIyQsRkhYjJCwJEYnRHEegoH85ezp2FXVNZr94j90+9DJri6hS8rKadSBc+dbyTl9QVkFJdp5sEzbcsq1ObMi5jDJONk9MELCYoSExQgJixESFiMkLEZIWBKjE6I4D8FA/nI2Qnp9/f1m15fgq1gjZCyMk90XIyQsRkhYjJCwGCFhMULCYoSEJTE6IYrzEAzkL+cjZKoX7wjJONl9MULCYoSExQgJixESFiMkLEZIWBKjE6I4D8FA/kr6CLl+406VllUk+8smrK4aIbt6nNx3uEYHj1Ur/6T7f0h7CkZIWIyQsBghYTFCwmKEhMUICUtidEIU5yEYyF9JHyFv7D9G1/UeoWFj5mrRe59q6+79qrlSm+zL6LISPUJ21TiZsYtxMtEYIWExQsJihITFCAmLERIWIyQsidEJUZyHYCB/JX2ErKur1559h7X0X98oPOVl9eo3Wr36jdbIqa/og1VrlZN3XI2NTcm+rE7naoRknAweRkhYjJCwGCFhMULCYoSExQgJS2J0QhTnIRjIX87vCXmltk47Mg/qrWVfavjEF9Wr7yjdOni868vqcEEbIRkn3WGEhMUICYsREhYjJCxGSFiMkLAkRidEcR6CgfzlfISsqq5Rxo4cLXrvUz08dr5u6DNSdz/6pOvL6nCpMkK2N05mHruknQfLtCW7XJt2VzJOdhIjJCxGSFiMkLAYIWExQsJihIQlMTohivMQDOSvpI+Q1TW12rY7V2988Ln+On6BbugzUv0fnq45ry7T2g3bdKGkLNmX5KtUHyFjyT55kXEyToyQsBghYTFCwmKEhMUICYsREpbE6IQozkMwkL+SPkL26jtKv+kzStOfe1dffrdJZ85eSPYldGnddYRknIwfIyQsRkhYjJCwGCFhMULCYoSEJTE6IYrzEAzkr6SPkHNeXaY7/zJDv7trrMY/uVjLP/leBw6fcP4wmssVVRo2Zq7Wb9zV4vWMHTlKC4Vb6NVvdOTjPW2ETMQ4eaigWsfOuP+XiV+MkLAYIWExQsJihITFCAmLERKWxOiEKM5DMJC/nN0TsujcRX39/WY99dIH6vvgNP1h8Hg9/tQSrfhsvfLyC5N6LUv++YV6D5uiG/qMbHOE7P/wdF2uqIoor6iKfJwRsmvHya17a7RrX2qOk4yQsBghYTFCwmKEhMUICYsREpbE6IQozkMwkL+cP5jGq/D0OS1YvEK9+o5SWijs5BoeSJ/X5gg54K8zY/41jJCMkx5GSFiMkLAYIWExQsJihITFCAlLYnRCFOchGMhfTkfIoydOa9WXGzRpzpu6dfB4pYXCGhJ+Wi+//bGT64k1Qt7QZ6RuHzpZdz/6pBYsXqHLfCck42QbGCFhMULCYoSExQgJixESFiMkLInRCVGch2AgfyV9hPx87a+aseBd3T50stJCYf35vkma9fz7+mb9FhVfLE325bSorRHyUlm58vILVXS+RJn7j+hvj7+gqXPfjnw8K6dRpyouIcGOX7yso+fKlVdYpX2Hrigzt067sxtiyj3UoCPHGnSyqFHFJY0qq2hUbX1iNTQ2qaGhKeFfB6mhsbFJ9Q2J/RrVVxqc/zwRnPOA1NHY1KS6FD8P1bX8+6erNDVJdQG4DgRDkzgP19KT/v0jyfk1IDg4D8FA/kr6CHlj/zFKn7FQyz/5XoeOnkz2l2+3tkbIq8vcf0Q39BmppqbmB+lk5TTq8IViOLL/5CVlHy/T3sOXtT27Ulv2VGvTjisx7cy+oqwDV5R3rFbHT13RqfNXdKGsa9TUNqiypr7L/n5IbVfqGlVendjzUByAnyc6pra+UZcr65xfB4KhvqFJZZwH/E9DY5MuVdQ6vw4EQ1OTdPGy++sIsp703z+SnF8DgoPzEAzkr6SPkLW1dcn+kh2uIyPkll379cd7JkZ+zNuxgymroER78pvf1r05s0K/7qxK+Nu6eTs2LN6ODYu3Y8Oqreft2Iji7diweDs2LIm33yKK8xAM5C9n94RsbGzSzsw8rfnmZ6355mdl7j/i6lIitTVCfrj6O63fuFNni0uUl1+oYWPm6qW3VkU+zgiZWvyOkwVnY/8hgRESFiMkLEZIWIyQsBghYTFCwpIYnRDFeQgG8peTETIvv1CDhs9WWiis6+8Yoet6j1BaKKwH0ufpeGFR0q/nhSUrdduQCbqhz0jdMnCsbhsyIXJ/yq/WZeieEU/rpjvTFbp/il58c5Wqa2ojfy0jZPcQ7zi5PatGe3Kbx8n8k83jJCMkLEZIWIyQsBghYTFCwmKEhCUxOiGK8xAM5K+kj5DnL5TqT/dO1N8nv6ScvONqaGi+Yf+Bwyc0cuor6j1sii5eupzsy+p0jJDdWzzj5Pa9dco+UNdqnHT9L0m4wQgJixESFiMkLEZIWIyQsCRGJ0RxHoKB/JX0EfKVd1Zr0PDZutLGvSGv1Nbprr/N0tLlXyf7sjodI7Yus0gAACAASURBVGTP1NY4uX1vnbburuvwd066/pcnEosREhYjJCxGSFiMkLAYIWFJjE6I4jwEA/kr6SPkkL8/pRWfrY/58Y/WrNOQ8NNJvCJ/MULCU1B2UXlnSn2/rdv1v1TRNRghYTFCwmKEhMUICYsREpbE6IQozkMwkL+SPkLePCBdWbn5MT++Zdd+3TwgPYlX5C9GSHgKyi7q2KWLrV7PLSpWVkGJdh8u1Y7cy9qcWaFftlczTnZzjJCwGCFhMULCYoSExQgJS2J0QhTnIRjIX0kfIdNCYR0+dirmxw8cPqG0UDh5F+QzRkh4Yo2QsTBOdm+MkLAYIWExQsJihITFCAlLYnRCFOchGMhfTkbItRu2KfvA0TZ9+8MWRkikpHhHyFgYJ7sHRkhYjJCwGCFhMULCYoSEJTE6IYrzEAzkLycjZEekSoyQ8HTVCBlLvOPkxh3N42TWwRrlHm0eJ08V84ecZGGEhMUICYsREhYjJCxGSFgSoxOiOA/BQP5K+gh5+uyFDkmVGCHhSfQIGQvjZDAxQsJihITFCAmLERIWIyQsidEJUZyHYCB/JX2E7G4xQsLjaoSMhXHSLUZIWIyQsBghYTFCwmKEhCUxOiGK8xAM5C9nI+RHa9bFtOrLDdq0PVu1tXWuLq/DMULCE7QRMhbGyeRghITFCAmLERIWIyQsRkhYEqMTojgPwUD+cjZC/t89j+vWweM1JPx0C2mhsAY+Mkt/GDxefR+cFvi3ZjNCwpMqI2Qs3ji582CZtuWUa3NmRbsPw2GcbB8jJCxGSFiMkLAYIWExQsKSGJ0QxXkIBvKXsxHy/tHP6qt1Ga1en7HgXS355xeqq2/QjAXvavKctxxcXcdjhIQn1UfIWHJOX2Cc7ARGSFiMkLAYIWExQsJihIQlMTohivMQDOQvZyPkb/uP0f5Dx1u9/vOWTP35vkmSpG27c3X70MnJvrS4YoSEp7uOkLEwTraPERIWIyQsRkhYjJCwGCFhSYxOiOI8BAP5y9kI+ef7JrX5nZA/b96rXv1GS5J2ZeXpt/3HJPvS4ooREp6eNkLG0tlxct/hGh08Vq38k+7/h6UrMELCYoSExQgJixESFiMkLInRCVGch2AgfzkbIf/58X9125AJ+mpdhk4VFev8hVJl7Ninu/42S2NnvS5J+mztRg3460xXl9ihGCHhYYRsX7zjZMau1B4nGSFhMULCYoSExQgJixESlsTohCjOQzCQv5yNkI2NTXp3xTe66c50pYXCEROeekPFF0slSRk7crR+4y5Xl9ihGCHhYYTsnO46TjJCwmKEhMUICYsREhYjJCyJ0QlRnIdgIH85GyG9qqprlH3gqPYdOKqLly67vpy4Y4SEhxGya6X6OMkICYsREhYjJCxGSFiMkLAkRidEcR6CgfzlfIRM9Rgh4WGETA47Tm7JLtem3ZWBHCcZIWExQsJihITFCAmLERKWxOiEKM5DMJC/GCF9xggJDyOkW9knLyrz2KXAjJOMkLAYIWExQsJihITFCAlLYnRCFOchGMhfjJA+Y4SEhxEymPyOk8fOdO5/nBghYTFCwmKEhMUICYsREpbE6IQozkMwkL8YIX3GCAkPI2RqiXec3Lq3Rrv2NY+ThwquPU4yQsJihITFCAmLERIWIyQsidEJUZyHYCB/MUL6jBESHkbI7qGrxklGSFiMkLAYIWExQsJihIQlMTohivMQDOQvRkifMULCwwjZvcU7Tu47WK99ebUd/s5JdG+MkLAYIWExQsJihIQlMTohivMQDOSvpI6QL721qsNSJUZIeBghe6asghLtyW8eJzdnVujXnVX6eWuNdu6t1+Zdte1+52T+yWoVnOUPnj0BIyQsRkhYjJCwGCFhSYxOiOI8BAP5K6kjZPqMhR2WKjFCwsMICSv/3GXlFJS1Gidj2Z5Voz25jJPdFSMkLEZIWIyQsBghYUmMTojiPAQD+Yu3Y/uMERIeRkhYhZdLdLTkQqvXY33nJONk98YICYsREhYjJCxGSFgSoxOiOA/BQP5yPkJWVdfo4JECHTxSoKrqGteXE3eMkPAwQsKKNULGwjjZvTFCwmKEhMUICYsREpbE6IQozkMwkL+cjZB1dfV65Z3V+k2fUUoLhZUWCqtX31F6beka1dU3uLqsuGOEhIcREla8I2RbcouKlVVQot2HS7Uj97I2Z1bol+3VHRonc482j5OnivmDbRAwQsJihITFCAmLERKWxOiEKM5DMJC/nI2Qr76zWn0fnKYfN+3WueJLOld8Ses37tIdD0zVa++ucXVZcccICQ8jJKyuGCFjiXec3LijeZzMOsg46QojJCxGSFiMkLAYIWFJjE6I4jwEA/nL2Qj5p3snauPWrFav/7wlU3++b5KDK+pcjJDwMELCSuQIGQvjZHAxQsJihITFCAmLERKWxOiEKM5DMJC/nI2QvfqOUl5+YavX8/IL1avfaAdX1LkYIeFhhITlYoSMJd5xMmNX8zi573CNDh6rVv5J9/9jn+oYIWExQsJihITFCAlLYnRCFOchGMhfzkbIh8c9p5ff/rjV68+/sVJ/Gb/AwRV1LkZIeBghYQVphIwl5/QFZRWUaOfBMm3LKdfmzIp27zd59Th57Iz7/whIFYyQsBghYTFCwmKEhCUxOiGK8xAM5C9nI+Tu7EO6sf8YDR01R3MXfqS5Cz/SvSOe0U13pmtvzhFXlxV3jJDwMELCSoURMparx8lNuyvbHSe37q3Rrn3N4+ShAsbJtjBCwmKEhMUICYsREpbE6IQozkMwkL+cjZCSdLa4RK+9u0ZjZy3S2FmLtPC9T3Su+JLLS4o7Rkh4GCFhpfIIGUv2yYvKPHZJOw+WaUt2fONk/slqFZztuX/IZoSExQgJixESFiMkLInRCVGch2AgfzkdIbtDjJDwMELC6o4jZCx2nNycWaFfd1a1O05uz6rRntyeNU4yQsJihITFCAmLERKWxOiEKM5DMJC/nI6QxwuLtHbDNn3yzc+tpEqMkPAwQsLqSSNkLFkFJdqTH9842V2f1M0ICYsREhYjJCxGSFgSoxOiOA/BQP5yNkIu/+R7pYXCumXgY+r74DT1feiJFlIlRkh4GCFhMULGFs+Tujfu6B7jJCMkLEZIWIyQsBghYUmMTojiPAQD+cvZCPnn+yZp4XufqLGxydUldEmMkPAwQsJihIxPblFxXOPk1U/qzj/p/j9I2sMICYsREhYjJCxGSFgSoxOiOA/BQP5yNkLeMnCsMnbkuPryXRYjJDyMkLAYIbvG1U/q3pxZ0e5buq8eJ4PypG5GSFiMkLAYIWExQsKSGJ0QxXkIBvKXsxFy1gvv69V3Vrv68l0WIyQ8jJCwGCET6+pxMp4ndR8qSP44yQgJixESFiMkLEZIWBKjE6I4D8FA/nI2Ql4oKdMf75mopf/6Ru+v/I/eW/FtC6kSIyQ8jJCwGCHdsE/q3pId3ziZyCd1M0LCYoSExQgJixESlsTohCjOQzCQv5yNkO989JXSQmGF7p+i+0Y+o6Gj5rSQKjFCwsMICYsRMljsONnRJ3Xvye26cZIREhYjJCxGSFiMkLAkRidEcR6CgfzlbIT8w+Dx+se//+Pqy3dZjJDwMELCYoRMDVkFJdqTH9842ZkndTNCwmKEhMUICYsREpbE6IQozkMwkL+cjZC3Dh7Pg2nQrTBCwmKETG3xPKl7445rj5OMkLAYIWExQsJihIQlMTohivMQDOQvZyPkS2+t0tyFH7n68l0WIyQ8jJCwGCG7n9yi4rjGSfuk7qLz9Tp5ltEJzRghYTFCwmKEhCUxOiGK8xAM5C9nI+T8Rct1Xe8RemHJSr22dE0rqRIjJDyMkLAYIXuOq5/UvTmzotUguTOzXhk7r7QaJw8eq1b+Sff/MYXkYoSExQgJixESlsTohCjOQzCQv5yNkKOmvdquVIkREh5GSFiMkLDjZE5+tbZnxf6uyauf1H2ooFrHzrj/jywkBiMkLEZIWIyQsCRGJ0RxHoKB/OVshOwuMULCwwgJixES1snLJcq/WKwD51o+qXtLdrk27a6Ma5z0+6RuuMcICYsREhYjJCyJ0QlRnIdgIH8xQvqMERIeRkhYjJCw7AgZix0nO/qk7j25zeNk/knGyVTCCAmLERIWIyQsidEJUZyHYCB/ORshz5y9oJkL3tPAR2bpz/dNaiVVYoSEhxESFiMkrI6MkLFkFZRoT37842SsJ3XDPUZIWIyQsBghYUmMTojiPAQD+cvZCDl84osa+MhMfbj6O321LqOVVIkREh5GSFiMkLD8jJCxxPOk7o07msfJrIOMk0HACAmLERIWIyQsidEJUZyHYCB/ORshb+w/Rj9u2u3qy7fqckWVho2Zq/Ubd7V4vbGxSa++s1q3Dh6vmweka9YL76vmSm3k44yQ8DBCwmKEhJWIEbItuUXFnRoneVJ3cjFCwmKEhMUICUtidEIU5yEYyF/ORsi7H31Sa3/c5urLt2jJP79Q72FTdEOfka1GyDXf/KyBj8zUscIinS0u0fCJL2jhe59EPs4ICQ8jJCxGSFjJGiFjsU/q3pZTrs2ZFe2+pTtjF+NkIjFCwmKEhMUICUtidEIU5yEYyF/ORsi1G7bp3hHPqK6+wdUltOqB9HmtRsi/T35JKz5bH/nxr9uy1eeBaZEfM0LCwwgJixESlusRMparx8l4n9R97Iz7/xBMRYyQsBghYTFCwpIYnRDFeQgG8pezEfLEybO6dfB49eo3Wjf0Ganr7xjRgovaGiFvHzpZm7ZnR35ccOqc0kJhVdc0vyWbERIeRkhYjJCwgjpCxmKf1L0lO/5xkid1t48REhYjJCxGSFgSoxOiOA/BQP5yNkI+kD5P9418Rp+t3agNGXv0U8beFlxd09Uj5O/uGqsdmQcjPz5XfElpobBKSsslNY+QR0suADp5uUQFZRedXweC4VT5pSSdh4tIAacrvfOQ2vLOlurAyTJl55dr1/5qbc+8oi276mLKyq1T7pE6HTlRp9Pn61R8qU6lFWhobFJFdb3z6/DjUkWt82voLhobpfIq99eBYGiSVFbp/jqCrKSiVpd6CEnJ+3rlCDpJzq8BtSJ/ORshe/Ud1eI7DINQrO+E3LwzJ/Ljtr4TsuDyRUBnKi/pVMUl59eBJCqL7UzlJZ0sL2n3c3wrDcCvATrkTGVp8/9REYBrSYQj58qUd6pcuScqlHmwRrv3XdGOvXUx5R6q1+Hj9TpZVK/iknpVVNer6krP0dDYpJraBufX4UdP+z1LpMamJlWn+HlA12lqkqoDcB1BVlFdr6qankFS8r5eAH5v0T5Jzq8B9Ulap7pvzkbIYWPm6svvNrn68m3W1gj56KQXtfLzHyI//mVrpkL3T4n8mLdjw8PbsWHxdmxYqfZ27K7SmSd1Zx2sUe7RGuWfrNap4u75FlXejg2Lt2PD4u3YsCTefosozkMwkL+cjZAbt2ZpwF9nqrSsQlXVNa24qK0RctWXP2rgI7N04uRZnb9QquETX9RLb62KfJwREh5GSFiMkLB66gjZltyi4k6Nk93pSd2MkLAYIWExQsKSGJ0QxXkIBvKXsxEyLRRuVzJ7YclK3TZkgm7oM1K3DByr24ZMUPHFUklSQ0OjXnxzlf4weLxuHpCuGQvebTGSMkLCwwgJixESFiPktV39pO7NmRXtPgwnY1fqjpOMkLAYIWExQsKSGJ0QxXkIBvKXsxFyb87hdqVKjJDwMELCYoSExQjZeVePk/E+qfvYGff/sXo1RkhYjJCwGCFhSYxOiOI8BAP5y9kI2V1ihISHERIWIyQsRsiul33yojKPXdLOg2Xakh3/OFlw1t3owwgJixESFiMkLInRCVGch2AgfyV1hMzLL0zKX5PMGCHhYYSExQgJixEyeew4uTmzQr/urGp3nNyeVaM9uc3jZP7J5IyTjJCwGCFhMULCkhidEMV5CAbyV1JHyM7c6zHZ94eMN0ZIeBghYTFCwmKEdC+roER78uMfJxPxpG5GSFiMkLAYIWFJjE6I4jwEA/mLEdJnjJDwMELCYoSExQgZXJ15UnfWQX/jJCMkLEZIWIyQsCRGJ0RxHoKB/JX0EXJI+Om4MEIiVTBCwmKEhMUImVpyi4o7NU529EndjJCwGCFhMULCkhidEMV5CAbyV1JHyP/8sLVTghwjJDyMkLAYIWExQnYPVz+pe3NmRbtv6c7Y1fY4yQgJixESFiMkLInRCVGch2Agf/F0bJ8xQsLDCAmLERIWI2T3dvU4ea0ndecealDWgSuRJ3UfO+P+P6jhDiMkLEZIWBKjE6I4D8FA/mKE9BkjJDyMkLAYIWExQvZM9kndW7Kj4+Tu7AZt2nGl1Ti5dW+Ndu2riYyTyXhSN9xjhITFCAlLYnRCFOchGMhfjJA+Y4SEhxESFiMkLEZIWMcvXlZOQWncT+red7j5YTiMk90LIyQsRkhYEqMTojgPwUD+YoT0GSMkPIyQsBghYTFCwjpVcUmHL7Q+D1kFJdqTfynucdLPk7rhHiMkLEZIWBKjE6I4D8FA/mKE9BkjJDyMkLAYIWExQsKKNULG0pkndWcd7NiTuuEeIyQsRkhYEqMTojgPwUD+YoT0GSMkPIyQsBghYTFCwop3hGxLblFxq3Eynid18zCc4GCEhMUICUtidEIU5yEYyF+MkD5jhISHERIWIyQsRkhYXTFCxhLvk7p5GI57jJCwGCFhSYxOiOI8BAP5ixHSZ4yQ8DBCwmKEhMUICSuRI2QssZ7U3dGH4XC/ycRhhITFCAlLYnRCFOchGMhfjJA+Y4SEhxESFiMkLEZIWC5GyFjieRiOvd8kD8PpOoyQsBghYUmMTojiPAQD+YsR0meMkPAwQsJihITFCAkrSCNkW9q632R7D8O5+n6TPAwnPoyQsBghYUmMTojiPAQD+YsR0meMkPAwQsJihITFCAkr6CNkLN446d1v8loPw7n6fpM8DKdtjJCwGCFhSYxOiOI8BAP5ixHSZ4yQ8DBCwmKEhMUICStVR8hYck5fiOt+k1v3trzfZE9/GA4jJCxGSFgSoxOiOA/BQP5yPkJWVV/RR2vWqbGxyfWldCpGSHgYIWExQsJihITV3UbIWOzDcK51v0nvYTg98X6TjJCwGCFhSYxOiOI8BAP5y/kIuerLH5UWCmv9xl2uL6VTMULCwwgJixESFiMkrJ4yQsbiPQynI/eb9B6G053vN8kICYsREpbE6IQozkMwkL+cjpANDY268y8z9NjMRXp43HMuL6XTMULCwwgJixESFiMkrJ4+QralrYfhtPddk/ZhOKl+v0lGSFiMkLAkRidEcR6CgfzldIRcv3Gn+j88XZVVNbp18HjtzTns8nI6FSMkPIyQsBghYTFCwmKE7Lic0xdaPAynI/eb9B6Gkyr3m2SEhMUICUtidEIU5yEYyF9OR8iHxz2nlZ//IEl6/f1PNfHpJS4vp1MxQsLDCAmLERIWIyQsRkj/7P0mt2SXp/T9JhkhYTFCwpIYnRDFeQgG8pezEXLPvsP63V1jVVlVI0k6W1yiXn1H6cTJs64uqVMxQsLDCAmLERIWIyQsRsjE8e436T0MJxXuN8kICYsREpbE6IQozkMwkL+cjZCPP7VEr76zusVrMxa8q/mLlru5oE7GCAkPIyQsRkhYjJCwGCGTqzP3m/Te0p2M+00yQsJihIQlMTohivMQDOQvJyPk8cIi/abPKJ05e6HF6/sPHdeN/ceopLTcxWV1KkZIeBghYTFCwmKEhMUIGQydud/kntyuv98kIyQsRkhYEqMTojgPwUD+cjJCHjh8Qus37mzzY1+ty9CxwqIkX1HnY4SEhxESFiMkLEZIWIyQwZbs+00yQsJihIQlMTohivMQDOQvpw+m6Q4xQsLDCAmLERIWIyQsRsjUlKj7TTJCwmKEhCUxOiGK8xAM5C9GSJ8xQsLDCAmLERIWIyQsRsjuoyvuN8kICYsREpbE6IQozkMwkL+cj5A1V2o1dNQcHT52yvWldCpGSHgYIWExQsJihITFCNn92ftNbslu/36Te7IbtCPrSkLuN4nUwwgJS2J0QhTnIRjIX85HyKrqGqWFwso+cNT1pXQqRkh4GCFhMULCYoSExQjZc9n7TW7OrNCvO6u0J7tBv25PzP0mkXoYIWFJjE6I4jwEA/mLEdJnjJDwMELCYoSExQgJixES1pnKUuWcbL7fpPeW7q663yRSDyMkLInRCVGch2AgfzFC+owREh5GSFiMkLAYIWExQsI6U1mqQ22cB3u/yW055Z2636TrP6ghfoyQsCRGJ0RxHoKB/OV8hKyrb9BHa9ap+GKp60vpVIyQ8DBCwmKEhMUICYsRElasETKWeO43+fPWGm3dW8P9JlMIIyQsidEJUZyHYCB/OR8hUz1GSHgYIWExQsJihITFCAkr3hEylrbuN9neOMn9JoOJERKWxOiEKM5DMJC/nI+QDQ2N+iljr8rKK11fSqdihISHERIWIyQsRkhYjJCwumqEjCWrgPtNphJGSFgSoxOiOA/BQP5yPkJyT0h0F4yQsBghYTFCwmKEhJXoEbIt3G8yuBghYUmMTojiPAQD+YsR0meMkPAwQsJihITFCAmLERKWixEyFu436R4jJCyJ0QlRnIdgIH8xQvqMERIeRkhYjJCwGCFhMULCCtIIGQv3m0weRkhYEqMTojgPwUD+CsAIeUU3D0hXTt5x15fSqRgh4WGEhMUICYsREhYjJKxUGCFjied+kxm7uN9kRzBCwpIYnRDFeQgG8pfzETLVY4SEhxESFiMkLEZIWIyQsFJ5hGxLvPeb3LqX+01ajJCwJEYnRHEegoH8xQjpM0ZIeBghYTFCwmKEhMUICau7jZCxxHu/ye1ZPfN+k4yQsCRGJ0RxHoKB/MUI6TNGSHgYIWExQsJihITFCAmrp4yQscRzv8mNO7r//SYZIWFJjE6I4jwEA/mLEdJnjJDwMELCYoSExQgJixESVk8fIWPpqfebZISEJTE6IYrzEAzkL0ZInzFCwsMICYsREhYjJCxGSFiMkB3XE+43yQgJS2J0QhTnIRjIX4yQPmOEhIcREhYjJCxGSFiMkLAYIf3LOX0h8pbuVL/fJCMkLInRCVGch2AgfzkdIauqa5SxI0efrd0oSWpoaNSurDxdKClzeVlxxQgJDyMkLEZIWIyQsBghYTFCJk72yYvak59a95tkhIQlMTohivMQDOQvZyNkXn6hQvdP0W13T1BaKBx5/fGnlujZ15a5uqy4Y4SEhxESFiMkLEZIWIyQsBghk897S3cQ7zfJCAlLYnRCFOchGMhfzkbI8JSX9eo7q9XU1NRihPwpY6/u/MsMV5cVd4yQ8DBCwmKEhMUICYsREhYjZDB495vcebDM6f0mGSFhSYxOiOI8BAP5y9kIecvAsTpVVCxJLUbIvPxC9eo7ytFVxR8jJDyMkLAYIWExQsJihITFCBlsyb7fJCMkLInRCVGch2AgfzkbIW+7e4IOHzslqeUI+d+ftqv3sCmOrir+GCHhYYSExQgJixESFiMkLEbI1JSo+00yQsKSGJ0QxXkIBvKXsxHymVc+1LR576iuviEyQh4tOKM7/zJDzy1e4eqyWpWxI0dpoXALvfqNjnycERIeRkhYjJCwGCFhMULCYoTsXvzeb5IREpbE6IQozkMwkL+cjZBllyv10GPz9ad7JyotFNaAv87Udb1H6P7Rz6qsvNLVZbUqY0eO+j88XZcrqiLKK6oiH2eEhIcREhYjJCxGSFiMkLAYIbu/eO43mZnTqG0Jut8kUo/E6IQozkMwkL+cjZCS1NDQqE3bs/Xuim+0dPnX+iljr+obGlxeUqsyduRowF9nxvw4IyQ8jJCwGCFhMULCYoSExQjZc7V1v8nMnEb9si0x95tE6pEYnRDFeQgG8pfTEfJa/WX8AteXoIwdObqhz0jdPnSy7n70SS1YvEKX+U5ItIEREhYjJCxGSFiMkLAYIWEVVZUq55S/+026/gM6unbscH0NCA7OQzCQvwI9Qg4aPtv1JehSWbny8gtVdL5EmfuP6G+Pv6Cpc9+OfDwrp1FFlaWAzlaVqqjK/XUgGJJxHk6WX2r+Ggi8s9XmTKDH6w7n4RT//uky3eE8oOucrY79sYKL5co/U6VDhdXaf6RWWfvrlZnT0KbcvAYdOdaoovNNulDSqMrqJjU2dQ9X6hrU2NjUI0hyfg0IDs5DMJC/GCHjLHP/Ed3QZ6SampoPX1ZOow6ePw+ooOyijl+66Pw6EAzed0Im+uvknS9GCjj5v/Pg+joQDKcqLunIxRQ/D8X8+6ernKks1eEL7q8DwVBUVapDxR3//ANnL2hfQYl25ZVp+/5ybcms1C/bamLatveKdufUaP+RKzpSWKOT567obElNSikqqe4xJDm/BgQH5yEYyF+MkHG2Zdd+/fGeiZEf83ZseHg7Nizejg2Lt2PD4u3YsHg7NqyiqlLlnff/92nrfpOx3tLd1v0mTxVzv8kgkHj7LaI4D8FA/mKEvEYfrv5O6zfu1NniEuXlF2rYmLl66a1VkY8zQsLDCAmLERIWIyQsRkhYjJCwumqEjCX7JPebTCUSoxOiOA/BQP5ihLxGX63L0D0jntZNd6YrdP8UvfjmKlXX1EY+zggJDyMkLEZIWIyQsBghYTFCwkr0CBlLVkGJdh8u1Y7cy9qcWaFftlfHHCczdjWPk/sO1+jgsWodO+N+FOiuJEYnRHEegoH8xQjpM0ZIeBghYTFCwmKEhMUICYsREparEbItOacvKKugRDsPlmlbTrk2Z1a0+5burXtrtGtf8zh5qKBaBWd5S3dXjB2urwHBwXkIBvIXI6TPGCHhYYSExQgJixESFiMkLEZIWEEaIWPJPnmR+00mcexwfQ0IDs5DMJC/GCF9xggJDyMkLEZIWIyQsBghYTFCwkqFETKWrIIS7jeZgLHD9TUgODgPwUD+cjZCHjl+Sp+vGZYo+QAAIABJREFU/VWNjU2R1xobm/Tld5t09MRpSVJdXb2ry+twjJDwMELCYoSExQgJixESFiMkrFQeIduSW1TM/SZ9jh2urwHBwXkIBvKXsxFy3OzFmr9oeavXX3prlSbNeTP5F9TJGCHhYYSExQgJixESFiMkLEZIWN1thIyF+012fOxwfQ0IDs5DMJC/nI2Qtw2ZoJ2Zea1e35tzRH+8Z6KDK+pcjJDwMELCYoSExQgJixESFiMkrJ4yQsbC/SZbjx2urwHBwXkIBvKXsxHypjvTte/A0Vav7ztwVLcNmeDgijoXIyQ8jJCwGCFhMULCYoSExQgJq6ePkLH01PtNSoxOiOI8BAP5y9kIOWb6Qs16/v1W94R88sUPNHr6a64uK+4YIeFhhITFCAmLERIWIyQsRkhYjJAd1xPuNykxOiGK8xAM5C9nI2RefqF+d9dYDRo+W8++tkzPvrZMdz/6pG66M11ZufmuLivuGCHhYYSExQgJixESFiMkLEZIWIyQ/nWn+01KjE6I4jwEA/nL2QgpSafPXtCr76xWeMrLGjH1Zb2wZKWOFxa5vKS4Y4SEhxESFiMkLEZIWIyQsBghYTFCJk4q3m9SYnRCFOchGMhfzkbI8ooqHSssavV27GOFRaqoTJ3fWEZIeBghYTFCwmKEhMUICYsREhYjZPJ19n6TB49VJ/x+kxKjE6I4D8FA/nI2Qs55dZlGTXu11euPP7VEzy1e4eCKOhcjJDyMkLAYIWExQsJihITFCAmLETIY/Nxv8lBB191vUmJ0QhTnIRjIX85GyDsemKpftma2en3T9mz1fegJB1fUuRgh4WGEhMUICYsREhYjJCxGSFiMkMHm536T+Sfjv9+kxOiEKM5DMJC/nI2QvfqOUl5+YavX8/ILdcvAxxxcUedihISHERIWIyQsRkhYjJCwGCFhMUKmps7ebzL3aPv3m5QYnRDFeQgG8pezEXLoqDn6x7//0+r1D1at1YOPzUv+BXUyRkh4GCFhMULCYoSExQgJixESFiNk93L1/Sbbe0t3W/eblBidEMV5CAbyl7MR8odfd+uGPiM1d+FHWvvjNq39cZvmL1quG/qM1Nffb3Z1WXHHCAkPIyQsRkhYjJCwGCFhMULCYoTs/tq632SsYTIrpzFh95tE6pEYIYOA/OVshJSkHzft1rAxc3X9HSN0Xe8RumfE0/r2hy0uLynuGCHhYYSExQgJixESFiMkLEZIWIyQPdfV95vctLtSWTmNCbvfJFKPxAgZBOQvpyOkV2Njk5qamlxfRqdihISHERIWIyQsRkhYjJCwGCFhMULCOltd2up+k7/urOqS+00i9UiMkEFA/nI2Ql4oKWtXqsQICQ8jJCxGSFiMkLAYIWExQsJihIR1tro05sf83m/S9ZCDzo1frq8BjJB+czZCpoXC7UqVGCHhYYSExQgJixESFiMkLEZIWIyQsNobIdsSz/0mf95aw/0mU4zECBkE5C9nI2RefmEr+w4c1W1DJmjthm2uLivuGCHhYYSExQgJixESFiMkLEZIWIyQsOIdIWNp636T7Y2T3G8ymCRGyCAgfwXinpC2ZWu+09hZi1xfRodjhISHERIWIyQsRkhYjJCwGCFhMULC6qoRMhbuN5laJEbIICB/BW6E3LQ9Wzf2H+P6MjocIyQ8jJCwGCFhMULCYoSExQgJixESVqJHyFi432QwSYyQQUD+CswI2dTUpAslZZo2b6kGPjLT9eV0OEZIeBghYTFCwmKEhMUICYsREhYjJCxXI2RbuN+kexIjZBCQvwL3YJpbBo7Vz5v3urqsuGOEhIcREhYjJCxGSFiMkLAYIWExQsIK0ggZC/ebTO745foawAjpN2cj5MatWa3szTmsyxVVri6pUzFCwsMICYsREhYjJCxGSFiMkLAYIWGlwggZC/ebTMz45foawAjpt8C8HTtVY4SEhxESFiMkLEZIWIyQsBghYTFCwkrlETIW7jfpb/xyfQ1ghPSbsxHy2x+26LWlayI/XrB4hf7vnsc1bvZiXSord3VZcccICQ8jJCxGSFiMkLAYIWExQsJihITVHUfItnTmfpPeW7p70v0mJUbIICB/ORshh46ao4+/+kmSlLEjRzf0GakVn63XmOkL9cwrH7q6rLhjhISHERIWIyQsRkhYjJCwGCFhMULC6ikjZCydud/kntzue79JiREyCMhfzkbImwek6/CxU5Kkhe99ovQZCyVJB48U6Pahk11dVtwxQsLDCAmLERIWIyQsRkhYjJCwGCFh9fQRMpbO3G8y62Dq329SYoQMAvKXsxHy/+55XEcLzkiS/vb4C3p72VeSpIJT53TTnemuLivuGCHhYYSExQgJixESFiMkLEZIWIyQsBgh4+Pdb9J7S3dH7je573Dq3G9SYoQMAvKXsxFy2rylGv/kYv3z4/8qLRRWVm6+JOmrdRka8venXF1W3DFCwsMICYsREhYjJCxGSFiMkLAYIWExQvpn7ze5Lac8pe83KTFCBgH5y9kIef5CqUZMfVm3DZmg19//VJJU39CggY/MjPw4FWKEhIcREhYjJCxGSFiMkLAYIWExQsJihEwce7/JLdmpcb9JiREyCMhfzkbI7hIjJDyMkLAYIWExQsJihITFCAmLERIWI2Ty2ftNbs6sCNT9JiVGyCAgfzFC+owREh5GSFiMkLAYIWExQsJihITFCAmLETI4gnC/SYkRMgjIX4yQPmOEhIcREhYjJCxGSFiMkLAYIWExQsJihAy2ZN9vUmKEDALyFyOkzxgh4WGEhMUICYsREhYjJCxGSFiMkLAYIVNTou43KTFCBgH5ixHSZ4yQ8DBCwmKEhMUICYsREhYjJCxGSFiMkN2L3/tNNjQyQgYB+YsR0meMkPAwQsJihITFCAmLERIWIyQsRkhYjJA9Q0fvN5mV05iw+02CETJZMUL6jBESHkZIWIyQsBghYTFCwmKEhMUICYsRsudq636TWTmNCbvfJBghkxUjpM8YIeFhhITFCAmLERIWIyQsRkhYjJCwGCFhna0uTdj9JsEImaycjZB/n/ySRk59pdXrI6e+or9PfsnBFXUuRkh4GCFhMULCYoSExQgJixESFiMkLEZIWO2dB7/3mzxVzDjJCJmcnI2QvYdNUd+Hnmj1et+HnlDvYVMcXFHnYoSEhxESFiMkLEZIWIyQsBghYTFCwmKEhNWZ89DR+03+vLWG+00yQiYl3o7tM0ZIeBghYTFCwmKEhMUICYsREhYjJCxGyP/f3p2/RVU+bhz/f7Ci1Ba/piiIG7lmYqUOsogYKIqiuGvhAn7ctywzi9xSM8ldU8kFRRE1RRQX3BUXFNln5v7+QDOdERTGYeaM9H5d131dzZmFB3yYnNvzPIcY01Tzob79Jhta0s1+k5SQTYUS0kOUkMQRSkhiDCUkMYYSkhhDCUmMoYQkxlBCEmMoIYkx3p4P7u43eez0f3O/SXjG9BKyorJK12/eM3sYr40SkjhCCUmMoYQkxlBCEmMoIYkxlJDEGEpIYgwlJDHGrPmQd+OhTl5u3H6TjiXdxv0mzS4NKSH9i2kl5OOSZ/oqZYECgiwKCLJIkqxWm/pETtKC7zaZNSy3UUISRyghiTGUkMQYSkhiDCUkMYYSkhhDCUmMoYQkxvjbfHAs6W7MfpNZJ1z3m3yTl3TDM6aVkMnfrFDk6Dm6WHjDWUJK0qY/DqlXxESzhuU2SkjiCCUkMYYSkhhDCUmMoYQkxlBCEmMoIYkx/lY6EXPzJswHx36Txy+UuL3f5JuypBueMa2EbNszXgVXbkqSSwl54dJ1tewcY9Ko3EcJSRyhhCTGUEISYyghiTGUkMQYSkhiDCUkMeZNKJ2I7/Imz4e8Gw+Ve+Wxc7/JVy3pduw3aVzSffOB/5ST8IxpJeR7odG6cfu+JNcSMiv7rD4KizNpVO6jhCSOUEISYyghiTGUkMQYSkhiDCUkMYYSkhjzJpdOpOnTHOfD6euPXPabdGdJt1n7TcIzppWQ0UnpWvHTNkn/lpCVVdUakjhb8RMWmDUst1FCEkcoIYkxlJDEGEpIYgwlJDGGEpIYQwlJjGmOpRN5/fxX5oNjSXdOwRO3l3RfvO79Jd3wjGkl5MXCG/ooLE4JkxYqIMiiSXNWqfOAMXq/W6wuXb1p1rDqsNnsmrNkrT4Ki1OrLjFK/maFKiqrnPdTQhJHKCGJMZSQxBhKSGIMJSQxhhKSGEMJSYz5r5ROpHH5r8+Hs7eKXZZ0H8p53uCS7pN/N/2SbnjGtBJSku7ce6j05es1JHG2LCNnaeaiDOcSbX+xbut+df1irK4U3dHdB480MP4bpS9f77yfEpI4QglJjKGEJMZQQhJjKCGJMZSQxBhKSGLMf710Iq5hPtSfvBsPXZZ0v2q/yQPZ/+43+bpLuuEZU0vIhqQuzDB7CBo8YqbWbNzlvH3waJ6C+41y3qaEJI5QQhJjKCGJMZSQxBhKSGIMJSQxhhKSGEPpRIxhPjQ+xiXd2X8/bXBJd9YJ1yXdV25TQnqLX5eQ3QeOM3sIat9nhA4dy3Pevn7zngKCLCqvqF2STQlJHKGEJMZQQhJjKCGJMZSQxBhKSGIMJSQxhtKJGMN88DxnbxXr9PVHOn6hREfPNn5J95mC2iXd1+9WmFVNNRuUkA14v1ussnMvOG/fe/BYAUEWPXryTFJtCUkIIYQQQgghhBBCmnfgGUrIBrTvM0J/HT/rvF3fmZCEEEIIIYQQQgghpHkHnqGEbMCghFT99Otu5+0/j+QqqO9Il8d48/Lv5M3J84oalTyvNn0cxD9SXmnV49Iq08dB/CMVVVY9elpp+jiIf6SqxqbiEuYDqU2N1a77T5rmip3kzY/NLt19ZP44iH+Ez5rEGOaDfwSeoYRsQMbmPer6RbKu3bir+8VPNDA+VTMXuV4wx+xfAuIfoYQkxlBCEmMoIYkxlJDEGEpIYgwlJDGGz5rEGOaDfwSeoYRsgNVqU+rCDH0YFqdWXWI0ZsYylZW7bkZq9i8B8Y9QQhJjKCGJMZSQxBhKSGIMJSQxhhKSGMNnTWIM88E/As9QQjYBs38JiH+EEpIYQwlJjKGEJMZQQhJjKCGJMZSQxBg+axJjmA/+EXiGErIJmP1LQPwjlJDEGEpIYgwlJDGGEpIYQwlJjKGEJMbwWZMYw3zwj8AzppWQz0rLdKXojmw2u/OYzWbXlaI7Kn1e+wd76myBWcNzi9m/BMQ/QglJjKGEJMZQQhJjKCGJMZSQxBhKSGIMnzWJMcwH/wg8Y1oJOWnOKkWOmlPn+FcpCzR93hoTRvT6zP4lIP4RSkhiDCUkMYYSkhhDCUmMoYQkxlBCEmP4rEmMYT74R+AZ00rIjv0S9eeR3DrHDx3LU0j/0SaM6PWZ/UtA/COUkMQYSkhiDCUkMYYSkhhDCUmMoYQkxvBZkxjDfPCPwDOmlZCBIZHKv1xU53j+5SK17jrUhBG9PrN/CYh/hBKSGEMJSYyhhCTGUEISYyghiTGUkMQYPmsSY5gP/hF4xrQSsk/kJH3787Y6x1dmZOqzoVN9PyAPmP1LQPwjlJDEGEpIYgwlJDGGEpIYQwlJjKGEJMbwWZMYw3zwj8AzppWQuw/m6O3gIZqS/oMy9xxV5p6jmjZ3td4OHqLfdv5l1rAAAAAAAAAANDHTSkhJ2nMoR59GT9FbHSPUokOEekZM0O+7D5s5JAAAAAAAAABNzNQS0sFms8tut5s9DAAAAAAAAABeYFoJ+e3P2zQwPlU227/lo81mV1RSmlat227WsAAAAAAAAAA0MdNKyLBB47VlR1ad45l7jqqHZYIJI2o8m82uOUvW6qOwOLXqEqPkb1aoorLK7GHByx6XPNNXKQv0Xmi02vVO0MqMzAafU11j1dT01fp6/k8+GCF8yd35sGjVZnX9IlnvhkYrNDxJP2/a46ORwhfcnQ+r1+9Uz4gJei80Wm0+Ga7xqSv1vKzCR6OFt73O/y8kqbKqWoNHzNTgETO9PEL40rmLV/Vp9BQFhkSq25fJyso+88rHhw+broAgi0saO4fg/9ydD5K0dddhfRo9Re+GRqtnhH9/ToJ73JkP2/dl13lvCAiyKGXmSh+OGN7k7vtDVvYZ9Y2arPZ9RqjrF2O19rd9Phop8PpMKyEDO0Up/3JRneMXC2+oZecYE0bUeOu27lfXL8bqStEd3X3wSAPjv1H68vVmDwteNmrqUo2YtEjFj0p0Nv+q2vaM16FjeS99/P7DuQrpP1otO8dQQjZD7s6H9OXrdTw3Xw8fP1VW9hm9FxqtY6fO+3DE8CZ358OWHVk6e+GKnpSUqvDaLYUNGs8qgGbE3fkgSVarTQmTFio8dholZDNSXV2jjv0StTIjU09Ly7Rt9xG17hqrR0+evfQ54cOma8O2A3paWuZMVVW1D0cNb3md+bB6w05nGfHw8VPduH3fhyOGN7k7H6prrC7vC09LyzRmxjIt+G6Tj0cOb3B3Ptx78FitusTozPlCSdK1G3fVpsdwnTid78thA24zrYQM6T9amXuO1jmeufeowgaNN2FEjTd4xEyt2bjLefvg0TwF9xtl4ojgbVVV1XWK87Rl6zQu9dsGn5u2dB0lZDPjyXxwGBifqh8N7yN4c3kyH6xWmy5dvanuA8cpO/eCN4cJH3nd+TA5bZXSl6/Xlh1ZlJDNyPHcfH0UFuey/dBnQ6dq8/ZDL31O+LDp2rb7iC+GBx9zdz6UPi9X665Dde3GXV8NET70Ou8PRtdu3NWHYXGvLLHx5nB3PhzPzVfLzjEqr/h3RWa/mKn6bedfXh8r4AnTSshla7aqbc94Ze49qnsPHuveg8favi9b7XonaP7KX80aVqO07zPC5YyG6zfvKSDI4vIGgOblatGdOn/GGzMPKDx2WoPPpYRsfjyZD5JUXlGlNj2G62jO394aInzIk/kQEGRRYEikNv3RuA8c8H+vMx/mr/xVKTNXym63U0I2Mxt+/1P9h7r+2Y+dsVxzl2946XPCh03XR2FxCvlslAYlpOrAkdPeHiZ8xN35sP9wrtr2jFfc+Hl6v1us/q/nV1q0arMvhgofeJ33B6PEKYs5C7IZcXc+VFfX6MuvvlHfqMnafzhXB4/mqVfERLb3gd8zrYS0Wm1KXZiht4OHOPezaNEhQhNnf6/q6hqzhtUo73eLdTlj5d6DxwoIsvCvUM3YhUvXFRBkcbmKe+aeo+oxOKXB51JCNj+ezAdJmjDrO0Unpbs8H28uT+aDzWbXpas31Sk8SbsP5nhzmPARd+fD1l2HFTd+nmqsVkmihGxmfty4S4MSUl2OTU5bpRnz1rz0ORcuXdeN2/dVdOu+Mjbv0buh0bpw6bq3hwofcHc+ZGzeoy6fj9W5i1dVWVWt8wXX1L7PCK3but8Xw4WXvc77g8PfF6/po7A4PS0t89bw4GOvMx9+2bJPgxJSFTMmXS06RLB/MN4IppWQDiVPn+vE6Xzl5F18Y0q89n1G6K/jZ523OROy+XOc2WLck2lj5oE6/1pVH0rI5seT+TB7yVr1jZqsZ/ylsdnwZD44/G/FBo2autQbw4OPuTsf5i7foMCQSAV2ilJgpyjnP84GdopSWXmlr4YNL9nw+58aEDfD5djYGcvd2ks8KilNy9ZsbeqhwQTuzoeMzXvrlBLzvt2o4SnzvTZG+I4n7w9RSWla8dM2bw0NJnB3Pvx1/Ky6fpHsPIGr8PptdeyXqIzNe70+VsATppeQb6JBCan66dfdztt/HslVUN+RJo4I3lZVVa3AkEiXPb7mLFmrsTOWN/hcSsjm53XmQ43VqgmzvtPA+G8oIJsZT94fHKbPW8PVLZsJT+cDZ0I2L9m5F9Tmk+EuZ8Z+Gj1FGzMPNPo1Po+bwRVPmwl350NW9hm1652g6hqr81jasnWaMOs7bw8VPvC67w/Hc/PVtme8yspZdtucuDsflq/5XcPGzXM5NjltlSbO/t6bwwQ8ZloJ2X3guFfGn2Vs3qOuXyTr2o27ul/8RAPjUzVzUYbZw4KXjZi0SIlTFuvh46e6cOm6Pu4Vr31Zp5z3z1r8S72nwFNCNk/uzIcaq1VRSWkanjJfj0ueqay8QmXlFZw93Yy8aj7YbHZFjpqjrOwzkmq3I4mfuFDZuRdU/KhEuw/mqFWXGOf9ePO5Mx9eRAnZvFRVVat9nxH6fu12PSst0x/7jqll5xg9ePhEUt35cOtusdKWrdP5gmsqflSinzftUasuMbr74JGZ3waaiLvzobq6Rp3CkzR7yVqVV1TpwqXratc7gT2lm4mG5sOFS9cVOWqO7hc/cXle+LDpWr1hpxlDhhe5Ox9OnilQy84xOnmmQJJ0+26xOoUnaeefx037HoDGMK2EXLNxV51898sfeic40u83XHbsZ/lhWJxadYnRmBnL+Jeo/4DiRyUaNm6e3g2NVtue8XWWRn351TdK/nqF8/b2fdlq02O43guN1ruh0WrTY7h2HTjh62HDS9yZD09Ly5x73xrj7//ggsZ71XyosVoVEGTRr5kHJUl2u10TZ3+v0PAkvRsarU8Gjdfvuw+bNXR4gTvz4UWUkM1P3vlC9Y2arHeCI9X1i7HafzjXed+L8+FJSaniJyxQu94Jatk5RuGx05STd9GsocML3JkPUu22T0MSZ6t116Hq9mUyV05vZl41H47m/K2AIIuKbt1zHtuXdUpBfUeq0rDlB5oPd+fD77sPq8fgFLXrnaDuA8dpPfvF4g3gd8uxUxf+rOmN2IwXAAAAAAAAwJvB70rIrOwz+jAszuxhAAAAAAAAAGgifldCLl39mz7uFW/2MAAAAAAAAAA0EdNKyMEjZrpkUEKqwgaNV0CQRUt+2GLWsAAAAAAAAAA0MdNKyLSl6+pk+ZrfdezkebOGBAAAAAAAAMAL/G45NgAAAAAAAIDmxbQS8klJqQqv33bevn7znn7Zsk9nzheaNSQAAAAAAAAAXmBaCTlq6hKNnbFckvTg4RN90H2Y2vaMV2BIpHYfzDFrWAAAAAAAAACamGklZPs+I3Qk55wk6edNe9QpPEnVNVZl7jmqvlGTzRoWAAAAAAAAgCZmWgn5bmi0bt0tliQlf71Ck9NWSZJu3L6vlp1jzBoWAAAAAAAAgCZmWgn5yaDx2nXghMrKK/Vxr3ht+uOQJCn33CW17zPCrGEBAAAAAAAAaGKmlZC/7fxLgSGRatNjuNr1TtDT0jJJ0uS0VYpNnmvWsAAAAAAAAAA0MdNKSEnKybuojM17dfufZdk1VquWr/ldeVwhGwAAAAAAAGg2TC0hGzIgbobZQwAAAAAAAADgIb8uIbsPHGf2EAAAAAAAAAB4iBISAAAAAAAAgFdRQgIAAAAAAADwKkpIAAAAAAAAAF5FCQkAAPAGixmTrpmLMsweBjywbfcRBQRZnMm/XCRJOnW2QAFBFpWVVzT6tfZmndTBo3leGWcPywTnGGPGpHvlawAAgOaLEhIAALglpP9oZxHRokOE3u8Wq+Ep83W16I7zMX0iJ7mUKucuXnXeV15RpaWrf9Mng8YrsFOUPgyLU2zy/3TyTIHL19myI0sBQRbFJv+vzhh+3LhLAUEWxU9YUOe+A0dOKyDIos4Dxshut9e5v+DKTY2aukRBfUfq3dBohfQfrdjkudq2+4jzMTVWq9Zs3KVeERPVsnOM2nwyXL2HTNTsJWt1+25xo39Wju+hvuw6cKLRr/MqxhKyrLxCAUEW/Zp5sEleu7kqK6/UyoxM9Y2arFZdYtS661D1sEzQlPQfdOZ8YZ3Hpy1dp4AgiybO/r7OfbfuFisgyKK2PeNVXlHlcl/+5SIFBFn0QfdhrxzPtt1H1MMyoc5xd0vIs/lX9UH3YXq/W6wuXb3ZqOcY2Wx2nS+4prY94185P1dmZFJCAgAAt1FCAgAAt4T0H61Zi3/Rpas3lX+5SPsP56rrF2M1JHG28zEVlVVatGqzQj4bpbLyCtlstWVgVVW1BsTN0Afdh+mHdTt04nS+9hzK0fCU+WrRIULb92U7X2PLjix90H2YAjtF6XzBNefx6hqrOvZLVLveCfWWkIlTFity9BwFBFnqFJtXi+7ovdBoDRs3T4eO5enshSvadeCEJs7+Xr0iJjofN3H29/qg+zBlbN6j039f1onT+Vq9fqd6D5noUlY2ZMuOLAWGRCrvfGGdlDx93ujXeRVjCVljtWrb7iMqunW/SV67OXpSUqqeERPUKTxJ67fu17mLV3XqbIF+zTyoIYmz6/z902q1qX2fEYpKSlPrrkPrFI2OEvLjXvFavWGny32jpi5Vu94JPikhS5+Xq/OAMdqw7YC+++UP9RicoorKqgaf5/Dw8VO9FxqtwE5RDZbklJAAAOB1UEICAAC3hPQfrZUZmS7H5i7foE7hSS7Hlq/5XSH9R7scW71+p97qGKGLhTfqvO6YGcvU5pPhKiuvlFRb4IX0H61pc1crYdJC5+M2bz+k/kOnaeTkxXVKyJJnzxXYKUo5eRc1KCFVk+ascrl/ZUamWncdKqvVVufrPykpdf53YEikvvvljzqPsdvtKnnW+PJwy44sBXaKavBxe7NOqm/UZAV2ilKn8CRNTV+t0uflkqRTZy9pQNwMBXaKUrveCUpbtk7VNVbnc19cjh3YKUr7sk41eoxS7dmhMWPS1apLjNr2jFf8hAW6fPWWJKm6ukZzlqxVu94JCuwUpc/jZijPcLagoyg7npuvwSNm6p3gSC1atdn5utFJ6WrZOUbteidofOpKl59zQ8Jjp+l/KzY4b2ds3quPwuJUY639/svKK/VWxwidOJ3f6NdM/nqF2vaM16Mnz+q9/+yFKy63Dx3LU+uuQ1X8qEStusTUKaEdJeTKjEy17zNClVXVkqTrN++pdddYLV61xeMS8mjO3/pi+Nd6NzRavYdMdPn5OyRNX6YRkxZJqp2nkaPnaEr6D6/8ui/zVscISkgAANDk/LqErK6uMXsIAADgBS/K/tbhAAAMpElEQVSWkFeL7ihs0Hj99Otul8fVV0I6ysP6FF6/rYAgi/YfzpX0bwl5884DvRMcqUtXb8pms6v7wHHa/9epekvIdVv3q1N4kux2uzb9cUjvd4t1ORts2+4jatEhot4lt0ZdPh+rhEkLPf67iKOErLFa68SxVHxf1im9ExypFT9t04VL17U366QiEmfpStEd3XvwWC07xyh1YYbzrNOO/RKVtnSd82t4WkLevlusD8PilPz1CuWdL9Sxk+c1OW2VNmYekCTNXJShkM9G6c8jucq/XKSv5/+k1l2H6sHDJ5L+Lco6fDpSqzfs1L6sUzp74Yru3H+kj8LitHzN77px+74uFt5Q5Og59Z69+jJzl2/QwPhvnLejktLUokOEjub8LUnKybuod0OjVfVP8deQsvJKBYZEavX6nQ0/+B+jpi5VysyVkqTkb1YocvQcl/sdJeSN2/fV9Ytk/bJlnyRpctoqzVr8izI27/W4hAwNT9KGbQf055FcDR4xU6H/zHGHTX8cUsd+iS5n1xY/KlHbnvHa+efxRn+vDpSQAADAG3xeQsZPWNCoAAAA/xTSf7RadIjQWx0j1KJDhAKCLOoUnqRDx1wvhlFfCflhWJwWr9pS7+tarTa9ExzpXNLqKCElaVzqt0qavky7DpxQr4iJstvt9ZaQ4bHTnK9f+rxc74VGK3PvUef9NVarEqcsVkCQRT0jJihl5kqt3rDTedafw7GT59Wmx3B90H2YYsaka+7yDdp14IRby1sd30NDe0J+Gj2l3jPWbDa70pevV5/ISS7Hdx/MUWBIpHOJrqcl5Jwla9V94Djnknnj139WWqZ3giO1N+uk87jdbleviImav/JXSf8WZXcfPHJ5fvry9YqfuNDl2MXCGwoIsjT655iVfUbvhkarusaqZ6Vlat01VpPTVmnGvDWSpFXrtssyclajv1fH18/OvdCoxz8tLdO7odHOxx/JOacWHSJcvldHCXnrbrE2bz+kkM9G6dbdYrXuGqv7xU+apIQ0nj2ak3dRAUEWZwl8peiOWneNVU7exTrPP3QsTx+GxenmnQeSav9My8orXFIfSkgAAOANPi8hA4Isatc7QSMnL35lAACAfzLuCVlw5aZyz13SyoxMvRMc6VL4vW4J+ePGXZJcS8jC67f1dvAQdfl8rHPfyBdLyCtFdxQQZHHZDzFp+rJ6y5IrRXe09rd9mjFvjT4bOlUBQRZN/6fYciivqNL+w7lauvo3JUxaqDafDFf7PiOcVy5uDMeekOcuXq2TkmfPZbfb9VbHiJfuMxmbPLdOQXnvwWMFBFl04dJ1SZ6XkBGJs+osW3c4e+GKAoIsul/8xOX45LRV+iql9mf/sn0LhyTOfmkBe/3mvUaNray8Qm8HD9GZ84Xavi9bw8bN08kzBQrqO9JZRC9d/Vujv1dHCXnq7CXnsXMXr7qMrcvnY533rd+6XyGfjXKedWiz2dXh05EuZwIbS8jqGqtCPhulbl8ma9rc1ZLUJCWk8Wd77cZdBQRZdMVwIajGevF7DQiy6GlpWZ3HUUICAABv8HkJuWbjLnXsl6j+Q6dp/1+n6r1qJQAA8F/17QkpScPGzVNUUprzdn0l5GdDpza4HPvAkdOSXEtIqbZ07PZlsvOMvRdLyHnfbqy38GrRIaJOifYixxmLhddvv/QxZeWV6jE4RaOmLn3la734uq/aE9Jms6tFhwhl7jla7/31lZB3Hzxq0hLSMnKWJqe5V0JOmtNwCfmq13VHeOw0/bhxl5KmL9OmPw7JZrPr417xOnO+UB37Jbq1H6Sj1DReQKaiskr5l4ucS82NJeSAuBn1zqmwQeOdjzGWkFJt6fh28BDn2YdNXUIW3br32iVkWXllnQsk1bc/KiUkAADwBlP2hKyusWrz9kPqPnCcekVM1PZ92XWWAAEAAP/0shIyctQcl2KivhJy1brteqtjRL1nEyZNX6Y2PYY7rz78YglZ/KjE5SxHYwlps9nVsV+iFn6/yVkoOfJ/Pb/S92u3S6q9UEp9pcvJM7Vlz/mCa6qsqn5pwTMkcbaGp8x/6c/mRY25ME0PywTnWXNGNptdacvWqfeQiS7Hdx04oXeCm2459tfzf9Ing8bXuxz7aWmZ3g4eoj2HcpzHX7Yc+8UScuaiDHUeMMZ5oRaH+n7+rzJ3+QZ9lbJAbT4Z7ryYzOS0VRqX+q1b+0E6JE5ZrKC+I53LmY2Wr/ndWUI6zjjc/9cpl/m0L+uUAoIszgvYvFhCVlVV6++L15yv6U8lZGNRQgIAAG8w9cI0Nptduw6cUJfPx6rbl8kuHywAAIB/enE59skzBZoxb41adIhw7gv5vKxCC77bpOB+o/S0tMxZPFVVVWtA3Ax90H2Yfli3Qzl5F7U366Tixs/TWx0jXC6i8WIJ+SJjCXkk55ze6hih4kcldR43Y94a9RicIqm2cOs+cJxWZmQqO/eC8s4XasuOLHUeMEbRSemy2+16XPJMgSGRSv56hXbsz9bZC1d0JOecpqav1tvBQ9w6886xHPvU2YI6cezzt3XXYb0dPEQrMzKVf7lIh0+c01cpC3Sl6I7uPnik90KjlbrwZ+VfLtKfR2ovTDN7yVrn1/C0hLxSdEfvhUZrfOpKnfnn7LhZi3/R5u2HJEmpC3+uc2GaVl1inGdHvqyEvHH7vt7vFqvI0XN0NOdv5V8u0pYdWeobNbnRY5Nq94Vs0SFCgxJSnccOHctTiw4Rbu0H6fDg4RN1/WKs2vcZodXrd+rU2UvKO1+oHfuzNSgh1VlCzl/5a73FoCT1GJzi3JfyxRLyRW9KCWm315bOT0vL9FbHCG3ddVhPS8vqvTgTJSQAAHgdppaQl67eVNL0ZXo7eIhGTl5c7wcHAADgX0L6j3ZZmtqqS4wGxn+jrOwzzsf0iZzk8phzF6867yuvqNLiVVvU7ctkBYZEOi/+cvJMgcvXcaeETP56haKT6i9Fcs9dco7hStEdzVmyVr0iJuqjsDi9Exyp0PAkLf1xq8rKKyXVnqm3ZUeWho6dq5DPRikwJFLvd4tVzJh05Z67VO/XeJnGXJhGqr3YTN+oyWrVJUYd+yVqavpqlT4vl1R7IZLwYdMVGBKpj3vFa86StS7FkKclpCTlXy5SVFKaPgqLU5tPhmt4ynwVXqu9WE9VVbVmLf5FH/eKV2BIpAbEzXDZU/FlJaRU+3e9r1IW6MOwOL0XGu1yBmVjOZZQr/lnr1BJqq6u0fvdYt3aD9Ko9Hm5lvywRX0iJ6ll5xgFdopSj8EpSl++XoXXbslmsyu436h6z/iVpG9/3qY2nwxXdXVNsykhix+V1DtP12/dX+exlJAAAOB1mFJCni+4poRJC/VOcKTGpX77yv2XAAAAgObsZSWkv6KEBAAAr8PnJWRs8lwFdorStLmrnRt2AwAAvElWb9ip7gPHvTIvu9iMr5w6W9DgGFMXZjT8Qm/A2Kqraxp8vQFxM7z4HXmGEhIAAPwX+LyEDAiy6IPuwxTSf/QrAwAAAPwXbNt9xGUJdH0XbvIHPSwTnGOkhAQAAO7yeQm5ZUdWowIAAAAAAACgeTD1wjQAAAAAAAAAmj/TSkibza7juflat3W/1m3d7/bVJgEAAAAAAAC8GUwpIfMvF6n7wHEKCLLorY4RatEhQgFBFvWLmaqrRXfMGBIAAAAAAAAAL/F5CXm/+Ik+7hWvwSNm6mz+VVmtNtlsdp0vuKYhibPV4dORevj4qa+HBQAAAAAAAMBLfF5Czl6yVt0HjlNlVXWd+yqrqtXty2QtXf2br4cFAAAAAAAAwEt8XkL2GJyiNRt3vfT+H9btUA/LBB+OCAAAAAAAAIA3+byEbNUlRqf/vvzS+w+fOKdWXWJ8OCIAAAAAAAAA3uTzEjIgyKKCKzdfev/5gmsKCLL4bkAAAAAAAAAAvMqUEjJz71HlnS+sN7/vPkwJCQAAAAAAADQjppSQjQkAAAAAAACA5sHnJeStu8WNCgAAAAAAAIDmweclJAAAAAAAAID/FkpIAAAAAAAAAF5FCQkAAAAAAADAqyghAQAAAAAAAHgVJSQAAAAAAAAAr6KEBAAAAAAAAOBVlJAAAAAAAAAAvIoSEgAAAAAAAIBXUUICAAAAAAAA8CpKSAAAAAAAAABeRQkJAAAAAAAAwKsoIQEAAAAAAAB4FSUkAAAAAAAAAK+ihAQAAAAAAADgVZSQAAAAAAAAALyKEhIAAAAAAACAV1FCAgAAAAAAAPAqSkgAAAAAAAAAXkUJCQAAAAAAAMCrKCEBAAAAAAAAeBUlJAAAAAAAAACvooQEAAAAAAAA4FWUkAAAAAAAAAC86v8BbPiR3oH4ZdIAAAAASUVORK5CYII=", "text/html": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "production_envelope = phenotypic_phase_plane(model, \n", " variables=[model.reactions.BIOMASS_Ecoli_core_w_GAM],\n", " objective=model.metabolites.succ_e)\n", "production_envelope.plot(plotter, height=400)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Set up a model that represents a reference state (in this case a model with a constrained growth rate)." ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "model.reactions.EX_o2_e.lower_bound = 0\n", "reference_model = model.copy()\n", "biomass_rxn = reference_model.reactions.BIOMASS_Ecoli_core_w_GAM\n", "biomass_rxn.lower_bound = 0.\n", "target = reference_model.metabolites.succ_e" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Set up the differential flux variability analysis strain design method." ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "diffFVA = DifferentialFVA(design_space_model=model,\n", " reference_model=reference_model,\n", " objective=target,\n", " variables=[biomass_rxn],\n", " normalize_ranges_by=biomass_rxn,\n", " points=10)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Run differential flux variability analysis (only on the surface of the production envelope)" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "58822f28111843bc9f281496f9b3acba", "version_major": 2, "version_minor": 0 }, "text/plain": [ "HBox()" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "result = diffFVA.run(surface_only=True)" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
lower_boundupper_boundgapsnormalized_gapsbiomassproductionKOflux_reversalsuddenly_essentialfree_fluxreactionexcluded
reaction
ACALD0.0000000.000000-8.279455NaN0.00000013.905778TrueFalseFalseFalseACALDFalse
ACALDt0.0000000.0000000.000000NaN0.00000013.905778TrueFalseFalseFalseACALDtFalse
ACKr-5.664889-5.664889-2.838696NaN0.00000013.905778FalseFalseFalseFalseACKrFalse
ACONTa0.4293330.4293330.200970NaN0.00000013.905778FalseFalseFalseFalseACONTaFalse
ACONTb0.4293330.4293330.200970NaN0.00000013.905778FalseFalseFalseFalseACONTbFalse
.......................................
TALA-0.033659-0.033659-0.0042070.0000000.1881452.123333FalseFalseFalseFalseTALAFalse
THD23.2259503.225950-0.4032440.0000000.1881452.123333FalseFalseFalseFalseTHD2False
TKT1-0.033659-0.033659-0.0042070.0000000.1881452.123333FalseFalseFalseFalseTKT1False
TKT2-0.101579-0.101579-0.0126970.0000000.1881452.123333FalseFalseFalseFalseTKT2False
TPI9.8128529.8128520.0233935.9056160.1881452.123333FalseFalseFalseFalseTPIFalse
\n", "

684 rows × 12 columns

\n", "
" ], "text/plain": [ " lower_bound upper_bound gaps normalized_gaps biomass \\\n", "reaction \n", "ACALD 0.000000 0.000000 -8.279455 NaN 0.000000 \n", "ACALDt 0.000000 0.000000 0.000000 NaN 0.000000 \n", "ACKr -5.664889 -5.664889 -2.838696 NaN 0.000000 \n", "ACONTa 0.429333 0.429333 0.200970 NaN 0.000000 \n", "ACONTb 0.429333 0.429333 0.200970 NaN 0.000000 \n", "... ... ... ... ... ... \n", "TALA -0.033659 -0.033659 -0.004207 0.000000 0.188145 \n", "THD2 3.225950 3.225950 -0.403244 0.000000 0.188145 \n", "TKT1 -0.033659 -0.033659 -0.004207 0.000000 0.188145 \n", "TKT2 -0.101579 -0.101579 -0.012697 0.000000 0.188145 \n", "TPI 9.812852 9.812852 0.023393 5.905616 0.188145 \n", "\n", " production KO flux_reversal suddenly_essential free_flux \\\n", "reaction \n", "ACALD 13.905778 True False False False \n", "ACALDt 13.905778 True False False False \n", "ACKr 13.905778 False False False False \n", "ACONTa 13.905778 False False False False \n", "ACONTb 13.905778 False False False False \n", "... ... ... ... ... ... \n", "TALA 2.123333 False False False False \n", "THD2 2.123333 False False False False \n", "TKT1 2.123333 False False False False \n", "TKT2 2.123333 False False False False \n", "TPI 2.123333 False False False False \n", "\n", " reaction excluded \n", "reaction \n", "ACALD ACALD False \n", "ACALDt ACALDt False \n", "ACKr ACKr False \n", "ACONTa ACONTa False \n", "ACONTb ACONTb False \n", "... ... ... \n", "TALA TALA False \n", "THD2 THD2 False \n", "TKT1 TKT1 False \n", "TKT2 TKT2 False \n", "TPI TPI False \n", "\n", "[684 rows x 12 columns]" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "result.solutions" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "application/vnd.plotly.v1+json": { "config": { "linkText": "Export to plot.ly", "plotlyServerURL": "https://plot.ly", "showLink": false }, "data": [ { "fillcolor": "#B3E2CD", "hoverinfo": "none", "mode": "markers", "name": "WT", "opacity": 0.3, "type": "scatter", "x": [ 0, 0, 0, -19.437336380718584 ], "y": [ 1, 2, 3, 4 ] }, { "fillcolor": "#FDCDAC", "hoverinfo": "none", "mode": "markers", "name": "Strain 5", "opacity": 0.3, "type": "scatter", "x": [ 0, 0, 0, -19.68740910039923 ], "y": [ 1.5, 2.5, 3.5, 4.5 ] } ], "layout": { "height": 432, "shapes": [ { "fillcolor": "#B3E2CD", "line": { "color": "#B3E2CD", "width": 0 }, "opacity": 0.3, "type": "rect", "x0": 0, "x1": 0, "y0": 0.75, "y1": 1.25 }, { "fillcolor": "#B3E2CD", "line": { "color": "#B3E2CD", "width": 0 }, "opacity": 0.3, "type": "rect", "x0": 0, "x1": 0, "y0": 1.75, "y1": 2.25 }, { "fillcolor": "#B3E2CD", "line": { "color": "#B3E2CD", "width": 0 }, "opacity": 0.3, "type": "rect", "x0": 0, "x1": 0, "y0": 2.75, "y1": 3.25 }, { "fillcolor": "#B3E2CD", "line": { "color": "#B3E2CD", "width": 0 }, "opacity": 0.3, "type": "rect", "x0": -19.43733638071854, "x1": -19.437336380718623, "y0": 3.75, "y1": 4.25 }, { "fillcolor": "#FDCDAC", "line": { "color": "#FDCDAC", "width": 0 }, "opacity": 0.3, "type": "rect", "x0": 0, "x1": 0, "y0": 1.25, "y1": 1.75 }, { "fillcolor": "#FDCDAC", "line": { "color": "#FDCDAC", "width": 0 }, "opacity": 0.3, "type": "rect", "x0": 0, "x1": 0, "y0": 2.25, "y1": 2.75 }, { "fillcolor": "#FDCDAC", "line": { "color": "#FDCDAC", "width": 0 }, "opacity": 0.3, "type": "rect", "x0": 0, "x1": 0, "y0": 3.25, "y1": 3.75 }, { "fillcolor": "#FDCDAC", "line": { "color": "#FDCDAC", "width": 0 }, "opacity": 0.3, "type": "rect", "x0": -19.68740910039923, "x1": -19.687409100399233, "y0": 4.25, "y1": 4.75 } ], "template": { "data": { "bar": [ { "error_x": { "color": "#2a3f5f" }, "error_y": { "color": "#2a3f5f" }, "marker": { "line": { "color": "#E5ECF6", "width": 0.5 } }, "type": "bar" } ], "barpolar": [ { "marker": { "line": { "color": "#E5ECF6", "width": 0.5 } }, "type": "barpolar" } ], "carpet": [ { "aaxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "baxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "type": "carpet" } ], "choropleth": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "choropleth" } ], "contour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "contour" } ], "contourcarpet": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "contourcarpet" } ], "heatmap": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmap" } ], "heatmapgl": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmapgl" } ], "histogram": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "histogram" } ], "histogram2d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2d" } ], "histogram2dcontour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2dcontour" } ], "mesh3d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "mesh3d" } ], "parcoords": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "parcoords" } ], "pie": [ { "automargin": true, "type": "pie" } ], "scatter": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatter" } ], "scatter3d": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatter3d" } ], "scattercarpet": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattercarpet" } ], "scattergeo": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergeo" } ], "scattergl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergl" } ], "scattermapbox": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattermapbox" } ], "scatterpolar": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolar" } ], "scatterpolargl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolargl" } ], "scatterternary": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterternary" } ], "surface": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "surface" } ], "table": [ { "cells": { "fill": { "color": "#EBF0F8" }, "line": { "color": "white" } }, "header": { "fill": { "color": "#C8D4E3" }, "line": { "color": "white" } }, "type": "table" } ] }, "layout": { "annotationdefaults": { "arrowcolor": "#2a3f5f", "arrowhead": 0, "arrowwidth": 1 }, "autotypenumbers": "strict", "coloraxis": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "colorscale": { "diverging": [ [ 0, "#8e0152" ], [ 0.1, "#c51b7d" ], [ 0.2, "#de77ae" ], [ 0.3, "#f1b6da" ], [ 0.4, "#fde0ef" ], [ 0.5, "#f7f7f7" ], [ 0.6, "#e6f5d0" ], [ 0.7, "#b8e186" ], [ 0.8, "#7fbc41" ], [ 0.9, "#4d9221" ], [ 1, "#276419" ] ], "sequential": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "sequentialminus": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ] }, "colorway": [ "#636efa", "#EF553B", "#00cc96", "#ab63fa", "#FFA15A", "#19d3f3", "#FF6692", "#B6E880", "#FF97FF", "#FECB52" ], "font": { "color": "#2a3f5f" }, "geo": { "bgcolor": "white", "lakecolor": "white", "landcolor": "#E5ECF6", "showlakes": true, "showland": true, "subunitcolor": "white" }, "hoverlabel": { "align": "left" }, "hovermode": "closest", "mapbox": { "style": "light" }, "paper_bgcolor": "white", "plot_bgcolor": "#E5ECF6", "polar": { "angularaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "radialaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "scene": { "xaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "yaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "zaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" } }, "shapedefaults": { "line": { "color": "#2a3f5f" } }, "ternary": { "aaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "baxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "caxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "title": { "x": 0.05 }, "xaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 }, "yaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 } } }, "title": { "text": "Compare WT solution 5" }, "width": 700, "xaxis": { "autorange": true, "range": [ -20.97077696131291, 1.283367860913681 ], "title": { "text": "Reactions" }, "type": "linear" }, "yaxis": { "autorange": true, "range": [ 0.7390560875512996, 4.760943912448701 ], "ticktext": [ "", "FBP", "G6PDH2r", "PGL", "PGK" ], "tickvals": [ 0, 1, 2, 3, 4 ], "title": { "text": "Flux limits" }, "type": "linear" } } }, "image/png": "iVBORw0KGgoAAAANSUhEUgAABSEAAAGwCAYAAABW/qB/AAAgAElEQVR4nOzd93eVZb7+8e/fsx3FwjCDc3QjGEHaiDIyKroxJoQmVaoFFFGRZgFFFBkQdVBBRyyIIE2KRFqQXgSBJCQkhIRUkn19f8BsCWm3IZs7z+d5v9d6rXUISXzwnHvJfZ3s5P+JiIiIiIiIiIiIKIn9P98PQERERERERERERLZjhCQiIiIiIiIiIqKkxghJRERERERERERESY0RkoiIiIiIiIiIiJIaIyQRERERERERERElNUZIIiIiIiIiIiIiSmqMkERERERERERERJTUGCGJiIiIiIiIiIgoqTFCEhERERERERERUVJjhCQiIiIiIiIiIqKkxghJRERERERERERESY0RkoiIiIiIiIiIiJIaIyQRERERERERERElNUZIIiIiIiIiIiIiSmqMkERERERERERERJTUGCGJiIiIiIiIiIgoqTFCEhERERERERERUVJjhCQiIiIiIiIiIqKkxghJRERERERERERESY0RkoiIiIiIiIiIiJIaIyQRERERERERERElNUZIIiIiIiIiIiIiSmqMkERERERERERERJTUGCGJiIiIiIiIiIgoqTFCEhERERERERERUVJjhCQiIiIiIiIiIqKkxghJRERERERERERESY0RkoiIiIiIiIiIiJIaIyQRERERERERERElNUZIIiIiIiIiIiIiSmqMkERERERERERERJTUGCGJiIiIiIiIiIgoqTFCEhERERERERERUVJjhCQiIiIiIiIiIqKkxghJRERERERERERESY0Rksih9z/+RpFoTDl5Bb4fhYiIiIiIiIgocDFCNtDh46f09Cvv6e5+Y9Tu7id0x31P6qHBL2jeos+Vd+6878dr8z0zfaEi0ZjOF5XUeXtFZZXadUlVl3+NqvcxB46cVCQa04cr1uix4S8pEo01K/2pGdfrj3TNI+SCpSu1ZuOOem9/dd7HikRj2r3v6LU+Yos6fPxUk/+Ot2T+4uW5iIiIiIiIiMhWjJBXtex/P+jGzo+rfbcMTZi2QPP+84Wmz/1IDw6cokg0pr/3GOz7Edt8X363WZFoTGs31R3dtu86kBi3snPP1fm9D1esUSQa0+Hjp7Tsy3WaNX9ZwrCnX1ckGtPwZ96o8/YV32y8bn+max0hO3QfpIkvLaj39uVfb9Cg8bN1/GT2tT5ii6odIVNHTa/z77bWydO5Xp6LiIiIiIiIiGzFCHlFP+85pBs6DdADac8pv6Co3u/v3HtYDw6c4uHJglVuXqEi0Zimz/2oztvfXvw/dX1orP5y1wCtXL2lzu+NnjJPt/ceqng8Xu/zrVq3XZForMGvJLxeJWuE9F3tCLn4k1W+H4WIiIiIiIiIDMcIeUWPDX9JN3VO1ansvEbfp+pSdeJ/jsfjWrp8tXo+OkHtuqTqzj7D9fzsxSoqvljnY9KfmqHu/cdr/+ETyhg3U7fek6477ntS8/7zheLxuHb9ckSpo6br1nvS1eXB0Xrngy/rjHG79x1VJBrTF6t+1PwlXyql3xi165KqB9KeU+bug3X+WQeP/qbHR76iu/qOVLsuqbqta7r6D5um9Vt213m/2s/59ZptWv71BvV8dIJu6pyqmW8vS7zPqnU/qV/GFN2Skqa/3puhQeNn68jx007/Lrs9PLbeYJs2ZoZemL1YfdMna8rMRXV+766+IzX8mTca/FwtHSELi4r14mtL1PWhsbolJU1dHxqrUZPnKuvA8Trvtzlzrx4Z+qJuvSddHboP0pCJc3T017p/zoZGyH6Dnle/Qc/X++dOmblI7e5+IvHrv9w1oN7LnLv3H9/o53V9pgVLVyoSjenQsVOaOmeJ7rjvSXXoPkiDxs9W9lm3sZQRkoiIiIiIiIiuR4yQv1dUfFE3dBqgwRNmO3/MK3M/UiQa08CxM/XeR1/r+dmL1a5LqvqkPqPyisrE+6U/NUO33pOudl1S9eiT0/T87MXqFZukSDSmIRPn6C93DdBjw1/S87MXq+ejExSJxvTd+szEx9cOhrd1Tdf//XOYJkxboFGT5+qWlDTdek+6Tpz64yWzG7dlqf+wFzXt9Q/0zgdfau77K/TQ4Bf0l7sGaEfWoXqf8+5+Y3RLSpoGjZ+tMc+/lRijasexYU+/ro8+X6MFS1eqc99R6thziM7mN/99MZ99daFu7Py4SsvKJUnV1TX6670Z+mrNVr385ofq9djExPueys5TJBrTB5+tbvBztWSEvFRdrb7pk/X3HoM1Z8Gn+nDFGk17fan+0XuYXp33ceL91mzcoRs6DVDvAZP01uIvNGfBp7qzz3B16D5IR0+cSbzftYyQX6z6UbekpOmRoS9qxTcbteKbjfp+48+Nfl7XZ6odIf/WY7BS+o3Rs68u1MCxMxWJxvTQ4Bec/j3VjpBvLf5CZ3LyVXKxzOnjiIiIiIiIiIj+TIyQv7frlyOKRGOat+hzp/c/djJbkWis3kts12zcoUg0poUff514W/pTM9S+20D9sHln4m0XS8t1W9d0/b3H4DpfpVhcUqpbUtI08rk3E2+rHQzHT3snMepJl79aLhKN6fnZi5t81srKKt3ZZ7jGT3un3ufsP2yazuTk13n/7LMFurHz43W+KlKSzuTkq93dT+i1dz9r8p8nSStXb1EkGtPmzL2SpH2Hfk2Mbd9v/LnOD67533c/KhKN6cCRkw1+rpaMkFkHjisSjdX7vpGlZRXaf/iEpMvD6F19R6rnoxPqjMY5eQW6rWu6MsbNSrztWkZIqfGXY1/9ef/MM9WOkPMWfV7nK3Snvb5UkWisya/ora2hH0zztx6DNXXOEpWWVTT78URERERERERELjFC/t6mn7L+1MtS3/voa0WiMR08+lu930vpN0YPD5ma+HXty7Gv7sGBU3Tf40/Xe/v9TzxbZ9y68qXTV3fvI+PUKzapztt+2nlAo6fMU0q/Mbr1nvTEuPTEmFedPufiT1YlfnhMRWVVHX3TJ9f5PI11Nv+8ItGY3li4XJK05NPv1OXB0ZKkgvPFdX5wzXMz3tffewxu8PtBSi0bIX85+GtiuC0rr2zwfbL2H2v0KzAnvfyubuz8eOJjr9cI+WeeqXaEvPr7l9b++/pp54EG/9xXdjb/vOYt+lyff7tJq9b9VOenk7v875mIiIiIiIiIyCVGyN/buffwn/pKyCkzFykSjTU4cD0x5lXd2Wd44teNjZCPj3yl3oAoSY8+OU19Up9J/LqpwXDg2Jnq2HNI4tdLl3+vSDSmLg+O1tz3V2jl91u0YetudX1orFJHTXf6nM/PXlzvq+Ou5PrDebo9PE6PDX9JkjTyuTc17sX5id/r+eiExMuiez46QcMmvdbo52np94QcOvE1RaIx3XpPuh59cppmzV9W56stv16zTZFoTBu27q73se988KUi0Vjip1ZfrxHyzzxTYyNk7VfIbtyW1ei/m+YaNum1Jr86lYiIiIiIiIjoz8QI+XuFRZe/Om/QeLfvCdnUCJk6errTCJk6anqDI+Rjw19yHiFTR03XHfc9Keny90H8W4/BenDglDov25akPqnPOI+QtX+27bsOaPe+o/UcPn6qoX8l9Xpuxvu6OSVNVVWXdGef4Vr2vx8Svzd5xiL1y5ii80UlzX4FaktHyJqauDZs3a0ps/6j+594Vjd0GqAbOg3QF6t+lPTH4Ldx2556Hzt/yeXB79ffciQ1MUJm1B9kW2OEdHmmxkbILZm/XPMIufzrDYpEY1q1bnuLPwcRERERERERUW2MkFfUb9DzurHz44mRp6Gqqi5Jkt798KsGv1IsHo+ry4OjnV6Ofa0jZE1NXHf1HalHn5wmScrNK1QkGtObC1fU+5x/ZoSs/bP9cvDXhv4VOPfVmq2J78sYicbq/GTtlau36MbOj2vl95e/d2Tt92lsqJaOkFeXk1egLv8apX8OuPwS+D2//ztoaACdMG1Bsy/Hfmz4Sw3+76+hEbJjzyF1vidnbVd/3j/zTMkcIWu/3cCWzF9a/DmIiIiIiIiIiGpjhLyijdv2KBKN6Z8Dnq4zNtW29+DxxE8dPnriTOJ7Dl7Ztz/8pEg0pvc+qvuDaZIxQn721eWvVlu6/HtJlwfSGzs/ruHPvFHn/dZt3qX23TKcR8jjJ7N1Q6cBSh09PTG61lZeUVnnp2w3Vd65y98XssejE9Sx55A63/MxO/ecItGY7n/iWf2tx2DV1DT8/SCllo2Q+w79qqz9x+q9/eEhU9Xj0QmSLn/laKcHRqh7//F1vqL1TE6+bklJU8a4mYm3NTRC1o6Ceef++GnhR389rV6PTaw3QvZ6bGKDL92++vP+mWdqjRHyh807Ez8gqLaiCxeV0m+M/t5jcJ0fjkNERERERERE1NIYIa+q9ivAbr0nXWOef0tvLFyul9/8UP8e/IIi0Zhu7z008b61P4U4bcwMvffR15oy6z+6qXOqeg+YVGdAaq0R8vbeQ/Xsqws19/0VGjV5buL7M145FE58aYEi0Zgyxs3Sq299rNTR0/WXuwbo5pQ05xFSkua+v0KRaExdHxqrGW//V4v++62efXWhbu89VM9MX+j87/PeR8YpEo1pyMQ59X6vy4OjG/29K2vJCFn7suaMcbO08OOvtXT59xo9ZZ4i0Zj+c8VXGX63PlM3dBqgXo9N1FuLv9CcBZ/qjvueVPtuGXW+crOhEbJ27Ov60FhNe/0DDXv6dd3UOVUdug+qN0LW/t/KpJff1YKlKxM/ebyhz+v6TK0xQg4aP1u3dU3XiGff1Ox3PtGUmYt0e++hikRj+vzbTY7/tomIiIiIiIiImo4RsoF2/XJEo6fM0119RyZGpYeHTNXCj79WYVFx4v3i8biWfPqdejw6QTd1TtUd9z2pKTMX1fvKstYaIcc8/5b+NXCybr0nXf/oPUxT5yxRcUlpnY8tLSvXczPeV8eeQ/S3HoOVMW6W9uw7+qdejl3b6g2Zio14WX+9N0O3pKSp56MT9PzsxQ3+RPDGqv3+kld+ZWht416cr0g0pkX//bbJz9GSEbLgfLHmLPhU/xzwtG5OSVOH7oPUb9Dz+mrN1nrvu3Fblh4eMlW3pKTpr/dmaPCE2Tp0rO73vWxoLJSkZV+uU7eHx+nmlDT9a+BkrVq3vcGXY18oKdW4qfPVsecQ3dY1XcOefr3Jz+vyTK0xQv64fa9GPvemuj08Trfek6723TIUG/Gyfty+t9mPJSIiIiIiIiJyjREyALkMhkRERERERERERG01RsgAxAhJRERERERERERBjhEyADFCEhERERERERFRkGOEDECMkEREREREREREFOQYIYmIiIiIiIiIiCipMUISERERERERERFRUmOEJCIiIiIiIiIioqTGCElERERERERERERJjRGSiIiIiIiIiIiIkhojJBERERERERERESU1RkgiIiIiIiIiIiJKaoyQRERERERERERElNQYIYmIiIiIiIiIiCipMUISERERERERERFRUmOEJCIiIiIiIiIioqTGCElERERERERERERJjRGSiIiIiIiIiIiIkhojJBERERERERERESU1RkgiIiIiIiIiIiJKaoyQRERERERERERElNQYIYmIiIiIiIiIiCipMUISERERERERERFRUmOEJCIiIiIiIiIioqTGCElERERERERERERJjRGSiIiIiIiIiIiIkhojJBERERERERERESU1RkgiIiIiIiIiIiJKaoyQRERERERERERElNQYIYmIiIiIiIiIiCipMUISERERERERERFRUmOEJCIiIiIiIiIioqTGCElERERERERERERJjRGSiIiIiIiIiIiIkhojJBERERERERERESU1RkgiIiIiIiIiIiJKaoyQRERERERERERElNQYIYmIiIiIiIiIiCipMUISERERERERERFRUmOEJCIiIiIiIiIioqTGCElERERERERERERJjRGSiIiIiIiIiIiIkhojJBERERERERERESU1RkgiIiIiIiIiIiJKaoyQRERERERERERElNQYIYmIiIiIiIiIiCipMUISERERERERERFRUmOEJCIiIiIiIiIioqTGCElERERERERERERJjRGSiIiIiIiIiIiIkhojJBERERERERERESU1RkhPnTlXFggVVTU6d6HC+3MAaF0XSqtUXHbJ+3MAaF35RRWqvFTj/TkAtL4g3SEAuDl3oUKSVFEVjv92EzFCesr34XfFCAnYxAgJ2MQICdgVpDsEADf5RZdHyEpGSApJjJCe8n34XTFCAjYxQgI2MUICdgXpDgHADSMkhS1GSE/5PvyuGCEBmxghAZsYIQG7gnSHAOCGEZLCFiOkp3wffleMkIBNjJCATYyQgF1BukMAcMMISWGLEdJTvg+/K0ZIwCZGSMAmRkjAriDdIQC4YYSksMUI6Snfh98VIyRgEyMkYBMjJGBXkO4QANwwQlLYYoT0lO/D74oRErCJERKwiRESsCtIdwgAbhghKWwxQnrK9+F3xQgJ2MQICdjECAnYFaQ7BAA3jJAUthghPeX78LtihARsYoQEbGKEBOwK0h0CgBtGSApbjJCe8n34XTFCAjYxQgI2MUICdgXpDgHADSNkOJo6Z4mem/F+4teVlVW6OSVNH32+JvG2A0dOqkP3Qdr6835ForFGnTyd6+OP0GoxQnrK9+F3xQgJ2MQICdjECAnYFaQ7BAA3jJBtq8oqKf+cVHheqq5pvc+7ekOm7n1kXOLXP+85pPbdMjR6yrzE25Yu/16DJ8zWpepqXSgp1YWSUh08+psi0ZjyC4oSb6upibfeg3mIEdJTvg+/K0ZIwCZGSMAmRkjAriDdIQC4YYRsOxUUxrV77x/2H4q32hBZdOGibug0QPkFRZKk+Uu+1Kz5yxS9f0TifUZNnqtFy76t83Gnc/IVicZUUVnVOg/SBmKE9JTvw++qxSPk2RLl7DusnH2HlX3stPc/B4C6GCEBmxghAbuCdIcA4IYRsu2Uta+mzgi5e29cZ7Jb76sO+6Q+o29/+EmSlDZmhvYd+lX3PjIu8fLqO/sM175Dv9b5GEZIarV8H35XLRkhs7PPK2f9euWuXZuQnbXP+58FwB8YIQGbGCEBu4J0hwDghhGybVRyUfUGyN174zpyrPVGyFfmfqQXZi9WdXWNOj0wQtXVNXpuxvta/vUGnTydq7/1GFzvpdaMkNRq+T78rlo0QmbtqzNA1jpztsT7nwfAZYyQgE2MkIBdQbpDAHDDCNl2SvZXQv6wead6D5ikrAPHNXDsTEnSytVbNPGlBVr+9QYNnfhavY9hhKRWy/fhd9WSETJne2aDI2T2iVzvfx4AlzFCAjYxQgJ2BekOAcANI2Tb6ervCbnvYE2r/nCa4pJS3dj5cc1Z8KkWLF0pScrJK9A9/35Kk15+V4s/WVXvYxghqdXyffhdtegrIQ8cqzdA5mxY7/3PAuAPjJCATYyQgF1BukMAcMMI2baqrJLy8uMqKGy9H0pzZQ8OnKK/3puhHVmHEm9L6TdGHXsN0YEjJ+u9PyMktVq+D7+rlv5gmiu/GjJnw3p+OA3QxjBCAjYxQgJ2BekOAcANI2S4mvn2Mt2ckqaqqkuJt016+V39vYHvBykxQlIr5vvwu2rxT8c+V6bsU/mXX4LN94IE2hxGSMAmRkjAriDdIQC4YYSksMUI6Snfh9/VtYyQANouRkjAJkZIwK4g3SEAuGGEpLDFCOkp34ffFSMkYBMjJGATIyRgV5DuEADcMEJS2GKE9JTvw++KERKwiRESsIkRErArSHcIAG4YISlsMUJ6yvfhd8UICdjECAnYxAgJ2BWkOwQAN4yQFLYYIT3l+/C7YoQEbGKEBGxihATsCtIdAoAbRkgKW4yQnvJ9+F0xQgI2MUICNjFCAnYF6Q4BwA0jJIUtRkhP+T78rhghAZsYIQGbGCEBu4J0hwDghhGSwhYjpKd8H35XjJCATYyQgE2MkIBdQbpDAHDDCElhixHSU74PvytGSMAmRkjAJkZIwK4g3SEAuGGEpLDFCOkp34ffFSMkYBMjJGATIyRgV5DuEADcMEKGp592HlD/YS+qfbcMdew5RP8aOFnvffR14vc7PTBCew8eb/Hn373vqCLRmC5VV7fo42fNX6ZINFZH+lMzWvw8jcUI6Snfh98VIyRgEyMkYBMjJGBXkO4QANwwQraxKitUc/aM4ufOStWXWu3Tnjydq9u6puvjL9Yqv6BIOXkF+mbtNg2d+Frifa51hCwrr9SBIydb/PGz5i/T+Gnv6EJJaUJpWUWLP19jMUJ6yvfhd8UICdjECAnYxAgJ2BWkOwQAN4yQbaea/FxVZW5KuJSV2WpD5JffbVb3/uMb/f1npi9UJBrTLSlpuq1rut776GsdPn5KHboP0rL//aBuD4/VnX2Gq7SsQn1Sn1H7bhm6rWu6+g+bpiPHT0uSjp3M1m1d0xOfc8Szb2rKrP8oY9xMte+WoX8NnKyTp882+gyz5i/TM9MXtsqft6kYIT3l+/C7YoQEbGKEBGxihATsCtIdAoAbRsi2U9XOLXVGyKrMTar5reVfmXhlx05m68bOj2vu+yu0//AJVVXVHzev/krIw8dPKRKN6elX3lPW/mM6cvy0qi5V6+c9h1R04aLKyiu1+JNVGjDy5cQ/4+oRsnv/8Vq7aYdOZedpyMQ5Gj/tnUafcdb8Zbqta7qi94/Q/U88q4Uff93o+15LjJCe8n34XTFCAjYxQgI2MUICdgXpDgHADSNk2yheXFRvgKzK3KRLB7Na7Z+x9ed9yhg3S3+9N0PtuqSq/7Bp2vrzvsTvNzRCdug+qN7nWbl6iwZPmK27+43R7b2H6o77npTU8Ai5+JNViV9/tWarHhw4pdHnO3k6V8dOZutMTr7W/rhTd/YZrg9XrLmmP3NDMUJ6yvfhd8UICdjECAnYxAgJ2BWkOwQAN4yQbadkfiXklcXjcZ08nas3Fi7XzSlpyi8okuQ2Qq78fovu7DNc67fsVll5pQ4fP6WOvYZIan6EXLd5l/qkPuP8nO988KUeH/lKi/6MTcUI6Snfh98VIyRgEyMkYBMjJGBXkO4QANwwQradrv6ekFV7trfa94QsK6+s/8+rieumzqmJ4bHLv0Ypa/+xxO83NEI+++pCzXx7WZ33SdYIOWfBpxo1ea7z+7vGCOkp34ffFSMkYBMjJGATIyRgV5DuEADcMEK2sSorVJ17WjX5ua3607Hf/fArDXv6df2855DOF5XoVHaeps/9SN0eHquqS9WSpEeGvqi3F/9PF0vLVXC+uMERcsHSleqXMUWnsvN08OhvSn9qRquNkC/MXqxtO/crv6BIG7dlqWPPIVqzcUer/TuojRHSU74PvytGSMAmRkjAJkZIwK4g3SEAuGGEDEdHT5zR06+8p64PjdXNKWmK3j9CT73wlk5l5yXeZ+uOfery4Gi1u/sJffzF2gZHyIul5Ro8YbZuSUlTr9gkvblwRauNkC+9sVR39xujm1PS1O3hcfrvF2tb649fJ0ZIT/k+/K4YIQGbGCEBmxghAbuCdIcA4IYRksIWI6SnfB9+V4yQgE2MkIBNjJCAXUG6QwBwwwhJYYsR0lO+D78rRkjAJkZIwCZGSMCuIN0hALhhhKSwxQjpKd+H3xUjJGATIyRgEyMkYFeQ7hAA3DBCUthihPSU78PvihESsIkRErCJERKwK0h3CABuGCEpbDFCNtDUOUsUicYUicbUsdcQjXj2TeWdO5/4/dKycs1+5xN1e3ic2nVJ1R33PamMcTP17Q8/JX4/Eo0pJ68g8TGHj5/S7b2H6ofNOyUF5y8QjJCATYyQgE2MkIBdQbpDAHDDCElhixGygabOWaIps/6j8opK/XbmrIZOfE39h70oSSqvqFTf9MnqP+xF7cg6rKLiizqdk6/Pv92kh4dMlVR/hMw+W6DOfUdpxTcbE/8M34ff1fUcIU/llevIb2U6kVPu/c8NWMcICdjECAnYFaQ7BAA3jJAUthghG2jqnCWaOmdJ4tdHT5xRJBpTYVGxFi37Vnf3G6Pyisp6HxePxyXVHSGLLlxUr8cm6v2Pv6nzvr4Pv6vrNULuP1amddvKtXbzZbv2M0QCycQICdjECAnYFaQ7BAA3jJAUthghG+jqEfLI8dO6odMAlVdUKnX0dL327mdNfnztCHniVK4eHjJVs+Yvq/c+vg+/q+sxQp7I+WN8vNLBEwyRQLIwQgI2MUICdgXpDgHADSMkhS1GyAa6coQsuVimYU+/roxxsyRJ9z/xrD7+Ym3ifect+lwdug9Sh+6D1OmBEZL+GCH/PfgF3dQ5VcdPZtf7Z1RdqgmEmnhcl6qT+884X1SjzN3V9fx2ptr7nx9IiuprU1ZRrUvVNdekpiaumpr4NX8eAG1LdU1c8bi8PweA1if5P9tlFdXX/PcYNMD3303hTe3ZjseDsxFcCyJGyAaaOmeJbk5JU8eeQ3RT51QNmThH5wovSJJSR0/X7Hc+Sbxv7feE/GHzTt3WNV3SHyPknAWf6vnZi9X1obEqLCqu8884e748ECqranS+pDKp/4yTueVav7W+A8f8//mBpChsBecrrklJ+SVdrKi+5s8DoG0pLKlUVXWN9+cA0Pou3yH8PkNua/09Blf9vQ5hVVhy+du8VV2q8f4s1wMRI2QDTZ2zRJNefle5eYUqLauo83vvffS17uo7st7b9x48Xm+EzMkrUHV1jTLGzdIjQ19UZWVV4v19fxm0q+v1PSG3Z9V9KfbG7eU6lcfLsYFk4eXYgE28HBuwK0h3CABueDk2hS1GyAa6+ntCXllpWbl6D5ikfhlTtG3nfp0vKtGZnHzNX/JlgyNk7a/vf+JZjX3h7cTn8X34XV3Pn46998jlMXLPwTIGSCDJGCEBmxghAbuCdIcA4IYRksIWI2QDNTVCSlJxSalefvNDpfQbo5s6p+quviM1dOJrWr0hU1L9EVKScvMKdVffkXpj4XJJwfkLxPUcIQFcP4yQgE2MkEILXocAACAASURBVIBdQbpDAHDDCElhixHSU74PvytGSMAmRkjAJkZIwK4g3SEAuGGEpLDFCOkp34ffFSMkYBMjJGATIyRgV5DuEADcMEJS2GKE9JTvw++KERKwiRESsIkRErArSHcIAG4YISlsMUJ6yvfhd8UICdjECAnYxAgJ2BWkOwQAN4yQFLYYIT3l+/C7YoQEbGKEBGxihATsCtIdAoAbRkgKW4yQnvJ9+F0xQgI2MUICNjFCAnYF6Q4BwA0jJIUtRkhP+T78rhghAZsYIQGbGCEBu4J0hwDghhGSwhYjpKd8H35XjJCATYyQgE2MkIBdQbpDAHDDCElhixHSU74PvytGSMAmRkjAJkZIwK4g3SEAuGGEpLDFCOkp34ffFSMkYBMjJGATIyRgV5DuEADcMEJS2GKE9JTvw++KERKwiRESsIkRErArSHcIAG4YISlsMUJ6yvfhd8UICdjECAnYxAgJ2BWkOwQAN4yQFLYYIT3l+/C7YoQEbGKEBGxihATsCtIdAoAbRkgKW4yQnvJ9+F0xQgI2MUICNjFCAnYF6Q4BwA0jJIUtRkhP+T78rhghAZsYIQGbGCEBu4J0hwDghhGSwhYjpKd8H35XjJCATYyQgE2MkIBdQbpDAHDDCElhixHSU74PvytGSMAmRkjAJkZIwK4g3SEAuGGEpLDFCOkp34ffFSMkYBMjJGATIyRgV5DuEADcMEJS2GKE9JTvw++KERKwiRESsIkRErArSHcIAG4YISlsMUJ6yvfhd8UICdjECAnYxAgJ2BWkOwQAN4yQFLYYIT3l+/C7YoQEbGKEBGxihATsCtIdAoAbRkgKW4yQnvJ9+F0xQgI2MUICNjFCAnYF6Q4BwA0jJIUtRkhP+T78rhghAZsYIQGbGCEBu4J0hwDghhGSwhYjpKd8H35XjJCATYyQgE2MkIBdQbpDAHDDCElhixHSU74PvytGSMAmRkjAJkZIwK4g3SEAODhbonN7f1HVzq0q27ND2cdO+3+mJCNihPSU78PvihESsIkRErCJERKwK0h3CADNy9meqXPrf1BV5iaVbt2g3LVrlX0i1/tzJRMRI6SnfB9+V4yQgE2MkIBNjJCAXUG6QwBoWnb2eeWuXVtvhMzZtdf7syUTESOkp3wffleMkIBNjJCATYyQgF1BukMAaFr2qXxGSApljJCe8n34XTFCAjYxQgI2MUICdgXpDgGgeblbttV/Obbx7wtJxAjpKd+H3xUjJGATIyRgEyMkYFeQ7hAAmpedfV7ndu1WVeYmle38SdlHTnh/pmQjYoT0lO/D74oRErCJERKwiRESsCtIdwgAbvKLKiRJlVXh+G83ESOkp3wffleMkIBNjJCATYyQgF1BukMAcMMISWGLEdJTvg+/K0ZIwCZGSMAmRkjAriDdIQC4YYSksMUI6Snfh98VIyRgEyMkYBMjJGBXkO4QANwwQlLYYoT0lO/D74oRErCJERKwiRESsCtIdwgAbhghKWwxQnrK9+F3xQgJ2MQICdjECAnYFaQ7BAA3jJAUthghPeX78LtihARsYoQEbGKEBOwK0h0CgBtGSApbjJCe8n34XTFCAjYxQgI2MUICdgXpDgHADSMkha1mR8jzRSU6euJM4tcnTuXqwxVrtGff0aQ+mPV8H35XjJCATYyQgE2MkIBdQbpDAHDDCElhq9kRctTkuRr7wtuSpLxz59Wh+yD9o/cwteuSqu/WZyb9Aa3m+/C7YoQEbGKEBGxihATsCtIdAoAbRkgKW82OkHf2Ga7NmXslSR98tlp39xujqkvVWrl6i+5/4tmkP6DVfB9+V4yQgE2MkIBNjJCAXUG6QwBwwwhJYavZEfLmlDSdzsmXJI2bOl/PvrpQkvTbmbO69Z705D6d4XwffleMkIBNjJCATYyQgF1BukMAcMMISWGr2RGy56MTtGrddpWWVej//jlMn321QZK0c+9h3dlneNIf0Gq+D78rRkjAJkZIwCZGSMCuIN0hALhhhKSw1ewI+fm3m9SuS6o69hqiO+57UhdKSiVJz766UBnjZib9Aa3m+/C7YoQEbGKEBGxihATsCtIdAoAbRkgKW82OkJKUufugli7/Xmd+f1n2pepqvb34f9rNT8hucb4PvytGSMAmRkjAJkZIwK4g3SEAuGGEpLDV7Aj59uL/qbSsot7bjxw/rWVfrkvGM4Ui34ffFSMkYBMjJGATIyRgV5DuEADcMEJS2Gp2hIxEY8ovKKr39n2HflX7bhlJeagw5Pvwu2KEBGxihARsYoQE7ArSHQKAG0ZIClstHiG/WrNVf+8xOCkPFYZ8H35XjJCATYyQgE2MkIBdQbpDAHDDCElhq9ERsldsknrFJikSjal7//GJX/eKTVK3h8cqEo1p+tyPruezmsr34XfFCAnYxAgJ2MQICdgVpDsEADeMkBS2Gh0h3//4G73/8TeKRGN6c+GKxK/f//gbfbhijbbvOnA9n9Ncvg+/K0ZIwCZGSMAmRkjAriDdIQC4YYSksNXsy7FXrduuysqq6/Esocr34XfFCAnYxAgJ2MQICdgVpDsEADeMkBS2mh0hKTn5PvyuGCEBmxghAZsYIQG7gnSHAOCGEZLCVqMjZKcHRuj7jT+rT+ozTaKW5fvwu2KEBGxihARsYoQE7ArSHQKAG0ZICluNjpCvzvtY+w79qgVLVzaJWpbvw++KERKwiRESsIkRErArSHcIAG4YISls8XLsRiosKlYkGlMkGlO7u59Q9/7jteKbjXXe56edB5T+1Az9rcdg3ZKSpp6PTtDMt5cp+2yBJGnqnCWaOmdJg5/f9+F3xQgJ2MQICdjECAnYFaQ7BAA3jJAUtpxGyHg8rrP553XsZHY9VqsdIXPzCnWhpFRrNu7QLSlp2py5V5L09Zptuq1ruhYt+1ansvN0vqhEu345oqlzluijz9dIYoQE0HYxQgI2MUICdgXpDgHAzfFTFTp5Kq7fztToVF659+dJNqJmR8itP+/TXX1HJr4q8GpWqx0hC4uKE28bNXmeps5ZosrKKv2j9zD994u1DX5sPB6XxAgJoO1ihARsYoQE7ArSHQJA83buL9OGbRXalVWjzF3VWretXCdybA+RRM2OkPc+Mk6Dxs9W1v5jOpWdp9M5+XVYraERcuRzb2r2O59oR9Zh3dBpgCorq5r8HE2NkGfPlwdCZVWNzpdUen8OAHXlFVVck4vll1RaUX3NnwdA23K+pFJV1XHvzwGg9Uny/gy+//4DWJFTUKH1W8u1aXvl5RFyd7XWby1X1kH/z5ZMRM2OkDenpGnD1t3X41naVFePkFsyf9HNKWnas++o1m7aof/757DE+1ZWVqlD90EJS5d/L6npEbLqUk0g1MTjulTt/zkA/KGsovqaP0d1TVzVNXHvfxYAretSdVzxOGcbsKgt3CFa4+8gAGp0vujy8Pjz7mrtyqrRzt+HyANHbJ8xomZHyNiIl7X4k1XX41naVLUjZMeeQ9S+20Dd3nuoVq7eIknakXVYkWhM5RWVki6//Lr2K0MfH/mK3v/4G0m8HBtA28XLsQGbeDk2YFeQ7hAAmrd+W3mdl2Ov3VyuvUf8P1cyETU7QmYdOK47+wxX5u6D2pF1WDuyDtVhtdoR8tCxUzpXeEE1NfHE71VUVunvPQZryaff1fu4QeNnM0ICaPMYIQGbGCEBu4J0hwDQvIMnyvXj9j9GyK27bH8/yDPnGCHJYYRcu2mHbuuazg+muaoV32zUzSlpeueDL3XiVK6Kii9q597D6hWbVGeEfGb6QhUWFScUXbgoKTh/gWCEBGxihARsYoQE7ArSHQKAm7OFlSq5KBUVh+O/3UTNjpBd/jVKY55/KzG0XSgprcNqzY2QkrRx2x7FRrysv96bodu6puv+J57VnAWfKievQNLlEfLq0bZD90GSgvMXCEZIwCZGSMAmRkjAriDdIQC4yf/9h05VVoXjv91ETj+Y5sfte6/Hs4Qq34ffFSMkYBMjJGATIyRgV5DuEADcMEJS2Gp2hBw68TW988GX1+NZQpXvw++KERKwiRESsIkRErArSHcIAG4YISlsNTtCLl2+Wu27Zei79Zlat3lXPdSyfB9+V4yQgE2MkIBNjJCAXUG6QwBwwwhJYavZEfIvdw1oErUs34ffFSMkYBMjJGATIyRgV5DuEADcMEJS2Gp2hKTk5Pvwu2KEBGxihARsYoQE7ArSHQKAG0ZICluMkJ7yffhdMUICNjFCAjYxQgJ2BekOAcANIySFrUZHyE9XrtfJ02e1ZuOOJlHL8n34XTFCAjYxQgI2MUICdgXpDgHADSMkha1GR8hINKYvv9usdnc/0SRqWb4PvytGSMAmRkjAJkZIwK4g3SEAuGGEpLDV6AhZWlauS9XV1/NZQpXvw++KERKwiRESsIkRErArSHcIAG4YISls8T0hPeX78LtihARsYoQEbGKEBOwK0h0CgBtGSApbDY6QL72x1Bm1LN+H3xUjJGATIyRgEyMkYFeQ7hAA3DBCUthqcIRMf2qGM2pZvg+/K0ZIwCZGSMAmRkjAriDdIQC4YYSksMXLsT3l+/C7YoQEbGKEBGxihATsCtIdAoAbRkgKW4yQnvJ9+F0xQgI2MUICNjFCAnYF6Q4BwA0jJIUtRkhP+T78rhghAZsYIQGbGCEBu4J0hwDghhGSwhYjpKd8H35XjJCATYyQgE2MkIBdQbpDAHDDCElhixHSU74PvytGSMAmRkjAJkZIwK4g3SEAuGGEpLDFCOkp34ffFSMkYBMjJGATIyRgV5DuEADcMEJS2Gp2hNy2c3+jv7f1532t+jBhyvfhd8UICdjECAnYxAgJ2BWkOwQAN4yQFLaaHSHbdUnVOx98qXg8nnhbdXWNXn/vM93QaUBSH85yvg+/K0ZIwCZGSMAmRkjAriDdIQC4YYSksNXsCJm5+6A69x2ljHGzVHThorLPFqj/sBf1tx6D9dWardfjGU3m+/C7YoQEbGKEBGxihATsCtIdAoAbRkgKW07fE/JCSalGTZ6rzn1HqWPPIRow8mVlny1I9rOZzvfhd8UICdjECAnYxAgJ2BWkOwQAN4yQFLacRsiy8kpNnrFIN3QaoJs6p2rJp9/VeXk2/fl8H35XjJCATYyQgE2MkIBdQbpDAHDDCElhq9kRcv/hE+ref7yi94/QTzsPaEfWIXV6YIQGjp2p/IKi6/GMJvN9+F0xQgI2MUICNjFCAnYF6Q4BwE1+QanixUWqvFDs/VmuByKnH0wzaPxsFRYVJ95WWFSsgWNn6h+9hyX14Szn+/C7YoQEbGKEBGxihATsCtIdAkDzso+d1rlNG1WVuUmlWzco56dM78+UbETNjpCLln3b4Euv4/G4FixdmZSHCkO+D78rRkjAJkZIwCZGSMCuIN0hADQvZ8N6nVv/Q2KEzF27VtkHjnl/rmQicvqekNT6+T78rhghAZsYIQGbGCEBu4J0hwDQtOwTucpdu7beCJmz3fZXQxI1O0JOmbmoSdSyfB9+V4yQgE2MkIBNjJCAXUG6QwBoxtmSBkfI7Kx9/p8tiYiaHSGHTXqtnoFjZyoSjWnQ+NnX4xlN5vvwu2KEBGxihARsYoQE7ArSHQJA87Kz9tUZIXM2rFd29nnvz5VMRC1+OfbgCbO1aNm3rfksocr34XfFCAnYxAgJ2MQICdgVpDsEADf5J7NV/ethlZ/4VWfOlnh/nmQjavEI+f3Gn9W576jWfJZQ5fvwu2KEBGxihARsYoQE7ArSHQKAm/yiCklSZVU4/ttN1OIRcsU3G9W+W0ZrPkuo8n34XTFCAjYxQgI2MUICdgXpDgHADSMkha1mR8i576+o482FKzRl1n/UofsgjZ4y73o8o8l8H35XjJCATYyQgE2MkIBdQbpDAHDDCElhq9kRsm/65HpSR03X6+99poul5dfjGU3m+/C7YoQEbGKEBGxihATsCtIdAoAbRkgKWy1+OTZdW74PvytGSMAmRkjAJkZIwK4g3SEAuGGEpLDFCOkp34ffFSMkYBMjJGATIyRgV5DuEADcMEJS2GpwhHzpjaXOqGX5PvyuGCEBmxghAZsYIQG7gnSHAOCGEZLCVoMjZPpTM5xRy/J9+F0xQgI2MUICNjFCAnYF6Q4BwA0jJIUtXo7tKd+H3xUjJGATIyRgEyMkYFeQ7hAA3DBCUthqdITs8uBo7T98IvHrk6dzVVlZdV0eKgz5PvyuGCEBmxghAZsYIQG7gnSHAOCGEZLCVqMjZCQa0+59RyVJl6qrFYnGtPfg8ev2YNbzffhdMUICNjFCAjYxQgJ2BekOAcANIySFLUZIT/k+/K4YIQGbGCEBmxghAbuCdIcA4IYRksIWI6SnfB9+V4yQgE2MkIBNjJCAXUG6QwBwwwhJYYsR0lO+D78rRkjAJkZIwCZGSMCuIN0hALhhhKSw1eQI+e/BLyj9qRlKGzNDkWhMD/3+6ytRy/J9+F0xQgI2MUICNjFCAnYF6Q4BwA0jJIWtRkfIiS8tcEIty/fhd8UICdjECAnYxAgJ2BWkOwQAN4yQFLYaHSEpufk+/K4YIQGbGCEBmxghAbuCdIcA4IYRksIWI6SnfB9+V4yQgE2MkIBNjJCAXUG6QwBwwwhJYYsR0lO+D78rRkjAJkZIwCZGSMCuIN0hALhhhKSwxQjpKd+H3xUjJGATIyRgEyMkYFeQ7hAA3DBCUthihPSU78PvihESsIkRErCJERKwK0h3CABuGCEpbDFCesr34XfFCAnYxAgJ2MQICdgVpDsEADeMkBS2GCE95fvwu2KEBGxihARsYoQE7ArSHQKAG0ZIClttfoT8ZOU69U2frFvvSVeH7oPUb9Dz+s8nq1RWXpl4n937jiptzAy175ahm1PSlJtXKEkqLCpWJBpTJBpTu7ufUPf+47Xim42Jj+sVm6RINKYbOg3Q33oM1uMjX9H3G3+u88+fOmeJps5ZUudtF0pKFYnGlF9QJElate4nPZD2nNp3G6j/++cwTZm5SBWVVU3+uXwffleMkIBNjJCATYyQgF1BukMAcMMISWGrTY+Qz89erM59R2n1hkzlFxQpv6BIG7buVsa4mTpw5KQkadcvR/SP3sP0ycp1yskrUHbuOVX+PgDWjpC5eYW6UFKqNRt36JaUNG3O3Cvp8gj51ZqtKi0r129nzmrp8u91e++hemvxF4lncBkhP/p8jVat+0m5eYU6ePQ39Xx0gt5YuLzJP5vvw++KERKwiRESsIkRErArSHcIAG4YISlstdkRct+hX3Vj58d1+PipBn8/Ho9LkvoPe1GffbWhwfepHSELi4oTbxs1eV5iVOwVm6RV67bX+Zj9h0/oxs6P6+Tps5LcRsire2Phcg2d+Fri1x17DdEHn63Wvwe/oHZdUnX0xBnvh98VIyRgEyMkYBMjJGCXxAgJWHIip1y79lVoV1aN9vxSrV+Olnt/pmQjarMj5PwlX+qRoS82+T7Fv4+Bo6fM0x33Pan23QZq5HNvqrikVFLDI+TI597U7Hc+kdTwCClJ/xo4WR98tlrS5RGyXZdUdeg+KOGv92Y0OUIOHDtTr737WeLXHXsN0cNDpurH7Xt17GS2ysorvR9+V4yQgE2MkIBNjJCAXRIjJGDJ5h3l2rDt8giZuataazeX6+AJ20MkUZsdIV9+80ONmjwv8eutP++rMwQeOHJSh46dUiQa0/otu1VaVq7cvELFRryc+LirR8gtmb/o5pQ07dl3VFLjI+SQiXMSI+LUOUs0YdoCnc7JTzh49LdGR8gV32xU9P4RKjj/x/DZsdeQxMvHaysurQqES9VxlVZc8v4c8KgMbU3RxWv/HBVV1aqoqvH+ZwHQukorLqm6Ju79OQC0Pknen6E1/g4CXBe+71DNyCus0pbMKm3bUaVdWTXasadaWzKrdOCo/2dLJqI2O0LOX/Kl/j34hcSvyysqEyPgrfeka+/B4zp8/PIIWVMTT7zfpp+y1L5bhqQ/RsiOPYeofbeBur33UK1cvSXxvk19JeSHK9ZI+nMvx161brv+0XtYvcGxoRHyQmlVIFyqrlFp+SXvzwGPLqKtKbpYpeLSS9ekoqpGlZdqrvnzAGhbSiuqL4+QbeBZALQuSd6f4fzFSu9/DwKc+L5DNeNcUZU2Z1Zp6xUj5ObMKu0/6v/ZkomozY6QWfuP6YZOA7Tv0K/1fu+2rpdHyNKycv3lrgE6dOyP7xu5fstu3d1vjKQ/RshDx07pXOGFOmOl1PAIeeDISbW7+wnl5BVIch8hl/3vB93Vd2SdZ6mtoRHS95dBu+Ll2IBNF0qrxMuxAXt4OTZgV5DuEACat3VX/ZdjH/nN/3MlE1GbHSEladLL7+rOPsO18vstys0rVH5BkdZt3qV2XVK19+BxSdJTL7ylx4a/pPyCIuWdO6/+w6bp9fcuv5S6oe8JeWV//HTsCp3JydfHX6zV//1zmBZ+/HXifVxGyHmLPte9j4zT0V9Pq7SsPKH2h+cwQgJoaxghAZsYIQG7gnSHANC8U3nl2nOwUln7avTLgWrz3w/yzDlGSGrjI2RNTVwffLZa9z3+tG5OSVPHXkP06JPTtHT5apWVV0qSLpaW6+lX3lPHnkN0Z5/hmv3OJ6q6VC3JbYSMRGOKRGPq0H2QHh/5itZu2lHnfVxGyD6pzyQ+z5Vqf58REkBbwwgJ2MQICdgVpDsEADf5RRWSpMqqcPy3m6hNj5CW8334XTFCAjYxQgI2MUICdgXpDgHADSMkhS1GSE/5PvyuGCEBmxghAZsYIQG7gnSHAOCGEZLCFiOkp3wffleMkIBNjJCATYyQgF1BukMAcMMISWGLEdJTvg+/K0ZIwCZGSMAmRkjAriDdIQC4YYSksMUI6Snfh98VIyRgEyMkYBMjJGBXkO4QANwwQlLYYoT0lO/D74oRErCJERKwiRESsCtIdwgAbhghKWwxQnrK9+F3xQgJ2MQICdjECAnYFaQ7BAA3jJAUthghPeX78LtihARsYoQEbGKEBOwK0h0CgBtGSApbjJCe8n34XTFCAjYxQgI2MUICdgXpDgHADSMkhS1GSE/5PvyuGCEBmxghAZsYIQG7gnSHAOCGEZLCFiOkp3wffleMkIBNjJCATYyQgF1BukMAcMMISWGLEdJTvg+/K0ZIwCZGSMAmRkjAriDdIQC4YYSksMUI6Snfh98VIyRgEyMkYBMjJGBXkO4QANwwQlLYYoT0lO/D74oRErCJERKwiRESsCtIdwgAbhghKWwxQnrK9+F3xQgJ2MQICdjECAnYFaQ7BAA3jJAUthghPeX78LtihARsYoQEbGKEBOwK0h0CgBtGSApbjJCe8n34XTFCAjYxQgI2MUICdgXpDgHADSMkhS1GSE/5PvyuGCEBmxghAZsYIQG7gnSHAOCGEZLCFiOkp3wffleMkIBNjJCATYyQgF1BukMAcMMISWGLEdJTvg+/K0ZIwCZGSMAmRkjAriDdIQC4YYSksMUI6Snfh98VIyRgEyMkYBMjJGBXkO4QANwwQlLYYoT0lO/D74oRErCJERKwiRESsCtIdwgAbhghKWwxQnrK9+F3xQgJ2MQICdjECAnYFaQ7BAA3jJAUthghPeX78LtihARsYoQEbGKEBOwK0h0CgBtGSApbjJCe8n34XTFCAjYxQgI2MUICdgXpDgHATf7x31T962FVnDyh7Ozz3p8n2YgYIT3l+/C7YoQEbGKEBGxihATsCtIdAkDzsrP26dz6H1SVuUmlWzcoZ/1680MkESOkp3wffleMkIBNjJCATYyQgF1BukMAaMbZEuWuXVtnhMxdu1bZWfv8P1sSETFCesr34XfFCAnYxAgJ2MQICdgVpDsEgKZln8htcITM2Z7p/dmSiYgR0lO+D78rRkjAJkZIwCZGSMCuIN0hADQvZ8P6+l8JeeCY9+dKJiJGSE/5PvyuGCEBmxghAZsYIQG7gnSHANC87BO5Ordp4x/fE/LnXd6fKdmIGCE95fvwu2KEBGxihARsYoQE7ArSHQKAm/yCUsWLi1R5odj7s1wPRIyQnvJ9+F0xQgI2MUICNjFCAnYF6Q4BwE1+UYUkqbIqHP/tJmKE9JTvw++KERKwiRESsIkRErArSHcIAG4YISlsMUJ6yvfhd8UICdjECAnYxAgJ2BWkOwQAN4yQFLYYIT3l+/C7YoQEbGKEBGxihATsCtIdAoAbRkgKW4yQnvJ9+F0xQgI2MUICNjFCAnYF6Q4BwA0jJIUtRkhP+T78rhghAZsYIQGbGCEBu4J0hwDghhGSwhYjpKd8H35XjJCATYyQgE2MkIBdQbpDAHDDCElhixHSU74PvytGSMAmRkjAJkZIwK4g3SEAuGGEpLDFCOkp34ffFSMkYBMjJGATIyRgV5DuEADcMEJS2GKE9JTvw++KERKwiRESsIkRErArSHcIAG4YISlsMUJ6yvfhd8UICdjECAnYxAgJ2BWkOwQAN4yQFLYYIT3l+/C7YoQEbGKEBGxihATsCtIdAoAbRkgKW4yQnvJ9+F0xQgI2MUICNjFCAnYF6Q4BwA0jJIUtRkhP+T78rhghAZsYIQGbGCEBu4J0hwDghhGSwhYjpKd8H35XjJCATYyQgE2MkIBdQbpDAHDDCElhixHSU74PvytGSMAmRkjAJkZIwK4g3SEAuGGEpLDFCOkp34ffFSMkYBMjJGATIyRgV5DuEADcMEJS2GKE9JTvw++KERKwiRESsIkRErArSHcIAG4YISlsMUJ6yvfhd8UICdjECAnYxAgJ2BWkOwQAN4yQFLYYIT3l+/C7YoQEbGKEBGxihATsCtIdAoAbRkgKW4yQnvJ9+F0xQgI2MUICNjFCAnYF6Q4BwA0jJIUtRkhP+T78rhghAZsYIQGbGCEBu4J0hwDghhGSwhYjpKd8H35XjJCATYyQgE2MkIBdQbpDAHDDCElhixGyiQqLihWJxuqZ+/4KrVq3vc7bOnQfpBHPCVWwNAAAEbZJREFUvqmC88X1Prbd3U+oe//xWvHNxsTn9n34XTFCAjYxQgI2MUICdgXpDgHADSMkhS1GyCaqHRJ/O3NWF0pKEyoqq7Rq3Xb1ik3SpepqVVZW6cSpXN3z76e0+JNVdT42N69QF0pKtWbjDt2SkqbNmXslBecvEIyQgE2MkIBNjJCAXUG6QwBwk51fqZKLcRWcD8d/u4kYIZuodkgsLCqu93u1I+SV79u9/3ht27m/0Y8dNXmeps5ZIik4f4FghARsYoQEbGKEBOwK0h0CQPMOnijXj9srtCurRpm7qrV1V7n3Z0o2IkbIJmpuhLyx8+O6474ndXvvofrLXQM0YOTLys0rbPRjRz73pma/84mk4PwFghESsIkRErCJERKwK0h3CADNW7+tXBu2/TFCrt1crr1H/D9XMhExQjZR7ZD413sz1KH7IHXoPkj/HPC0pMsjZPf+45VfUKS8c+d19NfTemH2YvV8dIKqq2vqjZBbMn/RzSlp2rPvqCTpUnVNIMTjcVXXxL0/B4A/lFdWq7omfk1q4pdd6+cB0PbE4/L+DABan+T/bJdXVnv/exBgwfkLNfp5d7V27K7Wrqwa7cy6/OuDR2yfMSJGyCaqHRL3Hz6h0zn5Op2Tr5y8Akn1X44tSWXllYpEYzpw5GTiYzv2HKL23Qbq9t5DtXL1lsT75haWB0JFVY0Kiyu8PweAq5y/NiVll3SxvPqaPw+AtqWwuFJV1TXenwNA65Pk/Rm8//0HMCL7XIXWbS3Xxp8qL38l5O5qrdtarj0H/D9bMhExQjbRn/mekJJUWlauSDSmk6dzEx976NgpnSu8oJrf/7+Xtfn+MmhXvBwbsImXYwM28XJswK4g3SEANG/n/rI6L8dev61cJ3Jsf19IIkbIJmp2hHxsoioqq1ReUakzOfka9+J89U2frHg83uTHSsH5CwQjJGATIyRgEyMkYFeQ7hAA3Bw/VaGc3LhOnq7WqTzbA+SZc4yQxAjZZM2NkJFoLKF9twwNnfiass8WNPuxUnD+AsEICdjECAnYxAgJ2BWkOwQAN/lFFZKkyqpw/LebiBHSU74PvytGSMAmRkjAJkZIwK4g3SEAuGGEpLDFCOkp34ffFSMkYBMjJGATIyRgV5DuEADcMEJS2GKE9JTvw++KERKwiRESsIkRErArSHcIAG4YISlsMUJ6yvfhd8UICdjECAnYxAgJ2BWkOwQAN4yQFLYYIT3l+/C7YoQEbGKEBGxihATsCtIdAoAbRkgKW4yQnvJ9+F0xQgI2MUICNjFCAnYF6Q4BwA0jJIUtRkhP+T78rhghAZsYIQGbGCEBu4J0hwDghhGSwhYjpKd8H35XjJCATYyQgE2MkIBdQbpDAHDDCElhixHSU74PvytGSMAmRkjAJkZIwK4g3SEAuGGEpLDFCOkp34ffFSMkYBMjJGATIyRgV5DuEADcMEJS2GKE9JTvw++KERKwiRESsIkRErArSHcIAG4YISlsMUJ6yvfhd8UICdjECAnYxAgJ2BWkOwQAN+cuXB4hKxghKSQxQhIREREREREREVFSY4QkIiIiIiIiIiKipMYISUREREREREREREmNEZKIiIiIiIiIiIiSGiMkNdobC5er28PjdHNKmlL6jdEHn62u8/t7Dx7XA2nPqV2XVN37yDht3LbH05MSUUv6YtWP6hWbVO/ts+YvUyQaqyP9qRkenpCIWtKxk9m6JSVNhUXFdd5eUxPX9Lkf6e89Buu2ruka9+J8lVdUenpKImppFZVV9f47HYnGtHvfUd+PRkR/osKiYg2d+JpuSUnTHfc9qQVLV/p+JKKkxwhJjTbj7f/qp50HdK7wgjZu26NbUtK0dcc+SVJV1SXd1XekFixdqQslpfryu81q3y1DBeeLm/msROS7s/nn1e3hsfpbj8GNjpDjp72jCyWlCaVlFR6elIj+bKmjp+uO+55UJBqrN0J+/MVadXt4rI6dzFZOXoH6D3tRM97+r6cnJaKWVjtCHv31dJ3/VldX1/h+NCL6E42aPE/Dn3lD+QVFyjpwXP/oPUwbtu72/VhESY0Rkpzr///bu/fvquo7j8P/T6YWW6dDqdZYEKMiXmgrQtXITZsYQaWoCKkUBYRWUcRry4x4wYLVKoxWqUCl4qUtkTIFRG4GFFDkTqCBJO/5QTn1QOIA405En2et80P2+easD2etsw7rlf3du25CHp3zUpLkrYZV+W6fmrS2tpWe//HVt+WZF17tqvGA47RoydsdRsixk2d0wUTAF6Fp/4F2I+SVwydm5qff40nyp9eXpeePbujs8YD/p8MRcstH27t6FOAENTcfTLezh2TVmsbSsSkPPJWbJzzchVNB8URIjsn+A83p0bc2r//1H0mS2c8vzKVXjytbM2r8g/nVg7O7YjzgBHxehDytalgq+41IvyH1mfHUC10wHXCiOoqQZ14yvOwMi/c2bU1FZXX2H7AlG04mhyPk9y++LlUDRuW6sfdm7Yb3u3os4Disb9x81HfwnHmL0v+acZ/zW3DyEyG/hvY1HciOXXvafXS05fLWSb/O0JG/TFvbJ2c+PjrnpVxx3YSyNfVTZmT83TMLnx9o3z+bD3b42d69t+mo9R1FyMb3t2Zd4+Z8sGVbFrzWkDMvGZ4nn32lM/4JQDtaWlo7/Gzv2LWn9N18WEcR8jvnXZM3G1aWft760Y5UVFa7lAp8SRzrZ721tS1v/2NNNm/9OO+u35T6KTPS68c3uMYrnERWvvteKiqry77D581/PX2vHN2FU0HxRMivoTvueSxVA0a1+3josblHrb9z+qz0G1KfPZ+JGLOfX5jLasaXrRs1/kHXloIuNPflJR1+tmtHTz1qfUcR8kgPPz43g66/s4iRgWPw3qatHX62qwaMyr6mA2XrP+9MyD+/tbzsdZ0JCV8ex/tZP6ylpTWn9h6av/19dSdPDJyow2dCNjcfLB2bM2/RUbsN4atGhKRDh1pacuukX+fyujvKAmSSvNmwMj0uqC37y80Ph/48c+Yt6uQpgRN1rBFy6iNP54bbpnfCRMAXoaMIecV1E/LY0y+Xfl64pCGV/UZ09njAF2z/geac0nOQLdlwEmluPphuvQaXXRNy8vRZGTX+wS6cCoonQtKuQy0tGTJySmpHT/10m/aBNO0/UDpborn5YM68ZHh+M+uF7NnblP9+5Y18+5xh+ejjnV08OXCsOoqQ4++emTcbVmbb9l1Z/Oby9LigNq8sXtoFEwInoqMI+cQz83PuwJuyYeOWfLhtZy6vm5CJ057ooimBE7Xgz0vz+O/mZ8PGLdm89ePcPOHh9BtSX3bDSODLb/jYabn+5/fl4x27s/Ld93LGRXX+z81XnghJu3bvbUpFZfVRj/Mvv7m0ZtmKtek3pD7f7Dk45w4clQWvNXThxMCx2vzh9vToW5vvnHdNvvGDq9Kjb20m3fdk6fmJ057I2f1H5tTeQ3PuwJvy2+cWdOG0wPHoN6Q+PfrWpqKyOj0uqC3b1tXS0poJ9z6R7n1qclrVsPxs/ANp2t/+9k7gy2vZirX5ybV3pHufmnTvU5NrR9+TD7Zs6+qxgOO0bfuu/PTmu3Nq76E5/cK6PDDzua4eCQonQgIAAAAAhRIhAQAAAIBCiZAAAAAAQKFESAAAAACgUCIkAAAAAFAoERIAAAAAKJQICQAAAAAUSoQEAAAAAAolQgIAAAAAhRIhAQAAAIBCiZAAAAAAQKFESAAAAACgUCIkAAAAAFAoERIAAAAAKJQICQAAAAAUSoQEAAAAAAolQgIAAAAAhRIhAQAAAIBCiZAAAAAAQKFESACATvLsi4tTUVldenTvU5O6MfdmfePmTp3jtKpheWnRX8qOTbn/qZz1wxGdOgcAAF8fIiQAQCd59sXF6dZrcN5dvynvrN2YhUsa8qNht+XMS4Zn7779nTZHexFyxeoN+ePiv3XaDAAAfL2IkAAAneTZFxen29lDyo4tX7U+FZXVeWPpitKxXbv3Zcydv0mPvrX5znnXZND1d2b1uk2l5x//3fycO/CmfPucYenepybVIyZl+cp1Za+7et2mDPvZL3Na1bCcfmFd6m69J2vWv5++1beWnY35jR9clSR55Il5uWTw2NLvHzx4KJOnz8r3L74u3c4ekgE147NsxdrS80uXr05FZXUWv7k8A2tvz6m9h+biQWPK1hxqacnUR55Or0tvzLd6D815P7kpdz00+4t5MwEAOKmIkAAAnaS9CHk45q1a05gkaWlpzaVXj8sNt03P6nWbsvGDDzN5+qyc9cMRadp/IEny2+cWZPbcRVm+an3e/seaDB87LWdcVJeWltYkyQdbtqV7n5rcdPtDWbZibd7424rUT5mROfMWZeeuvfn2OcPy9Lw/Zdv2Xdm2fVeSoyPkxGlPpNePb8jCJQ1ZtaYxt099LP9+7tX56OOdZXP37j8ys+cuysIlDbly+MT07j8ybW1tSZIHZj6X3v1HZunyd/Pepq35w8K3Ujt6arFvMgAAX0oiJABAJ/lshGxpac26xs0ZWHt7fnLtHaU1r76xLKdfWJeDh1pKx9ra2tKjb20Wv7m83ddtfH9rKiqrs+7Ta0tOnj4r519+c1pb28rWHf65ve3Yn42Qe/Y25Zs9B5dtz25ra8tFV43J1EeeTvKvCLlz197Smr8ueycVldWlUDmi/r5cc9OvjuMdAgDgq0qEBADoJIdvTPONH1xV2g49ov6+7Nqzr7TmgZnPlW2X/uxjzrxFSZIPt+3Mrx6cnYG1t6dqwKh878JrU1FZnf95Z32S5KrrJ2Xs5BkdzvF/RcjlK9elorI6H27bWbamfsqMXDv6niT/ipCHz85Mkg0bt5TF0AWvNeSUnoPS54pb8ou7Z2buy0uy/0DzCb57AACczERIAIBOcuSNaS69elwuHjQme/Y2ldbc/1+/z7kDR3X4Gk37D6Sy34hcXndHXlr0Vlas3lCKhocjZPWISamf8sVHyLGTPz9CHnlGZpJs/ODD/OdTL+aG2+5P9z41uWTw2LKzPAEA+HoQIQEAOsmR14TcvnNPqgaMypXDJ6a5+WCS5JXFS/NvZ11VFvIOa2lpLW15/njH7tLxHbv2lEXI26c+lguuuKXD7djfu/DazJv/etlzn42Qu/c25ZSegzL/1b+Wnu9oO/bnRchDLeWxcdWaxlRUVmfN+veP5e0CAOArRIQEAOgk7d2YpvH9rTn9wroMHzstra1tOdTSkn5D6lM1YFT+sPCtrFrTmAWvNaRuzL35y9ursvWjHTml56DcOX1W3mpYlWdeeDX9rxlXFiHXNW7Ot3oPzS0THsnfV6zNshVrM+m+J/PMC68mSa4cPjGDrr8zq9Y0lq4zeeSNaSbc+/hRN6Y5rWpY6ezIY4mQdbfek4cem5t31m7Me5u25q6HZqdH39oc+Kct2QAAXzciJABAJ2kvQibJ8lXrc1rVsIy769Ekya7d+zLurkdT2W9EuvUanN79R2bs5BmlO1nPm/96el16Y07tPTSX1YzPvD++XhYhk0/OOhwyckq+26cmPS6oTe3oqVm74ZMzEFes3pCLB41Jt16DUzXgk63fR0bI5uaDmXTfkznjorp06zU4l9WMz9Ll75aeP5YI+fzLr6X/T3+R7n1q8h/n/zSDb5ycFas3fEHvJgAAJxMREgAAAAAolAgJAAAAABRKhAQAAAAACiVCAgAAAACFEiEBAAAAgEKJkAAAAABAoURIAAAAAKBQIiQAAAAAUCgREgAAAAAolAgJAAAAABRKhAQAAAAACiVCAgAAAACFEiEBAAAAgEKJkAAAAABAoURIAAAAAKBQIiQAAAAAUCgREgAAAAAolAgJAAAAABRKhAQAAAAACiVCAgAAAACFEiEBAAAAgEKJkAAAAABAoURIAAAAAKBQIiQAAAAAUCgREgAAAAAolAgJAAAAABRKhAQAAAAACiVCAgAAAACFEiEBAAAAgEKJkAAAAABAoURIAAAAAKBQIiQAAAAAUCgREgAAAAAolAgJAAAAABRKhAQAAAAACiVCAgAAAACFEiEBAAAAgEKJkAAAAABAof4XXFANDQ+h4CAAAAAASUVORK5CYII=", "text/html": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "result.plot(plotter, 5, variables=['FBP', 'G6PDH2r', 'PGL', 'PGK'])" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Downloading Map from https://escher.github.io/1-0-0/6/maps/Escherichia%20coli/iJO1366.Central%20metabolism.json\n" ] }, { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "8c09229a8e6e4a9a8dceb50e96eecd32", "version_major": 2, "version_minor": 0 }, "text/plain": [ "Builder(reaction_data={'ACALD': -39.11622412349536, 'ACKr': 115.52882712367142, 'ACt2r': 115.52882712367142, '…" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "result.display_on_map(2, map_name=\"iJO1366.Central metabolism\")" ] } ], "metadata": { "anaconda-cloud": {}, "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.9.2" } }, "nbformat": 4, "nbformat_minor": 4 }