Skip to content

Instantly share code, notes, and snippets.

@mr-eyes
Last active September 23, 2018 13:21
Show Gist options
  • Save mr-eyes/1b22d47a44f747deab3f40ee38e5ab08 to your computer and use it in GitHub Desktop.
Save mr-eyes/1b22d47a44f747deab3f40ee38e5ab08 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import matplotlib.pyplot as plt\n",
"import tqdm\n",
"import gc\n",
"import json\n",
"from rna_graph import *"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Create Graph Object"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"RNA = RNA_Graph(\"test3.gfa\",\"t_test.tsv\")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Coloring low complex contigs"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"color_map = []\n",
"for node in RNA.G:\n",
" if RNA.maxDustWindow(RNA.nodes[node][-25:], 25, 25) > 1:\n",
" color_map.append('blue')\n",
" elif RNA.maxDustWindow(RNA.nodes[node][:25], 25, 25) > 1:\n",
" color_map.append('blue')\n",
" else:\n",
" color_map.append('red')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Drawing"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD8CAYAAAB+UHOxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xl8TPf6wPHPTGISk8SSSBBB7VGx01C1VO1RO7VW0Vv7rRbVlturuvvRXleraK2tWqpV+1JLxVakoYidxJJQskuzZ76/P065lGwyyZlknvfrNS9Zzjl5xsx8n/PdDUophRBCCLtj1DsAIYQQ+pAEIIQQdkoSgBBC2ClJAEIIYackAQghhJ2SBCCEEHZKEoAQQtgpSQBCCGGnJAEIIYSdkgQghBB2ShKAEELYKUkAQghhpyQBCCGEnZIEIIQQdkoSgBBC2ClJAEIIYackAQghhJ2SBCCEEHbKUe8AhCgMUlLg5EkIDoY//gCDAcqVg8aNoU4dMJn0jlCI3JMEIEQWzp2D2bNhxQpwdIT0dEhK0n5nNoPRCBYLDBsGr70GVavqG68QuWGQTeGFeFhqKrz7Lnz2GaSlaQV/VooV0xLE1KkwZYr2tRC2ThKAEH8TEwPPPgsXLkBiYu7OdXGBunVhxw5wc8uf+ISwFkkAQtznzz+hWTM4f16rBTwOJyeoVw/27dO+FsJWySggIe7z+utw8eLjF/6gdRiHhMC//mW9uITID1IDEOIvBw5A+/b/6+TNq+LF4eBBaNDAOtcTwtqkBiDEX955J+vC32SCr7+GsDCIj4djx6BTp8yPT06GGTOsHqYQViMJQAjg2jXtbj0rjo7aca1bQ8mSMG0arFkDlSs/+nilYMsWiIy0frxCWIMkACGAPXuyH7qZmKgNDb1yRSvcN2+G0FBtMlhmnJxg/37rxiqEtUgCEAKt/T8hIXfneHlBzZpah29mEhLg8OG8xSZEfpEEIATasM/ccHTUZgcvW6bNFs6MxZL174XQkyQAIYCMjJwfazDAN99oQ0XHjcv++OxmEQuhF0kAQqA15+TUokVQtiz07p2zwj019RozZsygc+fOnD9/noyMDFJSUh4/WCGsROYBCLsXGhrKiBGn+OWXDiiV9dTdL7/UxvW3a6fNGs5eAo6Ok8jIWIhSCl9fXy5fvkynTp1Yv369VeIX4nFJDUDYnejoaGbOnEnz5s1xcXGhatWqBAd/gdGYdTtQpUowapSWAG7ehDt3tMfAgVmdZSA9fR9377POnj1Leno69evXt94TEuIxSQ1AFHmpqamsWLGClStXEhQURExMDG5ubjRo0IC+ffsybNgwXFxcqVQJrl+37t/28blDVFRZkpKSMBgM3P9xc3d3x9/fn4EDB9KvXz9MsqmAKGBSAxBFjsViYdOmTfTt25fy5cvj7OzMmDFjiIyMZNy4cVy7do34+HgCAwMZP348rq6uGAwwaVIGjo7JVovDxQU+/NCNwMBASpcuzbBhw3B2dsZsNnPp0iVGjhxJeHg4I0aMwNnZGR8fHwYMGMD27duxWCxWi0OIzEgNQBQJR48eZeHChezevZuwsDAMBgPVqlWjQ4cOjBw5Ej8/vyzP37JlC/37DyYh4RBQA6Xydm/k4AANG2pzAIxGSE5OxsnJiY0bNxIcHMz06dMfOD44OJiFCxeyc+dOwsLCAKhSpQrt2rVj5MiRNJAFhUQ+kAQgCqXQ0FDmz5/Pli1bOH/+PGlpaVSoUIHWrVszfPhw2rRpg9GYfSGekJBAt27d+OWXX+jVqxdTp66kZctiOezgzZybm7ZWULVquT/XYrGwfft2li1bxr59+7hx4wYmk4latWrRtWtXRo4cSaVKlfIW4F9iY2MpVaqUVa4lCiElRAH4448/1LRp09S4cePU3r17lVJKWSyWHJ8fFRWlPvnkE9WsWTNlNpsVoMqUKaMCAgLUd999p9LS0nId0/z585WTk5Py8vJSv/76672f79qllNmslLbgQ+4eBoNSrq5KHTyY63AylZKSopYtW6Y6duyo3N3dFaBcXV1VixYt1OzZs1VcXNxjXffDDz9UjRo1UrVr11abN2+2XsCi0JAEIPJdWlqaevfdd9WLL76oFixYoBo1aqROnjyZ5TkpKSlq8eLFqn379qp06dIKUG5ubqply5bqv//9r7pz585jx3PlyhXl5+enjEajeu2111RGRsZDxwQFKVWpUu4SgdmsVI0aSmXz1PIsJiZGzZo1Sz399NPKxcVFAcrd3V116tRJffPNNyolJSXba6xbt07Vq1dPKaXUjh07VNeuXfM3aGGTJAGIAlGpUiWVkJCglFJq9uzZatKkSermzZv3fm+xWNT+/ftVnz59VLly5ZTBYFDOzs6qYcOG6l//+pe6du2aVeJ48803lYODg/L19VWXL1/O8tikJKWmTVOqRAml3Nwyv+N3c1OqVCml3n9fqdRUq4SZK2FhYeqtt95SdevWVSaTSRkMBuXt7a369++vgoKCHqppxcTEqNdff13NmzdPKaVUcHCw6tu37wOvh7APMgpI5Lv4+Hj8/f05efIkAK1btyYpKYnjx4/fOyYtLY29e/dy4sQJ+vTpw4kTJ0hKSiI4OJgZM2bg4+OTpxiCg4Px9vZm9uzZzJ49mzNnzlClSpUsz3F2hvfeg9u3YfFiGD1a2++3fHnw9ob69WHsWG1ZiFu3tA3hixXLU5iPpXLlynz44YecOHGClJQUjhw5wvPPP8/Ro0c5fPjwQyOKQkNDiY+Pp1WrVgCkpKRQvnx5rl69WvDBC11lswCuEHlnMBioUqUKZ86coVmzZnh7e+Pl5cXZs2fp2LEjACaTibfeeou3337bqn87PT2dF198kVWrVvHMM89w+vTpXHd6mkzQp4/2KAyaNGlCkyZNAFBKYTAYHvj9zZs3SUlJoXbt2ve+T0xMpGrVqgUeq9CXJACR71xcXKhUqRK///47AKVKlSI9Pf2hgvjvBVVebdiwgcGDB6OUYt26dXTv3t2q1y8MHvV/GhsbS2RkJEajkcjISC5evIi3tzceHh46RCj0JE1Atu70afj8c+jXD5o00QaXd+igtU3s3p27ZSzzWUpKCosXL6ZDhw5MnDiRP/8aS2k0GmndujWHDh0iIiKC4sWLc/DgQby9vfMljvj4eNq0aUOPHj3o0qULMTExdln4Z6Zs2bKUKFGCpKQkfvrpJ4KCgujZsyeg/d+NGjWKzp078+2335KamqpztPdJT4cdO2D6dG0xpoYNtc/EgAHaIk25XdNbyDBQm7V5s1INGihVvLhSzs4P9z46Omq9j56eSn3yiVI5GPlhbRkZGWrjxo0PdNw6OTmphg0bqs8+++yh0TUTJkxQ/fv3V/7+/qpbt27q6tWrVo9p7ty5ymQyqXLlyqkjR45Y/fpFQVpampoyZYry8fFRvXv3Vvv377/3O4vFot5///1Hdihv27btkSOm8l1SklIzZijl7q695x0cHv48FC+uPZ56SqmdOws+xkJKEoCtiY5Wqlev3I8/rFZNqWPH8j28I0eOqJdffllVrVpVGY1G5eDgoGrUqKHGjRuX7dDO1NRUFRgYqDZs2PBY4/azEhYWpnx9fZXRaFSTJ0/O1RwDexYfH5/l748ePapGjhypqlWrdu/1rl69uho1apQ6VgDvN/Xrr0pVrKgV7rn5PAwerFQ2z01IArAtERFKVa6slMn0eLOQzGalduywakiXL19WkydPVn5+fvfuCH18fNTAgQPVrl279LkjvE9GRoaaNGmSMhqN6sknn1RhYWG6xlOUZWRkqC1btqh+/fqp8uXL36vx1atXT7399tvqypUr1v2D69c//ow8JyelatVSKjLSujEVMZIAbEVCglJVq2pNO4/zhr8/CQQFPXYYd2fc+vv735tx6+HhkacZt/nlyJEjqly5cspkMqkvvvhC73DsTkpKilq6dOkDM5RdXFzuzVCOjY19/IsHBj5+4X/3UayYUnXqKJWcbL0nXcRIArAVr7zy6Lb+x3lUrqy1m97n8uXLasaMGQ/92aSkJPX1119nOuP2cZcZyE8pKSmqT58+ymAwqDZt2thkjPbIGjOUlVJK3bmjVNmy1vksmM1KvfFG/j7xQkwSgC347bfctXFm9yheXKl//1tFR0erZs2aqVdffVWVKlVKOTg4qPDwcLVx40bVu3fvhzpurTnjNr/8+OOPytXVVbm5uamNGzfqHY7IQlhYmHrzzTdz36H86qvWuxm6+3k4c6Zgn3whIQnAFvTtq5TRmPWb+JtvtD6CuDilzp1TasSILI9PL1FCVa1YURkMBmUwGBRw7+u7HbdjxoxRv//+u97PPkdiYmJUixYtlMFgUAMHDrSppiiRMznqUE5IyL7pZ+xYpY4e1Zp2lizJPgE4Omo1bPEQWQ5abwkJ4OkJydlsRPLkk3DxIqSmQq1a8MsvEBAAwcGPPDweGAJs+NvP27Rpw65du3K0VLKtmDNnDpMnT8bT05ONGzfSqFEjvUMSeXR3yeslS5awf/9+bt68iclkYqKXF+/cvIlTWlrmJ/fsCRYLdOwIxYvDsGHZ/0GzGWJitGnd4p7CUwoUVcHBOXtTnj6tFf7wv3ubLBabdwHa/jUL1GAw4OXlhdFoJCwsrNAU/qGhodSqVYuJEycyceJEwsPDpfAvIoxGI507d2bNmjVERESQnJzMggULaJ2RkXXhD7BuHaxfD1FROf+DDg4QEpK3oIugwlESFGXBwZCSkrNjv/gC/vwTzp2DGzdgy5ZMD3UABlepwtixY6lbty4LFy4kLi6OvXv3WifufGSxWHjttdeoXr06Tk5OhIWF8dFHH+kdlshHJpOJoUOH0sHNLX/+gMWSaW3ZnslaQHr744+cJ4CxY2H8eGjeHNq0yfY8D4uFzz///IGfubq6PmagBePQoUP06NGDuLg45s2bx8iRI/UOSRSk6Oj8uW5ysrasq3iA1AD0ltsF0CwWOHAAfHy09YmteW0dpaam0rt3b1q0aIGfnx+RkZFS+Nuj/HrPGgza5sziAfI/ojcfH60jK7ccHbPfcDafFluztrVr1+Lu7s7OnTvZvHkzu3btsvmaisgnZcvmz3WdnQvN56EgSQLQW+PG2e8i4ukJL7wALi7aXUyHDtoKiLt2ZX6O0QitW1s3ViuLjo7m6aefpl+/fvTq1YuoqCg6d+6sd1hCT888k/0xDg7g5KT9e//X2WncOO/xFTGSAPRWv762zG1WlNKae65f14ayzZoFEybAxo2Zn+PiAm3bWjdWK5o9ezblypXj2rVrHDt2jOXLl+PoKF1Sdq9DB8iuI3jaNK1N/623YMgQ7etp07I+x8EBata0XpxFhMwDsAWjR8PXX2efCHKjfHktYdhYu+eFCxcICAjg8uXLvPXWW7z33nt6hyRsSVoaeHlBbKz1runkBK+9BjKS7CG2VTrYq4kTrbuZrIuLdkdkQ4W/xWJh/Pjx+Pr6YjabuXr1qhT+4mHFinG2e3f+tOY1HRxg3DhrXrHIsJ0Swp5Vr67tKG425/1aRiP4+sKoUXm/lpXs37+fcuXK8dVXX7Fw4UKOHz+eb7uBicIrISGBtm3bUnfZMmJcXVHWGBHk4gKffAIVKuT9WkWQJABbMWUKNGigjVbIixIl4PvvbeLuPyUlhR49etCqVSsaNGhAdHQ0I0aM0DssYYO+/PJLypQpQ0hICPt//RWfI0cw5HUkWPHi0KIFjBljnSCLIP1LCaFxdNT2O23Y8PFqAsWKgbs77NsHVapYP75cWr16Ne7u7uzdu5dt27axY8cOzNao4Ygi5erVq/j5+TFu3DjGjBnDjRs38Pf3h9q1Yc8eKFlS+2zkUpLRSFLjxtqSETZwM2Sr5H/Glri4wN69Wp9A8eI5G9p297w2bbS1Tvz88jXE7ERGRuLv78+AAQPo168fUVFRdOjQQdeYhG2aMmUKVapUISMjg4sXL/Lpp58+uE5V48Zw8iT4+2vv8RxQDg4kAtMtFsocPcq0998nKSkpf55AUaDvYqQiUyEh2r6mzs5KlSyp7W50dxk4g0HbHNvZWalmzZTasEEpG9gD9+OPP1aOjo6qUqVK2e4PLOxXUFCQKl++vCpWrJj6z3/+k/0JFotSa9cq1bixtra/m5v2Gbj7eTCZlCpRQqnixdWfAwaoOiaTAu49Vq5cmf9PqpCSYaC2Lj4eDh+GoCA4f15bCqJsWXjqKe3OqGJFvSPk3LlzdOnShStXrjBt2jSmT5+ud0jCBqWnpzNkyBBWr15NixYt2LhxI6VKlcrdRUJD4cgROHpUW9vHwUFbHr1JE/D3x2I24+TkRPpfQ6qLFStG3759+eyzz/Dy8sqHZ1W4SQIQj81isTB27FgWLlxI/fr12bJlC+XKldM7LGGDNmzYwODBg1FK8c0339CjR498+1u+vr5ERkYS9ddy0Q4ODjz99NMEBgbm298srKQPQGTrUfcIgYGBeHl5sXTpUhYvXkxwcLAU/uIh8fHxtG7dmh49etCpUydiYmLytfAHOHPmDO3atbv3va+vL2vXrs3Xv1lYSQIQDzlz5gz//Oc/mTlzJnFxcRj+Nh47KSmJL7/8kqZNmxIVFcXQoUN1ilTYss8//xxPT0/Onz/P4cOHWbNmTYEs92EwGPDz86NUqVL06dOHc+fOFZpNkAqaNAGJB1y+fJmOHTsybtw4Dh8+jIeHB23btqVnz54PHJeeni5r94hHunLlCp06deL8+fNMnDiRjz/+uMAL4OTkZCwWC2azmSeeeIIyZcoQFBRUoDEUBpIWxQPOnj1Ls2bNePXVV5k3bx6NGjVi69atnDp16oHjpPAXf2exWJg8eTJVq1bFYDBw+fJlZs6cqcvdt7Oz8715J1u3buXYsWMsWrSowOOwdZIAxAP8/PwICQnh0qVLlCpVilatWlGxYkVpQxVZOnr0KBUqVOC///0vc+fO5fTp01SuXFnvsACoXbs2Y8eOZezYscTHx+sdjk2RBCAeULZsWbp27cq3334LQLVq1ahfvz5xcXHcuXNH5+iErUlNTaVv3774+/vj6+vL7du3GWODSy/85z//oXTp0nTt2lXvUGyKJADxACcnJ9q3b8/169f58ccfAWjZsiV79uwhOTlZ5+iELVm3bh0eHh5s376dDRs2sGfPHkqUKKF3WI9kNBrZvHkz+/fvZ82aNXqHYzOkIVfck5GRgYODA82bN+f69et88MEHVKpUiYsXL1KuXDnS0tL0DlHYgNjYWLp27crBgwcZMGAAy5YtKxR9Qo0aNWLIkCEMGzaMbt264ZzXhReLABkFJEhMTGTw4MHMmzfvgbH88+fPJzg4mGPHjvHFF1/w1FNP6RilsAVz5sxh8uTJeHp6snHjRho1aqR3SLlisVgoU6YM/v7+bN26Ve9wdCcJwM4tW7aMUaNGYTab2b17N/Xr13/g92lpaRSz5mY1olAKDQ2lU6dOXLp0icmTJ/NRId5dKzAwkDZt2rB161Y6duyodzi6kj4AO3Xz5k0aNWrE8OHDGTp0KLdv336o8Aek8LdzFouFCRMmUL16dUwmE6GhoYW68Ado1aoVPXr0oF+/fvfWDLJXkgDs0IwZM/Dx8SEuLo7Tp08zf/58mSkpHnLo0CHKly/P/PnzmTdvHidPnqSiDSw+aA2rVq3CYrEwcOBAvUPRle333AirOXXqFAEBAURERPD+++/z5ptv6h2SsEGpqan079+fn376iWeffZb169fjmtfduWyMyWTi+++/JyIi4t7gB3skfQB2wGKx8PLLL7N06VKaNm3K5s2bKVOmjN5hCRu0du1aXnrpJRwcHFi5ciVdunTRO6R8lZKSgpOTk95h6Ebq/UXczz//jIeHB6tXr2bFihUcPnxYCn/xkOjoaJ5++mn69etHr169iIqKKvKFP/DIwj81NZXY2FgCAwO5ceOGDlEVHGkCKqISExPp0aMHO3fu5Pnnn2f16tUy7lk80uzZs3nrrbcoW7YswcHBNGjQQO+QdJGRkUFQUBA///wz69ato3Llypw/f54NGzZQtWpVvcPLF1IDKIIWLVqEh4cHx44dIzAwkPXr10vhLx5y4cIFatasyZQpU5gyZQrXrl2z28IftJnNixYtIikpiSVLlvDjjz8yZMgQPvvsM71DyzdSAyhCIiIi6Ny5M6dOnWLMmDHMmTNHRveIh1gslnurvfr5+XH16lW8vb31DktXixYtYt68eUyaNIkBAwYA2pLSFy9e5IknnkAp9dC+GEWBJIAi4p133uHDDz+kSpUqnDt3jurVq+sdkrBB+/fvp1evXsTHx7NgwQJefvllvUPSXVpaGr/++iuff/45zZs3B7Q9DQIDAzEajYwcObJIFv4gCaDQO3HiBAEBAfzxxx988sknTJw4Ue+QhA1KTk6mf//+bNiwgXbt2vHTTz/dWy/f3hUrVozffvuN559/noSEBPbt28fx48cJDQ2lXbt2lClTpsjWAGQYaCFlsVgYNmwY33zzDc2aNWPTpk24u7vrHZawQatXr2b48OGYTCZWr15Nhw4d9A7J5uzbt48ZM2Zw+/ZtWrdujZubGz179qRx48Z6h5avJAEUQlu3bmXAgAGkp6ezdOlS+vTpo3dIwgZFRkYSEBDA0aNHGTp0KIsWLZI+oSxcvXoVR0dHnJ2dKVmypF1MDpN3QyGSkJBAu3btCAgI4LnnniM6OloKf/FIn3zyCeXLl+fmzZucOHGCJUuWSOGfjUqVKuHt7Y27uzsODg6kpqbqHVK+k3dEIbFw4ULKlCnDyZMnOXDgAD/88AMmk0nvsISNOXfuHNWqVWPq1KlMnTqVK1eu4Ofnp3dYhY7FYmHHjh0sXLhQ71DylSQAG3f9+nXq1avH6NGjGTVqFDdu3Lg3UkGIuywWC6NHj+bJJ5+kZMmSXL9+nenTp+sdVqFlNBo5dOgQY8eOJTIyUu9w8o30Adiwt99+m5kzZ1KtWjW2bt1aZGcjirzZu3cvvXv35s8//2T+/PkMHTpU75CKjMqVK+Pp6UlQUJDeoeQLqQHYoODgYHx8fJg1axazZs3i3LlzUviLhyQnJ9O1a1eeffZZmjZtSlRUlBT+VrZt2zaOHTvGkiVL9A4lX0gCsCHp6ekMGjSIJk2a8MQTT3Dr1i0mTJigd1jCBq1YsQJ3d3cOHjzIzz//zNatW2Vcfz6oXbs2Y8aMYfTo0cTHx+sdjtVJE5CN2LRpE4MGDcJisbBs2TJ69eqld0jCBt26dYsuXboQHBzMiBEjWLBggYzuyWcWiwVvb29q1arF3r179Q7HquSdo7P4+HieffZZunXrRseOHYmKipLCXzzSBx98QIUKFYiKiuLkyZN89dVXUvgXAKPRyKZNm9i3bx9r167VOxyrknePjubNm4enpydnz57l119/Zc2aNTK0Uzzk9OnTVKlShenTp/Pvf/+b0NBQ6tSpo3dYdqVJkyYMHjyYoUOHkpycrHc4ViMJQAdXr16lTp06jB8/nvHjxxMeHs5TTz2ld1jCxlgsFv7xj3/g5+eHh4cH4eHhTJs2Te+w7NbSpUtxcnKid+/eeodiNZIACtiUKVOoUqUKFouFy5cvM2vWLKnGi4fs3r0bT09PVqxYwfLlywkKCsLLy0vvsOya0Wjkxx9/ZOvWrWzfvl3vcKxCOoELSFBQEN26dSMqKopZs2Yxfvx4vUMSNigxMZE+ffqwbds2unTpwtq1a2UzHxvTs2dPdu/eTVRUFI6OhXtBZbn1zGfp6en079+fp556iho1anD79m0p/MUjLVu2DA8PDw4fPszu3bvZtGmTFP42aPXq1VgsFgYPHqx3KHkmCSAfrV+/Hnd3d7Zs2cJPP/3E3r17KVGihN5hCRtz8+ZNGjVqxPDhwxk6dCi3b9+mTZs2eoclMmEymVi5ciVr1qzhwIEDeoeTJ5IA8kFsbCytWrWiZ8+eBAQEEB0dTbdu3fQOS9igGTNm4OPjQ1xcHKdPn2b+/PnSJ1QIdO3aleeee47u3btjsVj0Duex2XcfgFJw+jQcPQrHjkFsLDg7g58fNG4MTZpALodlzp07l0mTJuHh4cGGDRto0qRJPgUvCrNTp04REBBAREQE7733Hm+++abeIYlcSkxMpEyZMgwaNIivvvpK+2FKilae/PYbhIRo35curZUnTZtCrVpgSzuLKXuUkqLUl18qVaWKUi4uSrm6KqWlA+3h7KxUiRJKlSyp1JQpSt269cjLHD9+XPXo0UOlpaWpsLAw5evrq4xGo3rjjTcK+AmJwiIjI0MNGzZMGQwG9dRTT6nbt2/rHZLIg++++04ZDAYVsnOnUq+/rpSbm1Z2ODs/WKa4uiplNitVs6ZSixcrlZamd+hKKaXsLwEcO6ZU9epawX//C5TZw8lJe0FXr753ifDwcLVt2zbl5+enHB0d1TPPPKOMRqOqU6eOunLlio5PTtiyHTt2qFKlSimz2ay+++47vcMR1mCxqBnVqqk7BoNSJlPOyhQXF6Xq1FEqJETv6O0sAfzwg5aFDYacvVD3P8xmpcaPV8piUc8++6wyGo3KyclJAQpQ7733nt7PTtiohIQE1b59e2UwGFS3bt1UUlKS3iEJa8jIUOqll5TFbM59eWIwaGXK1q26PgX76W3avh0GD4bERO0lyK3ERFi0iAu9e3Po0CEsFgspKSkAGAwGdu3aZeWARVGwaNEiPDw8OHbsGIGBgaxfv16GdhYVY8bAmjUYEhNzf65SWpnSuzfs22f92HLIPjqBIyOhenWIi8vzpRKBTsDdl+yJJ56gcePGNG/enIkTJ+b5+qJoiIiIoHPnzpw6dYoxY8YwZ84cGd1TlGzeDP36aYV4Xnl6wsWLoMMQcft4R44eDUlJVrmUGVgNmIsVw9HREX9/f9auXSuFv7jnnXfeoVKlSiQmJnL27Fnmzp0rhX9RkpgIL75oncIf4M4d0Kn8KPrvyogI2LQJUlOzP7Z6dS1RfPNNlod5mc38sWABqamprFq1ykqBisLm75XnEydOULFiRT7++GM+/vhjLly4QI0aNXSKTuSblSu14Z1ZKV0afvwREhIgLAwGDMj82ORk+PZbiImxapg5UfQTwFdf5bzN/4svtDG82XBITMR13jwMtjSeV+S7M2dccjEGAAAW/UlEQVTO8M9//pOZM2cSFxf3wOufkpLC3r17qVixIjdv3mTSpEk6Riry1axZ8OefWR/zxRfaTWfZsjBoEHz5JTz5ZObHG42wbJl148yBot8H0LQp5GRD5xdegF69tIlh1avDkCFZH1+smPYmKFbMOnEKm3b58mU6duzIuHHjOHz4MB4eHrRt25aePXveOyY9Pb3QLw4msnHnDnh4QFpa5seYzdrdvJ8fXLig/Wz5cggPh7feyvy8tm2hgAeTFO0agFLabLzsuLnBjBnw+us5v7azc86uLYqEs2fP0qxZM1599VXmzZtHo0aN2Lp1K6dOnbp3jBT+duD4ca2Az0rNmpCe/r/CH+D33yG7TXyOH897fLlUtBNAcnLO2v7few8WLdIydE4ZjVr/grALfn5+hISEcOnSJUqVKkWrVq2oWLFikdsiUGQjPByyW/vH1RX+voF8XJx2o5mV2Ni8xfYYinYCUCr7dTfq14d27eCzzx7v+sIulC1blq5du/Ltt98CUK1aNerXr09cXBx37tzROTpRYHLymU9IeHhIZ4kSWvNRdtcu4DKlaCcAZ2ftTj0rbdrAE0/A1atw4wZMmqRNzvjtt6zPUwrKlLFWpMLGOTk50b59e65fv86PP/4IQMuWLdmzZ0+R2iNWZMPTM/sy5fx5cHTU+hLvql8/+yZjV9cCXyiuaDdaGo1Qo0bW//ELF8L9QzknTdISwujRWV87MRHq1bNKmMJ2ZWRk4ODgAEDz5s25fv06H3zwAZUqVeLixYuUK1eOtKw6BEXR0rBh9nOKEhO1IaAzZsDLL0ODBtC9Ozz9dNbn1a1rvThzqGgnANCad86d0zplHiUp6cEXNCFB6zuIjMz6ulWrQvHi1otT2JT09HSGDh3KP/7xD1q0aEGxvyb+DRgwgLi4OBYuXMixY8f44osv8Pb21jtcUVA8PLRaQHb9hWPGwOLFcOsWREVpN5SnT2d+vJMTdOhg3VhzoOgPAz13LmdZOzdcXGD2bBg50nrXFDZj06ZNDBo0CIvFwpo1a+jcufNDx6SlpVFMhgDbp5kzYfp065YpTk7achA+Pta7Zg4U7T4A0DZgaNwY/qrGW4WDgza5QxQp8fHxtGnThm7dutGhQweioqIeWfgDUvjbsxEjsu8HyA1HR62looALf7CHBACwdKmWYa3BbIYFC7QOG1FkzJs3D09PT86dO8ehQ4f4/vvvMeVyNzhhJzw8tBYAFxfrXM/ZWeuL1IF9JIBq1eDTT7OfwJGd4sWhY0dt1rAoEq5cuUKdOnUYP34848aNIzw8HH9/f73DErbulVegeXOt8M6L4sW1wl+nfiT7SACgtddPnfr4ScBshmee0UYMyRpARcIbb7xB1apVsVgsXL58mdmzZ8uqnSJnDAZYv15rXn7cwSBms9afkNVCcfnMvt7tb7+trfRZsmTON3s3GLQX+PXXYcuWXG8SL2xPUFAQ3t7ezJkzh//85z+cOXOGypUr6x2WKGzMZm3tntGjc5cEnJy01UK//x7Gjcu/+HLAvhIAaAu+XbigdeSYzdr07Efd9bm4aC9Up05w6JC2XISs9VKopaen079/f5566ilq1KjBrVu3GD9+vN5hicLMyUnrD9i7F557Tvv+UX0DRqM2G9jVFUaNgkuXoEuXgo/3b4r+MNCs/PknbN0Kv/6qPe7c0V7AunWhRQut8Jcx3kXC+vXrGfLXCq/ffvst3bp10zkiUSRdvQo7dsCBA3DqlLYWWcmSWn9B8+ZaH6INzR+y7wQgirzY2Fi6devG/v37eeGFF/jmm29k1U4h/mJ/TUDCbsydO5eyZcty8eJFjhw5wsqVK6XwF+I+kgBEkRMWFoavry8TJkxgwoQJRERE0KRJE73DEsLmSAIQRYbFYuH111+nWrVqODo6EhoayieffKJ3WELYLKkPiyLh8OHDdO/enZiYGD7//HNGZ7eaqxBCagCicEtNTaVPnz40b96cJ598ktu3b0vhL0QOSQIQhdYPP/yAu7s7P//8Mxs3bmT37t2U+PtOTEKITEkCEIVObGwsLVq0oG/fvvTs2ZOoqCgCAgL0DkuIQkcSgChUPv30U7y8vLhy5QrBwcEyrl+IPJAEIAqFS5cuUbNmTd544w0mT57M9evXadCggd5hCVGoSQIQNs1isfDPf/6TmjVr4uzsTFhYGB988IHeYQlRJEjdWdisAwcO0LNnT+Lj4/nyyy955ZVX9A5JiCJFagDC5qSmptKzZ09atmxJvXr1iIyMlMJfiHwgCUDYlDVr1lC6dGn27NnDli1b2LlzJ66y/aYQ+UISgLAJ0dHRNGvWjP79+9OnTx+ioqLo1KmT3mEJUaRJAhC6mzlzJmXLliUiIoLjx4+zbNkyHBwc9A5LiCJPOoGFbi5cuEDnzp0JCwtj6tSpvPvuu3qHJIRdkQQgCpzFYmH8+PHMnz+fevXqce3aNcqXL693WELYHUkAokDt27ePXr16kZCQwNdff82wYcP0DkkIuyV9AKJAJCcn8/zzz9O6dWsaN25MVFSUFP5C6EwSgMh3K1euxMPDgwMHDrBjxw62bduG2WzWOywh7J4kAJFvIiMjadq0KYMGDWLAgAFERkbSrl07vcMSQvxFEoDIFx999BHly5fn9u3bnDx5kq+//hqjUd5uQtgS6QQWVnXmzBkCAgK4evUq77zzDu+8847eIQkhMiG3ZMIqLBYLr7zyCnXq1MHd3Z3w8HAp/IWwcVIDEHn2yy+/0Lt3b5KSkli6dCkvvvii3iEJIXJAagDisSUnJ9OlSxfatm1Ls2bNiIyMlMJfiEJEEoB4LMuXL8fd3Z3Dhw+zc+dONm/eLEM7hShkJAGIXLl58yaNGzfmpZdeYsiQIdy+fZu2bdvqHZYQ4jFIAhA5NmPGDHx8fIiJiSEkJIQFCxbI0E4hCjHpBBbZCgkJISAggPDwcGbMmMHbb7+td0hCCCuQ2zeRKYvFwogRI6hbty5eXl6Eh4dL4S9EESI1APFIu3btom/fviQnJ/Ptt98ycOBAvUMSQliZ1ADEAxITE+nUqRPt27fnmWeeITo6Wgp/IYooSQDiniVLluDh4UFQUBC//PILGzZswNnZWe+whBD5RBKAICIigoYNG/Lyyy8zfPhwbt26RatWrfQOSwiRzyQB2Lnp06dTqVIl7ty5w9mzZ/niiy9kaKcQdkI6ge3UiRMn6Nq1Kzdu3ODDDz/kjTfe0DskIUQBkwRg65SC69fht9/g2jXIyAAPD2jYEHx9wTF3L6HFYmH48OEsX74cf39/jh07hoeHRz4FL4SVpaXB6dNw/DjExICDA1SuDE2agLe33tEVOpIAbFVMDCxaBJ9+CrGxUKwYpKRoCcFkAoNBSwaDB8Nrr2nJ4BHmzJlDixYtaNKkCdu3b+eFF14gPT2dVatW0a9fvwJ+UkI8ppMntc/CqlXaZ0EpSE3VPgdOTtrXnp4waRK89BKUKKF3xIWCQSml9A5C/M3338PLL0N6OiQmZn2so6P2gRg9Gj74AP4atXPixAlSU1N5+umn8fLyolatWuzZs4fu3buzevVqTCZTATwRIfLozz9h4kRYvlwr5DMysj7ebNYSwvLl0LVrwcRYiEkCsCUZGVrBv2ZN9gX/35nNUKECBAYSlpxMjRo1KF26NLdv3wagePHi/Pzzz7Ro0SIfAhciH1y7Bi1bwq1bkJSUu3PNZhg+HP77X62WIB5JmoBshVIwbBj88EPuC3/QzgkNhebNeb1yZdLT0+8V/gBpaWm4u7tbMWAh8tEff4C/v1b4Z3fX/yiJibB4sfa5+vxz68dXRMh4P1vx7bfw44+PV/jflZ5OxvXrDN+3796PDAYDFSpUoEGDBkRHR1shUCHymVLwwgtw+/bjFf53JSbCkiWwfr31YitipAnIFty6BdWrw507VrlcAjDCZGKLyURCQgIHDx6kefPmVrm2EPlu2TIYO1Zr/7eGUqW02nGpUta5XhEiNQBbMHeu1sGVlT17tHbQO3e0x9mzmR7qCizy8uLgwYNERUVJ4S8KD6Vg6tSsC/+7n4G7j/R0ra0/Mykp2og68RBJAHpLT9faKFNSsj923Dhwc9MemQz7vMs1Joa6KSnS7m9l27Ztw9fXl5o1a/LJJ5/oHU7Rs2cPxMVlfczdz4CbG5Qrp90Yff995scnJcHs2VpyEQ+QBKC3kJC8tXNmJiUFtm+3/nXtmMViYdy4cWzfvp2QkBBWrlzJ2SxqYiJr586dw2KxEBcXx40bN7QfbtmSu6af3r21JtT7+r0eKS4OwsIeO9aiShKA3n77DSyWnB370Udax9j+/dC6ddbHpqfDL7/kOTzxP0eOHKFGjRpUrlyZYsWK0b9/f9ZLB+Nje/LJJ6lWrRpjxoyhQoUKDBs2jJSdO3N3pz50qDbmPzsODtpnTTxAhoHq7fLlnN3xTJmiTYFPTYX+/WHjRmjQQDs/EzHBwXw1c6YVg7UvvXr1onr16ve+Dw8Pp2LFive+9/Hx4ciRIw+cExISwubNmwssxsLMYrEQFhbGtWvXUEqxdOlSPgbK5vQClSppN0IjRmR/bFKS1hEsHiAJQG9paTk77v6CZvlyGDAAunTJcozznZgYPv744zwGaL/q1av3QALIiZCQEPk/z6WM+5pAc9UkMWSIVhvOSdOOUlqtWDxAEoDePDy05Rxy++ZUKtsZjpXq1iX699/zEJy4X4UKFbh69eq9769fv06FChUeOKZfv36yxlIOGQwGXF1dcXZ2Jjo6GpPJhMFkgvj4nF3gxRchp8nWZJJhoI8gfQB6q19fm7aelZIloUMHbY0TBwcYOBBatYJt27I87aiDA9evX7disPatadOmXLx4kStXrpCamsqqVavo1q2b3mEVWkOHDmXVqlWsWLGCcePGERoaSpkOHXJ2cvPm2tInWY3+uV+xYlqTqXiATATTW1SU9kbOahhomTLa6AhfX23E0Nmz8K9/wc6dmZ6S6ODABCcnvkpMxNXVlQYNGtCnTx+GDRtGCVkp8bFt27aNV199FYvFwogRI3jzzTf1Dqlo+fJLbUXP7GbEz5+v3Ti9+GLOrmsyaSvsZnezZWckAdiC556D3bute02zGf74g+jUVL7++mvWrVvHiRMnSExMxMPDA39/fwYOHMgLL7yAYy73FBAi3+Tkhii3DAbo2VNbZ0s8QBKALfj5Z+0Naq2p705O2vC4BQse+lVoaCjz589ny5YtnD9/nrS0NLy9vWndujXDhw/n2WeflS0hhb7699fWxcrpAInsmM2waxc0a2ad6xUhkgBsgVIQEKC9SbNbEiInSpeGixchB7OAjx49ysKFC9m9ezdhYWEYDAaqVq1K+/btGTlyJPXq1ct7PELkxs2bULOmddbGcnKCHj20jWTEQyQB2Ipbt6B2bcjrip1ms7afQEBArk+1WCxs2bKFZcuWsX//fv744w9MJhO1a9fm+eef55VXXsHHxydv8QmRE6tWaeP787I6rsGg7RJ27pyMAMqEJABbcvKktgFGfPzjrVtSvLi2bd6oUVYJJzk5mRUrVrB69WqCgoKIiYmxWodybGwspeRDKbLyf/8H06c/XhIwGrVC/9AhrTYhHkkSgK25cEHrDwgNzfkb32TSCv9ly6B793wLLTo6OtMO5UGDBtGvX78cdSh/9NFHrF27lqSkJGbNmkWXLl3yLWZRyH33HYwcCcnJOZ8r4+KijZj74Qdtw3iROSVsT1qaUh9/rFSpUkq5uSml1QcefpjNSjk7KzVokFK3bhV4mJcuXVJvvPGG8vPzUyaTSRkMBlWhQgU1cOBAFRgYqCwWy0PnrFu3TtWrV08ppdSOHTtU165dCzpsUdhERCjVu7f2Xnd2zvzz4OamlIeHUnPmKJWRoXfUhYLUAGxZWpq25s+WLXDwINy4oS0cV6IENGkC7dppIyY8PPSOFHiwQ7lly5bMnTsXNze3e7+PjY3lvffeo3r16owePZpjx47x0UcfMXfuXMqWzfEKMMJe3bql1Qh274bgYK2T2GjUho22aKH1ewUEaJMlRY7IAHBbVqwY9OqlPQqBpk2b0rRpUwCUUhj+tlRFaGgo8fHxtGrVCoCUlBTKly/P1atXJQGI7Hl5wYQJ2kNYhQz4Fvni74U/wM2bN0lJSaF27dr3vk9MTKRq1aoFHZ4QAqkBiAIUGxtLZGQkRqORyMhILl68iLe3Nx4eHqSnp7N48WLWr1/P4MGD6du3r23NUA4J0TYdOXAArl3ThhhWrqw1PbRqBbVq6R2hELlmQ58wUdSVLVuWEiVKkJSUxE8//URQUNC9tXQMBgMODg5cvXqVl156iUGDBuk/Q1kpWLsWZsz4374Lfx+ZtXat1i/j66sNWZTF4UQhIp3AosCkp6czbdo0VqxYgb+/P6+99hotWrR45LFZzVAePXo0fn5++RtsRAQMHqztw5DTJTpcXLQNSpYs0dqrhbBxkgCELu7cufPACKGsZDdDedSoUXh7e1svuLNn4ZlntH1kc7tPg8mkLcVx6BBUqWK9mITIB5IARKHzqBnKbm5uNGjQgN69e+dtyeuICKhXT1uS43E/GkajVgM4eVJbylsIGyUJQBR6mc1QbtasGYMGDXpkh/KSJUto3rw5vr6+//uhUtC2rbbNYF63DzSZoHNn+OmnvF1HiHwkCUAUOY9a8rpChQq0atWK4cOH07JlS0qVKoWjoyNTp07lhx9+YOPGjZTduxeGD7festwuLloncadO1rmeEFYmCUAUeX/vUAZtoppSCqPRiMVioVy5coS5uuJ08WL2F3zhBfj3v6FSJW3p4pde0moNj9K8uTaLWwgbJAlA2BWLxcKQIUNYuXIl97/16wCHAZfsLtCuHXz9tZYEjhyB8uW1n0dEPPp4Z2c4cwaeeCLvwQthZTITWNgVo9FIeHg4BoPhgVFIrR8xc/mR3n1Xmxdw+LDWZxARkXnhD9pyHgcO5DFqIfKHJABhdxYvXsyVK1fYunUrTzzxBL169WJ8o0bZ3/0bjdoifJ6e2rLd167B3LnaXX5mEhLg11+tGb4QViNNQEKANoErMDDrY8qX1+72g4Lg+ee11VrXr4dffoFp0zI/r1s37TghbIzUAIQAbW2f7CQlaf/Onat1/kZFaTuwZbehTUEvYSFEDsk7UwjQRvRkJzZWa/a5v9KcXQXaYMjZtYXQgSQAIUBb1dNszv64JUtg/HitH6BUKXjtNdi0KfPjXV2hWTPrxSmEFclqoEKAtqRzTrz3nra8w/nz2j61a9bABx9kfnxampZchLBB0gksxF0NG8Lx49a9ZqtWsHevda8phJVIE5AQd/3739ryDdZiNsM771jvekJYmSQAIe7q3h1attQmb+WVk5M2VPS55/J+LSHyiTQBCXG/W7egbl2IjNR2+nocDg7g7Q0nTmgdxULYKKkBCHE/Ly9tMxdPT21J59xyctIK/4MHpfAXNk8SgBB/V7WqtplLx4656xMwm6FHD/j9d/Dxyb/4hLASSQBCPIqnJ2zYoA3z9PfX1vv5+xaWBoP2M2dnbbTPhg2wapW2JaQQhYD0AQiRE5cva806hw5ps4ENBqhcWVvvv0ULme0rCiVJAEIIYaekCUgIIeyUJAAhhLBTkgCEEMJOSQIQQgg7JQlACCHslCQAIYSwU5IAhBDCTkkCEEIIOyUJQAgh7JQkACGEsFOSAIQQwk5JAhBCCDslCUAIIeyUJAAhhLBTkgCEEMJOSQIQQgg7JQlACCHslCQAIYSwU5IAhBDCTkkCEEIIOyUJQAgh7JQkACGEsFOSAIQQwk5JAhBCCDslCUAIIeyUJAAhhLBTkgCEEMJOSQIQQgg7JQlACCHslCQAIYSwU5IAhBDCTkkCEEIIO/X/VTSYiNSqML8AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"length = nx.get_edge_attributes(RNA.G, 'length')\n",
"pos = nx.circular_layout(RNA.G, scale=1.5)\n",
"nx.draw_networkx(RNA.G, pos, font_color=\"white\")\n",
"nx.draw_networkx_edge_labels(RNA.G, pos, edge_labels=length)\n",
"nx.draw(RNA.G, pos, node_color=color_map, node_size=500)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Unlinking low complex contigs"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"RNA.unlink_dusted()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Drawing"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD8CAYAAAB+UHOxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xd4VFX6wPHvTCoTApEQSujVAKFIMSBSREAkGDpSRcBVaYsKiAqrLnZ+oLIoUqUISxFBQhUFpIuEgEDokNACQhKSENOZ8/vjKCuQSia5k8z7eZ55SLn35B1m5rz3nmpSSimEEEI4HLPRAQghhDCGJAAhhHBQkgCEEMJBSQIQQggHJQlACCEclCQAIYRwUJIAhBDCQUkCEEIIByUJQAghHJQkACGEcFCSAIQQwkFJAhBCCAclCUAIIRyUJAAhhHBQkgCEEMJBSQIQQggHJQlACCEclCQAIYRwUM5GByBEYZCSAkePQmgo/P47mExQrhw0aQL16oGrq9ERCpF7kgCEyMKpUzBtGixdCs7OkJ4OSUn6dxYLmM1gtcKQIfDqq1C9urHxCpEbJtkUXoj7pabCv/8Nn30GaWm64s+Ki4tOEBMnwoQJ+msh7J0kACHucfMmPPEEnDkDiYm5O9fDA+rXhy1bwNMzf+ITwlYkAQjxN3/8Ac2bw+nT+i7gQbi5QYMGsGuX/loIeyWjgIT4m9deg7NnH7zyB91hHBYG//qX7eISIj/IHYAQf9qzBzp0+F8nb14VKwZ790KjRrYpTwhbkzsAIf709ttZV/6urjBvHkREQHw8HDoEnTplfnxyMkyebPMwhbAZSQBCAJcu6av1rDg76+PatIGSJWHSJFi5EqpUyfh4pWDjRoiKsn28QtiCJAAhgO3bsx+6mZioh4ZeuKAr9w0bIDxcTwbLjJsb7N5t21iFsBVJAEKg2/8TEnJ3TpkyULu27vDNTEIC7N+ft9iEyC+SAIRAD/vMDWdnPTt40SI9WzgzVmvWvxfCSJIAhABu3875sSYTfPONHio6alT2x2c3i1gIo8iEdSHQzTk5NX8+lC0LnTvnrHIvW/bB4xIiP8kdgHB44eHhxMSsw2RKyfbYr76COnXgmWf0MM/seHjomcVC2CNJAMLhxMTEMGXKFFq0aIGHhwfVq1cnNPRLzOas24EqV4aXX9YTu65dg1u39KN//8zPUUoSgLBfMhNYFHmpqaksXbqUZcuWERISws2bN/H09KRRo0b07t2bIUOG4OFRnMqV4fJl2/5tPz84ccK2ZQphK9IHIIocq9XKxo0bWbRoEbt37+b333/Hzc2NOnXqMGrUKF588UUqVqx433njxt1m3Lg00tPdbRKHhwe89ZZNihIiX8gdgCgSDhw4wJw5c9i2bRsRERGYTCZq1KhBx44deemll/D398/y/I0bN9K370ASEvYBtVAqb62jTk7wyCN6DoBZGlqFnZI7AFEohYeHM2vWLDZu3Mjp06dJS0ujQoUKtGnThrlz59K2bVvMOah5ExISCAoK4ueff6ZHjx5MnFidVq3M/PFH3uKzWGD5cqn8hX2TOwBRIK5fv86MGTOIjY2ld+/etG7dGqUUJpMpR+fHxMQwb9481qxZw5EjR0hMTKR06dIEBAQwYMAAevfujXMut+GaPXs2Y8aMoWTJkgQHBxMQEADAtm16lE9uN4MBPUfAw0NvCNOiRe7PF6IgSQIQ+S49PZ0PP/yQc+fO0bJlS2bPns2iRYuybJbJScdt8eLFHyieixcvEhgYyPHjxxkzZgxTp069727h4EHo0UMv5JbTRGCxQIUKsHo1ZNPiJIR9UEIUgMqVK6uEhASllFLTpk1T48aNU9euXbvze6vVqnbv3q169eqlypUrp0wmk3J3d1ePPPKI+te//qUuXbpkkzjeeOMN5eTkpPz8/NT58+ezPDYpSalJk5QqUUIpT0+l9KDOux8mk/6dl5dS77+vVGqqTcIUokBIC6XId/Hx8QQEBHD06FEA2rRpQ1JSEocPH75zTFpaGjt27ODIkSP06tWLI0eOkJSURGhoKJMnT85w1E5uhIaG4uvry7Rp05g2bRonTpygWrVqWZ7j7g7vvQc3bsDXX8Pw4Xq/3/LlwdcXGjaEkSP1shDXr+sN4V1c8hSmEAVKmoBEvrt16xbvv/8+fn5+DBkyhKtXrzJ37lxKlizJmDFj7hynctEnkFPp6ek899xzLF++nMcff5zg4GC8vLxs+jeEKKxkFJDIdx4eHlSuXJnffvsNAC8vL9LT0++riG1d+QcHBzNw4ECUUqxZs4auXbvatHwhCjtJAPbu+HE9LGXnTjh/Xi9b6eMDrVpBy5Z6eyonJ6OjBCAlJYWlS5eyfPly6tevz+TJk/Hw8MBsNtOmTRsWL15MZGQkvr6+7N27l1atWuVLHPHx8QQFBbFz50769OnDkiVLcj1CSNih9HT9Wdi7V++yEx2t3/u1akHr1vDkk3qDBpFj8qmwVxs36kblU6d0f+O9K49t3653HXd3h3Hj4JVX9Ka1BSijGbeurq7UrVuXSpUqUaxYsTvH+vv789hjjzF27FjCw8MpW7Ysfn5+No/piy++YOzYsZQqVYr9+/fTrFkzm/8NUcCSk+H//g8+/xzS0vSwrL+v333wIKxdq7+uXx8+/FAnA5E9Y/ugxX1iYpTq0UMpiyXjYScZPSwWpWrUUOrQoXwP79dff1UvvPCCql69ujKbzcrJyUnVqlVLjRo1Sh09ejTLc1NTU9XOnTtVcHCwSktLs2lcERERys/PT5nNZjV+/HhltVptWr4wyC+/KFWpklLFiuXu8zBwoFLx8UZHb/ckAdiTyEilqlRRytU152/2e9/4W7bYNKTz58+r8ePHK39/f+Xq6qpMJpOqWLGi6t+/v9q6dau6ffu2Tf9ebt2+fVuNGzdOmc1mVbduXRUREWFoPMKG1q7N3YXQ3x9ubko9/LBSUVFGPwu7JqOA7MUff0CDBnDxYt62kLJYdH9BVjuVZ+GvGberV6/m6NGjJCYm4u3tTfPmzR94xm1+OXDgAEFBQcTExPDZZ58xYsQIo0MStrJrF3Tq9GDTsf/i4qL7BA4eBDc328VWlBidgcSfXnxRKXf3B7vaufdRpYqexfQ358+fV5MnT77vzyYlJal58+apDh06qIceekgBytPTU7Vq1Ur95z//UXFxcQX0H5BzKSkpqlevXspkMqm2bdvaZYwiD27dUqpsWdt8FiwWpV5/3ehnZLckAdiDgwdz18aZ3aNYMaXeeUfFxMSo5s2bqzFjxigvLy/l5OSkrly5otatW6d69ux5Z8atm5ubzWfc5pfVq1er4sWLK09PT7Vu3TqjwxH5YcwY210M/fV5OHHC6GdllyQB2IPevZUym7N+E3/zje4jiItT6tQppYYNy/L49BIlVPVKlZTJZFImk0kBd77+q+N2xIgR6rfffjP62efIzZs3VcuWLZXJZFL9+/e3eSeysBMJCdm3+48cqdSBA0olJyu1YEH2CcDZWd9hi/tIH4DREhL0uP7sNpitWxfOnoXUVHj4Yfj5ZwgMhNDQDA+PBwYBwff8vG3btmzdujVHSyXbi+nTpzN+/Hh8fHxYt24djRs3NjokkV+++QZGjNCfi8x07w5WKzz1lB4KPWRI9uVaLHDzZoEPlbZ3hacWKKpCQ3P2pjx+XFf+8L9rmxo1Mj3cA2j358xak8lEmTJlMJvNREREFJrKPzw8nIcffpixY8cyduxYrly5IpV/Ubd1a9aVP8CaNXrcf3R0zst1coKwsLzFVgQVjpqgKAsNhZSUnB375Zd6tNCpU3D1qp4slgknYGC1aowcOZL69eszZ84c4uLi2LFjh23izkdWq5VXX32VmjVr4ubmRkREBB999JHRYYmC8Msv+VOu1Zrp3bIjs4/xfI7s999zngBGjoTRo/VOI23bZnuet9XKF198cdfPHnQN/YKyb98+unXrRlxcHDNnzuSll14yOiRRkGJi8qfc5GS9rKu4i9wBGC23C6BZrbBnD1SsqNcntmXZBkpNTaVnz560bNkSf39/oqKipPJ3RPn1njWZZH/ODMj/iNEqVtQdWbnl7JxlHwCgF60vBFatWkWpUqX46aef2LBhA1u3brX7OxWRT8qWzZ9y3d0LzeehIEkCMFqTJtnvIuLjA88+qzebNZuhY0fo1093mGXGbNYrhdqxmJgYHnvsMfr06UOPHj2Ijo7m6aefNjosYaTHH8/+GCcnPbPXyenur7PzgLPjizJJAEZr2DD7pR+U0s09ly/roWxTp+rVP9ety/wcDw9o1862sdrQtGnTKFeuHJcuXeLQoUMsXrzYbpaYEAbq2BE8PbM+ZtIk3ab/5pswaJD+etKkrM9xcpKlojMg8wDswfDhMG9e3tYAulf58jph2Fm755kzZwgMDOT8+fO8+eabvPfee0aHJOxJWhqUKQOxsbYr080NXn0VZCTZfeyrdnBUY8fadjNZDw99RWRHlb/VamX06NH4+flhsVi4ePGiVP7ifi4unOzalT9sWaaTE4waZcsSiwz7qSEcWc2aevMXiyXvZZnN4OcHL7+c97JsZPfu3ZQrV465c+cyZ84cDh8+jK90yIl7JCQk0K5dO+ovWsTN4sVRthgR5OEBn3wCFSrkvawiSBKAvZgwARo10qMV8qJECfj2W7u4+k9JSaFbt260bt2aRo0aERMTw7Bhw4wOS9ihr776itKlSxMWFsbuX36h4q+/YsrrSLBixfS2qbJMeKaMryWE5uwMW7bAI4882J2AiwuUKqXXUa9Wzfbx5dKKFSsoVaoUO3bsYPPmzWzZsgWLLe5wRJFy8eJF/P39GTVqFCNGjODq1asEBARAnTp629OSJfVnI5eSzGaSmjTRS0bYwcWQvZL/GXvi4QE7dug+gWLFcr7Zu4eHnhkcFgb+/vkaYnaioqIICAigX79+9OnTh+joaDp27GhoTMI+TZgwgWrVqnH79m3Onj3Lp59+evc6VU2awNGjEBCg3+M5oJycSATetVopfeAAk95/n6SkpPx5AkWBsYuRikyFhel9Td3dlSpZUikXl/8tb2syKeXpqX/XvLlSwcFK2cEeuB9//LFydnZWlStXznZ/YOG4QkJCVPny5ZWLi4v6/PPPsz/BalVq1SqlmjTRa/t7eurPwF+fB1dXpUqUUKpYMfVHv36qnqurAu48li1blv9PqpCSYaD2Lj4e9u+HkBA4fVovBVG2LDz6qL4yqlTJ6Ag5deoUnTt35sKFC0yaNIl3333X6JCEHUpPT2fQoEGsWLGCli1bsm7dOry8vHJXSHg4/PorHDig1/ZxctLLozdtCgEBWC0W3NzcSP9zSLWLiwu9e/fms88+o0yZMvnwrAo3SQDigVmtVkaOHMmcOXNo2LAhGzdupFy5ckaHJexQcHAwAwcORCnFN998Q7du3fLtb/n5+REVFUX0n8tFOzk58dhjj7Fz5858+5uFlfQBiGxldI2wc+dOypQpw8KFC/n6668JDQ2Vyl/cJz4+njZt2tCtWzc6derEzZs387XyBzhx4gTt27e/872fnx+rVq3K179ZWEkCEPc5ceIE//znP5kyZQpxcXGY7hmPnZSUxFdffUWzZs2Ijo5m8ODBBkUq7NkXX3yBj48Pp0+fZv/+/axcubJAlvswmUz4+/vj5eVFr169OHXqVKHZBKmgSROQuMv58+d56qmnGDVqFPv378fb25t27drRvXv3u45LT0+XtXtEhi5cuECnTp04ffo0Y8eO5eOPPy7wCjg5ORmr1YrFYqFq1aqULl2akJCQAo2hMJC0KO5y8uRJmjdvzpgxY5g5cyaNGzdm06ZNHDt27K7jpPIX97JarYwfP57q1atjMpk4f/48U6ZMMeTq293d/c68k02bNnHo0CHmz59f4HHYO0kA4i7+/v6EhYVx7tw5vLy8aN26NZUqVZI2VJGlAwcOUKFCBf7zn/8wY8YMjh8/TpUqVYwOC4A6deowcuRIRo4cSXx8vNHh2BVJAOIuZcuWpUuXLixZsgSAGjVq0LBhQ+Li4rh165bB0Ql7k5qaSu/evQkICMDPz48bN24wwg6XXvj888956KGH6NKli9Gh2BVJAOIubm5udOjQgcuXL7N69WoAWrVqxfbt20lOTjY4OmFP1qxZg7e3Nz/88APBwcFs376dEiVKGB1WhsxmMxs2bGD37t2sXLnS6HDshjTkijtu376Nk5MTLVq04PLly3zwwQdUrlyZs2fPUq5cOdLS0owOUdiB2NhYunTpwt69e+nXrx+LFi0qFH1CjRs3ZtCgQQwZMoSgoCDc87rwYhEgo4AEiYmJDBw4kJkzZ941ln/WrFmEhoZy6NAhvvzySx599FEDoxT2YPr06YwfPx4fHx/WrVtH48aNjQ4pV6xWK6VLlyYgIIBNmzYZHY7hJAE4uEWLFvHyyy9jsVjYtm0bDRs2vOv3aWlpuNhysxpRKIWHh9OpUyfOnTvH+PHj+agQ7661c+dO2rZty6ZNm3jqqaeMDsdQ0gfgoK5du0bjxo0ZOnQogwcP5saNG/dV/oBU/g7OarXyyiuvULNmTVxdXQkPDy/UlT9A69at6datG3369LmzZpCjkgTggCZPnkzFihWJi4vj+PHjzJo1S2ZKivvs27eP8uXLM2vWLGbOnMnRo0epZAeLD9rC8uXLsVqt9O/f3+hQDGX/PTfCZo4dO0ZgYCCRkZG8//77vPHGG0aHJOxQamoqffv25fvvv+eJJ55g7dq1FM/r7lx2xtXVlW+//ZbIyMg7gx8ckfQBOACr1coLL7zAwoULadasGRs2bKB06dJGhyXs0KpVq3j++edxcnJi2bJldO7c2eiQ8lVKSgpubm5Gh2EYue8v4n788Ue8vb1ZsWIFS5cuZf/+/VL5i/vExMTw2GOP0adPH3r06EF0dHSRr/yBDCv/1NRUYmNj2blzJ1evXjUgqoIjTUBFVGJiIt26deOnn37imWeeYcWKFTLuWWRo2rRpvPnmm5QtW5bQ0FAaNWpkdEiGuH37NiEhIfz444+sWbOGKlWqcPr0aYKDg6levbrR4eULuQMogubPn4+3tzeHDh1i586drF27Vip/cZ8zZ85Qu3ZtJkyYwIQJE7h06ZLDVv6gZzbPnz+fpKQkFixYwOrVqxk0aBCfffaZ0aHlG7kDKEIiIyN5+umnOXbsGCNGjGD69Okyukfcx2q13lnt1d/fn4sXL+Lr62t0WIaaP38+M2fOZNy4cfTr1w/QS0qfPXuWqlWropS6b1+MokASQBHx9ttv8+GHH1KtWjVOnTpFzZo1jQ5J2KHdu3fTo0cP4uPjmT17Ni+88ILRIRkuLS2NX375hS+++IIWLVoAek+DnTt3Yjabeemll4pk5Q+SAAq9I0eOEBgYyO+//84nn3zC2LFjjQ5J2KHk5GT69u1LcHAw7du35/vvv7+zXr6jc3Fx4eDBgzzzzDMkJCSwa9cuDh8+THh4OO3bt6d06dJF9g5AhoEWUlarlSFDhvDNN9/QvHlz1q9fT6lSpYwOS9ihFStWMHToUFxdXVmxYgUdO3Y0OiS7s2vXLiZPnsyNGzdo06YNnp6edO/enSZNmhgdWr6SBFAIbdq0iX79+pGens7ChQvp1auX0SEJOxQVFUVgYCAHDhxg8ODBzJ8/X/qEsnDx4kWcnZ1xd3enZMmSDjE5TN4NhUhCQgLt27cnMDCQJ598kpiYGKn8RYY++eQTypcvz7Vr1zhy5AgLFiyQyj8blStXxtfXl1KlSuHk5ERqaqrRIeU7eUcUEnPmzKF06dIcPXqUPXv28N133+Hq6mp0WMLOnDp1iho1ajBx4kQmTpzIhQsX8Pf3NzqsQsdqtbJlyxbmzJljdCj5ShKAnbt8+TINGjRg+PDhvPzyy1y9evXOSAUh/mK1Whk+fDh169alZMmSXL58mXfffdfosAots9nMvn37GDlyJFFRUUaHk2+kD8COvfXWW0yZMoUaNWqwadOmIjsbUeTNjh076NmzJ3/88QezZs1i8ODBRodUZFSpUgUfHx9CQkKMDiVfyB2AHQoNDaVixYpMnTqVqVOncurUKan8xX2Sk5Pp0qULTzzxBM2aNSM6OloqfxvbvHkzhw4dYsGCBUaHki8kAdiR9PR0BgwYQNOmTalatSrXr1/nlVdeMTosYYeWLl1KqVKl2Lt3Lz/++CObNm2Scf35oE6dOowYMYLhw4cTHx9vdDg2J01AdmL9+vUMGDAAq9XKokWL6NGjh9EhCTt0/fp1OnfuTGhoKMOGDWP27NkyuiefWa1WfH19efjhh9mxY4fR4diUvHMMFh8fzxNPPEFQUBBPPfUU0dHRUvmLDH3wwQdUqFCB6Ohojh49yty5c6XyLwBms5n169eza9cuVq1aZXQ4NiXvHgPNnDkTHx8fTp48yS+//MLKlStlaKe4z/Hjx6lWrRrvvvsu77zzDuHh4dSrV8/osBxK06ZNGThwIIMHDyY5OdnocGxGEoABLl68SL169Rg9ejSjR4/mypUrPProo0aHJeyM1WrlH//4B/7+/nh7e3PlyhUmTZpkdFgOa+HChbi5udGzZ0+jQ7EZSQAFbMKECVSrVg2r1cr58+eZOnWq3MaL+2zbtg0fHx+WLl3K4sWLCQkJoUyZMkaH5dDMZjOrV69m06ZN/PDDD0aHYxPSCVxAQkJCCAoKIjo6mqlTpzJ69GijQxJ2KDExkV69erF582Y6d+7MqlWrZDMfO9O9e3e2bdtGdHQ0zs6Fe0FlufTMZ+np6fTt25dHH32UWrVqcePGDan8RYYWLVqEt7c3+/fvZ9u2baxfv14qfzu0YsUKrFYrAwcONDqUPJMEkI/Wrl1LqVKl2LhxI99//z07duygRIkSRocl7My1a9do3LgxQ4cOZfDgwdy4cYO2bdsaHZbIhKurK8uWLWPlypXs2bPH6HDyRBJAPoiNjaV169Z0796dwMBAYmJiCAoKMjosYYcmT55MxYoViYuL4/jx48yaNUv6hAqBLl268OSTT9K1a1esVqvR4Twwx+4DUAqOH4cDB+DQIYiNBXd38PeHJk2gaVPI5bDMGTNmMG7cOLy9vQkODqZp06b5FLwozI4dO0ZgYCCRkZG89957vPHGG0aHJHIpMTGR0qVLM2DAAObOnat/mJKi65ODByEsTH//0EO6PmnWDB5+GOxpZzHliFJSlPrqK6WqVVPKw0Op4sWV0ulAP9zdlSpRQqmSJZWaMEGp69czLObw4cOqW7duKi0tTUVERCg/Pz9lNpvV66+/XsBPSBQWt2/fVkOGDFEmk0k9+uij6saNG0aHJPLgv//9rzKZTCrsp5+Ueu01pTw9dd3h7n53nVK8uFIWi1K1ayv19ddKpaUZHbpSSinHSwCHDilVs6au+P/+AmX2cHPTL+iKFXeKuHLlitq8ebPy9/dXzs7O6vHHH1dms1nVq1dPXbhwwcAnJ+zZli1blJeXl7JYLOq///2v0eEIW7Ba1eQaNdQtk0kpV9ec1SkeHkrVq6dUWJjR0TtYAvjuO52FTaacvVB/f1gsSo0erZTVqp544gllNpuVm5ubAhSg3nvvPaOfnbBTCQkJqkOHDspkMqmgoCCVlJRkdEjCFm7fVur555XVYsl9fWIy6Tpl0yZDn4Lj9Db98AMMHAiJifolyK3ERJg/nzM9e7Jv3z6sVispKSkAmEwmtm7dauOARVEwf/58vL29OXToEDt37mTt2rUytLOoGDECVq7ElJiY+3OV0nVKz56wa5ftY8shx+gEjoqCmjUhLi7PRSUCnYC/XrKqVavSpEkTWrRowdixY/NcvigaIiMjefrppzl27BgjRoxg+vTpMrqnKNmwAfr00ZV4Xvn4wNmzYMAQccd4Rw4fDklJNinKAqwALC4uODs7ExAQwKpVq6TyF3e8/fbbVK5cmcTERE6ePMmMGTOk8i9KEhPhuedsU/kD3LoFBtUfRf9dGRkJ69dDamr2x9asqRPFN99keVgZi4XfZ88mNTWV5cuX2yhQUdjce/N85MgRKlWqxMcff8zHH3/MmTNnqFWrlkHRiXyzbJke3pmVhx6C1ashIQEiIqBfv8yPTU6GJUvg5k2bhpkTRT8BzJ2b8zb/L7/UY3iz4ZSYSPGZMzHZ03heke9OnDjBP//5T6ZMmUJcXNxdr39KSgo7duygUqVKXLt2jXHjxhkYqchXU6fCH39kfcyXX+qLzrJlYcAA+OorqFs38+PNZli0yLZx5kDR7wNo1gxysqHzs89Cjx56YljNmjBoUNbHu7joN4GLi23iFHbt/PnzPPXUU4waNYr9+/fj7e1Nu3bt6N69+51j0tPTC/3iYCIbt26BtzekpWV+jMWir+b9/eHMGf2zxYvhyhV4883Mz2vXDgp4MEnRvgNQSs/Gy46nJ0yeDK+9lvOy3d1zVrYoEk6ePEnz5s0ZM2YMM2fOpHHjxmzatIljx47dOUYqfwdw+LCu4LNSuzakp/+v8gf47TfIbhOfw4fzHl8uFe0EkJycs7b/996D+fN1hs4ps1n3LwiH4O/vT1hYGOfOncPLy4vWrVtTqVKlIrdFoMjGlSuQ3do/xYvDvRvIx8XpC82sxMbmLbYHULQTgFLZr7vRsCG0bw+fffZg5QuHULZsWbp06cKSJUsAqFGjBg0bNiQuLo5bt24ZHJ0oMDn5zCck3D+ks0QJ3XyUXdkFXKcU7QTg7q6v1LPSti1UrQoXL8LVqzBunJ6ccfBg1ucpBaVL2ypSYefc3Nzo0KEDly9fZvXq1QC0atWK7du3F6k9YkU2fHyyr1NOnwZnZ92X+JeGDbNvMi5evMAXiivajZZmM9SqlfV//Jw58PehnOPG6YQwfHjWZScmQoMGNglT2K/bt2/j5OQEQIsWLbh8+TIffPABlStX5uzZs5QrV460rDoERdHyyCPZzylKTNRDQCdPhhdegEaNoGtXeOyxrM+rX992ceZQ0U4AoJt3Tp3SnTIZSUq6+wVNSNB9B1FRWZdbvToUK2a7OIVdSU9PZ/DgwfzjH/+gZcuWuPw58a9fv37ExcUxZ84cDh06xJdffomvr6/R4YqC4u2t7wKy6y8cMQK+/hquX4foaH1Befx45se7uUFpSc1dAAAUI0lEQVTHjraNNQeK/jDQU6dylrVzw8MDpk2Dl16yXZnCbqxfv54BAwZgtVpZuXIlTz/99H3HpKWl4SJDgB3TlCnw7ru2rVPc3PRyEBUr2q7MHCjafQCgN2Bo0gT+vI23CScnPblDFCnx8fG0bduWoKAgOnbsSHR0dIaVPyCVvyMbNiz7foDccHbWLRUFXPmDIyQAgIULdYa1BYsFZs/WHTaiyJg5cyY+Pj6cOnWKffv28e233+Kay93ghIPw9tYtAB4etinP3V33RRrAMRJAjRrw6afZT+DITrFi8NRTetawKBIuXLhAvXr1GD16NKNGjeLKlSsEBAQYHZawdy++CC1a6Mo7L4oV05W/Qf1IjpEAQLfXT5z44EnAYoHHH9cjhmQNoCLh9ddfp3r16litVs6fP8+0adNk1U6RMyYTrF2rm5cfdDCIxaL7E7JaKC6fOda7/a239EqfJUvmfLN3k0m/wK+9Bhs35nqTeGF/QkJC8PX1Zfr06Xz++eecOHGCKlWqGB2WKGwsFr12z/DhuUsCbm56tdBvv4VRo/IvvhxwrAQAesG3M2d0R47FoqdnZ3TV5+GhX6hOnWDfPr1chKz1Uqilp6fTt29fHn30UWrVqsX169cZPXq00WGJwszNTfcH7NgBTz6pv8+ob8Bs1rOBixeHl1+Gc+egc+eCj/ceRX8YaFb++AM2bYJfftGPW7f0C1i/PrRsqSt/GeNdJKxdu5ZBf67wumTJEoKCggyOSBRJFy/Cli2wZw8cO6bXIitZUvcXtGih+xDtaP6QYycAUeTFxsYSFBTE7t27efbZZ/nmm29k1U4h/uR4TUDCYcyYMYOyZcty9uxZfv31V5YtWyaVvxB/IwlAFDkRERH4+fnxyiuv8MorrxAZGUnTpk2NDksIuyMJQBQZVquV1157jRo1auDs7Ex4eDiffPKJ0WEJYbfkflgUCfv376dr167cvHmTL774guHZreYqhJA7AFG4paam0qtXL1q0aEHdunW5ceOGVP5C5JAkAFFofffdd5QqVYoff/yRdevWsW3bNkrcuxOTECJTkgBEoRMbG0vLli3p3bs33bt3Jzo6msDAQKPDEqLQkQQgCpVPP/2UMmXKcOHCBUJDQ2VcvxB5IAlAFArnzp2jdu3avP7664wfP57Lly/TqFEjo8MSolCTBCDsmtVq5Z///Ce1a9fG3d2diIgIPvjgA6PDEqJIkHtnYbf27NlD9+7diY+P56uvvuLFF180OiQhihS5AxB2JzU1le7du9OqVSsaNGhAVFSUVP5C5ANJAMKurFy5koceeojt27ezceNGfvrpJ4rL9ptC5AtJAMIuxMTE0Lx5c/r27UuvXr2Ijo6mU6dORoclRJEmCUAYbsqUKZQtW5bIyEgOHz7MokWLcHJyMjosIYo86QQWhjlz5gxPP/00ERERTJw4kX//+99GhySEQ5EEIAqc1Wpl9OjRzJo1iwYNGnDp0iXKly9vdFhCOBxJAKJA7dq1ix49epCQkMC8efMYMmSI0SEJ4bCkD0AUiOTkZJ555hnatGlDkyZNiI6OlspfCINJAhD5btmyZXh7e7Nnzx62bNnC5s2bsVgsRoclhMOTBCDyTVRUFM2aNWPAgAH069ePqKgo2rdvb3RYQog/SQIQ+eKjjz6ifPny3Lhxg6NHjzJv3jzMZnm7CWFPpBNY2NSJEycIDAzk4sWLvP3227z99ttGhySEyIRckgmbsFqtvPjii9SrV49SpUpx5coVqfyFsHNyByDy7Oeff6Znz54kJSWxcOFCnnvuOaNDEkLkgNwBiAeWnJxM586dadeuHc2bNycqKkoqfyEKEUkA4oEsXryYUqVKsX//fn766Sc2bNggQzuFKGQkAYhcuXbtGk2aNOH5559n0KBB3Lhxg3bt2hkdlhDiAUgCEDk2efJkKlasyM2bNwkLC2P27NkytFOIQkw6gUW2wsLCCAwM5MqVK0yePJm33nrL6JCEEDYgl28iU1arlWHDhlG/fn3KlCnDlStXpPIXogiROwCRoa1bt9K7d2+Sk5NZsmQJ/fv3NzokIYSNyR2AuEtiYiKdOnWiQ4cOPP7448TExEjlL0QRJQlA3LFgwQK8vb0JCQnh559/Jjg4GHd3d6PDEkLkE0kAgsjISB555BFeeOEFhg4dyvXr12ndurXRYQkh8pkkAAf37rvvUrlyZW7dusXJkyf58ssvZWinEA5COoEd1JEjR+jSpQtXr17lww8/5PXXXzc6JCFEAZMEYO+UgsuX4eBBuHQJbt8Gb2945BHw8wPn3L2EVquVoUOHsnjxYgICAjh06BDe3t75FLwQNpaWBsePw+HDcPMmODlBlSrQtCn4+hodXaEjCcBe3bwJ8+fDp59CbCy4uEBKik4Irq5gMulkMHAgvPqqTgYZmD59Oi1btqRp06b88MMPPPvss6Snp7N8+XL69OlTwE9KiAd09Kj+LCxfrj8LSkFqqv4cuLnpr318YNw4eP55KFHC6IgLBZNSShkdhLjHt9/CCy9AejokJmZ9rLOz/kAMHw4ffAB/jto5cuQIqampPPbYY5QpU4aHH36Y7du307VrV1asWIGrq2sBPBEh8uiPP2DsWFi8WFfyt29nfbzFohPC4sXQpUvBxFiISQKwJ7dv64p/5crsK/57WSxQoQLs3ElEcjK1atXioYce4saNGwAUK1aMH3/8kZYtW+ZD4ELkg0uXoFUruH4dkpJyd67FAkOHwn/+o+8SRIakCcheKAVDhsB33+W+8gd9Tng4tGjBa1WqkJ6efqfyB0hLS6NUqVI2DFiIfPT77xAQoCv/7K76M5KYCF9/rT9XX3xh+/iKCBnvZy+WLIHVqx+s8v9Lejq3L19m6K5dd35kMpmoUKECjRo1IiYmxgaBCpHPlIJnn4UbNx6s8v9LYiIsWABr19outiJGmoDswfXrULMm3Lplk+ISgGGurmx0dSUhIYG9e/fSokULm5QtRL5btAhGjtTt/7bg5aXvjr28bFNeESJ3APZgxgzdwZWV7dt1O+itW/px8mSmhxYH5pcpw969e4mOjpbKXxQeSsHEiVlX/n99Bv56pKfrtv7MpKToEXXiPpIAjJaertsoU1KyP3bUKPD01I9Mhn3+pfjNm9RPSZF2fxvbvHkzfn5+1K5dm08++cTocIqe7dshLi7rY/76DHh6Qrly+sLo228zPz4pCaZN08lF3EUSgNHCwvLWzpmZlBT44Qfbl+vArFYro0aN4ocffiAsLIxly5ZxMos7MZG1U6dOYbVaiYuL4+rVq/qHGzfmrumnZ0/dhPq3fq8MxcVBRMQDx1pUSQIw2sGDYLXm7NiPPtIdY7t3Q5s2WR+bng4//5zn8MT//Prrr9SqVYsqVarg4uJC3759WSsdjA+sbt261KhRgxEjRlChQgWGDBlCyk8/5e5KffBgPeY/O05O+rMm7iLDQI12/nzOrngmTNBT4FNToW9fWLcOGjXS52fiZmgoc6dMsWGwjqVHjx7UrFnzzvdXrlyhUqVKd76vWLEiv/76613nhIWFsWHDhgKLsTCzWq1ERERw6dIllFIsXLiQj4GyOS2gcmV9ITRsWPbHJiXpjmBxF0kARktLy9lxf69oFi+Gfv2gc+csxzjfunmTjz/+OI8BOq4GDRrclQByIiwsTP7Pc+n235pAc9UkMWiQvhvOSdOOUvquWNxFEoDRvL31cg65fXMqle0Mx8r16xPz2295CE78XYUKFbh48eKd7y9fvkyFChXuOqZPnz6yxlIOmUwmihcvjru7OzExMbi6umJydYX4+JwV8NxzkNNk6+oqw0AzIH0ARmvYUE9bz0rJktCxo17jxMkJ+veH1q1h8+YsTzvg5MTly5dtGKxja9asGWfPnuXChQukpqayfPlygoKCjA6r0Bo8eDDLly9n6dKljBo1ivDwcEp37Jizk1u00EufZDX65+9cXHSTqbiLTAQzWnS0fiNnNQy0dGk9OsLPT48YOnkS/vUv+OmnTE9JdHLiFTc35iYmUrx4cRo1akSvXr0YMmQIJWSlxAe2efNmxowZg9VqZdiwYbzxxhtGh1S0fPWVXtEzuxnxs2bpC6fnnstZua6ueoXd7C62HIwkAHvw5JOwbZtty7RY4PffiUlNZd68eaxZs4YjR46QmJiIt7c3AQEB9O/fn2effRbnXO4pIES+yckFUW6ZTNC9u15nS9xFEoA9+PFH/Qa11dR3Nzc9PG727Pt+FR4ezqxZs9i4cSOnT58mLS0NX19f2rRpw9ChQ3niiSdkS0hhrL599bpYOR0gkR2LBbZuhebNbVNeESIJwB4oBYGB+k2a3ZIQOfHQQ3D2LORgFvCBAweYM2cO27ZtIyIiApPJRPXq1enQoQMvvfQSDRo0yHs8QuTGtWtQu7Zt1sZyc4Nu3fRGMuI+kgDsxfXrUKcO5HXFTotF7ycQGJjrU61WKxs3bmTRokXs3r2b33//HVdXV+rUqcMzzzzDiy++SMWKFfMWnxA5sXy5Ht+fl9VxTSa9S9ipUzICKBOSAOzJ0aN6A4z4+Adbt6RYMb1t3ssv2ySc5ORkli5dyooVKwgJCeHmzZs261COjY3FSz6UIiv/93/w7rsPlgTMZl3p79un7yZEhiQB2JszZ3R/QHh4zt/4rq668l+0CLp2zbfQYmJiMu1QHjBgAH369MlRh/JHH33EqlWrSEpKYurUqXTu3DnfYhaF3H//Cy+9BMnJOZ8r4+GhR8x9953eMF5kTgn7k5am1McfK+XlpZSnp1L6fuD+h8WilLu7UgMGKHX9eoGHee7cOfX6668rf39/5erqqkwmk6pQoYLq37+/2rlzp7Jarfeds2bNGtWgQQOllFJbtmxRXbp0KeiwRWETGalUz576ve7unvnnwdNTKW9vpaZPV+r2baOjLhTkDsCepaXpNX82boS9e+HqVb1wXIkS0LQptG+vR0x4exsdKXB3h3KrVq2YMWMGnp6ed34fGxvLe++9R82aNRk+fDiHDh3io48+YsaMGZQtm+MVYISjun5d3xFs2wahobqT2GzWw0ZbttT9XoGBerKkyBEZAG7PXFygRw/9KASaNWtGs2bNAFBKYbpnqYrw8HDi4+Np3bo1ACkpKZQvX56LFy9KAhDZK1MGXnlFP4RNyIBvkS/urfwBrl27RkpKCnXq1LnzfWJiItWrVy/o8IQQyB2AKECxsbFERUVhNpuJiori7Nmz+Pr64u3tTXp6Ol9//TVr165l4MCB9O7d275mKIeF6U1H9uyBS5f0EMMqVXTTQ+vW8PDDRkcoRK7Z0SdMFHVly5alRIkSJCUl8f333xMSEnJnLR2TyYSTkxMXL17k+eefZ8CAAcbPUFYKVq2CyZP/t+/CvSOzVq3S/TJ+fnrIoiwOJwoR6QQWBSY9PZ1JkyaxdOlSAgICePXVV2nZsmWGx2Y1Q3n48OH4+/vnb7CRkTBwoN6HIadLdHh46A1KFizQ7dVC2DlJAMIQt27dumuEUFaym6H88ssv4+vra7vgTp6Exx/X+8jmdp8GV1e9FMe+fVCtmu1iEiIfSAIQhU5GM5Q9PT1p1KgRPXv2zNuS15GR0KCBXpLjQT8aZrO+Azh6VC/lLYSdkgQgCr3MZig3b96cAQMGZNihvGDBAlq0aIGfn9//fqgUtGuntxnM6/aBrq7w9NPw/fd5K0eIfCQJQBQ5GS15XaFCBVq3bs3QoUNp1aoVXl5eODs7M3HiRL777jvWrVtH2R07YOhQ2y3L7eGhO4k7dbJNeULYmCQAUeTd26EMeqKaUgqz2YzVaqVcuXJEFC+O29mz2Rf47LPwzjtQubJeuvj55/VdQ0ZatNCzuIWwQ5IAhEOxWq0MGjSIZcuW8fe3fj1gP+CRXQHt28O8eToJ/PorlC+vfx4ZmfHx7u5w4gRUrZr34IWwMZkJLByK2WzmypUrmEymu0Yhtclg5nKG/v1vPS9g/37dZxAZmXnlD3o5jz178hi1EPlDEoBwOF9//TUXLlxg06ZNVK1alR49ejC6cePsr/7NZr0In4+PXrb70iWYMUNf5WcmIQF++cWW4QthM9IEJAToCVw7d2Z9TPny+mo/JASeeUav1rp2Lfz8M0yalPl5QUH6OCHsjNwBCAF6bZ/sJCXpf2fM0J2/0dF6B7bsNrQp6CUshMgheWcKAXpET3ZiY3Wzz99vmrO7gTaZcla2EAaQBCAE6FU9LZbsj1uwAEaP1v0AXl7w6quwfn3mxxcvDs2b2y5OIWxIVgMVAvSSzjnx3nt6eYfTp/U+tStXwgcfZH58WppOLkLYIekEFuIvjzwChw/btszWrWHHDtuWKYSNSBOQEH955x29fIOtWCzw9tu2K08IG5MEIMRfunaFVq305K28cnPTQ0WffDLvZQmRT6QJSIi/u34d6teHqCi909eDcHICX184ckR3FAthp+QOQIi/K1NGb+bi46OXdM4tNzdd+e/dK5W/sHuSAIS4V/XqejOXp57KXZ+AxQLdusFvv0HFivkXnxA2IglAiIz4+EBwsB7mGRCg1/u5dwtLk0n/zN1dj/YJDobly/WWkEIUAtIHIEROnD+vm3X27dOzgU0mqFJFr/ffsqXM9hWFkiQAIYRwUNIEJIQQDkoSgBBCOChJAEII4aAkAQghhIOSBCCEEA5KEoAQQjgoSQBCCOGgJAEIIYSDkgQghBAOShKAEEI4KEkAQgjhoCQBCCGEg5IEIIQQDkoSgBBCOChJAEII4aAkAQghhIOSBCCEEA5KEoAQQjgoSQBCCOGgJAEIIYSDkgQghBAOShKAEEI4KEkAQgjhoCQBCCGEg5IEIIQQDkoSgBBCOChJAEII4aAkAQghhIOSBCCEEA5KEoAQQjgoSQBCCOGgJAEIIYSD+n/dZm0IAEHEMgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"length = nx.get_edge_attributes(RNA.G, 'length')\n",
"pos = nx.circular_layout(RNA.G, scale=1.5)\n",
"nx.draw_networkx(RNA.G, pos, font_color=\"white\")\n",
"nx.draw_networkx_edge_labels(RNA.G, pos, edge_labels=length)\n",
"nx.draw(RNA.G, pos, node_color=color_map, node_size=500)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Test 2"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [],
"source": [
"RNA2 = RNA_Graph(\"test4.gfa\",\"t_test.tsv\")"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD8CAYAAAB+UHOxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xl8TPf6wPHPTGISk8SSSBBB7VGx01C1VO1RO7VW0Vv7rRbVlturuvvRXleraK2tWqpV+1JLxVakoYidxJJQskuzZ76/P065lGwyyZlknvfrNS9Zzjl5xpzzfc75rgallEIIIYTdMeodgBBCCH1IAhBCCDslCUAIIeyUJAAhhLBTkgCEEMJOSQIQQgg7JQlACCHslCQAIYSwU5IAhBDCTkkCEEIIOyUJQAgh7JQkACGEsFOSAIQQwk5JAhBCCDslCUAIIeyUJAAhhLBTkgCEEMJOSQIQQgg75ah3AEIUCikpcPIkBAfDH3+AwQDlykHjxlCnDphMekcoRK5JAhAiK+fOwezZsGIFODpCejokJWm/M5vBaASLBYYNg9deg6pV9Y1XiFwwyKLwQjxCaiq8+y589hmkpWkFf1aKFdMSxNSpMGWK9rUQNk4SgBB/FxMDzz4LFy5AYmLu9nVxgbp1YccOcHPLn/iEsBJJAELc788/oVkzOH9eewp4HE5OUK8e7NunfS2EjZJeQELc7/XX4eLFxy/8QWswDgmBf/3LenEJkQ/kCUCIuw4cgPbt/9fIm1fFi8PBg9CggXWOJ4SVyROAEHe9807Whb/JBF9/DWFhEB8Px45Bp06Zb5+cDDNmWD1MIaxFEoAQANeuaXfrWXF01LZr3RpKloRp02DNGqhc+dHbKwVbtkBkpPXjFcIKJAEIAbBnT/ZdNxMTta6hV65ohfvmzRAaqg0Gy4yTE+zfb91YhbASSQBCgFb/n5CQu328vKBmTa3BNzMJCXD4cN5iEyKfSAIQArRun7nh6KiNDl62TBstnBmLJevfC6EjSQBCAGRk5HxbgwG++UbrKjpuXPbbZzeKWAidSAIQArTqnJxatAjKloXevXNUuF9LTWXGjBl07tyZ8+fPk5GRQUpKSh6CFcI6ZByAsHuhoaGcGjGCDr/8glN2l8OXX2r9+tu100YNZyMBmOToyMKMDJRS+Pr6cvnyZTp16sT69eut8waEeEzyBCDsTnR0NDNnzqR58+a4uLhQtWpVvggOJsOYzeVQqRKMGqUlgJs34c4d7TVwYKa7GIB96encvc86e/Ys6enp1K9f34rvSIjHI08AoshLTU1lxYoVrFy5kqCgIGJiYnBzc6NBgwb07duXYcOG4eriohXw169b9W/f8fGhbFQUSUlJGAwG7r/c3N3d8ff3Z+DAgfTr1w+TrCkgCpg8AYgix2KxsGnTJvr27Uv58uVxdnZmzJgxREZGMm7cOK5du0Z8fDyBgYGMHz8eV1dXMBjImDSJZGtO4+zigtuHHxIYGEjp0qUZNmwYzs7OmM1mLl26xMiRIwkPD2fEiBE4Ozvj4+PDgAED2L59OxaLxXpxCJEJeQIQRcLRo0dZuHAhu3fvJiwsDIPBQLVq1ejQoQMjR47Ez88vy/23bNnC4P79OZSQQA3AmNfLwsEBGjbUxgAYjSQnJ+Pk5MTGjRsJDg5m+vTpD2weHBzMwoUL2blzJ2FhYQBUqVKFdu3aMXLkSBrIfEIiH0gCEIVSaGgo8+fPZ8uWLZw/f560tDQqVKhA69atGT58OG3atMGYXZ0+kJCQQLdu3fjll1/o1asXK6dOpVjLljlq4M2Sm5s2V1C1arne1WKxsH37dpYtW8a+ffu4ceMGJpOJWrVq0bVrV0aOHEmlSpXyFt9fYmNjKVWqlFWOJQohJUQB+OOPP9S0adPUuHHj1N69e5VSSlkslhzvHxUVpT755BPVrFkzZTabFaDKlCmjAgIC1HfffafS0tJyHdP8+fOVk5OT8vLyUr/++uv/frFrl1Jms1LahA+5exkMSrm6KnXwYK7jyUxKSopatmyZ6tixo3J3d1eAcnV1VS1atFCzZ89WcXFxj3XcDz/8UDVq1EjVrl1bbd682WrxisJDEoDId2lpaerdd99VL774olqwYIFq1KiROnnyZJb7pKSkqMWLF6v27dur0qVLK0C5ubmpli1bqv/+97/qzp07jx3PlStXlJ+fnzIajeq1115TGRkZD28UFKRUpUq5SwRms1I1aiiVzXvLq5iYGDVr1iz19NNPKxcXFwUod3d31alTJ/XNN9+olJSUbI+xbt06Va9ePaWUUjt27FBdu3bN15iFbZIEIApEpUqVVEJCglJKqdmzZ6tJkyapmzdv3vu9xWJR+/fvV3369FHlypVTBoNBOTs7q4YNG6p//etf6tq1a1aJ480331QODg7K19dXXb58OeuNk5KUmjZNqRIllHJzy/yO381NqVKllHr/faVSU60SZ26EhYWpt956S9WtW1eZTCZlMBiUt7e36t+/vwoKCnroSSsmJka9/vrrat68eUoppYKDg1Xfvn0f+DyEfZBeQCLfxcfH4+/vz8mTJwFo3bo1SUlJHD9+/N42aWlp7N27lxMnTtCnTx9OnDhBUlISwcHBzJgxAx8fnzzFEBwcjLe3N7Nnz2b27NmcOXOGKlWqZL2TszO89x7cvg2LF8Po0dp6v+XLg7c31K8PY8dq00LcuqUtCF+sWJ7ifByVK1fmww8/5MSJE6SkpHDkyBGef/55jh49yuHDhx/qURQaGkp8fDytWrUCICUlhfLly3P16tUCj13oy4p93oR4NIPBQJUqVThz5gzNmjXD29sbLy8vzp49S8eOHQEwmUy89dZbvP3221b92+np6bz44ousWrWKZ555htOnT+e+0dNkgj59tFch0KRJE5o0aQKAUgqDwfDA72/evElKSgq1a9e+931iYiJVq1Yt8FiFviQBiHzn4uJCpUqV+P333wEoVaoU6enpDxXEfy+o8mrDhg0MHjwYpRTr1q2je/fuVj1+YfCo/9PY2FgiIyMxGo1ERkZy8eJFvL298fDw0CFCoSepArJ1p0/D559Dv37QpInWt7xDB61qYvfu3M1imc9SUlJYvHgxHTp0YOLEifz5V1dKo9FI69atOXToEBERERQvXpyDBw/i7e2dL3HEx8fTpk0bevToQZcuXYiJibHLwj8zZcuWpUSJEiQlJfHTTz8RFBREz549Ae3/btSoUXTu3Jlvv/2W1NRUnaP9n/R02LEDpk/XpmJq2FC7JAYM0KZoyu2M3gLpBmqzNm9WqkEDpYoXV8rZ+eHGR0dHrfHR01OpTz5RKgc9P6wtIyNDbdy48YGGWycnJ9WwYUP12WefPdS7ZsKECap///7K399fdevWTV29etXqMc2dO1eZTCZVrlw5deTIEasfvyhIS0tTU6ZMUT4+Pqp3795q//79935nsVjU+++//8gG5W3btj26x1Q+S0pSasYMpdzdtVPeweHhy6F4ce311FNK7dxZ4CEWWpIAbE10tFK9euW++2G1akodO5bv4R05ckS9/PLLqmrVqspoNCoHBwdVo0YNNW7cuGy7dqampqrAwEC1YcOGx+q3n5WwsDDl6+urjEajmjx5cq7GGNiz+Pj4LH9/9OhRNXLkSFWtWrV7n3f16tXVqFGj1LECON9+/VWpihW1wj03l8PgwUpl89aEkgRgWyIilKpcWSmT6fEGIZnNSu3YYdWQLl++rCZPnqz8/Pzu3RH6+PiogQMHql27dulyR3i/jIwMNWnSJGU0GtWTTz6pwsLCdI2nKMvIyFBbtmxR/fr1U+XLl7/3xFevXj319ttvqytXrlj1761f//jj8ZyclKpVS6nISKuGVORIArAVCQlKVa2qVe08zhl/fxIICnrsMO6OuPX397834tbDwyNPI27zy5EjR1S5cuWUyWRSX3zxhd7h2J2UlBS1dOnSB0You7i43BuhHBsb+9jHDgx8/ML/7qtYMaXq1FEqOdmKb7qIkQRgK1555dF1/Y/zqlxZqzi9z+XLl9WMGTMe+rNJSUnq66+/znTE7eNOM5CfUlJSVJ8+fZTBYFBt2rSxyRjtkTVGKCul1J07SpUta51LwWxW6o038vmNF2KSAGzBb7/lrpIzu1fx4kr9+98qOjpaNWvWTL366quqVKlSysHBQYWHh6uNGzeq3r17P9Rwa80Rt/nlxx9/VK6ursrNzU1t3LhR73BEFsLCwtSbb76Z6wblV1+13r3Q3cvhzJkCfvOFhCQAW9C3r1JGY9Zn8TffaG0EcXFKnTun1IgRWW6fXqKEqlqxojIYDMpgMCjg3td3G27HjBmjfv/9d73ffY7ExMSoFi1aKIPBoAYOHGhTVVEiZ3LSoJyQkH3Vz9ixSh09qlXtLFmSfQJwdNQesMXDZDpovSUkgKcnJCdnvd2TT8LFi5CaCrVqwS+/QEAABAc/cvN4YAiw4W8/b9OmDbt27crRVMm2Ys6cOUyePBlPT082btxIo0aN9A5J5NHdKa+XLFnC/v37uXnzJiaTCS+vidy8+Q5paU6Z7tuzJ1gs0LEjFC8Ow4Zl//fMZoiJ0QZ1i/8pPKVAURUcnLOz8vRprfCH/93cZDHXvAvQ9q9RoAaDAS8vL4xGI2FhYYWm8A8NDaVWrVpMnDiRiRMnEh4eLoV/EWE0GuncuTNr1qwhIiKC5ORkFixYQEZG6ywLf4B162D9eoiKyvnfc3CAkJA8Bl0EFY6SoCgLDoaUlJxt+8UX2kIl587BjRuwZUummzoAg6tUYezYsdStW5eFCxcSFxfH3r17rRN3PrJYLLz22mtUr14dJycnwsLC+Oijj/QOS+Qjk8nE0KFDcXPrkC/Ht1gyfVi2azIXkN7++CPnCWDsWBg/Hpo3hzZtst3Pw2Lh888/f+Bnrq6ujxlowTh06BA9evQgLi6OefPmMXLkSL1DEgUoOjp/jpucrE3qKh4kTwB6y+0EaBYLHDgAPj7a9MTWPLaOUlNT6d27Ny1atMDPz4/IyEgp/O1Qfp2yBgMUkprPAiX/JXrz8dFasnLL0TH79WbzabI1a1u7di3u7u7s3LmTzZs3s2vXLpt/UhH5o2zZ/Dmus3OhuRwKlCQAvTVunP0iIp6e8MIL4OKi3cZ06KBNgbhrV+b7GI3QurV1Y7Wy6Ohonn76afr160evXr2Iioqic+fOeocldPTMM9lv4+AATk7av/d/nZ3GjfMeX1EjCUBv9etr89xmRSmtuuf6da0v26xZMGECbNyY+T4uLtC2rXVjtaLZs2dTrlw5rl27xrFjx1i+fDmOjtIkZe86dAA3t6y3mTZNq9N/6y0YMkT7etq0rPdxcICaNa0XZ1Eh4wBswejR8PXX2SeC3ChfXksYNlbxeeHCBQICArh8+TJvvfUW7733nt4hCRuSlgZeXhAba71jOjnBa6+BdCR7mG2VDvZq4kTrriXr4qLdEtlQ4W+xWBg/fjy+vr6YzWauXr0qhb94SLFi0L37WeBPqx3TwQHGjbPa4YoU2ykh7Fn16tqC4mZz3o9lNIKvL4walfdjWcn+/fspV64cX331FQsXLuT48eP5thqYKLwSEhJo27Yty5bVxdU1Bm0Gk7xxcYFPPoEKFawQYBEkCcBWTJkCDRpo3RXyokQJ+P57m7j7T0lJoUePHrRq1YoGDRoQHR3NiBEj9A5L2KAvv/ySMmXKEBISwq+/7ufIER9cXfPWJ7R4cWjRAsaMsVKQRZD+pYTQODpqC542bPh4TwLFioG7O+zbB1WqWD++XFq9ejXu7u7s3buXbdu2sWPHDszWeMIRRcrVq1fx8/Nj3LhxjBkzhhs3buDv70/t2rBnD5QsqV0auWU0JtG4cRLr19vEvZDNkv8aW+LiAnv3am0CxYvnrG/b3f3atNEmO/Hzy9cQsxMZGYm/vz8DBgygX79+REVF0aFD/gzvF4XblClTqFKlChkZGVy8eJFPP/30gXmqGjeGkyfB3187xXPCwUEBiVgs0zl6tAzvvz+NpKSk/HkDRYG+k5GKTIWEaAubOjsrVbKktrzR3WngDAZtdWxnZ6WaNVNqwwalbGAN3I8//lg5OjqqSpUqZbs+sLBfQUFBqnz58qpYsWLqP//5T7bbWyxKrV2rVOPG2tz+bm7aJXD3cjCZlCpRQvvdgAF/KpOpjgLuvVauXFkA76pwkm6gti4+Hg4fhqAgOH9emwqibFl46int1qhiRb0j5Ny5c3Tp0oUrV64wbdo0pk+frndIwgalp6czZMgQVq9eTYsWLdi4cSOlSpXK1TFCQ+HIETh6VJvbx8FBmx29SRPtcjCbLTg5OZH+V5fqYsWK0bdvXz777DO8vLzy420VapIAxGOzWCyMHTuWhQsXUr9+fbZs2UK5cuX0DkvYoA0bNjB48GCUUnzzzTf06NEj3/6Wr68vkZGRRP01X7SDgwNPP/00gYGB+fY3CytpAxDZetQ9QmBgIF5eXixdupTFixcTHBwshb94SHx8PK1bt6ZHjx506tSJmJiYfC38Ac6cOUO7du3ufe/r68vatWvz9W8WVpIAxEPOnDnDP//5T2bOnElcXByGv03RmJSUxJdffknTpk2Jiopi6NChOkUqbNnnn3+Op6cn58+f5/Dhw6xZs6ZApvswGAz4+flRqlQp+vTpw7lz5wrNIkgFTaqAxAMuX75Mx44dGTduHIcPH8bDw4O2bdvSs2fPB7ZLT0+XuXvEI125coVOnTpx/vx5Jk6cyMcff1zgBXBycjIWiwWz2cwTTzxBmTJlCAoKKtAYCgNJi+IBZ8+epVmzZrz66qvMmzePRo0asXXrVk6dOvXAdlL4i7+zWCxMnjyZqlWrYjAYuHz5MjNnztTl7tvZ2fneuJOtW7dy7NgxFi1aVOBx2DpJAOIBfn5+hISEcOnSJUqVKkWrVq2oWLGi1KGKLB09epQKFSrw3//+l7lz53L69GkqV66sd1gA1K5dm7FjxzJ27Fji4+P1DsemSAIQDyhbtixdu3bl22+/BaBatWrUr1+fuLg47ty5o3N0wtakpqbSt29f/P398fX15fbt24yxwbkX/vOf/1C6dGm6du2qdyg2RRKAeICTkxPt27fn+vXr/PjjjwC0bNmSPXv2kJycrHN0wpasW7cODw8Ptm/fzoYNG9izZw8lSpTQO6xHMhqNbN68mf3797NmzRq9w7EZUpEr7snIyMDBwYHmzZtz/fp1PvjgAypVqsTFixcpV64caWlpeocobEBsbCxdu3bl4MGDDBgwgGXLlhWKNqFGjRoxZMgQhg0bRrdu3XDO68SLRYD0AhIkJiYyePBg5s2b90Bf/vnz5xMcHMyxY8f44osveOqpp3SMUtiCOXPmMHnyZDw9Pdm4cSONGjXSO6RcsVgslClTBn9/f7Zu3ap3OLqTBGDnli1bxqhRozCbzezevZv69es/8Pu0tDSKWXOxGlEohYaG0qlTJy5dusTkyZP5qBAvrxUYGEibNm3YunUrHTt21DscXUkbgJ26efMmjRo1Yvjw4QwdOpTbt28/VPgDUvjbOYvFwoQJE6hevTomk4nQ0NBCXfgDtGrVih49etCvX797cwbZK0kAdmjGjBn4+PgQFxfH6dOnmT9/voyUFA85dOgQ5cuXZ/78+cybN4+TJ09S0QYmH7SGVatWYbFYGDhwoN6h6Mr2W26E1Zw6dYqAgAAiIiJ4//33efPNN/UOSdig1NRU+vfvz08//cSzzz7L+vXrcXV11TssqzKZTHz//fdERETc6/xgj6QNwA5YLBZefvllli5dStOmTdm8eTNlypTROyxhg9auXctLL72Eg4MDK1eupEuXLnqHlK9SUlJwcnLSOwzdyHN/Effzzz/j4eHB6tWrWbFiBYcPH5bCXzwkOjqap59+mn79+tGrVy+ioqKKfOEPPLLwT01NJTY2lsDAQG7cuKFDVAVHqoCKqMTERHr06MHOnTt5/vnnWb16tfR7Fo80e/Zs3nrrLcqWLUtwcDANGjTQOyRdZGRkEBQUxM8//8y6deuoXLky58+fZ8OGDVStWlXv8PKFPAEUQYsWLcLDw4Njx44RGBjI+vXrpfAXD7lw4QI1a9ZkypQpTJkyhWvXrtlt4Q/ayOZFixaRlJTEkiVL+PHHHxkyZAifffaZ3qHlG3kCKEIiIiLo3Lkzp06dYsyYMcyZM0d694iHWCyWe7O9+vn5cfXqVby9vfUOS1eLFi1i3rx5TJo0iQEDBgDalNIXL17kiSeeQCn10LoYRYEkgCLinXfe4cMPP6RKlSqcO3eO6tWr6x2SsEH79++nV69exMfHs2DBAl5++WW9Q9JdWloav/76K59//jnNmzcHtDUNAgMDMRqNjBw5skgW/iAJoNA7ceIEAQEB/PHHH3zyySdMnDhR75CEDUpOTqZ///5s2LCBdu3a8dNPP92bL9/eFStWjN9++43nn3+ehIQE9u3bx/HjxwkNDaVdu3aUKVOmyD4BSDfQQspisTBs2DC++eYbmjVrxqZNm3B3d9c7LGGDVq9ezfDhwzGZTKxevZoOHTroHZLN2bdvHzNmzOD27du0bt0aNzc3evbsSePGjfUOLV9JAiiEtm7dyoABA0hPT2fp0qX06dNH75CEDYqMjCQgIICjR48ydOhQFi1aJG1CWbh69SqOjo44OztTsmRJuxgcJmdDIZKQkEC7du0ICAjgueeeIzo6Wgp/8UiffPIJ5cuX5+bNm5w4cYIlS5ZI4Z+NSpUq4e3tjbu7Ow4ODqSmpuodUr6TM6KQWLhwIWXKlOHkyZMcOHCAH374AZPJpHdYwsacO3eOatWqMXXqVKZOncqVK1fw8/PTO6xCx2KxsGPHDhYuXKh3KPlKEoCNu379OvXq1WP06NGMGjWKGzdu3OupIMRdFouF0aNH8+STT1KyZEmuX7/O9OnT9Q6r0DIajRw6dIixY8cSGRmpdzj5RtoAbNjbb7/NzJkzqVatGlu3bi2yoxFF3uzdu5fevXvz559/Mn/+fIYOHap3SEVG5cqV8fT0JCgoSO9Q8oU8Adig4OBgfHx8mDVrFrNmzeLcuXNS+IuHJCcn07VrV5599lmaNm1KVFSUFP5Wtm3bNo4dO8aSJUv0DiVfSAKwIenp6QwaNIgmTZrwxBNPcOvWLSZMmKB3WMIGrVixAnd3dw4ePMjPP//M1q1bpV9/PqhduzZjxoxh9OjRxMfH6x2O1UkVkI3YtGkTgwYNwmKxsGzZMnr16qV3SMIG3bp1iy5duhAcHMyIESNYsGCB9O7JZxaLBW9vb2rVqsXevXv1Dseq5MzRWXx8PM8++yzdunWjY8eOREVFSeEvHumDDz6gQoUKREVFcfLkSb766isp/AuA0Whk06ZN7Nu3j7Vr1+odjlXJ2aOjefPm4enpydmzZ/n1119Zs2aNdO0UDzl9+jRVqlRh+vTp/Pvf/yY0NJQ6deroHZZdadKkCYMHD2bo0KEkJyfrHY7VSALQwdWrV6lTpw7jx49n/PjxhIeH89RTT+kdlrAxFouFf/zjH/j5+eHh4UF4eDjTpk3TOyy7tXTpUpycnOjdu7feoViNJIACNmXKFKpUqYLFYuHy5cvMmjVLHuPFQ3bv3o2npycrVqxg+fLlBAUF4eXlpXdYds1oNPLjjz+ydetWtm/frnc4ViGNwAUkKCiIbt26ERUVxaxZsxg/frzeIQkblJiYSJ8+fdi2bRtdunRh7dq1spiPjenZsye7d+8mKioKR8fCPaGy3Hrms/T0dPr3789TTz1FjRo1uH37thT+4pGWLVuGh4cHhw8fZvfu3WzatEkKfxu0evVqLBYLgwcP1juUPJMEkI/Wr1+Pu7s7W7Zs4aeffmLv3r2UKFFC77CEjbl58yaNGjVi+PDhDB06lNu3b9OmTRu9wxKZMJlMrFy5kjVr1nDgwAG9w8kTSQD5IDY2llatWtGzZ08CAgKIjo6mW7dueoclbNCMGTPw8fEhLi6O06dPM3/+fGkTKgS6du3Kc889R/fu3bFYLHqH89jsuw1AKTh9Go4ehWPHIDYWnJ3Bzw8aN4YmTSCX3TLnzp3LpEmT8PDwYMOGDTRp0iSfgheF2alTpwgICCAiIoL33nuPN998U++QRC4lJiZSpkwZBg0axFdffaX9MCVFK09++w1CQrTvS5fWypOmTaFWLbCllcWUPUpJUerLL5WqUkUpFxelXF2V0tKB9nJ2VqpECaVKllRqyhSlbt165GGOHz+uevToodLS0lRYWJjy9fVVRqNRvfHGGwX8hkRhkZGRoYYNG6YMBoN66qmn1O3bt/UOSeTBd999pwwGgwrZuVOp119Xys1NKzucnR8sU1xdlTKblapZU6nFi5VKS9M7dKWUUvaXAI4dU6p6da3gv/8Dyuzl5KR9oKtX3ztEeHi42rZtm/Lz81OOjo7qmWeeUUajUdWpU0dduXJFxzcnbNmOHTtUqVKllNlsVt99953e4QhrsFjUjGrV1B2DQSmTKWdliouLUnXqKBUSonf0dpYAfvhBy8IGQ84+qPtfZrNS48crZbGoZ599VhmNRuXk5KQABaj33ntP73cnbFRCQoJq3769MhgMqlu3biopKUnvkIQ1ZGQo9dJLymI25748MRi0MmXrVl3fgv20Nm3fDoMHQ2Ki9hHkVmIiLFrEhd69OXToEBaLhZSUFAAMBgO7du2ycsCiKFi0aBEeHh4cO3aMwMBA1q9fL107i4oxY2DNGgyJibnfVymtTOndG/bts35sOWQfjcCRkVC9OsTF5flQiUAn4O5H9sQTT9C4cWOaN2/OxIkT83x8UTRERETQuXNnTp06xZgxY5gzZ4707ilKNm+Gfv20QjyvPD3h4kXQoYu4fZyRo0dDUpJVDmUGVgPmYsVwdHTE39+ftWvXSuEv7nnnnXeoVKkSiYmJnD17lrlz50rhX5QkJsKLL1qn8Ae4cwd0Kj+K/lkZEQGbNkFqavbbVq+uJYpvvslyMy+zmT8WLCA1NZVVq1ZZKVBR2Pz94fnEiRNUrFiRjz/+mI8//pgLFy5Qo0YNnaIT+WblSq17Z1ZKl4Yff4SEBAgLgwEDMt82ORm+/RZiYqwaZk4U/QTw1Vc5r/P/4gutD282HBITcZ03D4Mt9ecV+e6cwxSNAAAXAElEQVTMmTP885//ZObMmcTFxT3w+aekpLB3714qVqzIzZs3mTRpko6Rinw1axb8+WfW23zxhXbTWbYsDBoEX34JTz6Z+fZGIyxbZt04c6DotwE0bQo5WdD5hRegVy9tYFj16jBkSNbbFyumnQTFilknTmHTLl++TMeOHRk3bhyHDx/Gw8ODtm3b0rNnz3vbpKenF/rJwUQ27twBDw9IS8t8G7NZu5v384MLF7SfLV8O4eHw1luZ79e2LRRwZ5Ki/QSglDYaLztubjBjBrz+es6P7eycs2OLIuHs2bM0a9aMV199lXnz5tGoUSO2bt3KqVOn7m0jhb8dOH5cK+CzUrMmpKf/r/AH+P13yG4Rn+PH8x5fLhXtBJCcnLO6//feg0WLtAydU0aj1r4g7IKfnx8hISFcunSJUqVK0apVKypWrFjklggU2QgPh+zm/nF1hb8vIB8Xp91oZiU2Nm+xPYainQCUyn7ejfr1oV07+Oyzxzu+sAtly5ala9eufPvttwBUq1aN+vXrExcXx507d3SOThSYnFzzCQkPd+ksUUKrPsru2AVcphTtBODsrN2pZ6VNG3jiCbh6FW7cgEmTtMEZv/2W9X5KQZky1opU2DgnJyfat2/P9evX+fHHHwFo2bIle/bsKVJrxIpseHpmX6acPw+Ojlpb4l3162dfZezqWuATxRXtSkujEWrUyPo/fuFCuL8r56RJWkIYPTrrYycmQr16VglT2K6MjAwcHBwAaN68OdevX+eDDz6gUqVKXLx4kXLlypGWVYOgKFoaNsx+TFFiotYFdMYMePllaNAAuneHp5/Oer+6da0XZw4V7QQAWvXOuXNao8yjJCU9+IEmJGhtB5GRWR+3alUoXtx6cQqbkp6eztChQ/nHP/5BixYtKPbXwL8BAwYQFxfHwoULOXbsGF988QXe3t56hysKioeH9hSQXXvhmDGweDHcugVRUdoN5enTmW/v5AQdOlg31hwo+t1Az53LWdbODRcXmD0bRo603jGFzdi0aRODBg3CYrGwZs0aOnfu/NA2aWlpFJMuwPZp5kyYPt26ZYqTkzYdhI+P9Y6ZA0W7DQC0BRgaN4a/HuOtwsFBG9whipT4+HjatGlDt27d6NChA1FRUY8s/AEp/O3ZiBHZtwPkhqOjVlNRwIU/2EMCAFi6VMuw1mA2w4IFWoONKDLmzZuHp6cn586d49ChQ3z//feYcrkanLATHh5aDYCLi3WO5+ystUXqwD4SQLVq8Omn2Q/gyE7x4tCxozZqWBQJV65coU6dOowfP55x48YRHh6Ov7+/3mEJW/fKK9C8uVZ450Xx4lrhr1M7kn0kANDq66dOffwkYDbDM89oPYZkDqAi4Y033qBq1apYLBYuX77M7NmzZdZOkTMGA6xfr1UvP25nELNZa0/IaqK4fGZfZ/vbb2szfZYsmfPF3g0G7QN+/XXYsiXXi8QL2xMUFIS3tzdz5szhP//5D2fOnKFy5cp6hyUKG7NZm7tn9OjcJQEnJ2220O+/h3Hj8i++HLCvBADahG8XLmgNOWazNjz7UXd9Li7aB9WpExw6pE0XIXO9FGrp6en079+fp556iho1anDr1i3Gjx+vd1iiMHNy0toD9u6F557Tvn9U24DRqI0GdnWFUaPg0iXo0qXg4/2bot8NNCt//glbt8Kvv2qvO3e0D7BuXWjRQiv8pY93kbB+/XqG/DXD67fffku3bt10jkgUSVevwo4dcOAAnDqlzUVWsqTWXtC8udaGaEPjh+w7AYgiLzY2lm7durF//35eeOEFvvnmG5m1U4i/2F8VkLAbc+fOpWzZsly8eJEjR46wcuVKKfyFuI8kAFHkhIWF4evry4QJE5gwYQIRERE0adJE77CEsDmSAESRYbFYeP3116lWrRqOjo6EhobyySef6B2WEDZLnodFkXD48GG6d+9OTEwMn3/+OaOzm81VCCFPAKJwS01NpU+fPjRv3pwnn3yS27dvS+EvRA5JAhCF1g8//IC7uzs///wzGzduZPfu3ZT4+0pMQohMSQIQhU5sbCwtWrSgb9++9OzZk6ioKAICAvQOS4hCRxKAKFQ+/fRTvLy8uHLlCsHBwdKvX4g8kAQgCoVLly5Rs2ZN3njjDSZPnsz169dp0KCB3mEJUahJAhA2zWKx8M9//pOaNWvi7OxMWFgYH3zwgd5hCVEkyLOzsFkHDhygZ8+exMfH8+WXX/LKK6/oHZIQRYo8AQibk5qaSs+ePWnZsiX16tUjMjJSCn8h8oEkAGFT1qxZQ+nSpdmzZw9btmxh586duMrym0LkC0kAwiZER0fTrFkz+vfvT58+fYiKiqJTp056hyVEkSYJQOhu5syZlC1bloiICI4fP86yZctwcHDQOywhijxpBBa6uXDhAp07dyYsLIypU6fy7rvv6h2SEHZFEoAocBaLhfHjxzN//nzq1avHtWvXKF++vN5hCWF3JAGIArVv3z569epFQkICX3/9NcOGDdM7JCHslrQBiAKRnJzM888/T+vWrWncuDFRUVFS+AuhM0kAIt+tXLkSDw8PDhw4wI4dO9i2bRtms1nvsISwe5IARL6JjIykadOmDBo0iAEDBhAZGUm7du30DksI8RdJACJffPTRR5QvX57bt29z8uRJvv76a4xGOd2EsCXSCCys6syZMwQEBHD16lXeeecd3nnnHb1DEkJkQm7JhFVYLBZeeeUV6tSpg7u7O+Hh4VL4C2Hj5AlA5Nkvv/xC7969SUpKYunSpbz44ot6hySEyAF5AhCPLTk5mS5dutC2bVuaNWtGZGSkFP5CFCKSAMRjWb58Oe7u7hw+fJidO3eyefNm6dopRCEjCUDkys2bN2ncuDEvvfQSQ4YM4fbt27Rt21bvsIQQj0ESgMixGTNm4OPjQ0xMDCEhISxYsEC6dgpRiEkjsMhWSEgIAQEBhIeHM2PGDN5++229QxJCWIHcvolMWSwWRowYQd26dfHy8iI8PFwKfyGKEHkCEI+0a9cu+vbtS3JyMt9++y0DBw7UOyQhhJXJE4B4QGJiIp06daJ9+/Y888wzREdHS+EvRBElCUDcs2TJEjw8PAgKCuKXX35hw4YNODs76x2WECKfSAIQRERE0LBhQ15++WWGDx/OrVu3aNWqld5hCSHymSQAOzd9+nQqVarEnTt3OHv2LF988YV07RTCTkgjsJ06ceIEXbt25caNG3z44Ye88cYbeockhChgkgBsnVJw/Tr89htcuwYZGeDhAQ0bgq8vOObuI7RYLAwfPpzly5fj7+/PsWPH8PDwyKfghbCytDQ4fRqOH4eYGHBwgMqVoUkT8PbWO7pCRxKArYqJgUWL4NNPITYWihWDlBQtIZhMYDBoyWDwYHjtNS0ZPMKcOXNo0aIFTZo0Yfv27bzwwgukp6ezatUq+vXrV8BvSojHdPKkdi2sWqVdC0pBaqp2HTg5aV97esKkSfDSS1CihN4RFwoGpZTSOwjxN99/Dy+/DOnpkJiY9baOjtoFMXo0fPAB/NVr58SJE6SmpvL000/j5eVFrVq12LNnD927d2f16tWYTKYCeCNC5NGff8LEibB8uVbIZ2Rkvb3ZrCWE5cuha9eCibEQkwRgSzIytIJ/zZrsC/6/M5uhQgUIDCQsOZkaNWpQunRpbt++DUDx4sX5+eefadGiRT4ELkQ+uHYNWraEW7cgKSl3+5rNMHw4/Pe/2lOCeCSpArIVSsGwYfDDD7kv/EHbJzQUmjfn9cqVSU9Pv1f4A6SlpeHu7m7FgIXIR3/8Af7+WuGf3V3/oyQmwuLF2nX1+efWj6+IkP5+tuLbb+HHHx+v8L8rPZ2M69cZvm/fvR8ZDAYqVKhAgwYNiI6OtkKgQuQzpeCFF+D27ccr/O9KTIQlS2D9euvFVsRIFZAtuHULqleHO3escrgEYITJxBaTiYSEBA4ePEjz5s2tcmwh8t2yZTB2rFb/bw2lSmlPx6VKWed4RYg8AdiCuXO1Bq6s7Nmj1YPeuaO9zp7NdFNXYJGXFwcPHiQqKkoKf1F4KAVTp2Zd+N+9Bu6+0tO1uv7MpKRoPerEQyQB6C09XaujTEnJfttx48DNTXtl0u3zLteYGOqmpEi9v5Vt27YNX19fatasySeffKJ3OEXPnj0QF5f1NnevATc3KFdOuzH6/vvMt09KgtmzteQiHiAJQG8hIXmr58xMSgps327949oxi8XCuHHj2L59OyEhIaxcuZKzWTyJiaydO3cOi8VCXFwcN27c0H64ZUvuqn5699aqUO9r93qkuDgIC3vsWIsqSQB6++03sFhytu1HH2kNY/v3Q+vWWW+bng6//JLn8MT/HDlyhBo1alC5cmWKFStG//79WS8NjI/tySefpFq1aowZM4YKFSowbNgwUnbuzN2d+tChWp//7Dg4aNeaeIB0A9Xb5cs5u+OZMkUbAp+aCv37w8aN0KCBtn8mYoKD+WrmTCsGa1969epF9erV730fHh5OxYoV733v4+PDkSNHHtgnJCSEzZs3F1iMhZnFYiEsLIxr166hlGLp0qV8DJTN6QEqVdJuhEaMyH7bpCStIVg8QBKA3tLScrbd/QXN8uUwYAB06ZJlH+c7MTF8/PHHeQzQftWrV++BBJATISEh8n+eSxn3VYHmqkpiyBDtaTgnVTtKaU/F4gGSAPTm4aFN55Dbk1OpbEc4Vqpbl+jff89DcOJ+FSpU4OrVq/e+v379OhUqVHhgm379+skcSzlkMBhwdXXF2dmZ6OhoTCYTBpMJ4uNzdoAXX4ScJluTSbqBPoK0Aeitfn1t2HpWSpaEDh20OU4cHGDgQGjVCrZty3K3ow4OXL9+3YrB2remTZty8eJFrly5QmpqKqtWraJbt256h1VoDR06lFWrVrFixQrGjRtHaGgoZTp0yNnOzZtrU59k1fvnfsWKaVWm4gEyEExvUVHaiZxVN9AyZbTeEb6+Wo+hs2fhX/+CnTsz3SXRwYEJTk58lZiIq6srDRo0oE+fPgwbNowSMlPiY9u2bRuvvvoqFouFESNG8Oabb+odUtHy5ZfajJ7ZjYifP1+7cXrxxZwd12TSZtjN7mbLzkgCsAXPPQe7d1v3mGYz/PEH0ampfP3116xbt44TJ06QmJiIh4cH/v7+DBw4kBdeeAHHXK4pIES+yckNUW4ZDNCzpzbPlniAJABb8PPP2glqraHvTk5a97gFCx76VWhoKPPnz2fLli2cP3+etLQ0vL29ad26NcOHD+fZZ5+VJSGFvvr31+bFymkHieyYzbBrFzRrZp3jFSGSAGyBUhAQoJ2k2U0JkROlS8PFi5CDUcBHjx5l4cKF7N69m7CwMAwGA1WrVqV9+/aMHDmSevXq5T0eIXLj5k2oWdM6c2M5OUGPHtpCMuIhkgBsxa1bULs25HXGTrNZW08gICDXu1osFrZs2cKyZcvYv38/f/zxByaTidq1a/P888/zyiuv4OPjk7f4hMiJVau0/v15mR3XYNBWCTt3TnoAZUISgC05eVJbACM+/vHmLSleXFs2b9Qoq4STnJzMihUrWL16NUFBQcTExFitQTk2NpZSclGKrPzf/8H06Y+XBIxGrdA/dEh7mhCPJAnA1ly4oLUHhIbm/MQ3mbTCf9ky6N4930KLjo7OtEF50KBB9OvXL0cNyh999BFr164lKSmJWbNm0aVLl3yLWRRy330HI0dCcnLOx8q4uGg95n74QVswXmROCduTlqbUxx8rVaqUUm5uSmnPAw+/zGalnJ2VGjRIqVu3CjzMS5cuqTfeeEP5+fkpk8mkDAaDqlChgho4cKAKDAxUFovloX3WrVun6tWrp5RSaseOHapr164FHbYobCIilOrdWzvXnZ0zvx7c3JTy8FBqzhylMjL0jrpQkCcAW5aWps35s2ULHDwIN25oE8eVKAFNmkC7dlqPCQ8PvSMFHmxQbtmyJXPnzsXNze3e72NjY3nvvfeoXr06o0eP5tixY3z00UfMnTuXsmVzPAOMsFe3bmlPBLt3Q3Cw1khsNGrdRlu00Nq9AgK0wZIiR6QDuC0rVgx69dJehUDTpk1p2rQpAEopDH+bqiI0NJT4+HhatWoFQEpKCuXLl+fq1auSAET2vLxgwgTtJaxCOnyLfPH3wh/g5s2bpKSkULt27XvfJyYmUrVq1YIOTwiBPAGIAhQbG0tkZCRGo5HIyEguXryIt7c3Hh4epKens3jxYtavX8/gwYPp27evbY1QDgnRFh05cACuXdO6GFaurFU9tGoFtWrpHaEQuWZDV5go6sqWLUuJEiVISkrip59+Iigo6N5cOgaDAQcHB65evcpLL73EoEGD9B+hrBSsXQszZvxv3YW/98xau1Zrl/H11bosyuRwohCRRmBRYNLT05k2bRorVqzA39+f1157jRYtWjxy26xGKI8ePRo/P7/8DTYiAgYP1tZhyOkUHS4u2gIlS5Zo9dVC2DhJAEIXd+7ceaCHUFayG6E8atQovL29rRfc2bPwzDPaOrK5XafBZNKm4jh0CKpUsV5MQuQDSQCi0HnUCGU3NzcaNGhA79698zbldUQE1KunTcnxuJeG0ag9AZw8qU3lLYSNkgQgCr3MRig3a9aMQYMGPbJBecmSJTRv3hxfX9///VApaNtWW2Ywr8sHmkzQuTP89FPejiNEPpIEIIqcR015XaFCBVq1asXw4cNp2bIlpUqVwtHRkalTp/LDDz+wceNGyu7dC8OHW29abhcXrZG4UyfrHE8IK5MEIIq8vzcogzZQTSmF0WjEYrFQrlw5wlxdcbp4MfsDvvAC/PvfUKmSNnXxSy9pTw2P0ry5NopbCBskCUDYFYvFwpAhQ1i5ciX3n/p1gMOAS3YHaNcOvv5aSwJHjkD58trPIyIevb2zM5w5A088kffghbAyGQks7IrRaCQ8PByDwfBAL6TWjxi5/EjvvquNCzh8WGsziIjIvPAHbTqPAwfyGLUQ+UMSgLA7ixcv5sqVK2zdupUnnniCXr16Mb5Ro+zv/o1GbRI+T09t2u5r12DuXO0uPzMJCfDrr9YMXwirkSogIUAbwBUYmPU25ctrd/tBQfD889psrevXwy+/wLRpme/XrZu2nRA2Rp4AhABtbp/sJCVp/86dqzX+RkVpK7Blt6BNQU9hIUQOyZkpBGg9erITG6tV+9z/0JzdA7TBkLNjC6EDSQBCgDarp9mc/XZLlsD48Vo7QKlS8NprsGlT5tu7ukKzZtaLUwgrktlAhQBtSueceO89bXqH8+e1dWrXrIEPPsh8+7Q0LbkIYYOkEViIuxo2hOPHrXvMVq1g717rHlMIK5EqICHu+ve/tekbrMVshnfesd7xhLAySQBC3NW9O7RsqQ3eyisnJ62r6HPP5f1YQuQTqQIS4n63bkHduhAZqa309TgcHMDbG06c0BqKhbBR8gQgxP28vLTFXDw9tSmdc8vJSSv8Dx6Uwl/YPEkAQvxd1araYi4dO+auTcBshh494Pffwccn/+ITwkokAQjxKJ6esGGD1s3T31+b7+fvS1gaDNrPnJ213j4bNsCqVdqSkEIUAtIGIEROXL6sVescOqSNBjYYoHJlbb7/Fi1ktK8olCQBCCGEnZIqICGEsFOSAIQQwk5JAhBCCDslCUAIIeyUJAAhhLBTkgCEEMJOSQIQQgg7JQlACCHslCQAIYSwU5IAhBDCTkkCEEIIOyUJQAgh7JQkACGEsFOSAIQQwk5JAhBCCDslCUAIIeyUJAAhhLBTkgCEEMJOSQIQQgg7JQlACCHslCQAIYSwU5IAhBDCTkkCEEIIOyUJQAgh7JQkACGEsFOSAIQQwk5JAhBCCDslCUAIIeyUJAAhhLBTkgCEEMJOSQIQQgg7JQlACCHs1P8DQYWYiBd+/YkAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"color_map = []\n",
"for node in RNA2.G:\n",
" if RNA2.maxDustWindow(RNA2.nodes[node][-25:], 25, 25) > 1:\n",
" color_map.append('blue')\n",
" elif RNA2.maxDustWindow(RNA2.nodes[node][:25], 25, 25) > 1:\n",
" color_map.append('blue')\n",
" else:\n",
" color_map.append('red')\n",
"\n",
"length = nx.get_edge_attributes(RNA2.G, 'length')\n",
"pos = nx.circular_layout(RNA2.G, scale=1.5)\n",
"nx.draw_networkx(RNA2.G, pos, font_color=\"white\")\n",
"nx.draw_networkx_edge_labels(RNA2.G, pos, edge_labels=length)\n",
"nx.draw(RNA2.G, pos, node_color=color_map, node_size=500)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [],
"source": [
"RNA2.unlink_dusted()"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD8CAYAAAB+UHOxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3XlcVOX+B/DPzLA5gJLgEuKuiYZbaoLmUrmUJrmbgpnVzbWfmbbbcrWyvHq7XcvM3cwrmmmigllq4q4ILqFiBoiKSoCABLLN8/vjSROVmQHOcGY4n/frNa+AOefxC815vuc8q04IIUBERJqjVzsAIiJSBxMAEZFGMQEQEWkUEwARkUYxARARaRQTABGRRjEBEBFpFBMAEZFGMQEQEWkUEwARkUYxARARaRQTABGRRjEBEBFpFBMAEZFGMQEQEWkUEwARkUYxARARaRQTABGRRjmpHQCRQ8jPB06eBGJigKtXAZ0OqFsX6NABePBBwMVF7QiJyowJgMic+Hhg3jxg9WrAyQkoKgLy8uR7RiOg1wMmEzB2LDB1KtCkibrxEpWBjpvCE91DQQHwz38Cn30GFBbKit8cZ2eZIN55B3jjDfk1kZ1jAiC607VrwKOPAr/9BuTmlu1cd3egdWtg+3bA09M28REphAmA6HZ//gkEBgJnz8qngPJwdQXatAH27JFfE9kpjgIiut2rrwLnzpW/8gdkh3FcHPDuu8rFRWQDfAIgumnfPqB37787eSuqWjVg/36gXTtlyiNSGJ8AiG567z3zlb+LC7BkCZCUBGRnA7GxwBNPlH78jRvAzJmKh0mkFCYAIgC4cEHerZvj5CSP69EDqFEDmDEDWLcOaNjw3scLAUREAGlpysdLpAAmACIA2LXL8tDN3Fw5NPT8eVm5b90KJCbKyWClcXUF9u5VNlYihTABEAGy/T8np2zn1K4NPPCA7PAtTU4OcOhQxWIjshEmACJADvssCycnOTt45Uo5W7g0JpP594lUxARABADFxdYfq9MBq1bJoaKTJ1s+3tIsYiKVcL46ESCbc6y1dClQpw7Qr591lXudOuWPi8iG+ARAmpeYmIjNGRnI1+ksH/zVV0DLlsCAAXKYpyXu7nJmMZEdYgIgzcnIyMCcOXMQFBQEd3d3NGnSBF/GxKBYb+FyaNAAGD9eTuy6cgW4fl2+Ro0q/RwhmADIbnEmMFV5BQUFWL16NdasWYPo6Ghcu3YNnp6eaNeuHYYNG4axY8fCw91dVvAXLyr7j/v7A6dPK1smkULYB0BVjslkQkREBFauXIm9e/fi6tWrcHV1RcuWLTF58mS89NJL8PPzu+u84unTUTh9OtyU6rR1dwfefluZsohsgE8AVCUcOXIEixYtws6dO5GUlASdToemTZuiT58+GDduHAICAsyeHxERgdBnnsGBnBw0B6Cv6GVhMADt28s5AJaalohUwicAckiJiYlYuHAhIiIicPbsWRQWFqJevXro0aMHFi9ejJ49e0JvRcWbk5OD4OBg/PLLLxg8eDCavPMO9N26yWWhK8JoBMLCWPmTXeMTAFWK1NRUzJ8/H5mZmRg2bBi6d+8OIQR01oy8gey4XbJkCTZu3IgTJ04gNzcXPj4+6Ny5M0JCQjBs2DA4lXEXrq+//hpTpkxBjRo1EB4ejs6dO8s3du6Uo3zKuhkMIOcIuLvLDWGCgsp+PlElYgIgmysqKsLHH3+M33//HV27dsXXX3+NlStXmm2Wsarj1sOjXPEkJyejf//+OHXqFKZMmYK5c+fe/bRw9CgweLBcyM3aRGA0AvXqARs2ABaanIjsgiCqBA0aNBA5OTlCCCHmzZsnpk+fLq5cuXLrfZPJJPbu3SuGDh0q6tatK3Q6nXBzcxPt27cX7777rrhw4YIicbz55pvCYDAIf39/kZCQYP7gvDwhZswQonp1ITw9hZCDOku+dDr5npeXEB9+KERBgSJxElUGNlCSzWVnZ6Nz5844efIkAKBHjx7Iy8vDsWPHbh1TWFiI3bt348SJExg6dChOnDiBvLw8xMTEYObMmfcctVMWMTEx8PX1xbx58zBv3jycPn0ajRs3Nn+Smxswaxbwxx/AsmXAhAlyv9/77wd8fYG2bYFJk+SyEKmpckN4Z+cKxUlUmdgERDZ3/fp1fPjhh/D398fYsWNx+fJlLF68GDVq1MCUKVNuHSfK0CdgraKiIjz77LMICwvDI488gvDwcHh5eSn6bxA5Ko4CIptzd3dHgwYNcPz4cQCAl5cXioqK7qqIla78w8PDERoaCiEENm7ciKefflrR8okcHROAvTt1So5KiYoCEhLkqpW1agHdugFdu8rdqQwGtaMEAOTn52P16tUICwtD69atMXPmTLi7u0Ov16NHjx745ptvkJKSAl9fX+zfvx/dunWzSRzZ2dkIDg5GVFQUhg8fjm+//bbMI4TI/hQVyUth/365x056uvzoN28OdO8OPP643J6BrMerwl5FRMg25fh42d1458Jju3bJTcfd3IDp04FXXpF71laie824dXFxQatWrVC/fn1Uq1bt1rEBAQHo0qULpk2bhsTERNSpUwf+/v6Kx/TFF19g2rRpqFmzJg4dOoROnTop/m9Q5bpxA/jXv4D//AcoLJSDsm5fvfvoUWDTJvl169bAxx/LZEBWULcPmu6SkSHE4MFCGI33HnVyr5fRKETTpkLExto8vMOHD4sXX3xRNGnSROj1emEwGETz5s3F5MmTxcmTJ82eW1BQIKKiokR4eLgoLCxUNK6kpCTh7+8v9Hq9eO2114TJZFK0fFLHwYNC1K8vRLVqZbscQkOFyM5WO3r7xwRgT1JShGjYUAgXF+s/7Xd+8rdvVzSkhIQE8dprr4mAgADh4uIidDqd8PPzE6NGjRI7duwQxcXFiv57ZVVcXCymT58u9Hq9aNWqlUhKSlI1HlLOpk1luw+6/eXqKkSLFkKkpan9W9g3jgKyF3/+CbRpAyQnV2wHKaNR9heY26jcjJszbjds2ICTJ08iNzcX3t7eCAwMLPeMW1s5cuQIgoODkZGRgc8++wwTJ05UOyRSyJ49wBNPlG8y9k3OzrJP4OhRwNVVudiqFLUzEP3lpZeEcHMr3+3Ona+GDeUkptskJCSImTNn3vXP5uXliSVLlojevXuL++67TwAQnp6eolu3buK///2vyMrKqqQ/gPXy8/PF0KFDhU6nEz179rTLGKn8rl8Xok4dZS4Fo1GI119X+zeyX0wA9uDo0bI1clp6VasmxPvvi4yMDBEYGCimTJkivLy8hMFgEJcuXRKbN28WQ4YMuTXj1tXVVfEZt7ayYcMG4eHhITw9PcXmzZvVDodsYMoU5e6Fbl4Op0+r/VvZJyYAezBsmBB6vflP8apVso8gK0uI+HghXnjB7PFF1auLJvXrC51OJ3Q6nQBw6+ubHbcTJ04Ux48fV/u3t8q1a9dE165dhU6nE6NGjVK8E5nsQ06O5Xb/SZOEOHJEiBs3hFi+3HICcHKSD9h0N/YBqC0nR47rt7S/bKtWwLlzQEEB0KIF8MsvQP/+QEzMPQ/PBjAaQPgdP+/Zsyd27Nhh1VLJ9uLzzz/Ha6+9hlq1amHz5s146KGH1A6JbGTVKmDiRHlZlGbQIMBkAvr2lSOhx461XK7RCFy7Vukjpe2e49QCVVVMjHWfylOnZOUP/H1z07RpqYe7A3jsr5m1Op0OtWvXhl6vR1JSksNU/omJiWjRogWmTZuGadOm4dKlS6z8q7gdO8xX/gCwcaMc95+ebn25BgMQF1ex2Koix6gJqrKYGCA/37pjv/xSjhaKjwcuX5aTxUphABDauDEmTZqE1q1bY9GiRcjKysLu3buViduGTCYTpk6dimbNmsHV1RVJSUmYPXu22mFRJTh40DblmkylPixrmn2M59Oyq1etTwCTJgEvvyw3GunZ0+J53iYTvvjiixI/K+8a+pXlwIEDGDhwILKysrBgwQKMGzdO7ZCoEmVk2KbcGzfkoq5UEp8A1FbWBdBMJmDfPsDPTy5PrGTZKiooKMCQIUPQtWtXBAQEIC0tjZW/BtnqI6vTcXfOe+GfRG1+frInq6ycnMz2AQCQa9Y7gPXr16NmzZr4+eefsXXrVuzYscPun1TINurUsU25bm4OczlUKiYAtXXoYHkTkVq1gBEj5F6zej3Qpw8wcqTsMSuNXi9XCrVjGRkZ6NKlC4YPH47BgwcjPT0dTz75pNphkYoeecTyMQaDnNlrMJT82pJyTo6v0pgA1Na2reWlH4SQzT0XL8qxbHPnytU/N28u/Rx3d+Cxx5SNVUHz5s1D3bp1ceHCBcTGxuKbb76xmyUmSD19+gCenuaPmTFDtum/9RYwerT8esYM8+cYDFwq+l44D8AeTJgALFlSsTWA7nT//TJh2FnD52+//Yb+/fsjISEBb731FmbNmqV2SGRHCguB2rWBzEzlynR1BaZOBTiQ7G72VTto1bRpyu4l6+4ub4nsqPI3mUx4+eWX4e/vD6PRiOTkZFb+dBdnZ+D11+XELaUYDMDkycqVV5XYTw2hZc2ayc1flPjU6/WAvz8wfnzFy1LI3r17UbduXSxevBiLFi3CsWPH4MseOSrFa68BDRooMyLI3R349FOgXr2Kl1UVsQnIXhQVyU7bmBjLy0KY4+Uly2jcWLnYyik/Px8jRoxAeHg4evXqhR9++AFGJW/tqMo6fRro3Bm4fr38ZVSrJndOjYy0q4dhu8I/i71wcgK2bwfaty/fk4CzM1CzplxI3Q4q/7Vr16JmzZrYvXs3tm3bhu3bt7PyJ6u1bCl3Pa1RQ14aZaXX56FDhzxs2sTK3xz+aeyJuzuwe7fsE6hWzfrN3t3d5czguDggIMCmIVqSlpaGzp07Y+TIkRg+fDjS09PRp08fVWMix9ShA3DypHwScHe37hyDQQDIhcn0AY4c8cGHH85AXl6eTeN0aGouRUpmxMXJjU3d3ISoUUMIZ+e/17fV6YTw9JTvBQYKER4uhB3sgfvJJ58IJycn0aBBA4v7AxNZy2QSYv16ITp0kGv7e3rKS+Dm5eDiIkT16vK9kSP/FC4uDwoAt15r1qxR+1ewW+wDsHfZ2cChQ0B0NHD2rFwKok4d4OGH5a1R/fpqR4j4+Hj069cP58+fx4wZM/DBBx+oHRJVUYmJwOHDwJEjcm0fg0Gujt6xo7wcjEYTXF1dUfTXkGpnZ2cMGzYMn332GWrXrq1y9PaHCYDKzWQyYdKkSVi0aBHatm2LiIgI1K1bV+2wSOP8/f2RlpaG9L/WizYYDOjSpQuioqJUjsz+sA+ALLrXPUJUVBRq166NFStWYNmyZYiJiWHlT3bh9OnT6NWr163v/f39sX79ehUjsl9MAHSX06dP4//+7/8wZ84cZGVlQXfHgOy8vDx89dVX6NSpE9LT0zFmzBiVIiW6m06nQ0BAALy8vDB06FDEx8c7zCZIlY1NQFRCQkIC+vbti8mTJ+PQoUPw9vbGY489hkGDBpU4rqioiGv3kN26ceMGTCYTjEYjGjVqBB8fH0RHR6sdlt1hWqQSzpw5g8DAQEyZMgULFizAQw89hMjISPz6668ljmPlT/bMzc3t1ryTyMhIxMbGYunSpSpHZX+YAKiEgIAAxMXF4ffff4eXlxe6d++O+vXrsw2VHFbLli0xadIkTJo0CdnZ2WqHY1eYAKiEOnXq4KmnnsK3334LAGjatCnatm2LrKwsXK/IvHwiFf3nP//Bfffdh6eeekrtUOwKEwCV4Orqit69e+PixYvYsGEDAKBbt27YtWsXblRkjSIiFen1emzduhV79+7FunXr1A7HbrAhl24pLi6GwWBAUFAQLl68iI8++ggNGjTAuXPnULduXRQWFqodIlG5PfTQQxg9ejTGjh2L4OBguLm5qR2S6jgKiJCbm4vQ0FAsWLCgxFj+hQsXIiYmBrGxsfjyyy/x8MMPqxglUcWZTCb4+Pigc+fOiIyMVDsc1TEBaNzKlSsxfvx4GI1G7Ny5E23bti3xfmFhIZyV3KyGSGVRUVHo2bMnIiMj0bdvX7XDURUTgEZduXIF/fr1w/Hjx/GPf/wDCxYs4GQZ0ozBgwdjx44dSE9P1/SQZl7xGjRz5kz4+fkhKysLp06dwsKFC1n5k6aEhYXBZDJh1KhRaoeiKu2mPg369ddf0b9/f6SkpODDDz/Em2++qXZIRKpwcXHBd999h5SUlFuDH7SITUAaYDKZ8OKLL2LFihXo1KkTtm7dCh8fH7XDIlJdfn4+XF1d1Q5DNXzur+J++ukneHt7Y+3atVi9ejUOHTrEyp/oL/eq/AsKCpCZmYmoqChcvnxZhagqD5uAqqjc3FwMHDgQP//8MwYMGIC1a9dy3DORGcXFxYiOjsZPP/2EjRs3omHDhjh79izCw8PRpEkTtcOzCT4BVEFLly6Ft7c3YmNjERUVhU2bNrHyJ7Jg48aNWLp0KfLy8rB8+XJs2LABo0ePxmeffaZ2aDbDJ4AqJCUlBU8++SR+/fVXTJw4EZ9//jlH9xBZYenSpViwYAGmT5+OkSNHApBLSp87dw6NGjWCEOKufTGqAiaAKuK9997Dxx9/jMaNGyM+Ph7NmjVTOyQih1BYWIiDBw/iiy++QFBQEADg/PnziIqKgl6vx7hx46pk5Q8wATi8EydOoH///rh69So+/fRTTJs2Te2QiByKs7Mzjh49igEDBiAnJwd79uzBsWPHkJiYiF69esHHx6fKPgFwGKiDMplMGDt2LFatWoXAwEBs2bIFNWvWVDssIoe0Z88ezJw5E3/88Qd69OgBT09PDBo0CB06dFA7NJtiAnBAkZGRGDlyJIqKirBixQoMHTpU7ZCIHF5ycjKcnJzg5uaGGjVqaGJyGBOAA8nJycHAgQOxc+dODBo0CGvWrIGLi4vaYRFVSQUFBVX++uIQEQexaNEi+Pj44OTJk9i3bx++//77Kv/hJFKLyWTC9u3bsWjRIrVDsSkmADt38eJFtGnTBhMmTMD48eNx+fLlWyMViMg29Ho9Dhw4gEmTJiEtLU3tcGyGTUB27O2338acOXPQtGlTREZGVtnZiET2qmHDhqhVqxaio6PVDsUm+ARgh2JiYuDn54e5c+di7ty5iI+PZ+VPpIJt27YhNjYWy5cvVzsUm2ACsCNFRUUICQlBx44d0ahRI6SmpuKVV15ROywizWrZsiUmTpyICRMmIDs7W+1wFMcmIDuxZcsWhISEwGQyYeXKlRg8eLDaIRERZIewr68vWrRogd27d6sdjqL4BKCy7OxsPProowgODkbfvn2Rnp7Oyp/Ijuj1emzZsgV79uzB+vXr1Q5HUUwAKlqwYAFq1aqFM2fO4ODBg1i3bh2HdhLZoY4dOyI0NBRjxozBjRs31A5HMWwCUkFycjKefPJJnDlzBlOnTsWcOXO4aieRnTOZTPDx8UFQUBC2bt2qdjiKYK1Tyd544w00btwYJpMJCQkJmDt3Lit/Igeg1+uxYcMGREZG4scff1Q7HEXwCaCSREdHIzg4GOnp6Zg7dy5efvlltUMionIYNGgQdu7cifT0dDg5OfaCyrz1tLGioiI888wzePjhh9G8eXP88ccfrPyJHNjatWthMpkQGhqqdigVxgRgQ5s2bULNmjURERGBH374Abt370b16tXVDouIKsDFxQVr1qzBunXrsG/fPrXDqRA2AdlAZmYmgoODsXfvXowYMQKrVq1y+EdFIiqpd+/eiI2NRWpqqsP242k7AQgBnDoFHDkCxMYCmZmAmxsQEAB06AB07AiUcVjm/PnzMX36dHh7eyM8PBwdO3a0UfBEpKbc3Fz4+PggJCQEixcvlj/Mz5f1ydGjQFyc/P6++2R90qkT0KIFYE87iwktys8X4quvhGjcWAh3dyE8PISQ6UC+3NyEqF5diBo1hHjjDSFSU+9ZzLFjx8TAgQNFYWGhSEpKEv7+/kKv14vXX3+9kn8hIlLD//73P6HT6UTczz8L8eqrQnh6yrrDza1kneLhIYTRKMQDDwixbJkQhYVqhy6EEEJ7TwDHjgHDhgGXLwN//mn5eFdX+Vq8GBg+HACQkpKCkydPYvr06Thz5gwCAwOxf/9+tGzZEhEREWjQoIGNfwkisgtCYFbz5piakAAPZ2egoMDyOe7uQKNGwLp1QKtWNg/RHG0lgA0bgNGjgbw8mZfLwmgEXngB+PxzPPb449i9ezecnZ2Rn58PAJg1axZmzJhhg6CJyC6ZTMALL0CsWwddbm7ZztXpgGrVgO+/B554wjbxWROGZhLAjz8CgwbJyr+8jEb81rcv2kRGlpgOrtPp0KNHD+zatUuBQInIIYwfD6xaBZS18r+d0Qhs2wZ066ZcXGWgjQSQlgY0awZkZVW4qFwATwDY89f3jRo1QocOHRAUFIRp06ZVuHwicgBbt8om4YpU/jfVqgWcOweoMERcGwlg2DAgPNy69jkrXAbQzNkZBUJgyJAhCAsLU6RcInIAublA/fpARoYy5bm5AaGhsp+xkjnm4NWySEkBtmyxrvJv1kw2Ea1aZfaw2kYjrn79NQoKClj5E2nNmjVyeKc5990n+xxzcoCkJGDkyNKPvXED+PZb4No1RcO0RtVPAIsXW9/h++WXcgyvBYbcXHgsWACdPY3nJaLKMXeu5RGEX34pbzrr1AFCQoCvvjI/4kevB1auVDZOK1T9BLBli+VsDQAjRsiJYDt2WFfu8eNAYWHFYiMix3L9OvD77+aPMRqBIUOAd9+ViWLfPtkEPXp06efk5gKbNysbqxWqdgIQQs7Gs8TTE5g5E3j1VevLdnOzrmwiqjqOHZMVvDkPPAAUFQG//fb3z44fBx580HLZlaxqJ4AbN6xr+581C1i6FLh0yfqy9XrZv0BE2nHpkhz/b46HB3DnBvJZWfJG05zMzIrFVg5Ve4UyISyvu9G2LdCrF9C+ffnKJyLtsOaaz8m5e0hn9eqy+chS2dbUWQqq2gnAzU3eqZvTs6eclp2cLL/38AAMBtlh06FD6ecJAfj4KBUpETmCWrUs1ylnzwJOTnJU4blz8mdt21puMvbwqPSF4qr+PICAAPN/+GrVSmbr6dNlQpgwQU4gK42Tk3zMq1ZNsVCJyM6lpwO+vpablteskTeJL74ItGsHREQAXbrI1YdL06WL7DCuRFW7DwCQzTvm1uLPywOuXv37lZMj+w7MVf4A0KQJK38irfH2lk8BlkycKOuH1FSZDCZMMF/5u7oCffooF6eVqv4TQHy8bN+vyBpAd3J3B+bNA8aNU65MInIMc+YAH3ygbJ3i6iqbi/z8lCvTClU/AQByoaUDB4DiYmXKq15djgbw8FCmPCJyHOnpQMOG1i0nbw0nJ6BvXzlnqZJV/SYgAFixQmZYJRiNwNdfs/In0ipvb9kC4O6uTHlubsCiRcqUVUbaSABNmwL//rflCRyWVKsmM/WIEcrERUSO6aWXgKAgWXlXRLVqsvL39VUmrjLSRgIAZHv9O++UPwkYjcAjjwBhYfa1pycRVT6dDti0SQ4VL+9gEKNR9ieYWyjOxrSTAADg7bflSp81ali/2fvNnXtefVUO5SrjJvFEVEUZjXLtsAkTypYEXF3laqHffQdMnmy7+KygjU7gO/3xB/D++3L1PYNBdubcOb3b3V2u5/HYY8Ds2XIiBxHRvRw5Arz1FrB3r+zUvbODWK+X/YZ/bSOJ99+XSUBl2kwAN/35JxAZCRw8KF/Xr8vs3Lo10LWr3KtTpbY5InJAycnA9u1yQtevv8oJYzVqyP6CoCDZh2hH84e0nQCIiDRMW30ARER0CxMAEZFGMQEQEWkUEwARkUYxARARaRQTABGRRjEBEBFpFBMAEZFGMQEQEWkUEwARkUYxARARaRQTABGRRjEBEBFpFBMAEZFGMQEQEWkUEwARkUYxARARaRQTABGRRjEBEBFpFBMAEZFGMQEQEWkUEwARkUYxARARaRQTABGRRjEBEBFpFBMAEZFGMQEQEWkUEwARkUYxARARaRQTABGRRjEBEBFpFBMAEZFGMQEQEWkUEwARkUYxARARaRQTABGRRjEBEBFplJPaAZAFQgAXLwJHjwIXLgDFxYC3N9C+PeDvDzjxfyFpSGEhcOoUcOwYcO0aYDAADRsCHTsCvr5qR+dwWHvYq2vXgKVLgX//G8jMBJydgfx8mRBcXACdTiaD0FBg6lSZDIiqqpMn5bUQFiavBSGAggJ5Hbi6yq9r1QKmTweeew6oXl3tiB2CTggh1A6C7vDdd8CLLwJFRUBurvljnZzkBTFhAvDRR4CbW+XESFQZ/vwTmDYN+OYbWckXF5s/3miUCeGbb4CnnqqcGB0YE4A9KS6WFf+6dZYr/jsZjUC9ekBUFFC3rm3iI6pMFy4A3boBqalAXl7ZzjUageefB/77X/mUQPfEBGAvhADGjAG+/77slf9NTk6An5/sL6hZU9n4iCrT1auynys11fJdf2mMRmDsWOCLL5SNrQrhKCB78e23wIYN5a/8AdlklJIi20CJHJUQwIgRwB9/lL/yB+S1tHw5sGmTcrFVMXwCsAepqUCzZsD168qUZzQCq1cDAwcqUx5RZVq5Epg0Sbb/K8HLC0hMlP+lEvgEYA/mz5cdXObs2iXbQa9fl68zZ0o/NjcXeOstZWMkqgxCAO+8Y77yv3kN3HwVFcm2/tLk58sRdXQXJgC1FRXJNsr8fMvHTp4MeHrKl6VhnxcuANHRysRIt2zbtg3+/v544IEH8Omnn6odTtWzaxeQlWX+mJvXgKenHPCQlydHzpUmLw+YN08mFyqBCUBtcXEVa+csTX4+8OOPyperYSaTCZMnT8aPP/6IuLg4rFmzBmfMPYmRWfHx8TCZTMjKysLly5flDyMiytb0M2SIbELds8f8cVlZQFJSuWOtqpgA1Hb0KGAyWXfs7NmyY2zvXqBHD/PHFhUBv/xS4fDob4cPH0bz5s3RsGFDODs745lnnsEmdjCWW6tWrdC0aVNMnDgR9erVw9ixY5H/889lu1MfM0aO+bfEYJDXGpXAmcBqS0iw7o7njTfkFPiCAuCZZ4DNm4F27eT5pbgWE4PFc+YoGKy2DB48GM2aNbv1/aVLl1C/fv1b3/v5+eHw4cMlzokitCuGAAAJt0lEQVSLi8PWrVsrLUZHZjKZkJSUhAsXLkAIgRUrVuATAHWsLaBBA3kj9MILlo/Ny5MdwVQCE4DaCgutO+72iuabb4CRI4F+/cyOcb5+7Ro++eSTCgaoXW3atCmRAKwRFxfHv3kZFd/WBFqmJonRo+XTsDVNO0LIp2IqgQlAbd7ecgJXWT+cQlic4digdWtkHD9egeDodvXq1UNycvKt7y9evIh69eqVOGb48OEYPnx4ZYfmkHQ6HTw8PODm5oaMjAy4uLhA5+ICZGdbV8CzzwLWJlsXFw4DvQf2AaitbVs5bt+cGjWAPn3kGicGAzBqFNC9O7Btm9nTjhgMuHjxooLBalunTp1w7tw5nD9/HgUFBQgLC0NwcLDaYTmsMWPGICwsDKtXr8bkyZORmJgInz59rDs5KEgufWJu9M/tnJ1lkymVwIlgaktPlx9kc8NAfXzk6Ah/fzli6MwZ4N13gZ9/LvWUXIMBr7i6YnFuLjw8PNCuXTsMHToUY8eORXWulFhu27Ztw5QpU2AymfDCCy/gzTffVDukquWrr+SKnpZmxC9cKG+cnn3WunJdXOQKu5ZutjSGCcAePP44sHOnsmUajcDVq8goKMCSJUuwceNGnDhxArm5ufD29kbnzp0xatQojBgxAk7cU4DshTU3RGWl0wGDBsl1tqgEJgB78NNP8gOq1NR3V1c5PO7rr+96KzExEQsXLkRERATOnj2LwsJC+Pr6okePHnj++efx6KOPQq9nyyCp6Jln5LpY1g6QsMRoBHbsAAIDlSmvCmECsAdCAP37yw+ppSUhrHHffcC5c1atCHrkyBEsWrQIO3fuRFJSEnQ6HZo0aYLevXtj3LhxaNOmTcXjISqLK1eABx5QZm0sV1e5JlZYWMXLqoKYAOxFairQsiWQkVGxcoxGuZ9A//5lPtVkMiEiIgIrV67E3r17cfXqVbi4uKBly5YYMGAAXnrpJfj5+VUsPiJrhIXJ8f0VWR1Xp5O7hMXHcwRQKZgA7MnJk3IDjOzs8q1bUq2a3DZv/HhFwrlx4wZWr16NtWvXIjo6GteuXVOsQzkzMxNevCjJnH/9C/jgg/IlAb1eVvoHDsinCbonJgB789tvsj8gMdH6D76Li6z8V64Enn7aZqFlZGSU2qEcEhKC4cOHW9WhPHv2bKxfvx55eXmYO3cu+vXrZ7OYycH973/AuHHAjRvWz5Vxd5cj5r7/Xm4YT6UTZH8KC4X45BMhvLyE8PQUQj4P3P0yGoVwcxMiJESI1NRKD/P3338Xr7/+uggICBAuLi5Cp9OJevXqiVGjRomoqChhMpnuOmfjxo2iTZs2Qgghtm/fLp566qnKDpscTUqKEEOGyM+6m1vp14OnpxDe3kJ8/rkQxcVqR+0Q+ARgzwoL5Zo/ERHA/v3A5cty4bjq1YGOHYFeveSICW9vtSMFULJDuVu3bpg/fz48PT1vvZ+ZmYlZs2ahWbNmmDBhAmJjYzF79mzMnz8fdepYvQIMaVVqqnwi2LkTiImRncR6vRw22rWr7Pfq319OliSrcAC4PXN2BgYPli8H0KlTJ3Tq1AkAIISA7o6lKhITE5GdnY3u3bsDAPLz83H//fcjOTmZCYAsq10beOUV+SJFcMA32cSdlT8AXLlyBfn5+WjZsuWt73Nzc9GkSZPKDo+IwCcAqkSZmZlIS0uDXq9HWloazp07B19fX3h7e6OoqAjLli3Dpk2bEBoaimHDhtnXDOW4OLnpyL59crc1nU52MHbtKtdlatFC7QiJysyOrjCq6urUqYPq1asjLy8PP/zwA6Kjo2+tpaPT6WAwGJCcnIznnnsOISEh6s9QFgJYvx6YOfPvfRfuHJm1fr3sl/H3l0MWuTgcORB2AlOlKSoqwowZM7B69Wp07twZU6dORdeuXe95rLkZyhMmTEBAQIBtg01JAUJD5T4M1i7R4e4uNyhZvly2VxPZOSYAUsX169dLjBAyx9IM5fHjx8PX11e54M6cAR55RO4jW9Z9Glxc5FIcBw4AjRsrFxORDTABkMO51wxlT09PtGvXDkOGDKnYktcpKUCbNnJJjvJeGnq9fAI4eVIu5U1kp5gAyOGVNkM5MDAQISEh9+xQXr58OYKCguDv7//3D4UAHntMbjNY0e0DXVyAJ58EfvihYuUQ2RATAFU591ryul69eujevTuef/55dOvWDV5eXnBycsI777yD77//Hps3b0ad3buB559Xbllud3fZSfzEE8qUR6QwJgCq8u7sUAbkRDUhBPR6PUwmE+rWrYskDw+4njtnucARI4D33wcaNJBLFz/3nHxquJegIDmLm8gOMQGQpphMJowePRpr1qzB7R/9BwEcAuBuqYBevYAlS2QSOHwYuP9++fOUlHsf7+YGnD4NNGpU8eCJFMaZwKQper0ely5dgk6nKzEKqcc9Zi7f0z//KecFHDok+wxSUkqv/AG5nMe+fRWMmsg2mABIc5YtW4bz588jMjISjRo1wuDBg/HyQw9ZvvvX6+UifLVqyWW7L1wA5s+Xd/mlyckBDh5UMnwixbAJiAiQE7iioswfc//98m4/OhoYMECu1rppE/DLL8CMGaWfFxwsjyOyM3wCIALk2j6W5OXJ/86fLzt/09PlDmyWNrSp7CUsiKzETyYRIEf0WJKZKZt9bn9otvQArdNZVzaRCpgAiAC5qqfRaPm45cuBl1+W/QBeXsDUqcCWLaUf7+EBBAYqFyeRgrgaKBEgl3S2xqxZcnmHs2flPrXr1gEffVT68YWFMrkQ2SF2AhPd1L49cOyYsmV27w7s3q1smUQKYRMQ0U3vvy+Xb1CK0Qi8955y5REpjAmA6Kannwa6dZOTtyrK1VUOFX388YqXRWQjbAIiul1qKtC6NZCWJnf6Kg+DAfD1BU6ckB3FRHaKTwBEt6tdW27mUquWXNK5rFxdZeW/fz8rf7J7TABEd2rSRG7m0rdv2foEjEZg4EDg+HHAz8928REphAmA6F5q1QLCw+Uwz86d5Xo/d25hqdPJn7m5ydE+4eFAWJjcEpLIAbAPgMgaCQmyWefAATkbWKcDGjaU6/137crZvuSQmACIiDSKTUBERBrFBEBEpFFMAEREGsUEQESkUUwAREQaxQRARKRRTABERBrFBEBEpFFMAEREGsUEQESkUUwAREQaxQRARKRRTABERBrFBEBEpFFMAEREGsUEQESkUUwAREQaxQRARKRRTABERBrFBEBEpFFMAEREGsUEQESkUUwAREQaxQRARKRRTABERBrFBEBEpFFMAEREGsUEQESkUUwAREQaxQRARKRRTABERBrFBEBEpFH/DyzEspiw5LuqAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"length = nx.get_edge_attributes(RNA2.G, 'length')\n",
"pos = nx.circular_layout(RNA2.G, scale=1.5)\n",
"nx.draw_networkx(RNA2.G, pos, font_color=\"white\")\n",
"nx.draw_networkx_edge_labels(RNA2.G, pos, edge_labels=length)\n",
"nx.draw(RNA2.G, pos, node_color=color_map, node_size=500)"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 2",
"language": "python",
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.15rc1"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment