My_sample_notebook.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Uncomment the following lines (by deleting the leading `#`) if you are running this in Colab. (Thanks to Jake VanderPlas for the tip!)**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# !pip install biopython\n",
    "# !pip install folium\n",
    "# !curl -O https://raw.githubusercontent.com/jperkel/example_notebook/master/NC_005816.gb"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This simple notebook demonstrates how users can interleave text, code, and results in a single document. We start with a simple calculation -- computing the first 25 numbers in the Fibonacci sequence, where each value equals the sum of the two previous values. The Jupyter notebook allows us to express that mathematically, using the typesetting language $\\LaTeX{}$: $$F_n = F_{n-1} + F_{n-2}$$\n",
    "Thus, the sequence is: 0, 1, 1, 2, 3, 5, 8, ..."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The first cell contains an IPython 'magic' code, '%matplotlib', which allows the notebook to display plots inline, in the body of the notebook."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233, 377, 610, 987, 1597, 2584, 4181, 6765, 10946, 17711, 28657, 46368]\n"
     ]
    }
   ],
   "source": [
    "# calculate the first 25 Fibonacci numbers\n",
    "f1 = 0\n",
    "f2 = 1\n",
    "ar = [f1, f2] # a list to hold the computed values. We know the first two numbers\n",
    "\n",
    "# we only need to run our calculation 23 times, because positions 1 and 2 are known\n",
    "for i in range (23):\n",
    "    f3 = f1 + f2\n",
    "    ar.append (f3)\n",
    "    f1 = f2\n",
    "    f2 = f3\n",
    "    \n",
    "print (ar) # below, you see the output of the code itself."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Plot the data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEWCAYAAACufwpNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAfRklEQVR4nO3deZhcVZ3/8fcnCVuLGDABIVsjREdwBKGBII7DoD8IiwYVWZ7+SUA0s4jLzDiCMj7EJf4UN8ZB/U0UNEDLIohEBidGER0VgQ6yBUQayAYMCSSBYFgkfOePewpuiu7qqupbVanqz+t56qm655577/dWJfXtc86tcxURmJmZFWFMqwMwM7PO4aRiZmaFcVIxM7PCOKmYmVlhnFTMzKwwTipmZlYYJxWriaS5ki5uwnGWSjq00cdpJEk/kTQ7vT5F0q9bHVM1JPVK+mmr4yiRtEzS21odh1XHScU2I+nJ3ON5SU/llnsLPtb3JD1bdswTACJi74i4vo59dksKSeMq1JktaYmkJyStknROvr6k6yU9nYvpngr7mivpz2Xn8PF0DkdGxIJaz6HVIqIvIg4fbF3u/f3PsvKLJc1tSoC2RXNSsc1ExPalB7ACeHuurK8Bhzwnf8yIuGy4DSoljCp1AR8FJgAHAW8FPlZW5/RcTK8dZn+XlZ3DOSOMrx3MkHRIq4OoRQH/bqwKTipWj60lXShpQ+qm6imtkLSbpCslrZH0gKQP13OAfJdHag1ckf4afgI4RdKBkvpTa+MRSV9Nm/4qPa9PrYaDy/cdEd+KiP+OiGcj4kGgDyj8CzK1eN6/eZH+XdLjkv4g6a25FbtJWihpraQBSR/IrZsr6fIK7/mZku5L6+6S9M6yOD4g6e7c+v1S+RRJP0yf1WOSzkvl1XTVnQN8bojzfsn2qXWzZ3r9PUnfTN2DT0r6jaRXSTpX0rr03ryxbLcHpNjXSfqupG1z+z5G0q2S1kv6raQ35NYtk3SGpNuBP0kal5YfTO/HPfnPwUbOScXq8Q7gUmA8sBAofRmNAX4M3AZMImsBfFTSEQUccxZwRTpmH/BvwL9FxA7AHsDlqd5b0vP41Gq4oYp9vwVYWlb2/yQ9mr7wDh1x9JmDgPvJWkhnAz+UtFNadwmwCtgNOA74fNmX3aDveXIf8FfAK4BPAxdL2hVA0nuAucDJwA5pP49JGgtcAywHusk+r0trOJdvAK9R/WMdxwP/SvZePAPcANySlq8AvlpWvxc4guyzfk3alpQgLwD+Fngl8B/AQknb5LY9CTia7L3bAzgdOCAiXp72uazOc7BBOKlYPX4dEddGxCbgImCfVH4AMDEiPpNaAfcD3wZOrLCvj6W/MNdLerRCvRsi4kcR8XxEPAX8GdhT0oSIeDIiflfPiUg6FegBvpwrPgN4NdkX7Xzgx5L2qLCb43PnsF7SbkPUWw2cGxF/Tt189wBHS5oCvBk4IyKejohbge8A781tO9R7TkT8ICIeSu/NZcC9wIFp9fvJuhhvjsxARCxP63cD/iUi/pSOW8uFBE8D8xiitVKFqyJiSUQ8DVwFPB0RF6bzuwwob6mcFxErI2JtOu5JqfwDwH9ExI0RsSmNYT0DzMht+/W07VPAJmAbYC9JW0XEsoi4r85zsEE4qVg9/if3eiOwbeqvngbslv+CBT4J7FJhX1+OiPHpMaFCvZVly6eR/cX6B0k3Szqm1pOQdCzwBeDIiHghoaUvqA0R8Uz6kvoNcFSFXV2eO4fxEfHQEPUejM1ncF1O9sW+G7A2IjaUrZuUWx7qPUfSybnun/XA68n+4geYQtaSKTcFWB4Rz1U4r+F8G9hF0tvr2PaR3OunBlnevqx+/vMvvW+Q/Zv757J/c1Ny6zfbNiIGyMbT5gKrJV1a4Y8Aq4OTihVpJfBA2RfsyyOi0hdytTabTjsi7o2Ik4CdgS8CV0h6WXm9oUiaSfal+PaIuKOKY6v2kF9ikqT8fqYCD6XHTpJeXrbuweF2KGka2XmcDrwyIsYDd+biXUnW5VNuJTBVIxi8jog/k3W3fZbN358/kV0MUYrxVfUeI2dK7nXpfYPsPOaV/ZvriohL8qGWxf39iHgzWUIKsn8/VhAnFSvSTcATaSB0O0ljJb1e0gFFH0jS/5U0MSKeB9an4k3AGuB5su6robY9jGxc5t0RcVPZuvGSjpC0bRrU7SUbc1lUQNg7Ax+WtFUa63gdcG1ErAR+SzaOs20aaD4txTicUiJdk+I/laylUvIdsi7G/ZXZMyWim4CHgS9Ielk6bj0XK1xE1p00M1d2G7C3pH3TgPrcOvZb7oOSJqcxqE+SdZFBllD/TtJB6fxeJunosgT9AkmvlXRYGnN5mqxVtKmA+CxxUrHCpP7wtwP7Ag8Aj5J9qb2iAYebCSyV9CTZoP2JaVxgI1mf+29Sd8iMQbb9VIrpWr3425KfpHVbkY0TrEnxfwg4NiKG/K1KDW4Epqf9zgOOi4jH0rqTyAbMHyIbYzg7IhYPt8OIuAv4CtlA9yPAX5J115XW/yAd6/vABuBHwE65z2pPskvHVwEn1HpCaT9nAzvlyv4IfAb4Gdn4ThE/+vw+8FOyCx3uJ43lREQ/2bjKecA6YAA4pcJ+tiHr8nyUrEtxZ7IkZQWRb9JlZmZFcUvFzMwK46RiZmaFcVIxM7PCOKmYmVlhRt0EaxMmTIju7u5Wh2Fm1jaWLFnyaERMrKbuqEsq3d3d9Pf3tzoMM7O2IWl5tXXd/WVmZoVxUjEzs8I4qZiZWWGcVMzMrDBOKmZmVhgnFTOzTtbXB93dMGZM9txXzeTX9Rt1lxSbmY0afX0wZw5s3JgtL1+eLQP09jbkkG6pmJl1qrPOejGhlGzcmJU3iJOKmVmnWrGitvICOKmYmXWqqVNrKy+Ak4qZWaeaNw+6ujYv6+rKyhvEScXMrFP19sL8+TBtGkjZ8/z5DRukB1/9ZWbW2Xp7G5pEyrmlYmZmhXFSMTOzwjipmJlZYZxUzMysME4qZmZWGCcVMzMrjJOKmZkVxknFzMwK46RiZmaFcVIxM7PCOKmYmVlhnFTMzKwwTipmZlYYJxUzMyuMk4qZmRXGScXMzArjpGJmZoVxUjEzs8I4qZiZWWEanlQkjZX0e0nXpOXdJd0o6V5Jl0naOpVvk5YH0vru3D4+kcrvkXRErnxmKhuQdGajz8XMzCprRkvlI8DdueUvAl+LiOnAOuC0VH4asC4i9gS+luohaS/gRGBvYCbwzZSoxgLfAI4E9gJOSnXNzKxFGppUJE0Gjga+k5YFHAZckaosAI5Nr2elZdL6t6b6s4BLI+KZiHgAGAAOTI+BiLg/Ip4FLk11zcysRRrdUjkX+DjwfFp+JbA+Ip5Ly6uASen1JGAlQFr/eKr/QnnZNkOVv4SkOZL6JfWvWbNmpOdkZmZDaFhSkXQMsDoiluSLB6kaw6yrtfylhRHzI6InInomTpxYIWozMxuJcQ3c9yHAOyQdBWwL7EDWchkvaVxqjUwGHkr1VwFTgFWSxgGvANbmykvy2wxVbmZmLdCwlkpEfCIiJkdEN9lA+3UR0Qv8AjguVZsNXJ1eL0zLpPXXRUSk8hPT1WG7A9OBm4CbgenparKt0zEWNup8zMxseI1sqQzlDOBSSZ8Dfg+cn8rPBy6SNEDWQjkRICKWSrocuAt4DvhgRGwCkHQ6sAgYC1wQEUubeiZmZrYZZY2B0aOnpyf6+/tbHYaZWduQtCQieqqp61/Um5lZYZxUzMysME4qZmZWGCcVMzMrjJOKmZkVxknFzMwK46RiZmaFcVIxM7PCOKmYmVlhnFTMzKwwTipmZlYYJxUzMyuMk4qZmRXGScXMzArjpGJmZoVxUjEzs8I4qZiZWWGcVMzMrDBOKmZmVhgnFTMzK4yTipmZFcZJxczMCuOkYmZmhXFSMTOzwjipmJlZYZxUzMysME4qZmZWGCcVMzMrjJOKmZkVxknFzMwK46RiZmaFcVIxM7PCOKmYmVlhnFTMzKwwDUsqkraVdJOk2yQtlfTpVL67pBsl3SvpMklbp/Jt0vJAWt+d29cnUvk9ko7Ilc9MZQOSzmzUuZiZWXXqSiqSbqmi2jPAYRGxD7AvMFPSDOCLwNciYjqwDjgt1T8NWBcRewJfS/WQtBdwIrA3MBP4pqSxksYC3wCOBPYCTkp1zcysRepKKhGxXxV1IiKeTItbpUcAhwFXpPIFwLHp9ay0TFr/VklK5ZdGxDMR8QAwAByYHgMRcX9EPAtcmuqamVmLVJVUJO0i6Zj02LnanacWxa3AamAxcB+wPiKeS1VWAZPS60nASoC0/nHglfnysm2GKh8sjjmS+iX1r1mzptrwzcysRsMmFUnHAzcB7wGOB26UdFw1O4+ITRGxLzCZrGXxusGqlQ41xLpayweLY35E9EREz8SJE4cP3MzM6jKuijpnAQdExGoASROBn/FiF9awImK9pOuBGcB4SeNSa2Qy8FCqtgqYAqySNA54BbA2V16S32aocjMza4Fqur/GlBJK8lg120maKGl8er0d8DbgbuAXQKmlMxu4Or1emJZJ66+LiEjlJ6arw3YHppO1nG4GpqerybYmG8xfWMX5mJm1p74+6O6GMWOy576+Vkf0EtW0VP5L0iLgkrR8AnBtFdvtCixIV2mNAS6PiGsk3QVcKulzwO+B81P984GLJA2QtVBOBIiIpZIuB+4CngM+GBGbACSdDiwCxgIXRMTSKuIyM2s/fX0wZw5s3JgtL1+eLQP09rYurjLKGgPDVJLeDRxCNo7xq4i4qtGBNUpPT0/09/e3Ogwzs9p0d2eJpNy0abBsWUMPLWlJRPRUU7ealgoRcSVw5YiiMjOz+q1YUVt5iwyZVCQ9wBBXU5H9DGWPxoRkZmYvMXXq4C2VqVObH0sFlQbce4ADco+DgK+QdYHd2vjQzMzsBfPmQVfX5mVdXVn5FmTIpBIRj0XEY2RTqRxDdtXWwcDREfHuJsVnZmaQDcbPn5+NoUjZ8/z5W9QgPVTu/toKeB/wj8CvgVkRcV+zAjMzszK9vVtcEilXaaD+AbJLeM8FVgD7SNqntDIiftjg2MzMrM1USio/Ixuo3yc98gJwUjEzs80MmVQi4pQmxmFmZh3Ad340M7PCOKmYmVlhnFTMzKwwlS4pPiwirpP0rsHW++ovMzMrV+nqr78GrgPePsg6X/1lZmYvUenqr7PT86nNC8fMzNpZNTfb+nzpZltpecd0LxQzM7PNVDNQf2RErC8tRMQ64KjGhWRmZu2qmqQyVtI2pYV0a+BtKtQ3M7NRqpqbdF0M/FzSd8kG6N8HLGhoVGZm1paGTSoRcY6k24G3kd1L5bMRsajhkZmZWdsZNqlI2h24PiL+Ky1vJ6k7IpY1OjgzM2sv1Yyp/AB4Pre8KZWZmZltppqkMi4ini0tpNdbNy4kMzNrV9UklTWS3lFakDQLeLRxIZmZWbuq5uqvvwP6JJ1HNlC/Eji5oVGZmVlbqubqr/uAGZK2BxQRGxoflpmZtaNqWipIOhrYG9hWEgAR8ZkGxmVmZm2omrm//j9wAvAhsu6v9wDTGhyXmZm1oWoG6t8UEScD6yLi08DBwJTGhmVmZu2omqTyVHreKGk34M/A7o0LyczM2lU1YyrXpKnvvwTcQjb/13caGpWZmbWlaq7++mx6eaWka4BtI+LxxoZlZmbtqNqrv94EdJfqSyIiLmxgXGZm1oaqmVDyImAP4Fayeb8g6wJzUjEzs81UM1DfAxwSEf8QER9Kjw8Pt5GkKZJ+IeluSUslfSSV7yRpsaR70/OOqVySvi5pQNLtkvbL7Wt2qn+vpNm58v0l3ZG2+bpKP6IxM7OWqCap3Am8qo59Pwf8c0S8DpgBfFDSXsCZwM8jYjrw87QMcCQwPT3mAN+CLAkBZwMHAQcCZ5cSUaozJ7fdzDriNDOzglQzpjIBuEvSTcAzpcKIeMfQm0BEPAw8nF5vkHQ3MAmYBRyaqi0ArgfOSOUXRkQAv5M0XtKuqe7iiFgLIGkxMFPS9cAOEXFDKr8QOBb4SRXnZGZmDVBNUpk70oNI6gbeCNwI7JISDhHxsKSdU7VJZJNVlqxKZZXKVw1SbmZmLTJs91dE/BL4A/Dy9Lg7lVUlTUR5JfDRiHiiUtXBDl9H+WAxzJHUL6l/zZo1w4VsZmZ1qmbur+OBm8jm/DoeuFHScdXsXNJWZAmlLyJ+mIofSd1apOfVqXwVm0//Mhl4aJjyyYOUv0REzI+InojomThxYjWhm5lZHaoZqD8LOCAiZqc5wA4EPjXcRulKrPPJWjZfza1aCJSu4JoNXJ0rPzldBTYDeDx1ky0CDpe0YxqgPxxYlNZtkDQjHevk3L7MzKwFqhlTGRMRq3PLj1FdMjoEeC9wh6RbU9kngS8Al0s6DVhB1gICuBY4ChgANgKnAkTEWkmfBW5O9T5TGrQH/h74HrAd2QC9B+nNzFpI2cVWFSpIXwLeAFySik4Abo+IMxocW0P09PREf39/q8MwM2sbkpZERE81dauZ++tfJL0LeDPZ4Pj8iLhqhDGamVkHqqYbC+C3wC+B64AbGheOmdko0tcH3d0wZkz23NfX6ohGrJqrv95PdvXXO4HjyH6Y+L5GB2Zm1tH6+mDOHFi+HCKy5zlz2j6xVDOmcg/Z3R8fS8uvBH4bEa9tQnyF85iKmW0RuruzRFJu2jRYtqzZ0VRUy5hKNd1fq4ANueUNbP4LdzMzq9WKFbWVt4khB+ol/VN6+SDZDx6vJvvF+iyy7jAzM6vX1KmDt1SmTm1+LAWq1FIpTctyH/AjXpwC5WrSRJFmZlanefOgq2vzsq6urLyNDdlSiYhPNzMQM7NRpbc3ez7rrKzLa+rULKGUyttUpe6vcyPio5J+zCATNQ439b2ZmQ2jt7ftk0i5Sj9+vCg9f7kZgZiZWfurlFTWwAtT35uZmQ2r0kD9j0ovJF3ZhFjMzKzNVUoq+ZtgvbrRgZiZWfurlFRiiNdmZmaDqjSmso+kJ8haLNul16TliIgdGh6dmZm1lUq/UxnbzEDMzKz9VTv1vZmZ2bCcVMzMrDBOKmZmVhgnFTMzK4yTipmZFcZJxczMCuOkYmZmhXFSMTOzwjipmJlZYZxUzMysME4qZmZWGCcVMzMrjJOKmZkVxknFzMwK46RiZlaEvj7o7oYxY7Lnvr5WR9QSlW7SZWZm1ejrgzlzYOPGbHn58mwZoLe3dXG1gFsqZmYjddZZLyaUko0bs/JRxknFzGykVqyorbyDNSypSLpA0mpJd+bKdpK0WNK96XnHVC5JX5c0IOl2Sfvltpmd6t8raXaufH9Jd6Rtvi5JjToXM7OKpk6trbyDNbKl8j1gZlnZmcDPI2I68PO0DHAkMD095gDfgiwJAWcDBwEHAmeXElGqMye3XfmxzMyaY9486OravKyrKysfZRqWVCLiV8DasuJZwIL0egFwbK78wsj8DhgvaVfgCGBxRKyNiHXAYmBmWrdDRNwQEQFcmNuXmVlz9fbC/PkwbRpI2fP8+aNukB6af/XXLhHxMEBEPCxp51Q+CViZq7cqlVUqXzVI+aAkzSFr1TB1FDZHzawJentHZRIpt6UM1A82HhJ1lA8qIuZHRE9E9EycOLHOEM3MbDjNTiqPpK4r0vPqVL4KmJKrNxl4aJjyyYOUm5lZCzU7qSwESldwzQauzpWfnK4CmwE8nrrJFgGHS9oxDdAfDixK6zZImpGu+jo5ty8zM2uRho2pSLoEOBSYIGkV2VVcXwAul3QasAJ4T6p+LXAUMABsBE4FiIi1kj4L3JzqfSYiSoP/f092hdl2wE/Sw8zMWkjZxVOjR09PT/T397c6DDOztiFpSUT0VFN3SxmoNzOzDuCkYmZmhXFSMTOzwjipmJlZYZxUzMysME4qZmblfBfHuvnOj2Zmeb6L44i4pWJmlue7OI6Ik4qZWZ7v4jgiTipmZnm+i+OIOKmYmeX5Lo4j4qRiZpbnuziOiK/+MjMr57s41s0tFTMzK4yTipmZFcZJxczMCuOkYmadz9OuNI0H6s2ss3nalaZyS8XMOpunXWkqJxUz62yedqWpnFTMrLN52pWmclIxs87maVeayknFzDqbp11pKicVM2sv9Vwe3NsLy5bB889nz04oDeNLis2sffjy4C2eWypm1j58efAWz0nFzNqHLw/e4jmpmFnr1Do+4suDt3hOKmbWGqXxkeXLIeLF8ZFKicWXB2/xnFTMrDXqGR/x5cFbPCcVMytOLd1Z9Y6P+PLgLZqTipkNrtbxjlq7szw+0pGcVMxGg0YnCKi9O8vjIx3JScWs1er5hXgt2zQjQUDt3VkeH+lMEdHWD2AmcA8wAJw5XP39998/anbxxRHTpkVI2fPFFxdbf0vexnE1Nq6LL47o6orIvu6zR1dX5e1q3WbatM3rlh7Tpg19DGnwbaSht6nnONYWgP6o9ju52opb4gMYC9wHvBrYGrgN2KvSNjUnlVr/AzfjS6JZ2ziuxsdVzxdxrds0K0HUc/7WFkZTUjkYWJRb/gTwiUrb1JxUav3P1YwviWZt47gaH1c9X/i1btPMBFFPa9C2eKMpqRwHfCe3/F7gvEHqzQH6gf6pU6fW9m7W+h+4GV8SzdrGcTU+rmYkOycIG6Fakkq7D9RrkLJ4SUHE/IjoiYieiRMn1naEWi97rOcyyS11G8fV+LjquQKq1m3qHRD370GsHtVmny3xQTO6vzqlv7+TzqWT4ipt14yLDszqxCjq/hoH3A/szosD9XtX2sZXf3XIuXRSXGZbuFqSirL67UvSUcC5ZFeCXRARFX851dPTE/39/U2JzcysE0haEhE91dRt+zs/RsS1wLWtjsPMzPyLejMzK5CTipmZFcZJxczMCuOkYmZmhWn7q79qJWkNsLzOzScAjxYYTjsZzecOo/v8fe6jV+n8p0VEVb8cH3VJZSQk9Vd7WV2nGc3nDqP7/H3uo/Pcob7zd/eXmZkVxknFzMwK46RSm/mtDqCFRvO5w+g+f5/76FXz+XtMxczMCuOWipmZFcZJxczMCuOkUgVJMyXdI2lA0pmtjqfZJC2TdIekWyV19BTPki6QtFrSnbmynSQtlnRvet6xlTE20hDnP1fSg+nzvzXNDN5xJE2R9AtJd0taKukjqbzjP/8K517zZ+8xlWFIGgv8Efg/wCrgZuCkiLirpYE1kaRlQE9EdPyPwCS9BXgSuDAiXp/KzgHWRsQX0h8VO0bEGa2Ms1GGOP+5wJMR8eVWxtZoknYFdo2IWyS9HFgCHAucQod//hXO/Xhq/OzdUhnegcBARNwfEc8ClwKzWhyTNUhE/ApYW1Y8C1iQXi8g+8/WkYY4/1EhIh6OiFvS6w3A3cAkRsHnX+Hca+akMrxJwMrc8irqfLPbWAA/lbRE0pxWB9MCu0TEw5D95wN2bnE8rXC6pNtT91jHdf+Uk9QNvBG4kVH2+ZedO9T42TupDE+DlI22PsNDImI/4Ejgg6mLxEaPbwF7APsCDwNfaW04jSVpe+BK4KMR8USr42mmQc695s/eSWV4q4ApueXJwEMtiqUlIuKh9LwauIqsS3A0eST1OZf6nle3OJ6miohHImJTRDwPfJsO/vwlbUX2pdoXET9MxaPi8x/s3Ov57J1UhnczMF3S7pK2Bk4EFrY4pqaR9LI0cIeklwGHA3dW3qrjLARmp9ezgatbGEvTlb5Qk3fSoZ+/JAHnA3dHxFdzqzr+8x/q3Ov57H31VxXSZXTnAmOBCyJiXotDahpJryZrnQCMA77fyecv6RLgULIpvx8BzgZ+BFwOTAVWAO+JiI4czB7i/A8l6/4IYBnwt6Uxhk4i6c3AfwN3AM+n4k+SjS109Odf4dxPosbP3knFzMwK4+4vMzMrjJOKmZkVxknFzMwK46RiZmaFcVIxM7PCOKnYqCfprDQz6+1pJtaDWh3TSEj6XppZdpu0PCFNCmrWcONaHYBZK0k6GDgG2C8inpE0Adi6xWEVYRPwPrJpNsyaxi0VG+12BR6NiGcAIuLR0rQ0kvaX9Ms0keai3FQd+0u6TdINkr5UuveIpFMknVfasaRrJB2aXh+e6t8i6QdpjqXSvWo+ncrvkPQXqXx7Sd9NZbdLenel/QziXOAfJW32h6MyX5J0Z9r3CcW9lWZOKmY/BaZI+qOkb0r6a3hhHqR/B46LiP2BC4DSTALfBT4cEQdXc4DU+vlX4G1pYs5+4J9yVR5N5d8CPpbKPgU8HhF/GRFvAK6rYj95K4BfA+8tK38X2S+k9wHeBnypbCoOsxFx95eNahHxpKT9gb8C/ga4LN2IqR94PbA4mxaJscDDkl4BjI+IX6ZdXEQ2e3MlM4C9gN+kfW0N3JBbX5q4cAnZlz5kX/gn5uJcJ+mYYfZT7vNk81b9Z67szcAlEbGJbKLEXwIHMIrms7PGclKxUS99wV4PXC/pDrJJA5cAS8tbI5LGM/StD55j89b/tqXNgMURcdIQ2z2Tnjfx4v9JDXKc4fazmYgYkHQr2d378vswaxh3f9moJum1kqbnivYFlgP3ABPTQD6StpK0d0SsBx5PE/AB9Oa2XQbsK2mMpCm8OE3474BDJO2Z9tUl6TXDhPZT4PRcnDvWuZ95vNilBvAr4ARJYyVNBN4C3DTMPsyq5qRio932wAJJd0m6nax7aW66dfRxwBcl3QbcCrwpbXMq8A1JNwBP5fb1G+ABsplevwyUbs+6huw+55ekY/wO+Ith4vocsGMaUL8N+Jt69hMRS0txJFcBtwO3AdcBH4+I/wFIrRqzEfEsxWYjoOzWq9dExOtbHIrZFsEtFTMzK4xbKmZmVhi3VMzMrDBOKmZmVhgnFTMzK4yTipmZFcZJxczMCvO/ndIxBwv+BCsAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots()\n",
    "ax.plot (range(25), ar, \"ro\")\n",
    "## uncomment the following call to ax.plot() (by removing the leading '#') and select \n",
    "## 'Cell > Run All' (in Binder) or 'Run > Run All Cells' (Jupyter) to change the graph below\n",
    "# ax.plot (range(25), ar)\n",
    "ax.set (xlabel = \"Sequence No.\", ylabel = \"Fibonacci No.\", \n",
    "       title = \"The First 25 Fibonacci Numbers\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Anything you can do programmatically, can be documented in a notebook. Here we'll do some simple sequence analysis with Biopython. The following example is adapted from the [Biopython tutorial](http://biopython.org/DIST/docs/tutorial/Tutorial.html). "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "First, we'll read in a Genbank-formatted file, which represents a circular DNA, called a plasmid, from the bacterium, *Yersinia pestis*. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Seq('TGTAACGAACGGTGCAATAGTGATCCACACCCAACGCCTGAAATCAGATCCAGG...CTG', IUPACAmbiguousDNA())"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from Bio import SeqIO\n",
    "\n",
    "record = SeqIO.read(\"NC_005816.gb\", \"genbank\")\n",
    "record.seq"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "How long is the DNA, and how which genes does it encode?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Length: 9609 bp\n",
      "No. features: 41 \n",
      "\n",
      "Feat. 3: putative transposase\n",
      "Feat. 8: transposase/IS protein\n",
      "Feat. 15: putative replication regulatory protein\n",
      "Feat. 18: hypothetical protein\n",
      "Feat. 21: pesticin immunity protein\n",
      "Feat. 23: pesticin\n",
      "Feat. 29: hypothetical protein\n",
      "Feat. 32: outer membrane protease\n",
      "Feat. 35: putative transcriptional regulator\n",
      "Feat. 38: hypothetical protein\n"
     ]
    }
   ],
   "source": [
    "print(\"Length:\", len(record.seq), \"bp\")\n",
    "print(\"No. features:\", len(record.features),\"\\n\")\n",
    "\n",
    "for i,feat in enumerate (record.features):\n",
    "    if (feat.type == \"CDS\"):\n",
    "        product = feat.qualifiers['product'][0]\n",
    "        print (\"Feat. %d: %s\" % (i, product))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This plasmid encodes 41 features, including 10 genes. Let's focus on feature #18, which encodes a 'hypothetical protein'."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Gene location\n",
      "Start: 3485 \n",
      "End: 3857\n",
      "GeneID: ['GI:45478715', 'GeneID:2767720']\n"
     ]
    }
   ],
   "source": [
    "feat = record.features[18]\n",
    "print (\"Gene location\\nStart:\", feat.location.start, \"\\nEnd:\", feat.location.end)\n",
    "print (\"GeneID:\", feat.qualifiers.get(\"db_xref\"))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "These data show us that the gene is located between bases 3485 and 3857. Let's retrieve that segment and translate it."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "VSKKRRPQKRPRRRRFFHRLRPPDEHHKNRRSSQRWRNPTGLKDTRRFPPEAPSCALLFRPCRLPDTSPPFSLREAWRFLIAHAVGISVRCRSFAPSWAVCTNPPFSPTTAPYPVTIVLSPTR*\n"
     ]
    }
   ],
   "source": [
    "my_gene = record[3485:3857]\n",
    "print (my_gene.seq.translate(table=\"Bacterial\"))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now, let's find out more about this gene, using the NCBI Entrez database. We learned the accession number above: `GeneID: ['GI:45478715', 'GeneID:2767720']`.\n",
    "\n",
    "**20 May 2020**: Updated the following two cells to reflect changes in NCBI's Entrez database. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "from Bio import Entrez"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "LOCUS       AAS58761                 123 aa            linear   BCT 26-JUL-2016\n",
      "DEFINITION  conserved hypothetical protein (plasmid) [Yersinia pestis biovar\n",
      "            Microtus str. 91001].\n",
      "ACCESSION   AAS58761\n",
      "VERSION     AAS58761.1\n",
      "DBLINK      BioProject: PRJNA10638\n",
      "            BioSample: SAMN02602970\n",
      "DBSOURCE    accession AE017046.1\n",
      "KEYWORDS    .\n",
      "SOURCE      Yersinia pestis biovar Microtus str. 91001\n",
      "  ORGANISM  Yersinia pestis biovar Microtus str. 91001\n",
      "            Bacteria; Proteobacteria; Gammaproteobacteria; Enterobacterales;\n",
      "            Yersiniaceae; Yersinia.\n",
      "REFERENCE   1  (residues 1 to 123)\n",
      "  AUTHORS   Song,Y., Tong,Z., Wang,J., Wang,L., Guo,Z., Han,Y., Zhang,J.,\n",
      "            Pei,D., Zhou,D., Qin,H., Pang,X., Han,Y., Zhai,J., Li,M., Cui,B.,\n",
      "            Qi,Z., Jin,L., Dai,R., Chen,F., Li,S., Ye,C., Du,Z., Lin,W.,\n",
      "            Wang,J., Yu,J., Yang,H., Wang,J., Huang,P. and Yang,R.\n",
      "  TITLE     Complete genome sequence of Yersinia pestis strain 91001, an\n",
      "            isolate avirulent to humans\n",
      "  JOURNAL   DNA Res. 11 (3), 179-197 (2004)\n",
      "   PUBMED   15368893\n",
      "REFERENCE   2  (residues 1 to 123)\n",
      "  AUTHORS   Zhou,D., Tong,Z., Song,Y., Han,Y., Pei,D., Pang,X., Zhai,J., Li,M.,\n",
      "            Cui,B., Qi,Z., Jin,L., Dai,R., Du,Z., Wang,J., Guo,Z., Wang,J.,\n",
      "            Huang,P. and Yang,R.\n",
      "  TITLE     Genetics of metabolic variations between Yersinia pestis biovars\n",
      "            and the proposal of a new biovar, microtus\n",
      "  JOURNAL   J. Bacteriol. 186 (15), 5147-5152 (2004)\n",
      "   PUBMED   15262951\n",
      "REFERENCE   3  (residues 1 to 123)\n",
      "  AUTHORS   Song,Y., Tong,Z., Wang,L., Han,Y., Zhang,J., Pei,D., Wang,J.,\n",
      "            Zhou,D., Han,Y., Pang,X., Zhai,J., Chen,F., Qin,H., Wang,J., Li,S.,\n",
      "            Guo,Z., Ye,C., Du,Z., Lin,W., Wang,J., Yu,J., Yang,H., Wang,J.,\n",
      "            Huang,P. and Yang,R.\n",
      "  TITLE     Direct Submission\n",
      "  JOURNAL   Submitted (24-APR-2003) The Institute of Microbiology and\n",
      "            Epidemiology, Academy of Military Medical Sciences, No. 20,\n",
      "            Dongdajie Street, Fengtai District, Beijing 100071, People's\n",
      "            Republic of China\n",
      "COMMENT     Method: conceptual translation.\n",
      "FEATURES             Location/Qualifiers\n",
      "     source          1..123\n",
      "                     /organism=\"Yersinia pestis biovar Microtus str. 91001\"\n",
      "                     /strain=\"91001\"\n",
      "                     /sub_species=\"microtus\"\n",
      "                     /db_xref=\"taxon:229193\"\n",
      "                     /plasmid=\"pPCP1\"\n",
      "                     /note=\"biovar: Microtus\"\n",
      "     Protein         1..123\n",
      "                     /product=\"conserved hypothetical protein\"\n",
      "     CDS             1..123\n",
      "                     /locus_tag=\"YP_pPCP04\"\n",
      "                     /old_locus_tag=\"pPCP04\"\n",
      "                     /coded_by=\"AE017046.1:3486..3857\"\n",
      "                     /note=\"Best Blastp hit = gi|321919|pir||JQ1541\n",
      "                     hypothetical 16.9K protein - Salmonella typhi murium\n",
      "                     plasmid NTP16.\"\n",
      "                     /transl_table=11\n",
      "ORIGIN      \n",
      "        1 mskkrrpqkr prrrrffhrl rppdehhknr rssqrwrnpt glkdtrrfpp eapscallfr\n",
      "       61 pcrlpdtspp fslreawrfl iahavgisvr crsfapswav ctnppfsptt apypvtivls\n",
      "      121 ptr\n",
      "//\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "Entrez.email = \"A.N.Other@example.com\"  # Always tell NCBI who you are\n",
    "with Entrez.efetch (db=\"protein\",rettype=\"gb\",retmode=\"text\",id=\"AAS58761\") as handle:\n",
    "   print (handle.read())\n",
    "# handle = Entrez.efetch (db=\"nucleotide\",id=\"45478715\",rettype=\"gb\",retmode=\"text\")\n",
    "# print (handle.read())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Jupyter also supports interactive data exploration. Here we'll add an interactive map, just as in our June [mapping feature](https://www.nature.com/articles/d41586-018-05331-6) -- something you might do when working with geospatial data. This requires the Python Leaflet library, folium. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "import folium"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we create a simple map: a few points in London, Oxford and Cambridge, overlaid on either a street map, or on a map of geological data provided by the [Macrostrat Project](https://macrostrat.org/). "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "coords = { \n",
    "    0: { \"name\": \"Nature\", \"lat\": 51.533925, \"long\": -0.121553 },\n",
    "    1: { \"name\": \"Francis Crick Institite\", \"lat\": 51.531877, \"long\": -0.128767 },\n",
    "    2: { \"name\": \"MRC Laboratory for Molecular Cell Biology\", \"lat\": 51.524435, \"long\": -0.132495 },\n",
    "    3: { \"name\": \"Kings College London\", \"lat\": 51.511573, \"long\": -0.116083 },\n",
    "    4: { \"name\": \"Imperial College London\", \"lat\": 51.498780, \"long\": -0.174888 },\n",
    "    5: { \"name\": \"Cambridge University\", \"lat\": 52.206960, \"long\": 0.115034 },\n",
    "    6: { \"name\": \"Oxford University\", \"lat\": 51.754843, \"long\": -1.254302 },\n",
    "    7: { \"name\": \"Platform 9 3/4\", \"lat\": 51.532349, \"long\": -0.123806 }\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<folium.features.LatLngPopup at 0x120acde50>"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "m = folium.Map(location = [51.8561, -0.2966], tiles = 'CartoDB positron', zoom_start = 9)\n",
    "\n",
    "# add the locations to the map\n",
    "for key in coords.keys():\n",
    "    folium.CircleMarker(\n",
    "    location=[coords[key]['lat'], coords[key]['long']],\n",
    "    popup=coords[key]['name'],\n",
    "    color=('crimson' if coords[key]['name'] == 'Nature' else 'blue'),\n",
    "    fill=False,\n",
    "    ).add_to(m)\n",
    "\n",
    "# pull in the Macrostrat tile layer\n",
    "folium.TileLayer(tiles='https://tiles.macrostrat.org/carto/{z}/{x}/{y}.png', \n",
    "                 attr='Macrostrat', name='Macrostrat').add_to(m)\n",
    "folium.LayerControl().add_to(m) # allow user to switch between layers\n",
    "folium.LatLngPopup().add_to(m) # click on the map to get Lat/Long in a popup"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Draw the map. **Note that this map is interactive**: you can zoom, pan, click the points of interest, and alternate between the two layers (by clicking on the tiles icon in the upper-right corner of the map). If you click anywhere on the map, a popup will appear showing the latitude and longitude of that position."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe src=\"about:blank\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" data-html=PCFET0NUWVBFIGh0bWw+CjxoZWFkPiAgICAKICAgIDxtZXRhIGh0dHAtZXF1aXY9ImNvbnRlbnQtdHlwZSIgY29udGVudD0idGV4dC9odG1sOyBjaGFyc2V0PVVURi04IiAvPgogICAgCiAgICAgICAgPHNjcmlwdD4KICAgICAgICAgICAgTF9OT19UT1VDSCA9IGZhbHNlOwogICAgICAgICAgICBMX0RJU0FCTEVfM0QgPSBmYWxzZTsKICAgICAgICA8L3NjcmlwdD4KICAgIAogICAgPHNjcmlwdCBzcmM9Imh0dHBzOi8vY2RuLmpzZGVsaXZyLm5ldC9ucG0vbGVhZmxldEAxLjYuMC9kaXN0L2xlYWZsZXQuanMiPjwvc2NyaXB0PgogICAgPHNjcmlwdCBzcmM9Imh0dHBzOi8vY29kZS5qcXVlcnkuY29tL2pxdWVyeS0xLjEyLjQubWluLmpzIj48L3NjcmlwdD4KICAgIDxzY3JpcHQgc3JjPSJodHRwczovL21heGNkbi5ib290c3RyYXBjZG4uY29tL2Jvb3RzdHJhcC8zLjIuMC9qcy9ib290c3RyYXAubWluLmpzIj48L3NjcmlwdD4KICAgIDxzY3JpcHQgc3JjPSJodHRwczovL2NkbmpzLmNsb3VkZmxhcmUuY29tL2FqYXgvbGlicy9MZWFmbGV0LmF3ZXNvbWUtbWFya2Vycy8yLjAuMi9sZWFmbGV0LmF3ZXNvbWUtbWFya2Vycy5qcyI+PC9zY3JpcHQ+CiAgICA8bGluayByZWw9InN0eWxlc2hlZXQiIGhyZWY9Imh0dHBzOi8vY2RuLmpzZGVsaXZyLm5ldC9ucG0vbGVhZmxldEAxLjYuMC9kaXN0L2xlYWZsZXQuY3NzIi8+CiAgICA8bGluayByZWw9InN0eWxlc2hlZXQiIGhyZWY9Imh0dHBzOi8vbWF4Y2RuLmJvb3RzdHJhcGNkbi5jb20vYm9vdHN0cmFwLzMuMi4wL2Nzcy9ib290c3RyYXAubWluLmNzcyIvPgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL21heGNkbi5ib290c3RyYXBjZG4uY29tL2Jvb3RzdHJhcC8zLjIuMC9jc3MvYm9vdHN0cmFwLXRoZW1lLm1pbi5jc3MiLz4KICAgIDxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iaHR0cHM6Ly9tYXhjZG4uYm9vdHN0cmFwY2RuLmNvbS9mb250LWF3ZXNvbWUvNC42LjMvY3NzL2ZvbnQtYXdlc29tZS5taW4uY3NzIi8+CiAgICA8bGluayByZWw9InN0eWxlc2hlZXQiIGhyZWY9Imh0dHBzOi8vY2RuanMuY2xvdWRmbGFyZS5jb20vYWpheC9saWJzL0xlYWZsZXQuYXdlc29tZS1tYXJrZXJzLzIuMC4yL2xlYWZsZXQuYXdlc29tZS1tYXJrZXJzLmNzcyIvPgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL3Jhd2Nkbi5naXRoYWNrLmNvbS9weXRob24tdmlzdWFsaXphdGlvbi9mb2xpdW0vbWFzdGVyL2ZvbGl1bS90ZW1wbGF0ZXMvbGVhZmxldC5hd2Vzb21lLnJvdGF0ZS5jc3MiLz4KICAgIDxzdHlsZT5odG1sLCBib2R5IHt3aWR0aDogMTAwJTtoZWlnaHQ6IDEwMCU7bWFyZ2luOiAwO3BhZGRpbmc6IDA7fTwvc3R5bGU+CiAgICA8c3R5bGU+I21hcCB7cG9zaXRpb246YWJzb2x1dGU7dG9wOjA7Ym90dG9tOjA7cmlnaHQ6MDtsZWZ0OjA7fTwvc3R5bGU+CiAgICAKICAgICAgICAgICAgPG1ldGEgbmFtZT0idmlld3BvcnQiIGNvbnRlbnQ9IndpZHRoPWRldmljZS13aWR0aCwKICAgICAgICAgICAgICAgIGluaXRpYWwtc2NhbGU9MS4wLCBtYXhpbXVtLXNjYWxlPTEuMCwgdXNlci1zY2FsYWJsZT1ubyIgLz4KICAgICAgICAgICAgPHN0eWxlPgogICAgICAgICAgICAgICAgI21hcF9iYmM3NzJiNjc1NDg0ZjAwOWYyZjcxMTAwZGIzNjdmYiB7CiAgICAgICAgICAgICAgICAgICAgcG9zaXRpb246IHJlbGF0aXZlOwogICAgICAgICAgICAgICAgICAgIHdpZHRoOiAxMDAuMCU7CiAgICAgICAgICAgICAgICAgICAgaGVpZ2h0OiAxMDAuMCU7CiAgICAgICAgICAgICAgICAgICAgbGVmdDogMC4wJTsKICAgICAgICAgICAgICAgICAgICB0b3A6IDAuMCU7CiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgIDwvc3R5bGU+CiAgICAgICAgCjwvaGVhZD4KPGJvZHk+ICAgIAogICAgCiAgICAgICAgICAgIDxkaXYgY2xhc3M9ImZvbGl1bS1tYXAiIGlkPSJtYXBfYmJjNzcyYjY3NTQ4NGYwMDlmMmY3MTEwMGRiMzY3ZmIiID48L2Rpdj4KICAgICAgICAKPC9ib2R5Pgo8c2NyaXB0PiAgICAKICAgIAogICAgICAgICAgICB2YXIgbWFwX2JiYzc3MmI2NzU0ODRmMDA5ZjJmNzExMDBkYjM2N2ZiID0gTC5tYXAoCiAgICAgICAgICAgICAgICAibWFwX2JiYzc3MmI2NzU0ODRmMDA5ZjJmNzExMDBkYjM2N2ZiIiwKICAgICAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgICAgICBjZW50ZXI6IFs1MS44NTYxLCAtMC4yOTY2XSwKICAgICAgICAgICAgICAgICAgICBjcnM6IEwuQ1JTLkVQU0czODU3LAogICAgICAgICAgICAgICAgICAgIHpvb206IDksCiAgICAgICAgICAgICAgICAgICAgem9vbUNvbnRyb2w6IHRydWUsCiAgICAgICAgICAgICAgICAgICAgcHJlZmVyQ2FudmFzOiBmYWxzZSwKICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgKTsKCiAgICAgICAgICAgIAoKICAgICAgICAKICAgIAogICAgICAgICAgICB2YXIgdGlsZV9sYXllcl9iOTlmOTZjM2JlZTE0N2RiODQ2ZTJmMWI1ZTA5OWYxMiA9IEwudGlsZUxheWVyKAogICAgICAgICAgICAgICAgImh0dHBzOi8vY2FydG9kYi1iYXNlbWFwcy17c30uZ2xvYmFsLnNzbC5mYXN0bHkubmV0L2xpZ2h0X2FsbC97en0ve3h9L3t5fS5wbmciLAogICAgICAgICAgICAgICAgeyJhdHRyaWJ1dGlvbiI6ICJcdTAwMjZjb3B5OyBcdTAwM2NhIGhyZWY9XCJodHRwOi8vd3d3Lm9wZW5zdHJlZXRtYXAub3JnL2NvcHlyaWdodFwiXHUwMDNlT3BlblN0cmVldE1hcFx1MDAzYy9hXHUwMDNlIGNvbnRyaWJ1dG9ycyBcdTAwMjZjb3B5OyBcdTAwM2NhIGhyZWY9XCJodHRwOi8vY2FydG9kYi5jb20vYXR0cmlidXRpb25zXCJcdTAwM2VDYXJ0b0RCXHUwMDNjL2FcdTAwM2UsIENhcnRvREIgXHUwMDNjYSBocmVmID1cImh0dHA6Ly9jYXJ0b2RiLmNvbS9hdHRyaWJ1dGlvbnNcIlx1MDAzZWF0dHJpYnV0aW9uc1x1MDAzYy9hXHUwMDNlIiwgImRldGVjdFJldGluYSI6IGZhbHNlLCAibWF4TmF0aXZlWm9vbSI6IDE4LCAibWF4Wm9vbSI6IDE4LCAibWluWm9vbSI6IDAsICJub1dyYXAiOiBmYWxzZSwgIm9wYWNpdHkiOiAxLCAic3ViZG9tYWlucyI6ICJhYmMiLCAidG1zIjogZmFsc2V9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2JiYzc3MmI2NzU0ODRmMDA5ZjJmNzExMDBkYjM2N2ZiKTsKICAgICAgICAKICAgIAogICAgICAgICAgICB2YXIgY2lyY2xlX21hcmtlcl9mZTkwNDA5YWJiODQ0ZDViOGQ5YzYxYTg5ODgyODk0MCA9IEwuY2lyY2xlTWFya2VyKAogICAgICAgICAgICAgICAgWzUxLjUzMzkyNSwgLTAuMTIxNTUzXSwKICAgICAgICAgICAgICAgIHsiYnViYmxpbmdNb3VzZUV2ZW50cyI6IHRydWUsICJjb2xvciI6ICJjcmltc29uIiwgImRhc2hBcnJheSI6IG51bGwsICJkYXNoT2Zmc2V0IjogbnVsbCwgImZpbGwiOiBmYWxzZSwgImZpbGxDb2xvciI6ICJjcmltc29uIiwgImZpbGxPcGFjaXR5IjogMC4yLCAiZmlsbFJ1bGUiOiAiZXZlbm9kZCIsICJsaW5lQ2FwIjogInJvdW5kIiwgImxpbmVKb2luIjogInJvdW5kIiwgIm9wYWNpdHkiOiAxLjAsICJyYWRpdXMiOiAxMCwgInN0cm9rZSI6IHRydWUsICJ3ZWlnaHQiOiAzfQogICAgICAgICAgICApLmFkZFRvKG1hcF9iYmM3NzJiNjc1NDg0ZjAwOWYyZjcxMTAwZGIzNjdmYik7CiAgICAgICAgCiAgICAKICAgICAgICB2YXIgcG9wdXBfYzM5NGViMTlhMDNmNGYwOGIxNzFlNmJkMTRhYjljN2MgPSBMLnBvcHVwKHsibWF4V2lkdGgiOiAiMTAwJSJ9KTsKCiAgICAgICAgCiAgICAgICAgICAgIHZhciBodG1sX2JjZjI0NTA5MWQ1NDQ3NzU5ZTFmNTg2ZjA1YTBjYTlkID0gJChgPGRpdiBpZD0iaHRtbF9iY2YyNDUwOTFkNTQ0Nzc1OWUxZjU4NmYwNWEwY2E5ZCIgc3R5bGU9IndpZHRoOiAxMDAuMCU7IGhlaWdodDogMTAwLjAlOyI+TmF0dXJlPC9kaXY+YClbMF07CiAgICAgICAgICAgIHBvcHVwX2MzOTRlYjE5YTAzZjRmMDhiMTcxZTZiZDE0YWI5YzdjLnNldENvbnRlbnQoaHRtbF9iY2YyNDUwOTFkNTQ0Nzc1OWUxZjU4NmYwNWEwY2E5ZCk7CiAgICAgICAgCgogICAgICAgIGNpcmNsZV9tYXJrZXJfZmU5MDQwOWFiYjg0NGQ1YjhkOWM2MWE4OTg4Mjg5NDAuYmluZFBvcHVwKHBvcHVwX2MzOTRlYjE5YTAzZjRmMDhiMTcxZTZiZDE0YWI5YzdjKQogICAgICAgIDsKCiAgICAgICAgCiAgICAKICAgIAogICAgICAgICAgICB2YXIgY2lyY2xlX21hcmtlcl84MTFiNjdmOWIyYWY0NDUwOTlmYzdmMTIwMjgzM2Y2NiA9IEwuY2lyY2xlTWFya2VyKAogICAgICAgICAgICAgICAgWzUxLjUzMTg3NywgLTAuMTI4NzY3XSwKICAgICAgICAgICAgICAgIHsiYnViYmxpbmdNb3VzZUV2ZW50cyI6IHRydWUsICJjb2xvciI6ICJibHVlIiwgImRhc2hBcnJheSI6IG51bGwsICJkYXNoT2Zmc2V0IjogbnVsbCwgImZpbGwiOiBmYWxzZSwgImZpbGxDb2xvciI6ICJibHVlIiwgImZpbGxPcGFjaXR5IjogMC4yLCAiZmlsbFJ1bGUiOiAiZXZlbm9kZCIsICJsaW5lQ2FwIjogInJvdW5kIiwgImxpbmVKb2luIjogInJvdW5kIiwgIm9wYWNpdHkiOiAxLjAsICJyYWRpdXMiOiAxMCwgInN0cm9rZSI6IHRydWUsICJ3ZWlnaHQiOiAzfQogICAgICAgICAgICApLmFkZFRvKG1hcF9iYmM3NzJiNjc1NDg0ZjAwOWYyZjcxMTAwZGIzNjdmYik7CiAgICAgICAgCiAgICAKICAgICAgICB2YXIgcG9wdXBfMGVhY2Y5OTg0N2Q4NDY4ZGE3NTMxMGQ2ZmI2ODdjNDIgPSBMLnBvcHVwKHsibWF4V2lkdGgiOiAiMTAwJSJ9KTsKCiAgICAgICAgCiAgICAgICAgICAgIHZhciBodG1sX2Q4NGQ5ZmVjMTBkMTRhMmJiZGJiZDc2ODIxMTBkNDFjID0gJChgPGRpdiBpZD0iaHRtbF9kODRkOWZlYzEwZDE0YTJiYmRiYmQ3NjgyMTEwZDQxYyIgc3R5bGU9IndpZHRoOiAxMDAuMCU7IGhlaWdodDogMTAwLjAlOyI+RnJhbmNpcyBDcmljayBJbnN0aXRpdGU8L2Rpdj5gKVswXTsKICAgICAgICAgICAgcG9wdXBfMGVhY2Y5OTg0N2Q4NDY4ZGE3NTMxMGQ2ZmI2ODdjNDIuc2V0Q29udGVudChodG1sX2Q4NGQ5ZmVjMTBkMTRhMmJiZGJiZDc2ODIxMTBkNDFjKTsKICAgICAgICAKCiAgICAgICAgY2lyY2xlX21hcmtlcl84MTFiNjdmOWIyYWY0NDUwOTlmYzdmMTIwMjgzM2Y2Ni5iaW5kUG9wdXAocG9wdXBfMGVhY2Y5OTg0N2Q4NDY4ZGE3NTMxMGQ2ZmI2ODdjNDIpCiAgICAgICAgOwoKICAgICAgICAKICAgIAogICAgCiAgICAgICAgICAgIHZhciBjaXJjbGVfbWFya2VyXzAzMmIwYzcyN2UxOTQyY2I5MTAwNzhlMjg3ZDBiMTg5ID0gTC5jaXJjbGVNYXJrZXIoCiAgICAgICAgICAgICAgICBbNTEuNTI0NDM1LCAtMC4xMzI0OTVdLAogICAgICAgICAgICAgICAgeyJidWJibGluZ01vdXNlRXZlbnRzIjogdHJ1ZSwgImNvbG9yIjogImJsdWUiLCAiZGFzaEFycmF5IjogbnVsbCwgImRhc2hPZmZzZXQiOiBudWxsLCAiZmlsbCI6IGZhbHNlLCAiZmlsbENvbG9yIjogImJsdWUiLCAiZmlsbE9wYWNpdHkiOiAwLjIsICJmaWxsUnVsZSI6ICJldmVub2RkIiwgImxpbmVDYXAiOiAicm91bmQiLCAibGluZUpvaW4iOiAicm91bmQiLCAib3BhY2l0eSI6IDEuMCwgInJhZGl1cyI6IDEwLCAic3Ryb2tlIjogdHJ1ZSwgIndlaWdodCI6IDN9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2JiYzc3MmI2NzU0ODRmMDA5ZjJmNzExMDBkYjM2N2ZiKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBwb3B1cF9lZjczM2Q0MjY0YTI0MmQ3YjI2YTdkNmEzMDViZjQzOCA9IEwucG9wdXAoeyJtYXhXaWR0aCI6ICIxMDAlIn0pOwoKICAgICAgICAKICAgICAgICAgICAgdmFyIGh0bWxfOTIxMzgxMzY0NmQ1NDliYmI2MDY2MWM0MTk5OWE4MjIgPSAkKGA8ZGl2IGlkPSJodG1sXzkyMTM4MTM2NDZkNTQ5YmJiNjA2NjFjNDE5OTlhODIyIiBzdHlsZT0id2lkdGg6IDEwMC4wJTsgaGVpZ2h0OiAxMDAuMCU7Ij5NUkMgTGFib3JhdG9yeSBmb3IgTW9sZWN1bGFyIENlbGwgQmlvbG9neTwvZGl2PmApWzBdOwogICAgICAgICAgICBwb3B1cF9lZjczM2Q0MjY0YTI0MmQ3YjI2YTdkNmEzMDViZjQzOC5zZXRDb250ZW50KGh0bWxfOTIxMzgxMzY0NmQ1NDliYmI2MDY2MWM0MTk5OWE4MjIpOwogICAgICAgIAoKICAgICAgICBjaXJjbGVfbWFya2VyXzAzMmIwYzcyN2UxOTQyY2I5MTAwNzhlMjg3ZDBiMTg5LmJpbmRQb3B1cChwb3B1cF9lZjczM2Q0MjY0YTI0MmQ3YjI2YTdkNmEzMDViZjQzOCkKICAgICAgICA7CgogICAgICAgIAogICAgCiAgICAKICAgICAgICAgICAgdmFyIGNpcmNsZV9tYXJrZXJfYjA5M2EyOTFjYTUwNGVjMmFkYTAyNjhlNWRiMzhjMjcgPSBMLmNpcmNsZU1hcmtlcigKICAgICAgICAgICAgICAgIFs1MS41MTE1NzMsIC0wLjExNjA4M10sCiAgICAgICAgICAgICAgICB7ImJ1YmJsaW5nTW91c2VFdmVudHMiOiB0cnVlLCAiY29sb3IiOiAiYmx1ZSIsICJkYXNoQXJyYXkiOiBudWxsLCAiZGFzaE9mZnNldCI6IG51bGwsICJmaWxsIjogZmFsc2UsICJmaWxsQ29sb3IiOiAiYmx1ZSIsICJmaWxsT3BhY2l0eSI6IDAuMiwgImZpbGxSdWxlIjogImV2ZW5vZGQiLCAibGluZUNhcCI6ICJyb3VuZCIsICJsaW5lSm9pbiI6ICJyb3VuZCIsICJvcGFjaXR5IjogMS4wLCAicmFkaXVzIjogMTAsICJzdHJva2UiOiB0cnVlLCAid2VpZ2h0IjogM30KICAgICAgICAgICAgKS5hZGRUbyhtYXBfYmJjNzcyYjY3NTQ4NGYwMDlmMmY3MTEwMGRiMzY3ZmIpOwogICAgICAgIAogICAgCiAgICAgICAgdmFyIHBvcHVwXzcyYzY5ZTg0MjE3NzQ0M2NhYjlkZTJlZDExNTYyNTU3ID0gTC5wb3B1cCh7Im1heFdpZHRoIjogIjEwMCUifSk7CgogICAgICAgIAogICAgICAgICAgICB2YXIgaHRtbF8wOTAzMmYwZWIyMWY0NTNlOGExY2M3Njg0Yjg3NTc1NiA9ICQoYDxkaXYgaWQ9Imh0bWxfMDkwMzJmMGViMjFmNDUzZThhMWNjNzY4NGI4NzU3NTYiIHN0eWxlPSJ3aWR0aDogMTAwLjAlOyBoZWlnaHQ6IDEwMC4wJTsiPktpbmdzIENvbGxlZ2UgTG9uZG9uPC9kaXY+YClbMF07CiAgICAgICAgICAgIHBvcHVwXzcyYzY5ZTg0MjE3NzQ0M2NhYjlkZTJlZDExNTYyNTU3LnNldENvbnRlbnQoaHRtbF8wOTAzMmYwZWIyMWY0NTNlOGExY2M3Njg0Yjg3NTc1Nik7CiAgICAgICAgCgogICAgICAgIGNpcmNsZV9tYXJrZXJfYjA5M2EyOTFjYTUwNGVjMmFkYTAyNjhlNWRiMzhjMjcuYmluZFBvcHVwKHBvcHVwXzcyYzY5ZTg0MjE3NzQ0M2NhYjlkZTJlZDExNTYyNTU3KQogICAgICAgIDsKCiAgICAgICAgCiAgICAKICAgIAogICAgICAgICAgICB2YXIgY2lyY2xlX21hcmtlcl84MTk0N2MyOWU1YzU0Y2E4YjJkODU4ZTZjOWYwZjIxZiA9IEwuY2lyY2xlTWFya2VyKAogICAgICAgICAgICAgICAgWzUxLjQ5ODc4LCAtMC4xNzQ4ODhdLAogICAgICAgICAgICAgICAgeyJidWJibGluZ01vdXNlRXZlbnRzIjogdHJ1ZSwgImNvbG9yIjogImJsdWUiLCAiZGFzaEFycmF5IjogbnVsbCwgImRhc2hPZmZzZXQiOiBudWxsLCAiZmlsbCI6IGZhbHNlLCAiZmlsbENvbG9yIjogImJsdWUiLCAiZmlsbE9wYWNpdHkiOiAwLjIsICJmaWxsUnVsZSI6ICJldmVub2RkIiwgImxpbmVDYXAiOiAicm91bmQiLCAibGluZUpvaW4iOiAicm91bmQiLCAib3BhY2l0eSI6IDEuMCwgInJhZGl1cyI6IDEwLCAic3Ryb2tlIjogdHJ1ZSwgIndlaWdodCI6IDN9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2JiYzc3MmI2NzU0ODRmMDA5ZjJmNzExMDBkYjM2N2ZiKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBwb3B1cF82NjYwZDJlYTA5OTQ0YmM3OTNmZTY4ZTNhNjE3MGEwZSA9IEwucG9wdXAoeyJtYXhXaWR0aCI6ICIxMDAlIn0pOwoKICAgICAgICAKICAgICAgICAgICAgdmFyIGh0bWxfYTJjNDFlNzM4MDY0NGM0Yzk5MTUyNTM5MzQyOWU1ODkgPSAkKGA8ZGl2IGlkPSJodG1sX2EyYzQxZTczODA2NDRjNGM5OTE1MjUzOTM0MjllNTg5IiBzdHlsZT0id2lkdGg6IDEwMC4wJTsgaGVpZ2h0OiAxMDAuMCU7Ij5JbXBlcmlhbCBDb2xsZWdlIExvbmRvbjwvZGl2PmApWzBdOwogICAgICAgICAgICBwb3B1cF82NjYwZDJlYTA5OTQ0YmM3OTNmZTY4ZTNhNjE3MGEwZS5zZXRDb250ZW50KGh0bWxfYTJjNDFlNzM4MDY0NGM0Yzk5MTUyNTM5MzQyOWU1ODkpOwogICAgICAgIAoKICAgICAgICBjaXJjbGVfbWFya2VyXzgxOTQ3YzI5ZTVjNTRjYThiMmQ4NThlNmM5ZjBmMjFmLmJpbmRQb3B1cChwb3B1cF82NjYwZDJlYTA5OTQ0YmM3OTNmZTY4ZTNhNjE3MGEwZSkKICAgICAgICA7CgogICAgICAgIAogICAgCiAgICAKICAgICAgICAgICAgdmFyIGNpcmNsZV9tYXJrZXJfMTIyNjdkY2U2MmY4NGM2YWIzNGUyMzMzYTViOGVlNzYgPSBMLmNpcmNsZU1hcmtlcigKICAgICAgICAgICAgICAgIFs1Mi4yMDY5NiwgMC4xMTUwMzRdLAogICAgICAgICAgICAgICAgeyJidWJibGluZ01vdXNlRXZlbnRzIjogdHJ1ZSwgImNvbG9yIjogImJsdWUiLCAiZGFzaEFycmF5IjogbnVsbCwgImRhc2hPZmZzZXQiOiBudWxsLCAiZmlsbCI6IGZhbHNlLCAiZmlsbENvbG9yIjogImJsdWUiLCAiZmlsbE9wYWNpdHkiOiAwLjIsICJmaWxsUnVsZSI6ICJldmVub2RkIiwgImxpbmVDYXAiOiAicm91bmQiLCAibGluZUpvaW4iOiAicm91bmQiLCAib3BhY2l0eSI6IDEuMCwgInJhZGl1cyI6IDEwLCAic3Ryb2tlIjogdHJ1ZSwgIndlaWdodCI6IDN9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2JiYzc3MmI2NzU0ODRmMDA5ZjJmNzExMDBkYjM2N2ZiKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBwb3B1cF8wYzllNTg3YmU4Mjc0Nzg0OGM1ZDJiMmFiY2ViNzE3NyA9IEwucG9wdXAoeyJtYXhXaWR0aCI6ICIxMDAlIn0pOwoKICAgICAgICAKICAgICAgICAgICAgdmFyIGh0bWxfNjgyZDU4ZDViY2ExNDdkZGJiN2VkNTNmYmQzMDExYmUgPSAkKGA8ZGl2IGlkPSJodG1sXzY4MmQ1OGQ1YmNhMTQ3ZGRiYjdlZDUzZmJkMzAxMWJlIiBzdHlsZT0id2lkdGg6IDEwMC4wJTsgaGVpZ2h0OiAxMDAuMCU7Ij5DYW1icmlkZ2UgVW5pdmVyc2l0eTwvZGl2PmApWzBdOwogICAgICAgICAgICBwb3B1cF8wYzllNTg3YmU4Mjc0Nzg0OGM1ZDJiMmFiY2ViNzE3Ny5zZXRDb250ZW50KGh0bWxfNjgyZDU4ZDViY2ExNDdkZGJiN2VkNTNmYmQzMDExYmUpOwogICAgICAgIAoKICAgICAgICBjaXJjbGVfbWFya2VyXzEyMjY3ZGNlNjJmODRjNmFiMzRlMjMzM2E1YjhlZTc2LmJpbmRQb3B1cChwb3B1cF8wYzllNTg3YmU4Mjc0Nzg0OGM1ZDJiMmFiY2ViNzE3NykKICAgICAgICA7CgogICAgICAgIAogICAgCiAgICAKICAgICAgICAgICAgdmFyIGNpcmNsZV9tYXJrZXJfZTc2NDJjNGFhOTBjNGZlNjllNDJkMzg3OTgyY2U1MjUgPSBMLmNpcmNsZU1hcmtlcigKICAgICAgICAgICAgICAgIFs1MS43NTQ4NDMsIC0xLjI1NDMwMl0sCiAgICAgICAgICAgICAgICB7ImJ1YmJsaW5nTW91c2VFdmVudHMiOiB0cnVlLCAiY29sb3IiOiAiYmx1ZSIsICJkYXNoQXJyYXkiOiBudWxsLCAiZGFzaE9mZnNldCI6IG51bGwsICJmaWxsIjogZmFsc2UsICJmaWxsQ29sb3IiOiAiYmx1ZSIsICJmaWxsT3BhY2l0eSI6IDAuMiwgImZpbGxSdWxlIjogImV2ZW5vZGQiLCAibGluZUNhcCI6ICJyb3VuZCIsICJsaW5lSm9pbiI6ICJyb3VuZCIsICJvcGFjaXR5IjogMS4wLCAicmFkaXVzIjogMTAsICJzdHJva2UiOiB0cnVlLCAid2VpZ2h0IjogM30KICAgICAgICAgICAgKS5hZGRUbyhtYXBfYmJjNzcyYjY3NTQ4NGYwMDlmMmY3MTEwMGRiMzY3ZmIpOwogICAgICAgIAogICAgCiAgICAgICAgdmFyIHBvcHVwXzk1ZjEzNmQ2MzgzYjRkMTU4ZDM3NzgyOTY3N2M2YmU4ID0gTC5wb3B1cCh7Im1heFdpZHRoIjogIjEwMCUifSk7CgogICAgICAgIAogICAgICAgICAgICB2YXIgaHRtbF80NzRhY2IwZTVhNDQ0MWRlYTRlYjFlZTFmYzIzNjBhNCA9ICQoYDxkaXYgaWQ9Imh0bWxfNDc0YWNiMGU1YTQ0NDFkZWE0ZWIxZWUxZmMyMzYwYTQiIHN0eWxlPSJ3aWR0aDogMTAwLjAlOyBoZWlnaHQ6IDEwMC4wJTsiPk94Zm9yZCBVbml2ZXJzaXR5PC9kaXY+YClbMF07CiAgICAgICAgICAgIHBvcHVwXzk1ZjEzNmQ2MzgzYjRkMTU4ZDM3NzgyOTY3N2M2YmU4LnNldENvbnRlbnQoaHRtbF80NzRhY2IwZTVhNDQ0MWRlYTRlYjFlZTFmYzIzNjBhNCk7CiAgICAgICAgCgogICAgICAgIGNpcmNsZV9tYXJrZXJfZTc2NDJjNGFhOTBjNGZlNjllNDJkMzg3OTgyY2U1MjUuYmluZFBvcHVwKHBvcHVwXzk1ZjEzNmQ2MzgzYjRkMTU4ZDM3NzgyOTY3N2M2YmU4KQogICAgICAgIDsKCiAgICAgICAgCiAgICAKICAgIAogICAgICAgICAgICB2YXIgY2lyY2xlX21hcmtlcl81NjFkMjc3MjYyNGE0ZDlmYTgyYTYzZWQ0MTM5ODUwZiA9IEwuY2lyY2xlTWFya2VyKAogICAgICAgICAgICAgICAgWzUxLjUzMjM0OSwgLTAuMTIzODA2XSwKICAgICAgICAgICAgICAgIHsiYnViYmxpbmdNb3VzZUV2ZW50cyI6IHRydWUsICJjb2xvciI6ICJibHVlIiwgImRhc2hBcnJheSI6IG51bGwsICJkYXNoT2Zmc2V0IjogbnVsbCwgImZpbGwiOiBmYWxzZSwgImZpbGxDb2xvciI6ICJibHVlIiwgImZpbGxPcGFjaXR5IjogMC4yLCAiZmlsbFJ1bGUiOiAiZXZlbm9kZCIsICJsaW5lQ2FwIjogInJvdW5kIiwgImxpbmVKb2luIjogInJvdW5kIiwgIm9wYWNpdHkiOiAxLjAsICJyYWRpdXMiOiAxMCwgInN0cm9rZSI6IHRydWUsICJ3ZWlnaHQiOiAzfQogICAgICAgICAgICApLmFkZFRvKG1hcF9iYmM3NzJiNjc1NDg0ZjAwOWYyZjcxMTAwZGIzNjdmYik7CiAgICAgICAgCiAgICAKICAgICAgICB2YXIgcG9wdXBfZThlZDFiZmMwZWE3NDE5ZmI0MDliZTMxZmY4MTM5NWMgPSBMLnBvcHVwKHsibWF4V2lkdGgiOiAiMTAwJSJ9KTsKCiAgICAgICAgCiAgICAgICAgICAgIHZhciBodG1sXzAyMjRhZWI3OGE3ZjRiMTI4YzQ4YjM1NTgxNjk4NDkwID0gJChgPGRpdiBpZD0iaHRtbF8wMjI0YWViNzhhN2Y0YjEyOGM0OGIzNTU4MTY5ODQ5MCIgc3R5bGU9IndpZHRoOiAxMDAuMCU7IGhlaWdodDogMTAwLjAlOyI+UGxhdGZvcm0gOSAzLzQ8L2Rpdj5gKVswXTsKICAgICAgICAgICAgcG9wdXBfZThlZDFiZmMwZWE3NDE5ZmI0MDliZTMxZmY4MTM5NWMuc2V0Q29udGVudChodG1sXzAyMjRhZWI3OGE3ZjRiMTI4YzQ4YjM1NTgxNjk4NDkwKTsKICAgICAgICAKCiAgICAgICAgY2lyY2xlX21hcmtlcl81NjFkMjc3MjYyNGE0ZDlmYTgyYTYzZWQ0MTM5ODUwZi5iaW5kUG9wdXAocG9wdXBfZThlZDFiZmMwZWE3NDE5ZmI0MDliZTMxZmY4MTM5NWMpCiAgICAgICAgOwoKICAgICAgICAKICAgIAogICAgCiAgICAgICAgICAgIHZhciB0aWxlX2xheWVyXzlhZDBiMTg0NzJkMDRhYmM5YTZlMjY4NzgwZGQ2YzkyID0gTC50aWxlTGF5ZXIoCiAgICAgICAgICAgICAgICAiaHR0cHM6Ly90aWxlcy5tYWNyb3N0cmF0Lm9yZy9jYXJ0by97en0ve3h9L3t5fS5wbmciLAogICAgICAgICAgICAgICAgeyJhdHRyaWJ1dGlvbiI6ICJNYWNyb3N0cmF0IiwgImRldGVjdFJldGluYSI6IGZhbHNlLCAibWF4TmF0aXZlWm9vbSI6IDE4LCAibWF4Wm9vbSI6IDE4LCAibWluWm9vbSI6IDAsICJub1dyYXAiOiBmYWxzZSwgIm9wYWNpdHkiOiAxLCAic3ViZG9tYWlucyI6ICJhYmMiLCAidG1zIjogZmFsc2V9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2JiYzc3MmI2NzU0ODRmMDA5ZjJmNzExMDBkYjM2N2ZiKTsKICAgICAgICAKICAgIAogICAgICAgICAgICB2YXIgbGF5ZXJfY29udHJvbF9jMjM3YTEzYTI4OGI0MjViOTA3MzgwYzExZTE2ZjY2YyA9IHsKICAgICAgICAgICAgICAgIGJhc2VfbGF5ZXJzIDogewogICAgICAgICAgICAgICAgICAgICJjYXJ0b2RicG9zaXRyb24iIDogdGlsZV9sYXllcl9iOTlmOTZjM2JlZTE0N2RiODQ2ZTJmMWI1ZTA5OWYxMiwKICAgICAgICAgICAgICAgICAgICAiTWFjcm9zdHJhdCIgOiB0aWxlX2xheWVyXzlhZDBiMTg0NzJkMDRhYmM5YTZlMjY4NzgwZGQ2YzkyLAogICAgICAgICAgICAgICAgfSwKICAgICAgICAgICAgICAgIG92ZXJsYXlzIDogIHsKICAgICAgICAgICAgICAgIH0sCiAgICAgICAgICAgIH07CiAgICAgICAgICAgIEwuY29udHJvbC5sYXllcnMoCiAgICAgICAgICAgICAgICBsYXllcl9jb250cm9sX2MyMzdhMTNhMjg4YjQyNWI5MDczODBjMTFlMTZmNjZjLmJhc2VfbGF5ZXJzLAogICAgICAgICAgICAgICAgbGF5ZXJfY29udHJvbF9jMjM3YTEzYTI4OGI0MjViOTA3MzgwYzExZTE2ZjY2Yy5vdmVybGF5cywKICAgICAgICAgICAgICAgIHsiYXV0b1pJbmRleCI6IHRydWUsICJjb2xsYXBzZWQiOiB0cnVlLCAicG9zaXRpb24iOiAidG9wcmlnaHQifQogICAgICAgICAgICApLmFkZFRvKG1hcF9iYmM3NzJiNjc1NDg0ZjAwOWYyZjcxMTAwZGIzNjdmYik7CiAgICAgICAgICAgIHRpbGVfbGF5ZXJfOWFkMGIxODQ3MmQwNGFiYzlhNmUyNjg3ODBkZDZjOTIucmVtb3ZlKCk7CiAgICAgICAgCiAgICAKICAgICAgICAgICAgICAgIHZhciBsYXRfbG5nX3BvcHVwXzdhNWE0MmU4MmE1MjQ3MDRiZTdmMmU3NzIwOWI0NTI2ID0gTC5wb3B1cCgpOwogICAgICAgICAgICAgICAgZnVuY3Rpb24gbGF0TG5nUG9wKGUpIHsKICAgICAgICAgICAgICAgICAgICBsYXRfbG5nX3BvcHVwXzdhNWE0MmU4MmE1MjQ3MDRiZTdmMmU3NzIwOWI0NTI2CiAgICAgICAgICAgICAgICAgICAgICAgIC5zZXRMYXRMbmcoZS5sYXRsbmcpCiAgICAgICAgICAgICAgICAgICAgICAgIC5zZXRDb250ZW50KCJMYXRpdHVkZTogIiArIGUubGF0bG5nLmxhdC50b0ZpeGVkKDQpICsKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIjxicj5Mb25naXR1ZGU6ICIgKyBlLmxhdGxuZy5sbmcudG9GaXhlZCg0KSkKICAgICAgICAgICAgICAgICAgICAgICAgLm9wZW5PbihtYXBfYmJjNzcyYjY3NTQ4NGYwMDlmMmY3MTEwMGRiMzY3ZmIpOwogICAgICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgICAgIG1hcF9iYmM3NzJiNjc1NDg0ZjAwOWYyZjcxMTAwZGIzNjdmYi5vbignY2xpY2snLCBsYXRMbmdQb3ApOwogICAgICAgICAgICAKPC9zY3JpcHQ+ onload=\"this.contentDocument.open();this.contentDocument.write(atob(this.getAttribute('data-html')));this.contentDocument.close();\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
      ],
      "text/plain": [
       "<folium.folium.Map at 0x120ae1910>"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "m"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Document our session, for [computational reproducibility](https://www.nature.com/articles/d41586-018-05990-5)!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'commit_hash': '8bda98619',\n",
      " 'commit_source': 'installation',\n",
      " 'default_encoding': 'UTF-8',\n",
      " 'ipython_path': '/opt/anaconda3/lib/python3.7/site-packages/IPython',\n",
      " 'ipython_version': '7.12.0',\n",
      " 'os_name': 'posix',\n",
      " 'platform': 'Darwin-18.7.0-x86_64-i386-64bit',\n",
      " 'sys_executable': '/opt/anaconda3/bin/python',\n",
      " 'sys_platform': 'darwin',\n",
      " 'sys_version': '3.7.6 (default, Jan  8 2020, 13:42:34) \\n'\n",
      "                '[Clang 4.0.1 (tags/RELEASE_401/final)]'}\n"
     ]
    }
   ],
   "source": [
    "import IPython\n",
    "print(IPython.sys_info())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "biopython==1.76\n",
      "folium==0.11.0\n",
      "matplotlib==3.1.3\n"
     ]
    }
   ],
   "source": [
    "!pip freeze | grep -E 'folium|matplotlib|biopython'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}