{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Arcs coverage" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import coolbox\n", "from coolbox.api import *" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'0.3.5'" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "coolbox.__version__" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "arcs1 = \"../../../tests/test_data/bedpe_chr9_4000000_6000000.bedpe\"\n", "arcs2 = \"../../../tests/test_data/pairs_chr9_4000000_6000000.pairs\"" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "tags": [ "nbsphinx-thumbnail" ] }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA/QAAADUCAYAAAAyYWjJAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABQ90lEQVR4nO3dd5wkVbn/8U9tgmVhlyXnVEoWlWMEQUBMgFfEQFRAUEtEVH4i4DViRK9i9ogYURRQVBTkGhC5JMFCJCwSDpklLbvL5ty/P55T27VNz0xP6uqe/r5fr3rVdHV195mume566jznOUmtVkNEREREREREusu4qhsgIiIiIiIiIoOngF5ERERERESkCymgFxEREREREelCCuhFREREREREupACehEREREREZEupIBeREREREREpAspoBcRERERERHpQgroRURERERERLqQAnoRERERERGRLqSAXkRERERERKQLKaAXERERERER6UIK6EVERERERES6kAJ6ERERERERkS6kgF5ERERERESkCymgFxEREREREelCCuhFREREREREupACehEREREREZEupIBeREREREREpAspoBcRERERERHpQgroRURERERERLqQAnoRERERERGRLqSAXkRERERERKQLKaAXERERERER6UIK6EVERERERES6kAJ6ERERERERkS6kgF5ERERERESkCymgFxEREREREelCCuhFREREREREupACehEREREREZEupIBeREREREREpAspoBcRERERERHpQgroR0iSJO+uug3SHjrWvUPHunfoWPcGHefeoWPdO3Sse4eOdXMK6EeO/sB6h45179Cx7h061r1Bx7l36Fj3Dh3r3qFj3YQCehEREREREZEulNRqtZF9wiS5AthoRJ+0O2wEzKq6EdIWOta9Q8e6d+hY9wYd596hY907dKx7x7m1Wu3cqhvRaUY8oBcRERERERGR0aeUexEREREREZEupIBeREREREREpAspoAeSJNkySZIFSZLUkiRZt4X9D0uS5KYkSRYnSfJ0kiRXJEkypWGfNyZJcluSJEuSJJmRJMnhTZ5nWpIkP0qSZE6SJM8kSfLzJEk2HMnfTdY0mGOdJMmLkiT5UzzGs5Mk+UuSJC9tsp+OdQdIkuS4eFwbl2yAx7V0bHScO8dQjnWSJC+Ox+beJEkWJUlyV5Ikn0ySZO0m++pYd4ih/l+XHj8uSZI8PuaQJvfrWHeI4RxrnZdVa5jHrvJzrdjWkxu2vSQ+xxVJkqxV+h0HjBNkTUmSTEiS5IwkSe5JkmRpkiSPJElyTguP0/lZiyZU3YAO8WVgATBloB2TJDkR+BbwJeA0YDpwAKX3MkmSVwC/Br4DnAIcBPwiSZI5tVrtT6WnuxDYCTgRWAWcDfwW2GfYv5H0paVjnSTJ1sBfgJuBd8TNpwF/SpJkj1qt9mDcT8e68xwALC7dvm+A/Qc8NjrOHWswx/pwIMWOyT3AHsBn4vrNxU461h1rsP/XhROBLZvdoWPdsQZ1rHVe1lEGe+w68lwrSZIXAFcAOfCmWq22NEmSwTyFrOlHwKuATwP/AbYGdm3hcTo/a1WtVuvpBTuQs4EPAzVg3X723QiYD7xrgOf8X+DKhm2XA9eUbr88vt6+pW0vidsOrPp9GYvLII91BqwE1i9tmx63vVfHuvMW4LiBjmuTx7R0bHScO2sZ4rHeuMm2d8fn2VbHujOXoRzr0mOnA08BJ8TnOKThfh3rDlqG+H+t87IOWIb6f9op51pxn5Pjz7vFz41rgCnD/R17fQFeBywHdh3k43R+Noilp1PukyQZD3wTOIvWprt4W1z/pJ/nXAvYH7io4a5fAi9PkmRavP164IlarXZ1sUOtVrsRuD/eJyNoCMd6IrAC680vLIjbkvicOtbdb8Bjo+M8NtRqtaeabP5XXG8COtZj0GeAa4G/Nt6hYz1m6Lysu3XUuVaSJDtiGQP3AQfVarWFTXbbJUmS/4vDO+5OkuRNrTx3D3snFnDPGOTjdH42CD0d0GNXBtcGvt3i/i8F7gJOiOM/lidJ8o8kSfYq7ZNiH1D/aXjsndj7vWO8vXOTfYr9dm6xPdK6wR7rXwOLgK8kSbJJkiSbAOcAc4CL4z461p0pJEmyIrEx0u8ZYN9Wjo2Oc+cazLFuZi8s9e6ueFvHunMN6lgnSbIHcDyWkdWMjnXnGsyx1nlZZxnsZ3InnWttj138ewx4ba1Wm9fHfhcCvwMOA24DLk6S5PktPH+veilwd5Ik30qSZF5iNWwuSZJkiwEep/OzQejZgD4WQvgMcGqtVlve4sM2w8ZffAw4HXgDsBC4IkmSTeM+0+N6bsNj5zTcP73JPsV+05tslyEayrGu1Wozsat+bwaeiMth2Id80dOnY91ZHgM+Drwd+9/8B+CTJPlQP49p5djoOHeeoRzrNSRJshnw38D5pRM3HevOM9Rj/U3g27Va7d4+7tex7jxDOdY6L+sMQ/o/7bBzrVOBDbGe+WbPUzivVqv9T61W+yPwVuBe4MwWnr9XbYYNV3gBcAR2odUBv0n6L0yg87NB6OWieJ8D/lGr1S4fxGPGAesCb63ValcAJElyHfAgcDL2YVaoNTw2abK9cZ9iv2bbZegGfayTJNkc+BVWEOXEuPl9wGVJkuxVq9UeKu2uY90BarXa/2JjqQp/jOlYH0uS5Ou1Wm1VXw9tsq3ZsdFx7hDDONYAJEkyCUvRWwA0O+HUse4QQznWSZIcgQV5b2jlJRof3mS7jnUbDPH/WudlHWCon8kddq71F6yo3xeTJDm+FgdaN/Gb1S9Wq61KkuR3WGAvzSVxeWOtVnsaIEmSx4C/Y+/3s4ZElej8rEU92UOfJMlu2JiOTydJsn6SJOsD68S7pyVJMrmPh86O66uKDbFnJ6derbG4KrR+w2OL23NL+zXuU+w3t8l2GYJhHOvTsAteb6nValfEE4U3Y4VaihROHevO9ytgA2C7Pu5v5djoOHeHgY41ALFH4KdY4aODarXanNLdOtbdoc9jnSTJRGw2k7OBcfEzf2q8e0qSJOvFn3Wsu8NA/9c6L+tcrXwmd9K51u+A9wPHYp8ffXmyye3NW3j+XjUHuK0I5qNrgGX0X+le52eD0JMBPfBcbMzF9dhBnkN9bPUjWKpeM3diV3IaU0QSbBwmQMCqOTaOydg57nN3vP2fJvsU+zUb5yFDM9RjvTNwRzlFv1arLQPuwMbsgI51N+nrCmwrx0bHubsMdLX9HOCNWG9B43HRse4uzY71FGAr4KvUP/P/He/7JfVCiDrW3aWv/2udl3W+/j6TO+pcq1arfQcbonlakiSn9rHbJk1uP9bK8/eoO/vYXv4fbUbnZ4PQqwH9NdiYnfJSXI07CLu638wfsD/A/YsNsYKiI54w1Gq1pcDfeHb6zeHA9bVa7Zl4+4/AZnH+xOK5XgTsEO+TkTHUY/0gsHtMzQVWV9PcHXgAdKy7xJuxWQ0e7OP+AY+NjnPXGOhYkyTJmVgPzDG1Wu2axvt1rLtGf8d6Ac/+zD8y3vdR4GjQse4iA/1f67yscw34mUwHnmvVarVPAN8D/idJkmOa7LK6qn2SJOOwC8Q3tvr8PegPwB5JkmxU2rYv1tn27+YPAXR+NjhVz5vXKQtN5pfEUriuatjvt9iVuGOBg7ExIE8B00v7vAKbcuNrwH7Al7ArRa9peK4rsKkxDgMOxSq1/l/V78VYX1o51tjJwHLgsnicD8H+6ZcDz9ex7rwFq5Z7OjYFySHA+fE4v7+v49zqsdFx7qxlKMcaOCru8yPgZQ3LxqX9dKw7aBnq/3XDc2xH83nodaw7aBnGZ/hv0XlZxx+7ZsePDjnXojQPfbw9DrgESwt/Xdx2XNzvPmw4wOvi770CeGHVx6BTF2zI00NYpuwbsO/ih4E/9/e30erx1P92/P2qbkCnLDQP8m4ELmrYb13gu8DTwGKsiMbzmjzfocDtwFIsneOIJvusj51czgXmARcAG1X9Xoz1ZRDH+lXA1dgYvdnYScJ+OtaduQCfjx/Qi+L/Zg68vWGfZse5pWOj49w5y1CONfDj+H/fbDlOx7ozl6H+Xzfcvx1NAnod685ahvEZrvOyLjh2/Ry/ys+1aAjo47a1sCBzATb12nFxv5cA1wJLsAr3b676/e/0BXgOcDk2A8Uc7Pt4esM+Oj8bxpLEX1IaxJSfedgVnr9X3R4ZPTrWvUHHuXfoWPcOHeveoWPd3XT8pC/62xi+Xh1D34oXAbfrD6sn6Fj3Bh3n3qFj3Tt0rHuHjnV30/GTvuhvY5jUQy8iIiIiIiLShdRDLyIiIiIiItKFFNCLiIiIiIiIdCEF9CIiIiIiIiJdSAG9iIiIiIiISBdSQC8iIiIiIiLShRTQi4iIiIiIiHQhBfQiIiIiIiIiXUgBvYiIiIiIiEgXUkAvIiIiIiIi0oUmVN0A6W7BuwQ4Dngvz75A9Hvg7DTLl7S7XSLdxvkwGdgS2ArYHJgGTC0t68Vt68Xb46tp6ZDVgAXAfOCZuJ7XsDwBPBKXeXmW1qppqoiIiEh3UEAvw/UG4H3x5/OBC4F9gQ/G+xYDX6qkZSIdxPmQAOtjAXuzZcNBPuUSnh0Ql5fGwHkpsApYWVrXSrfL28ZjF+iKdePP44CJ1C8uNF5wmNZw31Rg3bhtPWCLFn6/Bc6HIrhvXJ7Ms3TVYN4sERERkbEoqdXUASJDE7zbDLgIWAc4K83yS0v37Qr8ELtodFKa5TdW00qR9nE+jAM2o3nAviUwpZ+HrwBmYgHro8Ac1gzKy8H5/DxLl43ObzE6nA/jsaB+apNlPexiR/m9W6efp1tO/b1qXB7ttvdGREREZKgU0MuQxFT7bwEvBf4GfCTN8lrDPidgqfiPA29Ls3xR2xsqMkqcDxsBuwK7ATsDW2M9z/1lPi2keRD6CPCEep1NzGaYzpoXQ1rNZqgBT2Hv6T3AHcAM4CG9vyIiIjLWKOVehuqFWDA/D/hCYzAf/RjYD9gFOBS4oE1tExlRzoep2N/xblgQvyuwSR+7F8Hk6h7j0s/PaFz4wOJ7NDsutzbe31BvoHHZHDs2mwB7lh620PkwAwvuiyD/CR0PERER6WbqoZchCd59BXglcG6a5ef2s9++wFex9Ng3pVm+sk1NFBkS58PawE6sGbxv02TXhVhQWASIDwAz8yxVEcgKxdT+zbBjtjP1LIpmF2BmUw/u7wBm5Fk6tz0tFRERERk+9dDLoAXvtsYK3y0DfjXA7tcAD2PpyPsBfx3VxokMgvNhApBSD953i7cbZ2xYBtxFPfC7A3hYKdydJ8/SlVhWxKPA9cX2hiESxXoDYJ+4FPvNpH6h5nbgP3mWariQiIiIdCT10MugBe8+ArwNuDTN8rNa2P9twEeAW9Msf+dot0+kL86HicAewN7YsJGdgEkNu60CAmsG7yHP0hVtbKqMsjhOf0vqwX1RC2Hthl1rwP1Y6v/1wI15ls5vY1NFRERE+qSAXgYleDcB+BNWmfrINMvvaeExk4ErsArfh6VZ/tDotlKkzvmwOfByYC/gJTy7evrDrBm836W0+d4U0/W3Z80g/7nYlH2FVVhwf11c7lamhoiIiFRFAb0MSvDuxcB3gfvTLH/rIB73GeD1wNfTLD9/tNon4nyYhPW+7xWX7Rt2uQ8LxP4B3JFn6bz2tlC6Sfx72glw2N/TC1hzSMZsrOf+OuCGPEufaXcbRUREpHdpDL0M1n5x/fdmdwbvtsDSVpcCt6RZvrC0/+vj4xXQy4hyPmxFPYB/EWumTS/CgvfrgOvzLH28/S2UbhXntL8tLj92PqyLZXq8HBu6sQlwcFxqzoc7qPfez1DvvYiIiIwm9dBLy+Lc87/HKkgfl2b57Q33vQX4EPUxyY8CZ6ZZPiN4tw5WEG8C8No0y2e3tfEypsRK9EWP6V5Y0cWyu6kHVbdq/LuMhjgOf3vqf4cvBCaWdnkGuAG4Fuu91+eeiIiIjCj10Mtg7IgF87OwMcdlp2GF8gBuxnqttgJ+GLw7Nc3y64J3N2I9WvsAv2tPk2WscD5sAhyIBU57smYxu/lY4FT0ws9qfwul18Q57O+Ly8+cD+tgGSJF7/0WwGvjgvPhP9jf6FXAnfHxIiIiIkOmHnppWfDu3cC7gUvSLP98aftewDeA5cAn0yz/U/BuEvBh4DDsAsDbsGDso8DVaZaf2u72S/eJAdIBWDrzi4Ak3lUD7qTeC39HnK5MpCPE3vutWXMoSPki1P3A5cAfNQxEREREhkoBvbQseHcB1kt/Sprl18Vt6wEXYj3y30yz/Cel/ccB3weej524fh2rdr8cODDNcs3tLM8SK42/BAvi9wfWinctA/4P6928Ic/SOZU0UGQInA9rYZkl+wCvBqbHu2pAjn1G/jXP0oXNn0FERETk2RTQS0tisbtLsQJjB6ZZvixuPwl4J1Yw6oQ0y1c1PG4b4BdYUPYOrNd+D+AjaZZf2b7fQDqd82FH4CCseOKGpbv+BVyGBTua/1u6nvNhAvAy7O99P+o990uxC1aXAf9Q1omIiIgMRAG9tCR49xbgDOAvaZafEbetjfUqTQXemWb5rX089gPA24H/xYqVvR/4bZrln21H26VzOR82Bl6H9cY/p3TXQ1hQ88c8S2dW0TaRdohV81+F/Q/sWbprNpbRdBk2172+rEVERORZVBRPWrVrXN9c2nYwFszfgfXQ9+VC4GhsDP2fG55PeozzYTKWSn8w8GLqc3rPwwKYy7Ex8QpgZMzLs3QBViT0d86HLahf4NoWOCouwflwGXBFnqVPVtZYERER6TjqoZeWBO9+ifWgHp9m+W1xfPzF2EnnR9Ms/9MAj/888Brg58CRcfM+aZYvHcVmS4dwPozDgvdiXPzkeNdybFz8ZcB1eZYur6aFIp0jFtTbFft/eS0wLd5VA27C/l/+lmep6pCIiIj0OPXQy4Biav0OwCosZR5gZyyYfwqbX34gv8QC+ldj1Z1T4LnA7f09SLqb82Ea8GbgLVjhxMK/saDkL3mWzquibSKdKman3AHc4Xw4B5sG72BgX6xg5EuAM50PlwMX5Fn6QFVtFRERkWopoJdW7IilRd9T6lHfL67/lmZ5K4WbbgOexIK6+7CAflcU0I9JzodtsEyM/6Jepf4RLJ3+8jxLH6mqbSLdJGatXA1c7XyYig1dOgh4ATYt6GHOh2uw7Kd/aqiKiIhIb1FAL63YJa7vLG3bL66vauUJ0iyvBe/+DryVeoC3Sz8PkS4T04RfiNVL2Jf6nPHXYcHGjQo2RIYuZrNcAlzifCjG2L8BeEVc7nY+/Az4s4aviIiI9AYF9NKKooDdDIDg3dZYCv4C1iySN5CrsIB+q4bnlS4Wp+B6FRbIF8d0OdYb//M8S++rqm0iY1WepQ8CX3A+fBcb0vI2LJvqLOD9zocLgUs0pEVERGRsU0AvrWjsoX9lXF+TZvmKQTzPzcBCLO1+PLB98G5ymuWLR6aZ0k5xuq03AUcAm8bNc7FiiRfnWTq7oqaJ9Iw8S+cC5zkffopVyD8Gu+B6MnCi8+FSbJy9hrmIiIiMQQropV/Bu3WA7YCVwD1x835xfdVgnivN8uXBu2uwqs2LgElYYbym89dLZ4pTax0BHAqsEzc/iKXVX5ZnqWYuEGmzPEuXAZc6H34PvAwL7F+K9dy/1flwFZYxc0tljRQREZERp4BeBlIUxLsrzfJlwbsNgD2wlOrrhvB8f8cC+knx9q4ooO8KMZA/ATgEy7AA+CfwM2zKuVVVtU1ETKxTcT1wvfPhOdhQmNdj00Xu73zIge/lWTqY4VIiIiLSoRTQy0CKMdFFuv1eWIB/Q5rlQ5kD+VpgBbAhNgZf4+g7nPNhM+CdwBuxQH4VcAXwszxL/1Nl20Skb3mW3gt82vnwbWJPPeCAc50PN2KB/b+rbKOIiIgMjwJ6GUgxfn5GXD8vrm8aypOlWb4weHcHlgq6Nqp037GcD5tggfyh2GfFKqzQ3Xl5lj5UYdNEZBDyLJ0FfMf5cD42neTRxPnsnQ83YIH9bVW2UURERIZGAb0MpLEgXnH7jmE85wxsDuX1ge2Cd+sMsbdfRoHzYWPgeKzg3USgBvwv8P08Sx+osGkiMgx5ls7Heud/iQX1R2Lj7V/mfLgO8HmWzujvOURERKSzKKCXPgXvpgDbYuPlQ/BuIlbEDuCuYTz1nViQuASbq3wn4F/DeD4ZAc6HKVggfxT1Ggd/xgJ5TT0nMkbEqey+63y4ACuedwQ2nGov58OVwDfzLH24yjaKiIhIaxTQS392wgLue2NBvJ2xHtsH0yxfOIznLXr7i7+/XVBAXxnnw3gsrT4DpsfNVwLnxjG4IjIG5Vn6DPDtGNi/HRtnfwCwr/PhImx4jeaxFxER6WAK6KU/jQXxitvDTcl8EJu2bn2syJoK41XA+ZBgvXIfwOatBptx4Kt5lt5eWcNEpK3yLJ0DfCOm4mfAG7BMnUOcD+cBF+dZurzKNoqIiEhz46pugHS0Isgr5p9vHE8/JGmWrwL+AyzFCuPt0P8jZKQ5H54LfBv4Ovb+zwROB05QMC/Sm/IsfTLP0rOw8fU3AVOBU4GLnA8HxIuAIiIi0kEU0Et/torroqJ5Y4/9cPwHWAZMBrYK3ulEsQ2cD9OdDx8DLsCqXC8Avga8Jc/Sv8Y5rEWkh+VZejdwEhbMPwhsDXwJK6i3Y5VtExERkTUp5V76s01cPxy8mwQ8BytmN5yCeIUZwErsotI6wAbA0yPwvNKE82Eclkb7AazXbSVwMVbw7pkq2yYinSde3Ls6Vr8/DHg38ELgZzE13+dZqtlJREREKqaAXpoK3q0DbIRVuH8c2Bkb737fCE0xV4zDL7JEtkEB/ahwPqTAmdhUgQA3AmfnWfpgZY0Ska6QZ+kKLOX+cmx8/duw8fUHOh++DFylzB4REZHqKOVe+lKk2z8ax7yPZLo9wCPAQuxvcAKW0ikjyPmwtvPh/Vh6/QuA2cDHgPcpmBeRwcizdEGepf8DHItdkN0E+DLwVefD5pU2TkREpIcpoJe+FOn2ozF+vlwYbxlWGG+b/h8hg+F8eAVwEXbyPQ74NTZO/gr1ponIUOVZeidwHDamfhGwD3Cx8+Edzgdl/YmIiLSZvnylL6vHz8d1UeF+uFPWlc0A9kcB/YhxPkzHqtUfGDfdDXwhz9LbqmuViIwleZauwtLw/wZ8CHgNcApwkPPh0zHoFxERkTZQD730pUiBfyh4txaQAquwAHGk3Em90r1S7ofJ+XAA1it/ILAYOAd4u4J5ERkNeZY+lWfpR7Fg/lGscOpPnA/vcT5MrLZ1IiIivUEBvfSlCLAfBnbE/lbuS7N8yQi+xgzqKfdba+q6oXE+THM+fA5LgZ0O/BM4PM/Sn+dZurLa1onIWJdn6XXA4cAvsO+Kd2GBvaa4ExERGWUK6KUv5YB++/jzvSP8Go9iYzATbOq6jUb4+cc858NewIXAa4ElWFB/Up6lMyttmIj0lDxLl+RZ+hVsertHsQvBP3U+HBunzRQREZFRoC9ZeZaYYr8hNlf5Ezy7QN6ISLO8hl0wWAZMAlQpuUWxgv3pwDewCyH/Bo7Is/SiOL5VRKTt8iy9GTgSK8Q5AXg/cK7zYYtKGyYiIjJGKaCXZorA+vFYjb7cWz/SHsbmup8E6ISvBc6H5wI/B94KrAC+Bbwrz9JHKm2YiAiQZ+miPEu/gI2tfxqbNvOXzofXVdowERGRMUhV7qWZIrB+NK5XF8gbhdd6CAX0LXE+JMAbgY9g79d9wMfyLB3JQoUiIiMiz9LrnA+HAx8FDgA+63zYE/ifPEuXVds6ERGRsUE99NJMEVg/FgvVjWYPvQL6FjgfJgOfBj6GvVe/xSrYK5gXkY6VZ+lcbCrNz2PDqw4Dfux80MwmIiIiI0ABvTRTBNYzsfHZk4F5aZbPG4XXKo+hV0DfhPNhB+CnwEFY4btP5Fn62TxLl1bbMhGRgeVZWsuz9BLgeOARrGDez50PB1bbMhERke6ngF6aKcbQz2R00+2L5y166FUUr4Hz4RAsmN8eS7F/R56ll1fbKhGRwcuz9C7gGOCv2MwmX3Q+fMT5MKnalomIiHQvBfTSzJZx3Y6AfjbwDDAe2Cp4p79JwPkwyfnwceBTwNrAZcCxeZbeV2nDRESGIc/SBcAZwJexop5vA85zPmxWacNERES6lIInaaboKX+M+pR1o1JBPU5d9xB2YjcZzUWP82ED4LtYAbxlwFnAp/IsXVxpw0RERkBMwb8QeCd24XhXbM7651XbMhERke6jgF7WELybDEzH0uBnUR/XPppToj2KCuMB4Hx4DvAT4PnAk8DxeZZemmdprdqWiYiMrDxLZ2Ap+DcBGwDf09R2IiIig6OAXhqt7p2Pc9CXC+SNlplYQD+RHg7onQ/7Aj/EjsHt2Hj5u6ptlYjI6MmzdB7wfuDX2EXdzzofTnI+6PxERESkBfrClEZFyvtTcV0ukDdaygF9z6XcOx8S58M7gK9ghaL+F3hPnqWzqm2ZiMjoy7N0BfBF4EvAKiwV/+w4XaeIiIj0QwG9NNowrp+O6fcbYOPbRzO4nAmsxHpnNhxg3zElVnf+JHAKkADfAT6mKelEpJfEcfUXYZ+FC4D9gR84HzattmUiIiKdTQG9NNogrmfz7PT70TITu2gwsfT6Y57zYR3ga8Ah2PzyH8mz9IcaLy8ivSrP0huA44CHsfnqf+R82KHSRomIiHQwBfTSqOghLwf0o5luD/A49YC+J3ronQ/rAx54CfZen5hn6ZWVNkpEpAPkWfoAcCzwL2ATbFq73SttlIiISIeaUHUDpOMUPeRPs+Z89KMmzfJFwbvZ2BR5mw+0f7dzPmwCfBvYHntvT8qzdDRnERAR6Sp5ls5zPpyMja3fB/iu8+HDwI1YrZGp2EXgcQ3LeKCGjcVfGX9eGW8vBublWbqsvb+NiIjI6FFAL43KKfd7xp9Hu4cebFq8FwBjeryk82FbLJjfDAjAyXmWPtX/o0RExibnw1rYxeMtgfWxQH0qsB4wDTtPmQikwJ+Ax7Ax9sN5zaXAPGA+8EzDel5pmYV9Nz2ZZ+loDjsTEREZMgX00mh1UTzaM2Vd4YHi9YN340Z5zH4lnA+7AN8ApgO3Ah+MUzaJiIxJzocEC8y3BLZqsmzcwtMsx4L4DYCtgQex74xl1HviV5WWhOY995OxiwVrxddt5bUBljsfZmLBfePyqHr8RUSkSgropVG5h77oLX+iDa/7GPVK9+vH1x8znA8vBL6OpYpeB5yeZ+nialslIjJy4jRzOwO7xmU7LGif0s/DVmIXjR/FPveb9ZQXy39hY+trwC/zLP3pENqYYIH9etSzAcpLkRkwFfsO3Ar7Xtw2Ls2e80ksuL8XuAOYATyoXn0REWkHBfSyWvBuHGsG9OXx9KNtNlYYbwKWJTBmAnrnw/OxnvnJWMroJ/MsXV5tq0REhs75MBF4Dha47xaX7WlebHcRzXu3HwGeyLN0ZYsv+43YU346cIrzoZZn6fmDaXecRWRRXFq6WB1nJCkyDBozDbbACvdtQn2YGsAi58OdWIBfBPmPaxYTEREZaQropWw9LC1xAZbiWE6/H21PYwH9eMbQ1HXOh92Ab2LB/OXAp9RrIyLdJs4H/yLqwftzsYyqslXA3dQD2HuxoH3uSAWyeZb+yvmwDPgE8AHnw/I8S385Es/dz2suAu6JyxqcD+OxmijbADtRv8CxKeDiUpjjfJiBvT+3AzfnWbpkNNsuIiJjnwJ6KStPWTcZG2e4FKsMPNpmY6mXRQ9913M+7IwVwFsH65n/tIJ5EekGzodJWKHSveLSbC74h6gH73cAd+VZunS025Zn6aXOhwnAR4EPx6D+16P9un20ZSU2XOBR4Ppiu/NhQyywL2cwTAf2jgvAMufDzdgwrOuBB9SDLyIig6WAXsrKKfare+fTLG/HCUbRQ78WY6CH3vmwI/AdYF3gSuATg0grFRFpO+fDltQD+BcDa5fuXgT8E7gNC97vzLN0ftsbGeVZeklM+z8NODMG9ZdW1Z5GeZY+DVwdl2Ls/hbUg/w9gV2Al8UF4DHnw/XAtcBNMTNARESkXwropazZ+Pl2jWUvxtB3fcq982EHLJifip3M/XeepSuqbZWIyJpiL/yLqAfx2zTscg/Wc3wtcGun1f7Is/TC2FP/IeDjzocVeZZeXnW7mok970VP/p8AnA8bAC/F3vuXA5sDh8VlhfPhFqz3/lrgPvXei4hIMwropaw8Zr6cft8Oc6gXxduoTa854mIPl8cq9V8HnNFpJ8Ei0rtiT/EewMHAq7HaKYUFwD+wAPKGPEufbH8LByfP0p/HnvqTgU85H5bkWXpl1e1qRZ6ls4E/An90PozDZggoLq7sjl1seRFwChCcD5cBV3TDcRERkfZJajVd8BUTvHsfcDwWkM4FzgB+k2b559r0+v/EUhAvSLP8Xe14zZHkfJgG/Ajr5boR+FA7xpOKiAzE+bA1cFBctizddTdwDRbE396tQ4OcD+8G3o3NTZ/lWXprxU0aFufDVKz3/uXAK7Gp9MCm7LsJuAz4m9LyRUREPfRSVu6V36j0c7vMiuvN2viaIyKmrn4VC+bvBk5TMC8iVYpB4WuwIH6P0l1PAlcAl+VZGqpo2yj4Pva9dRhwjvPh+DxLH6q4TUOWZ+k84M/An+Owgr2wrIp9gZfE5Uznw5XYDCo3quiqiEhvUkAvZeWieDuVfm6Xx+O6q1LuY6rkWcDzsRPlD+ZZurDaVolIL4rTqO2DBX/7UP+eX4wV6LwM+OdYC/7yLK05H87GpovbG/hmDOrbeVF6VMQaLFcDV8eLNAdiF2leQD3rYpbz4Qrg0jxL76uqrSIi0n4K6KWsyqJ4AI8V7QjejUuzvFtOON+PnWAtBE7R+EYRaTfnw2TgDcAxWDV1sHnhb8CC+KvyLG3HFKSVybN0pfPhTOB72PCtrzofsrE013vsub8EuCTWbCkC+q2xY3+M8+E64KdArkJ6IiJjnwJ6KauyKB7AU9gJ6ESsQvzcNr72kDgf3ga8HSvo9+E8S++tuEki0kNipfTDgbdin5sAjwC/Bv6YZ+msvh47FuVZusj58EHgx1hhuc86Hz4y1jISAPIsfRT4vvPhPOx3fQOWmVEU1rvT+XA+8NdurY0gIiIDU0AvAATvEpr30Lcz5b6Yum5CfP25bXztQXM+7A18ON48K8/Sm6psj4j0DufDtsDRwCHApLj5dqxn9qqxGMC2Ks/Sp50PpwA/BPYDPgCcU2mjRlHshb8NuM358B3s4s7bsCyFzwMznQ8XYOn4KqInIjLGjKu6AdIxpmA944vSLF9CNSn3jQF9x4qpjp/F/oe+16lzH4vI2OF8SJwPL3A+fAX4FVYAbhI2vvpE4Pg8S6/s5WC+kGfp/cD/A1YCRzsfXlNxk9oiz9K5eZZ+H7vQ83ngIWwIxoeBPzgfTnI+bNjfc4iISHdRQC+F1QF88G5tYB1gOTYvcbs8TRcE9M6HtYEvY/M3/x34QbUtEpGxzvmwJ/ZZcx42jdly4DfAm/MsPTXP0ls0XnpNeZbeDHwl3vy482GHKtvTTnmWLs2z9BLgLVgwfys2JOOdWGB/pgJ7EZGxQSn3UiiPmV8d3KdZ3s4TxKKHflKpPR3F+ZAAZwI7Ag8Dn1RvmIiMFufDNljK+CvjpnnAxcCFY6GCextcDDwPeD3wP86Ht/fSLCTx++kq4Crnwx7AO7C/pTcDr3c+/Bj4uaZZFRHpXuqhl0JRTGkuaxbHa6ensfTICaX2dJrDsKJDS7C55tuZwSAiPcL5MM358GEsIH0lNu2cBw7Js/S7CuZbE7MWPg/cC2wDfCpemO05eZbemmfph7Hx9VdjmXgnAb9xPhwUp2AVEZEuow9vKUyJ6wVUM34eYA5W5X4CsG6bX3tAzofdgdPizc+oor2IjDTnwyTnwzHAb4EjsO/p3wFvyrP0PBU1G7w4Xd9p2Pfb/tjMJD0rz9L78yw9FXgvcDewCXAW8JM4tENERLqIAnopFAH0QqqZso40y8tj9jsq5d75MA34Enax4Zd5lv5vxU0SkTEkFrx7FdYj/0GsRsc/gKPyLP1Mr00/N9LyLH0Y+ES8ebLzwVXZnk4QZ2Y5Bvg0Nm3sLsC5zocvx6EeIiLSBRTQS6Ec0FcxZV1hblxvVMFrN1UaN78JVljo69W2SETGEufDzlixu7OBLYH7gVOAk/MsvafKto0leZZejU1lNw44y/nQcZlg7ZZn6ao8S3+PDSf7HjacbH/gYufD/3M+TOn3CUREpHIK6KVQfGmXA/oqxmjOievpFbx2X14HHAgsAj6WZ+nyitsjImNATK9/HzZ3/POxz78vAEfkWXqdqtaPinOBGcCm1IdQ9bw8SxfH6e7eBFyKnR8eCVzkfHh5pY0TEZF+KaCXQnkMfRFMz+lj39FUpJVOq+C1n8X5sClwerz5lTxLZ1bZHhEZG2JNjp8DxwMJ8Avg0DxLf51n6cpKGzeG5Vm6Aku9XwYc7Hw4oOImdZQ8S5/Ks/Qs4GjqFz6+6Xz4hPNhvWpbJyIizSigl0I55b4I7udX0I65cV15KmSs+PsprC1XY70WIiJDFnvlT8FSv7cHHgJOyLP0K700nVqV8ix9APhavPnfzoeOGeLVKfIsvRu72PQN7OLHf2G99XtX2jAREXkWBfRSKPfQl4P7dpsb150wbu9w4MVYpsJnlf4qIsPhfNgB+DE2FzhYqv2ReZbeWlmjetevsKKD04CP9+pUdv3Js3RlnqU/BY4CbgM2Br7ufDjd+bB2ta0TEZGCAnoplMfQFwF9FXOsF+P216ngtVdzPmwPvD/e/JzmfBaRoYoV7A8HfgbsCDyC9cp/I8/SpdW2rjflWboKq+4+H9gbGzsuTcSMhhOw3voVwFuB850PO1XZLhERMQropVAO4svBfbvNjeu1g3eV/H2WqtpPAn6fZ+lVVbRDRLqf82F9bGaM07DPlN9hU9HdVmW7BPIsfRL4Yrz5AedDR02X2kliNfyfAsdiszBsj81bf4yyG0REqqWAXgrNxtBXEdDPB1YB46ku7f71wJ5Yqv05FbVBRLpcTLH/CbAXMA/4SJxTflG1LZOSPwHXYt83H6i4LR0vz9K7gLcDFwETgA8Cn3Y+TKqyXSIivUwBvRSK4HkRa46nb7cFwErsb7PtAX2s4vvBePPreZbOa3cbRKT7OR/2An6EzSt/JzYV3ZXVtkoaxdooX8YKvx3kfNiz4iZ1vDxLl+RZ+iUs62QJcBDgnQ8b9P9IEREZDQroheBdQj14XolNobQ4zfJVFTRnIdZDX0lAD7wX2AC4BbisgtcXkS4Wx8sfhVVRnwL8GXhXTO+WDpRn6SNYsUKAM5wPEypsTtfIs/RvwDuBJ4A9gJ86H3astlUiIr1HAb2AjeucgPVQrBW3VTV9Ujmgb+vUdc6HnYG3xNf/oqrai8hgOB8mAv8NnIp9hp0LfDTP0iWVNkxa8WOsWOEOwJHVNqV7xOnt3gHcDmwG/MD58MpqWyUi0lsU0AusOWa+yvHzxeuupM1j6OOc82dg/xO/yLP03na9toh0v1j87tvAodjF0TPzLD1XFwa7Q56ly4AvxZvvdj5sUmV7ukmepU8D7wb+CEwG/sf5cKyK5YmItIcC+g4SvNsteHdO8O6Twbt2jkUrF8Srcsq64nWr6KE/CNgdeAr4XhtfV0S6XKn43Z7ALODEPEv/XG2rZLDyLL0OuBILSk+puDldJV4Q+QTwHWzY3vuBT6lYnojI6FNA3yGCdzsA5wH7AG8Azg3eteuLsFMq3Bev29aAPp5wvCfe/KYqUItIq+JQnR9QL373jjxLZ1TbKhmGc4DlwGudD8+tujHdJM/SWp6lP6ReLO9g4BwF9SIio0sBfQeIRenOBCZi0+c8AmyHjUtrh3JV+yor3AMspl7lfr02veahwObAfcAVbXpNEelyMeD7DvZZ9XdU/K7r5Vn6GPBrrJf5vRU3pyuViuXNBl4KfFlBvYjI6FEl186wG/BCYC7wcWBHwANHBe9+nGb5ilF+/WZj6CsJ6NMsrwXvFgEbAtNH+/WcD5OBE+PN7+RZWkVl/2GJJ0pT47Jeab0u9j8+HrtAUsOyH1ZhY3znx+WZuJ4HzM+zdGWbfwWRrhPT7L+L/b9dDZyRZ+nyalslI+SHwBuBfZ0Pe+RZemvVDeo2eZbe7Xx4L1YYcm/gi86H0/U/IiIy8hTQd4ZXxfUVaZbPA/4ZvHsA66V3wD9G+fXL4+bL6fdVKVLe21FH4Ij4OndgPWwdx/kwHssg2KrJsgWwzgi/3lzgUSxTpHGZpSJf0uucD9tiF13XB65DwfyYkmfpbOfDL7Be5vc5HzJ97g1enqUhBvUe2Bf4nPPhTF00FhEZWQroKxbT7YuA/q+lu/4KnBDvG+2AvtwrXwT0S4J3WwCPpVne7hOZ4mLC+qP5Is6HqdSHNXy7E07YYlXgrYFd47I7sBP16QSbWcmavezFsgBYgfXIr8RSSIve+knYsZ6G9eYX6/Ww9319LHOk0VznwwzsAsgM4I48S2cP8dcV6TrOh62xAGUD4EbgtFgQTMaW84G3YhfVX4wdaxmk2FP/Piyb5QDgM86Hj3VjNpyISKdSQF+9HbFe1tnAv0vbi4B+v+DdF9MsH80vv8aU+7WADwDHA7cE796TZnk7r6gX6f5TR/l1jsEC2JvyLK3kZC0G8CmwF/ASLIhv9ns/SfMe80ewNPkRuRgRp+/bkObZAFtjgf5ecSke8zhwG3A9cIPGEMtY5XzYAgvmNwZuBk7Ns3Rpta2S0ZBn6Xznw0+Ak4GTnQ/HdsJF326UZ+mdzoeTsXoTrwGWOx8+raBeRGRkKKCvnovraxqC9nuAJ4BNsdT7+0axDeU0+w2xCwxz4rYXAG8HfjyKr99oflyPWkDvfFgHODze/O5ovU4fr70eFrzvBbwcaJzv+GmsB7zoBZ+RZ+kz7WhbPMF6Ki7/Kt8XLz5sjl102C0uOwObxeXVcb97sOD+WuBWpSLLWOB82BQL5jfFLr5+MM/SJdW2SkbZhcBR2GeeeumHIc/S250PpwDfxKrfr3A+fE5BvYjI8Cmgr94ucX17eWMsDncHdvK4C+0L6PcE1sYuJnwM+DLwnuDdb9MsnzuKbSgrgtfR7KF/I5aNcHM7Ch7F4nv7YfPdv5Q1Z5h4GhuHez1wC/BUJ/YExTbNjMtfYHWP/nbYhamXYxcqnhuXdwCLnA9XApcBuU7epBs5H9YGvo5d7LwDOEXTW459eZYudj5ciFW7PxoF9MOSZ+ktzocPAt/AvoNnYlM+iojIMCigr96ucd1s3uIZ2JizXbGAaLSUx9A/J/786zTL/xa8uwF4GfBK4Hej2IayIqCf0u9eQxSD0CPjzZ+PxmuUXudFWG/E/tSL163C0nWLXux7OjGAb0UM0O+Ly8Wx4v4LqKfl7wAcEpcnnQ+XA5fnWTqaF6hERkzMTPko9tn4EPD+PEurLBoq7fVrbPjb3s6H7fMsvb/qBnWzPEtz58PpwDlA5nyYkWfp9VW3S0Skmymgr1DwbgqwLbAcCE12uTOud2ly30gqAs3xWEZADQs0wcbyvwwrzteugL5I9x+VgB7rKd8CeBj4v5F+cufD+sBbgMNYM53+NuBy4E/tSqFvt1gc7Ma4fM35sA12QeP12Ht+HHCc8+FO4JfYe6GUfOlkb8Uya5YAH86zdF7F7ZE2yrN0rvPhD9jn+ZHA5ytuUtfLs/Qa58O5wHuwyvfH5Fk6s+p2iYh0q3ED7yKjaOe4vifN8mZBTRHQ7xy8Gz+K7SgqqO+OBfULgVlx21VYj/JLgnfrjWIbyooT5v4quw/HMXH9i5FMAXc+bOV8+AiWTZFhwfxM4PvAYXmWHp9n6cVjNZhvJs/Sh/Is/S5wKPAu4LdYJsguwKeB3zkf3uF8WLfPJxGpiPNhD+D/xZtnKbOkZ10Q14c4H6ZX2pKx4wdYx8FU4Esxu0tERIZAAX21ioD+P83ujHPSz8SmGNtuFNsxMa6fi/1NzCdOHZdm+RxsfP8ELOBvhyKgH/EveOfD84A94mv8foSeczfnw9nAJcDbsAsR12HjLt+YZ+n38ix9aCReq1vlWboqz9J/5Vn6WazK8VlYmv4mwCnA5c6HDzgfGgsEilTC+bABcDZ2kfOCPEv/VHGTpCJ5lj4AXIN9J7252taMDfFi+seBR7FzodOrbZGISPdSQF+tbeO6v16f4r5t+9mnZX309BeB87bYyeti6nPBgxWBgtFP/S8UPdgT+91raIqx85fkWbp4OE/kfNja+fAl4CfYkIRV2EWCw/MsPSXP0pu6dWz8aMqzdFmepZcCR2DTI+bYsI+3A791PpzsfBit4RYiA3I+jAe+iE1PdwtWxEt6W1Fv5W3OBw1XHAFx+MppwDLgjc6HQ6ttkYhId9KXUrW2juv+em+L+7buZ58BBe/GYRXr9wne5cAppTT/idjFnc3i7WXYeNFCu8byF4rq0eODd+MapvMbsjhd3H5YjYBfDeN5pgInYr3xE7D36iLgl5qDvXWxh+Za4Frnw67AsdiFkeOwkzsP/DbP0pXVtVJ61MnYjB9PA2fkWbqi4vZI9f6JXWDfASv4eXW1zRkb8iy92/nwOWwI1unOh7vzLG1WJFhERPqgHvpqbRPXD/ezT3HfNv3s04rDsUr147D5dI8v3TcJm6puFbAUC3jLJ7BFQL8r7bEstiFhZHvpD8R+15vyLH18sA92Pkx0PhyFjQM/CstmuBQbH/8NBfNDl2fpjDxLT8eC+VuB6cCZwC+cD3tV2TbpLfHv7e3ASuD0PEtnDfAQ6QEx26qYbebgKtsy1uRZehl2kX0i8MU4zauIiLRIAX1FgncTsfHDq4DH+tn1kbjefBivNQF4d7z5i7g+PnhXVLefBBRfoEuAZWmWl1PFH8R6zTcN3m0w1HYMQjmgH8lx9IfE9R8G+0Dnw07AT4FTsSI+NwFH51l6lgL5kZNn6e3YFFGnY/UjdgC+4Xz4vIpRyWhzPqwFnBFvfjvP0lsqbI50nj9i3037xkwtGTlfBe7CZkM5oeK2iIh0FQX01dkYC1hnpVneXzpn0ZO8WT/7DOT5wHrA/WmWfwXrAZ0IvCTeP5F64LwEm0ZvtZjyXkyrt90w2tGq5dhJ0zhGqIfe+bAV9j4sxir3t/q4Cc6HDDgfKxo4EwvqT8qz9O6RaJusKc/SWp6lf8Wm/vs69jf5GuAi58MBlTZOxrrjsIDiXupjpkUAiBdvb8S+l15dcXPGlDjl6Rew7/5jnA/bVdsiEZHuoYC+OkWAPlDq9+qAPo6DH4p94rqYc/2ahu3lgH5ZXBqNVOp/K5ZhmQsjmXJ/UFxfmWfpon73jOIc6j/ExsuPw+ZNPzzP0qtV7G70xeJ552PDRf6JpeF/yfnwSefDOv0/WmRwnA9bY3UcAL6o2g3SB6Xdj5KYofVbrDbN6c6HpNoWiYh0BxXFq04R0PeXbk+a5YuDd88A04D1gdlDeK0icC8C+f8DTgJeES8STGLNgH45z1YU52tHQF/00I9Iyn08KSjS7S/rb9/SY14DfAKrLfA48Ik8S28ebltk8PIsfdT5cBLWY/9B4A3Ans6HU/MsDf0+WKQF8TPiNOzz5g9KtZd+/A3L9NrD+bBNN05J6nxYF9iqtGyBfddNxM4Li0K5K7Dv42I9GxsGWCyPx+KmI+lbwAFYrZ9XA5ouUkRkAAroq9NqD32xz7T4mEEF9MG7rbDp6OYD/46b743PuRk2/+sk7G9hZVyaBfTFWP4tBvP6Q1QO6Eeih353rN1PYj29fXI+jAPeFReAK4Cz8yydPwLtkCGKJ40XOR/+CXwW2BH4kfPho3mWXtP/o0UGtD9WuXw+mqJO+pFn6WLnw1+xi8SvBb5fcZP6FOuOvASrRVIO4KeN0EuscD7MZM0g/9/AnUMN9PMsfcb58E3gY8Cpzodr8yxdONDjRER6mQL66hRF7loN6HfCAvDBTueyc1zfkmb5SoA0y2vBu5uxNPSdsDnAAWZhac3NUu5HYix/q4qieEX2wHAVGQp/6+8kw/mwNjZ1TjGn/NeAXyi9vnPkWXqf8+GdwCex3ptznA9fAy7QcZKhiMM3PhxvfifP0qFkQUlvuRIL6PehgwL6eEF6d+zi1F7YVLPN0taXAo+yZiC+AOuJL3rjV1HvrZ+AfRdvwpoXBjbCsvYaM/eecT5cD1wH3DCE/6lLgTcCzwPegxXMExGRPiigr85ge+hhaJXut43r+xu2F7d3wL6wa1jv/3Sa99C3M6Af0ZR74BVx3WdPrvNhE+Ar2AnQQuDMPEuvG4HXlhGWZ+kS58NHsb/hdwMfAlLnwxfyLG32tyvSnxOxQOVO4NcVt0W6w03YheddnQ8bVHkRyPmwIfUA/qXYLCyFZcDNwB2sGbzPGokLoPEi+JbUA/ztsYyALYDXxQXnwwzgeuBa4I6B6lPkWbrK+fBFrBjtEc6H3+dZes9w2ysiMlYpoK/OpnH9RAv7Fvts2u9ezW0X1w82bC9ub8+aAX1K8x76p7Ar9hsF7yalWd5sn5EyYvPQx0B9R6xSet7PPt/HTkweBT6UZ+l9w3ldGV3xZPRc58P9WFbFfwFTnQ9n5Fna36wRIqs5HzYAjog3vzgK44FlDIpp9/+kHkgPeirU4Yg1H14AvB3Yt+Huh6kHzzfnWbp4tNqRZ+kSbAac1bVMYtu2AfYGXg44YNe4nAA85Xz4BXBJnqUL+nnuu5wPF2NFUd9DPYtGREQaKKCvzrpx/UwL+xb7rNvvXs1tF9eNAf0Dcb099ndQA+bGbc8K1tMsXxm8ewq7qLAxFviOlpGch37vuL4xTouzBufDRoDHgvk7gffnWTp3mK8pbZJn6Z/jGM5vAfsBn3U+/LcqlEuL3oJ9xlydZ+kdVTdG1hS8S7AsssVpls+suj0NrsGC+VfQpoA+ptTvD7wD2C1uXoZlDFwLXJ9n6cN9PLwt4sXWB+NyQezF3xP7Ln4F9l17CnCi8+ESbFhbXx0bPwAOA17pfNi66t9NRKRTadq66hTBeZ9XqEuKfQYV0MeToSLlvjGgfwTrcd8cO6EtB/R9pS0XFflHO+1+JOeh7zPdPvbOeaw34W7gfQrmu08MxN6HDZU4EPhUPPEV6ZPzYS3grfHmz6psi/TpzcCFwKXBO1d1YxoU3ykvdz6MaueI82Gy8+FtwG+As7Fg/hkss+yQPEs/kGfpRZ0Y8OZZuiTP0uvyLP0ycCgWzP8Tq91zDHCp8+Es58OOTR47G/gjdnH/yPa1WkSku+iktwJxqriiEF0r6XBFhdfB9tCvF19nEfVgHYCYMv8kMJ56QD8v3t1XOn1xFX20A/oR6aF3PkzCxvOB9V6U75sGfBfLYAjASXmWzkO6Up6lM4D3Y3/rrwc+rqBeBvB6rGbIf4B/VdwWaRC8mwp8oLTpv+NF6o6QZ+lMrI7HFOD5o/EazodJzocTsOlWP0J9WNjZwMF5ln6vm4o45llai8F9hg0X+BP2PX8Q1pv/befDTg0P+3lc/5fzYSoiIvIsOuGtRhHML0qzvJUxm0VAP2WQr7O68F6a5c0K4DyOfZkWKfdL4va+eujnxvV6g2zHoMRq/MX7svYwnmp3YDJwbzmlz/kwHvgiVi/gfuC96pnvfnmW3ooFAEuwueqPrbZF0qniON+j482fa4aEjvQy7PP7dqy+yzbYZ3YnKS4Uv3Skn9j58ALgAuC9WKG7O4DTgTflWXpxHL/etfIsvTPP0o9ivfa/wDo3Xgqc73w4Jabqk2dpAG7AzgXeVFFzRUQ6mgL6ahSBeSvp9jACAX0f95cD+lXYHPTQdw/9UDMFhqIobDacgP55cd3Y+/Ze4MXYSeL7uqmHQ/qXZ+m/gDPizfc6H0b8RFvGhJdj9UOeBP5ScVukuWK60T9TT2/fp499q3JLXO8+Uk/ofFjH+XAGcB6WQfYgcBJwXJ6lfx1rhRvzLJ2ZZ+lXgIOxwD7BagT80vlQDLMoeumPcD4MdxieiMiYo4C+GoMN6Iv9RiugH4/10BdBdF899ENtx1AUbZk8jOcoAvrbig3Oh/2B47ALGGfkWfrkMJ5fOlCepdcA52Kfb593PgxlukcZ24re+Qs11WHnicPSioKm/xcX6LyA/va43n0khvg4H3bG6jm8BbvAfh5wZJ6lN471LJI8S+fFwP544F5sGjzvfDgJK/p3H1aQ98DqWiki0pkU0Fej6OFe2O9edUMqikc9oH+sj/sfw/4GxmEBfXHlvxMC+qINQwroY0ptEdDfHrdtg01xBvCNPEtvHlYLpZOdB1wHTAPOjvUURHA+PAdL7V2CFRmTzrM5lmb+VJrlD1GfcnSnGOx3hDxLZ2Hfo+tgGR9D4nxInA9HAT/GhhbcCxydZ6lvNjvLWJZn6e1YsbzvY+cl78Qu0F4edzkmfr+LiEjUMV+MPWawAf0S7Itt8iBPZjaK61l93P8U1kM/DusNKJ57jROI4N0ewbvjgbXipm5Iud8M2BAr9Pdw7D05Czvx+iv1FD4Zg2Ja6seAmdj8xydW2yLpIAfF9R9UCLNjFbOzPACQZvk8bIjUWsAmFbWpL0UG2PP63asP8bvpo8Cp2PC3i4Fj8yy9d2Sa133yLF2RZ+n3gAwbFrMHVkRvGbATw7h4IiIyFimgr0bRw91SQB8L5y1qeGwrisB7fh/3L8DS7cG+KCeWfgYgeLcJ8E1sWrCiyFg7Avph9dBTSrePqYpvxMY5zgLOGuvpi2IpnMDH4813OB+27W9/6RnFVJZXVtoK6c92cV2ebvWBhvs6RZF2P+iAPgbzH8OKvS0DPpJn6dl5li4dwfZ1rZhFdyTwDyzbanPsos4r+nuciEivUUBfjcHMQV8YStr9QGP1F2IBfYKdTBRpyeWU+/8qPc/WWIDdzpT7ofbQF0WKbnM+rI9NaQbw1TxLW82MkC6XZ+m/gUs3ZNba41h5ulI1e1usp7ADdoFUU9V1rq3j+qHStiK436bNbRnIkHroYzD/cew7dinwgTxLdZGpQZ6lzwAfwmYUWI5lbxxcaaNERDrMhKob0KMG1UPfsO9ge+jHAQcH754P/CLN8kWl+8s99Etp0kNPvQjRI9iUQdNoTw990Ya1g3fTgBOwNPwfpVneV8ZB2c5xPQM4GRuPeRNWMVl6yLm8a/4C1t1hFhttex4nXg7pH6puk1Sm6Nm7QcXwOlpRyHJmaVtRC2YzOstdWP2Z7Z0Pk1oZ8x6D+U9igekSLJjP+39U78qzdJnz4TTsfOV44PXOhz1VB0ek/YJ3E7HCnZsCP0mzfE7FTRLUQ1+VwVa5L+87mIB+CjAdeDM2VduhDfcvxP4GEuykYo0e+uDdVGA3LLj+KnbSMmWQbRiq1QE9cDhwFDaVTavz0BY9POOw33sFcLZS7XtL8G7C+jzzhmk88+SGPD0pJXxcBfJ6WlE5/Zp+95KqNZuh5fGG+zpCDOBnYt+jW7X4sAwL5hcDpyiYH1h8n08FAhbYf8/5sFH/jxKRUbA/8P+w4pXvqrgtEimgr0a7euinYD3TxfzyBzTcv4j630A5oC+C6e3i+j6syvBKbPxaO1Pu1wJeVdre+Ds8i/NhnU14YovncetaCateEzf/Ks/SB0a4jUMSvNs8eLdX8C6JtycG794avDs8eKdgc2Q5YOq6LJg7gRVLd2XGNsCrq26UtJ/zYS3gxfHmtVW2RQZU9NB3fEAfFUMDBhwO4Hx4JVa5fRXw/9TL3LoY1J+Dnbtsic1gokxTkfYqn4cf0Ekzj/QyHYRqtGsM/XSsh7uYju75wbsNizvTLC+nBpZT7otgeru4fiDN8oVYtdkEWL8IRkdR0bYNsFT/FdgFhd1job4+HcpvXvwRvvScE/jB1q/nj0djv/8Fo9raFgXvNgV+C3yD+pXN1wGnA6cBh1TTsjFrb4AE/rguC55+DveuO46VR2ssfU96EXaB8M48S5+uujHSXPCuuGi8HHimdNdTcd2JvbItBfTOh62w2VYAvpVn6Y2j2qqx6e/YEMCJwPOBU6ptjjQK3iUxLVvGmNjptHdp00YMcYYPGVkK6Ksx2Gnryvu2FNDHf7qp2HR3f8fGjyfALg27Fr33K3h2QL9lXD8c1w9hwfFaDL36fH9tfk7w7kXxYkHRhuICxDXY7wD18fFNPZd7Dk6oUSMZtzu3rw9cmWfpzP4e00YHUK9bUPQUlzMQDmxvc8a8NK7/PIWFD6zNklXTmbMr1nMvvaUYP690+862+oJ3muXlIVJDuajdLkVAv3VfO8SLiB/FLlZcCZzfhnaNRfdjAf2T2DnLUc6HPaptkhTiUM3zgUuDdwr0xp4dsPP/B4BL4rbdKmuNrKaAvhrtSLmfgqXQr8L+8e6P27dr2K/ovV9JPdAs5oAv0h6LYkSPx/0nDqIdLQne7Qj8EvDAa6lfaFgvru+nxWmLNuapvQFWMGHSpjyx9oH8uZMK4e1f+nn74N1OwEtL214UvGvHkIZesXo+63HU7p/Cwjmb8fhawNFVNkoq8cK4vr7SVshA+vp+XP0d2IYMscEqLnr3NzXmgcBLgHnA51XPZWji+3Y9Vn/glrj5DOfD+D4fJO20H9bpsjGt1zyS7rFdXN9H/ZxcUwJ3AAX01Rj1Hvq430QsMH6Y+gnHlg371UrrxpOkxsJE5YB+pHtJyr3U5THO68T1I/T9O6wWvBu/Dou2WcW48Q+w3dKJLF/8br4/wk1tXfBu2+DdV4N3R8eT0CJDoiiCtA/2fj4A3I39T24/iOd/fvDuoOCdTmYaxHFdxfCMx4HH1mP+nOnMGQfs7XyYGrwbp/FfY1+sKl6kQ4cq2yIDajokLQ4RW4ZdeO60WiP99tA7HyZjBd3AUu3ntqNRY1gorR8HdsSK/0r1yuOr9wveqcbB2LJFXD9C/XOvz3NyaR+dyFajHVXup2CB4irgCSw9DZ49/nBVaV0E9EWQv35cz47rp6n35I90L3J5TM7LSz8X89A/SX0M5cb9PM/mNZK1nmEaj7Ll4sksns8APfqj7FhgX2we3V2wVKW5wL/j/UVK2oMMcp7l4N3mwA+wMZlvHJnmjinTsc+4uTEYeGo8q1ZuyaMzgXG7cfvrgF8BF5VrS8iYtDEWBM7Os3QwF1Kl/fq74D2U4rDt8Bj2vblxH0XaDsH+Bu/EaqjI8BSBxBbYDDwAx6pAXkfYPa5r2LDPVmd+kO5QnH+X44r+zsmlTRTQV6P40hnMPMjFvq1+YU2J+64CZlEPhvsqKFfuoS8C+saekoXx+cYzgj30sYe06JVegJ14F4F8sX6K1ooibbqKceOfZsNxc5i+bDKLF2BzZbZd7DXfr7Tp9XE9k/r8ysWX3ePUMyFabW85fX/A6v89qPhbf6q83oH7ngB4HrcdiV082Y41M0Rk7ClSAh/qdy/pBP0NSRtsplpb5Fm6Cpgfb04t3xezQ46KN38U95XhKRchvArLcNsUfQ9WKo6fXx8bDnFD3Kx07LGlfF41K/6sgL4DKKCvRvG+D+aLvdi31WO2LhbQr6T/f7ziecvp9kVA35hJsCDeN456Ab2RsCFWaG8uMCNuWyuuyz30rVwN3GQl4yc9wzQWsc5TE1m+dID9R9MWrHlyVxTueYr671J8OLaagVC2Z+nnFyh1/FmK93GNgH577p8PsClPvLg0iLWcISJjT5EKrYC+8/WXwdbJhfHmxvXUhu2vwP7+ZmLBpwzf41itn02wDoBiFhvVRqlW0RnxGPXaS5V0qMioKc6rZmGfeSuAqXF2EqmQAoBqtCOgn4b1pK/E/umKgH6jhsCveN7xlHro4z7F+PVFcb0QC+jL+46E4gOiHNQWHw4TsHGT87GUf4AN+xkzvvFKxk+cz3q1FUy4Llnz+dutCCJWNdwuX2DZqLRtsAF9OTV/bTpzOqcqld/b1etpPLMW8MgGzF5vKWsVF4xarlsgXan4X3m4372kE/Q3rWunptxDfYq9aQ3bi5lLfqXe+ZGRZ+lK4NF4c2vgMuw8ZTfnw+Z9PlBGWznY6+RpJmXoiuP5ZJrlRQZwebtURAF9NdoR0Be9BIvSLK/FMcTz4+Ob9W4kpeeuYWO9k/j44rWLk6lx1Cvij4RmAf0k6hcN5sTfYQVWIXgcz+4FAWAZEzetkYybx9TaYiYXKV9VBfTFB9ztDbfLwfu00rbBjkcqevdnNtwWU7y3z5TXCUybxNI71mXBhOVMLC4cbaoMhzGtCOgf7Hcv6QSt9NB3YkA/L65XB/Qx3b7I/rm67S0a21bXnMmzdCn1FO9X9LG/jL7yuVxjFqKMDdPjek5cz27YLhXRCWw1yoFzq4p9Wz1mxVj7laVt/Y3DbxxD32wcY3EyNY6R7aEvj8kpAt21S/cvb/Jz01oCS1h70xrJuIVMWTmLje6Mmxt7TNql+HIrhhEU7ZhDPT2z+D2foR54Nr1YURa8m4xdmFlOveKvxjGtqXEaxmI9fmsengOwgHWXY+/7eOpFIGXsUQ999yi+exY1ua/beuh3j7cfRReTRlrxv1z8b18T1wroq1N0WsxC46vHqsYaYCsatktFFNBXox099MUY9/JrNPvHaxaY12heabhIuR83iHa0ojg5m0+9sFA5A2BFk5+bZggsY9KUGkmykvErayT3x80TKpq3uKic/gj23k3A3rcVaZbXqP8uSfx5MB+MxXOXU9v0xbmm4n1sDOgn7MjdCwDmsv4qBi4YKd2vPNWOdLbGC3Fl/X7+V6xZQL9rXN+keedHXOM0tjfF9S5N9pX2KDooFmGF8aA+fFK6XMxiTIBVpcxdBfQdQgF9NYrgcigBfauBabMe+v7+8Rp76It5fpeW9il+HumAvhx4lYNcStsLxe/T9MNjMZOnAdRIlj/NRnMZ/IWQkVR8uS0GllAf1lAOMBOGFtBPLj130ZO1dh/79qo+A/oduG8uwHzWW4t65sk6yJjjfBhPnMIzz9LFA+0vlevvgvdgM9XaqdlwgKL3+IH2NqUnNH7vPYHV29nA+dCJRRN7QbNzOQV6Y8dg4wppo078UuwF7eihL3owyv94/QXDRWBZaNbGVU3uHwnlHpnGgD6h+YdH0x6apay1LkCNZEnsEamyR6fxyy0p3Qb7vZKGfaC1tjZ7zzqx16pKfabcp4TZAEtYezL6QhrrimylZj2+0nn6+35c2bBPJynaVv4cLgJ6za4w8oqU34mweurA8nR20n4K6Me2xk6S8s86/6xYJ34p9oJOTLkv99Cvovk4/1ppv9HuoS9r9uHR9EuiRjIprosv+3579EdZ+WpmEbyXL1A09tAPpq3l96zK37GT9dlDvz5zV9qGCeWMCb1/Y1PxWbis0lZIq/rLYOvkHvpm39HFmOLH29yWXlD8P5en0H0irlVxuxrljiQFemOPeug7mA5ANe7A0nsH02M0F7iV1q/0P4kV4rmvtO0/WCXe8ontw9g8oQ9jaYHT4z6L4+vdX9p3ZXy+jagXdRsJT8TXepz677kF9h4tBO5u+B3ms2ahvNUWsc7DC1h3wTImFePrbqM+fKDdHsJ+l9lYYbyNsMJIRXXQ27Fx23dgwxmWY+2d/6xnerYlwL+x4zMzvs6T/T6i9zyGvS/FGPkV2Pu7ZBLL5s5io5kLmfIkVlRwMs2rasvYcAv2PyOdr/g8m9Xkvkeof6Z2msepf48V7sb+7vTZMvLmYp/nD5S2Baz+j97vahT/u09h5263oUKkY0kN+y4tD127H6sb0sp5q4yipFZTnRYRERERERGRbtOJaWsiIiIiIiIiMgAF9CIiIiIiIiJdSAG9iIiIiIiISBdSQC8iIiIiIiLShRTQi4iIiIiIiHQhBfQiIiIiIiIiXej/A9P7P8Z+tmDVAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "arcs1 = \"../../../tests/test_data/bedpe_chr9_4000000_6000000.bedpe\"\n", "\n", "frame = XAxis() + \\\n", " Arcs(arcs1, line_width=2, color=\"#3297dc\") + TrackHeight(5) + \\\n", " ArcsCoverage(arcs2, line_width=2, color=\"#dc9732\")\n", "\n", "frame.plot(\"chr9:4500000-6000000\")" ] } ], "metadata": { "@webio": { "lastCommId": null, "lastKernelId": null }, "celltoolbar": "编辑元数据", "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.7.10" } }, "nbformat": 4, "nbformat_minor": 4 }