Taller3-ManejoArchivosBioPythonSimon.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Taller 3 -Manejo de Archivos y BioPython\n",
    "## Fundamentos en Biología Computacional\n",
    "## 2019-1, Prof. Javier Correa\n",
    "\n",
    "Nota: Se calificará que el código sea lo más simple, claro, original y eficiente. No usar funciones no realizadas por ustedes mismos."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Competencias: 1. Resuelve problemas de ciencias de los datos utilizando razonamiento lógico y programación en python. 2. Manipula archivos con formatos utilizados en bioinformática para extaeer información relevante. 3. Utiliza efectivamente módulos y paquetes de python especializados para biología.\n",
    "\n",
    "Tener en cuenta: a. desarrollar y aplicar funciones en python, b. manipular archivos (entrada y salida), c. utilizar librerias especializadas (Biopython).<br>\n",
    "Ejercicio:<br>\n",
    "\n",
    "1. Un investigador desea desarrollar un Script en python, que le permita encontrar el grupo monifilético en común entre una lista de especies de plantas (archivo adjunto: lista_mercado.tsv). La información sobre la taxonomía de cada una de las especies fue descargada de la base de datos del GenBank, luego concatenada y se encuentra en formato .gb anexo <br> "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Especie\tid\n",
      "Quinua\t63459\n",
      "Tomate\t4081\n",
      "Oliva\t4146\n",
      "Yuca\t3983\n",
      "Papa\t4113\n",
      "Repollo\t3712\n",
      "Lechuga\t4236\n",
      "Durazno\t3760\n",
      "Piña\t4615\n",
      "Coco\t13894\n",
      "Lenteja\t3864\n",
      "Arveja\t3888\n",
      "Berenjena\t4111\n",
      "Garbanzo\t3827\n"
     ]
    }
   ],
   "source": [
    "%%bash\n",
    "cat lista_mercado.tsv"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "LOCUS       MG859922              122198 bp    DNA     circular PLN 30-APR-2018\n",
      "DEFINITION  Pisum sativum subsp. elatius isolate CE1 chloroplast, complete\n",
      "            genome.\n",
      "ACCESSION   MG859922\n",
      "VERSION     MG859922.1\n",
      "KEYWORDS    .\n",
      "SOURCE      chloroplast Pisum sativum subsp. elatius\n",
      "  ORGANISM  Pisum sativum subsp. elatius\n",
      "            Eukaryota; Viridiplantae; Streptophyta; Embryophyta; Tracheophyta;\n",
      "            Spermatophyta; Magnoliophyta; eudicotyledons; Gunneridae;\n",
      "            Pentapetalae; rosids; fabids; Fabales; Fabaceae; Papilionoideae;\n",
      "            Fabeae; Pisum.\n",
      "REFERENCE   1  (bases 1 to 122198)\n",
      "  AUTHORS   Bogdanova,V.S., Mglinets,A.V., Shatskaya,N.V., Kosterin,O.E. and\n",
      "            Vasiliev,G.V.\n",
      "  TITLE     Cryptic divergences in the genus Pisum L., as revealed by\n",
      "            comparative analysis of complete plastid genomes\n",
      "  JOURNAL   Unpublished\n",
      "REFERENCE   2  (bases 1 to 122198)\n",
      "  AUTHORS   Bogdanova,V.S., Mglinets,A.V., Shatskaya,N.V., Kosterin,O.E. and\n"
     ]
    }
   ],
   "source": [
    "%%bash\n",
    "head -20 Mercado.gb"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### $\\rightarrow$ Con el archivo \"Mercado.gb\";<br>\n",
    "- A. Desarrolle un script con el cual pueda extrer la información relacionada con la taxonomía de cada especie.<br>\n",
    "- B. Organice esta información en una objeto tipo lista o diccionario.<br>\n",
    "- C. Cree una función que reciba el objeto anterior y compare las clasificaciones entre las diferentes especies para retornar el ancestro monofilético común más reciente a todas ellas. <br> Ejemplo: Embryophyta"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Magnoliophyta\n"
     ]
    }
   ],
   "source": [
    "from Bio import SeqIO\n",
    "\n",
    "mercado = list(SeqIO.parse(\"Mercado.gb\",\"gb\")) \n",
    "#tax_list = [plant.annotations[\"taxonomy\"] for plant in mercado] \n",
    "tax_list = []\n",
    "for plant in mercado:\n",
    "    tax_list.append(plant.annotations[\"taxonomy\"])\n",
    "\n",
    "found = False\n",
    "for i,rank in enumerate(tax_list[0]):\n",
    "    for plant in tax_list:\n",
    "        if plant[i] != rank:\n",
    "            comun = plant[i-1]\n",
    "            found = True\n",
    "            break\n",
    "    if found:\n",
    "        break\n",
    "print(comun)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "2 . El mismo investigador está ahora interesado en hacer una caracterización funcional de los genomas de los cloroplastos que poseen estas especies. Para esto y continuando con el archivo del punto anterior (Mercado.gb), resuelva las siguientes preguntas:\n",
    "\n",
    "- A. Una forma de diferenciar una especie de otra es a través de la estimación del contenido de Guaninas y Citosinas en los genomas. Utilizando funciones de biopython responder; ¿Cuales son los valores de %GC para cada especie? retorne un diccionario con las especies como llave y su respectivo valor de %GC. Consulte el diccionario y muestre ¿Qué especie presenta el mayor %GC? ¿Y cual el más bajo? Represente los resultados en un dendograma donde se relacione las proximidades en los contenidos de %GC (puede usar Ploty, numpy, etc).\n",
    "\n",
    "- B. Otra forma de caracterizar los genomas es estimando la frecuencia con las que aparecen los codones en los genes codificantes. Usando biopython hacer una función que me extraiga los genes codificantes para cada genoma del cloroplasto del archivo .gb, definir una función que me permita estimar la frecuencia de cada uno de los codones usados por la célula. Con las herramientas anteriores y usando funciones accesorias grafique la frecuencia de estos codones versus el codón encontrado para cada una de las especies. Finalmente, responda: ¿Qué especie presenta el codón más inusual? ¿Qué aminoácido es el mas prevalente y menos prevalente?\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "#GC = lambda x: (x.count('G')+x.count('C'))/len(x)\n",
    "def GC(seq):\n",
    "    return (seq.count(\"G\")+seq.count(\"C\"))/len(seq)\n",
    "\n",
    "#gcs = {plant.annotations[\"organism\"]:GC(plant.seq) for plant in mercado}\n",
    "gcs = {}\n",
    "for plant in mercado:\n",
    "    sp = plant.annotations[\"organism\"]\n",
    "    gcs[sp] = GC(plant.seq) "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'Pisum sativum subsp. elatius': 0.3483117563298908, 'Solanum melongena var. insanum': 0.3771035546699235, 'Cocos nucifera': 0.374382835724441, 'Prunus persica': 0.3675961721275112, 'Cicer arietinum': 0.33906271195908044, 'Lactuca sativa': 0.37549420050794646, 'Lens culinaris': 0.3442549627135736, 'Olea europaea subsp. laperrinei': 0.37812247411569994, 'Solanum tuberosum': 0.3787344184609045, 'Ananas comosus': 0.3737314891377885, 'Chenopodium quinoa': 0.3724374706454366, 'Brassica oleracea var. capitata': 0.3636138387908663, 'Solanum lycopersicum': 0.3786184803026014, 'Manihot esculenta': 0.3587297851387091}\n"
     ]
    }
   ],
   "source": [
    "print(gcs)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.3483117563298908\n",
      "0.3771035546699235\n",
      "0.374382835724441\n",
      "0.3675961721275112\n",
      "0.33906271195908044\n",
      "0.37549420050794646\n",
      "0.3442549627135736\n",
      "0.37812247411569994\n",
      "0.3787344184609045\n",
      "0.3737314891377885\n",
      "0.3724374706454366\n",
      "0.3636138387908663\n",
      "0.3786184803026014\n",
      "0.3587297851387091\n"
     ]
    }
   ],
   "source": [
    "for i in gcs.values():\n",
    "    print(i)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "El mayor GC% lo tiene Solanum tuberosum y es 37.87344184609045\n",
      "El menor GC% lo tiene Cicer arietinum y es 33.90627119590805\n"
     ]
    }
   ],
   "source": [
    "def get_key(dicc,value):\n",
    "    index = list(dicc.values()).index(value)\n",
    "    return list(dicc.keys())[index]\n",
    "\n",
    "maxgc = max(gcs.values())\n",
    "mingc = min(gcs.values())\n",
    "\n",
    "print(f\"El mayor GC% lo tiene {get_key(gcs,maxgc)} y es {maxgc*100}\")\n",
    "print(f\"El menor GC% lo tiene {get_key(gcs,mingc)} y es {mingc*100}\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Aquí graficamos el dendograma, con ayuda de [dendrogram](https://docs.scipy.org/doc/scipy/reference/generated/scipy.cluster.hierarchy.dendrogram.html), [linkage](https://docs.scipy.org/doc/scipy/reference/generated/scipy.cluster.hierarchy.linkage.html#scipy.cluster.hierarchy.linkage), [pdist](https://docs.scipy.org/doc/scipy/reference/generated/scipy.spatial.distance.pdist.html#scipy.spatial.distance.pdist) y [expand_dims](https://docs.scipy.org/doc/numpy/reference/generated/numpy.expand_dims.html)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "from scipy.cluster import hierarchy\n",
    "from scipy.spatial.distance import pdist\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjYAAAD8CAYAAABkWu7OAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3Xe8XFW5//HPl94CAQKEQEiU0EFyJQjyo4QLFhCE60WE0BKFoNiQotI0gBS9UgULogQJoUmRFgmIobcEDyCC0hJCCaQ3AqZ8f3+sNbAzzpwz5+TkzDmT5/16zYuZvddae+090XnO2muvR7YJIYQQQmgEy9W7AyGEEEII7SUCmxBCCCE0jAhsQgghhNAwIrAJIYQQQsOIwCaEEEIIDSMCmxBCCCE0jAhsQgghhNAwIrAJIYQQQsOIwCaEEEIIDWOF1hTu0aOH+/btu5S6EkIIjWncuHFTbK9X736EsCxoVWDTt29fxo4du7T6EkIIDUnShHr3IYRlRdyKCiGEEELDiMAmhBBCq0kaL2nv/P5USVcuxWPNkfTxpdV+aCytuhUVQgihfUkaDxxt+75696WtbJ+7lNtfY2m2HxpLjNiEEELolCTFH9+h1SKwCSGETkrSfpKaJM2Q9KikTxT2jZd0kqRnJc2UdIOkVfK+HpLuzPWmSXpIUsX/v5e0jaR7c7l3JJ2atw+X9JNCuYGS3qjSxjBJI/L7vpIs6ShJr0uaIum0QtlPSXos9+1tSZdJWqmw35K+Kekl4KXCtn75/b6S/iFptqQ3JZ20BJc4NKCIhkPoSFdcASNH1rsXoQuQ9Eng98D+wFjgcOB2SVvY/iAXOxj4PPA+8AgwGPg1cCLwBlB6xHxnwBWO0Q24D/h5Ps6KwNbtdAq7AlsAmwNPSrrF9gvAQuB7+Zw2BkYBxwEXF+oeCOwEzKvQ7u+Ag20/JGlt4GPt1N/QIGLEJoSONHIkNDXVuxehazgG+I3tJ2wvtH018AEpSCm51PZbtqcBdwD98/b5wIZAH9vzbT9k+z8CG2A/YJLtC2y/b3u27Sfaqf9n2p5n+xngGWB7ANvjbD9ue4Ht8cBvgD3K6p5ne5rtSoHNfGBrSWvanm776Xbqb2gQMWITQkfr3x/GjKl3L0JHktpSqw9wlKRvF7atBPQqfJ5UeP9eYd//AcOA0UrHvsL2+RWO0Rt4pS2dq0F539YAkLQ5cCEwAFiN9Ds0rqzuxGba/V/gdOB8Sc8CP7T9WHt1OnR9MWITQgid00TgHNvdC6/VbF/XUsU88nKi7Y+TbjGdIGmvKsfYtEozc0mBR0nP1p5AFb8CXgQ2s70mcCpQHvlVGl1KO+ynbB8ArA/cBtzYTv0KDSICmxBCqL8VJa1SeK0A/Bb4uqSdlKwu6Qt5Xkyz8qTjfkrDNbNI81oWVih6J9BT0vGSVpbUTdJOeV8TsK+kdST1BI5vn1OlW+7THElbAt+otaKklSQdJmkt2/P56NxC+FAENiGEUH93kybKll7DbI8lzbO5DJgOvEyaHFyLzUiTgucAjwG/tD2mvJDt2cBnSKM6k0hPIe2Zd19DmhszHhgN3NDqs6rsJGAQMJsUvLW23SOA8ZJmAV8nTaoO4UOqPJ+ssgEDBjhyRYWwBAYOTP+NOTbLFEnjbA+odz9CWBbEiE0IIYQQGkY8FRWWDZ1l/ZimpvRUVAghhKUiRmzCsiHWjwkhhGVCjNiEZUdnWD+mNMcmhBDCUhEjNiGEEEJoGBHYhBBCJyVplKSj6t2P9iBpE0lzJC1f7760haTBkh6udz+KJI2RdHQb63bp76M5EdiEEEId5Szd8/KPzDuSrpK0BoDtfXKOqC4nn9fepc+2X7e9hu1YUK8OlqXvIwKbEEKov/1trwF8EtiRlAsphNAGEdiEEEInYftNYBSwLSx+qyGnSHhA0kxJUyTdkLf3leSchoEK9QZLekTSRZJmSHpV0i55+0RJ7zZ3uyuXe1XSbEmvSTosb99U0v2Spub+XCupe953DbAJcEceifp+sZ+SDpE0tuw435N0e3n/C314uPDZko6T9FLu19m5P49JmiXpRkkrVTmfNl3HjzbpF7nui8X8W81cp9L1r1ivJZJ2lvRo/u6ekTSwSrk2fx+5zGIjOpKGSRqR368iaURue4akpyRtUOs5dLR4KiqEjtbUFE9HhYok9Qb2BW6psPtsUmqDPUlZvluzkvFOwJXAusCZwPXAHUA/YA/gZkk3255T1p/VgUuBHW3/U9KGwDql3cB5wIPAmsDNpIzix9s+QtJuwNG278tt9S00fTvwW0mb2X4pbxsEXNCKc/o8sAMpQ/nTwC7AYcBUUhqJQ4FKt/GW9Dr+EegBfAm4RdLHgA+ofp2q1rM9rbmDSdoIuIuURuLPwF6k72pL25PLi9P276MlRwFrka71B0B/UuqPTilGbELoSIMGxQJ9oZLbJM0AHgYeAM6tUGY+0AfoZft9262ZyPqa7avyfIobSD9QZ9n+wPZo4N+kIKeSRcC2kla1/bbt5wFsv2z73tzGZOBCUpDUItvvAX8iBR9I2gzYkhTw1Oqntmfl/vwdGG37VdszSaNe/1Wl3pJcx3eBi23Pt30D8E/gC3lfxetUQ73mHA7cbftu24ts3wuMJQW/i1mS76MG80lBcT/bC22Psz2rndpudzFiE0JHGjo0vcKyRWqpxIGlv6Sb8X3SaMOTkqYDF9j+fY09eKfwfh6A7fJta5RXsj1X0ldIiSt/J+kR4ETbL0panzRKsRspY/dypGSdtRpJGqE5izRac1sOeGpV3v/yzz2r1FuS6/imF0+wOIEUIFW9Ts3Vq+F4fYAvS9q/sG1F4K/lBdvh+2jONaRg+Pp8e2sEcFrOsN7pxIhNCCF0AbYn2T7Gdi/gWOCXkvoBc3OR1QrFq/2ot+W499j+DLAh8CIpIzek2x4GPmF7TdLoQjGCaynD8migh6T+pJGbYs6TuSy981mS67iRtFiUugnwVm632nVqtl4LJgLX2O5eeK1u+/wKZZf0+6h6zfNI05m2tybd8tsPOLKG/tdFh43YdJZUPWEZ1XQxY/ofX+9ehNBmkr4MPGb7DdJf4gYW2p4s6U3gcEm/Ic2H2LSdjrkBaX7IX0ijIHOA0uPB3YCZwIw8F+TksurvAB+v1rbtBZL+CPwfaT7KvYXdTcCXJF1JGtn4GouPyLTZEl7H9YHvSPolcCCwFXB3C9epar0aujsCeErS54D7SKM1OwMv5/4XLdH3Qbrmh0gaBWwPHESa14OkPYEpwD+AWaRbU532MfEOG7GJVD0hhLBEdgSekDSHNBflu7Zfy/uOIf2QTQW2AR5tp2MuB5xIGl2YRpqzcVzedybp8fSZpAmu5ROezwNOz0/RnFSl/ZHA3sBNthcUtl9EmvfzDmkC8LVLfiofWpLr+ASwGelH/hzgINtTaf46NVcPSb+W9OtKHbU9ETgAOBWYTBrBOZnKv91L+n2cQQrkpue2ikMRPUmTn2cBL5DmgY2o1OfOQIvf9mvegAEDPHbs2JYLVlB6CKTeqXrCMir+AYY6kjTOdmuevgkNQtJg0tNIu9a7L8uKmGMTQgghhIYRgU0IIYQQGkYENiGEEMJSYnt43IbqWBHYhBBCCKFhRGATQgghhIYRgU0IIXRikk7N67l0OpKer5aUsYa6h0ka3c5dCiECmxBCqDdJgySNzZmX35Y0StKuALbPtX10S23Ug+1tbI9pqVylzNm2r7X92aXawbBMisAmhBDqSNIJwMWkxJcbkJbb/yVpYbaldcwlWnV+SeuHsDTFP86w7Ghq+mihvhA6AUlrkZJADrFdXCn2jvxC0jBSVuXD8+ddgZ8BWwOzgTNsD5e0MmlV24OBlYFbge/ZnpdvF40AfgF8j5S+4IiyvmxKym+0PSnNwD3AN23PyPvHA78CDgO2kLQ68DJp8bn7JC1HSjB5DNCdlF7g67anAQ/mw8zIKZM+A2xBYeE6SQa+QVrBtwdp5dtv2XaFa9AXeA1YMadmGEPKjP7fwCdISSIHk5JC7k/Kpv1l2+Nb+k5C1xcjNmHZMGgQ9O9f716EUO7TwCqkIKRFkjYBRpEClPWA/qQcPwA/BTbP2/oBGwE/KlTvScrJ1AeolGJepGX3e5FyGfUGhpWVORT4AtC9LAUCwHdIeZD2yG1MBy7P+3bP/+1uew3bj1U5xf1IKQ+2JwVon6tSrpJDSMHaRqTUAI8BV5HO+QXgx61oK3RhMWITlg1Dh6ZXCPWwWGLnxawLTKkQJFRzGHCf7evy56nA1Jw5+hhSZudp6ZA6lzTqcUouuwj4se0PKjVs+2XSCAzAZEkX8p/BwKU5f1Elx5JGWN7Ixx8GvC7piCrlKzk/jxDNkPRXUpD25xrrXmX7lXzsUcDWtu/Ln28Czm5FP0IXFoFNCCHUz1Sgh6QVagxuegOvVNi+HrAaME4fBVECli+UmWz7/WoNS1qfdOtmN1Km6OVIoy5F1YIaSCNBt0paVNi2kDRvqFaTCu/fA9ZoRd1i9u95FT63pq3QhcWtqBBCqJ/HgPdJt3BqMZF0m6XcFNKP9za2u+fXWraLP+YtZTw+L5f5hO01gcNJwVFRc21MBPYpHL+77VVsv1nDsVsylxS4lfRcwvZCA4vAJoQQ6sT2TNI8mMslHShpNUkrStpH0s8qVLkW2FvSwZJWkLSupP62F5Em/l6UR16QtJGk1sxR6QbMId0G2gg4uZWn82vgHEl98vHXk1R6smsy6VbYx1vZZkkTsLukTfKE61NaqhCWXRHYhBBCHdm+EDgBOJ0UAEwEvgXcVqHs68C+pCeHppF+8LfPu39AmiPzuKRZwH2kJ49qdSbwSWAmcBdwS/PF/8MlwO3AaEmzgceBnXK/3yM9sfWIpBmSdm5Nw7bvBW4AngXGAXe2sm9hGSK79hHCAQMGeOzYsW06UOkp2zFj2lQ9hBC6LEnjbA+odz9CWBbEiE0IIYQQGkYENiGEEEJoGBHYhBBCCKFhRGATQgghhIYRgU0IISzjJB0maXQHH3O8pL078phhcZLOkPTrevejvUVgE0IIdZR/4OdJmiNpuqS7JPXuyD7Yvtb2ZzvymKH+bJ9t++sAkvrlRKQ1k/RGTrDaqURgE0II9bd/XiV4Q1IqgF9UKyhp+Wr7lgWSGioVkJL4LW5HcTFDCKGTyLmc/ghsXdomabikX0m6W9JcYE9JX5D0N0mzJE3MCSdL5VeRNELS1LwY3lOSNsj7Bkt6VdJsSa9JOqyw/eFCG9tIulfSNEnvSDo1b/+UpMdyu29LukzSStXOR9IXJT2fy4+RtFWVcstJ+qGkV3K/b5S0Tt7XV5IlfU3S68D9eftNkiZJminpQUnbFNpbVdIFkibk/Q9LWjXv21nSo7lPzxRHHCQNkfRCvj6vSjq2Sn9XzvW3LWxbL4+8rS9pbUl3SpqcR+HulLRxoewYSedIeoSUE6vFFZkl7S7p8Xw+E5WTi+Zr3JT7/LqkMwp1+uVrd4ykt/Lre4X9P5E0PH98MG+bk187StpM0l/zdzJF0jVKKz8j6TpSFvdRufwJ+Xv8Y/5emv3Ol6aGinxD6EquGHcFI58bWe9uhE5E0mrAV0ir9hYNIq04vB+wErAzcCTwPLAtcK+kJtu3AUcBa5ESZn5AypA9T9LqpCSXO9r+p6QNgXUq9KEbadXinwP7AyvyUaC1EPgeMBbYGBgFHAdcXKGdzYHrSHmwxuR6d0ja2va/y4p/J5fbg7T68qXA5cChhTJ7AFuRUjOQj/1V4N/AT0npJvrnfT8HtgF2ISXW3AlYpJQq4i7gCFLW8L2AmyVtaXsy8G6+xq8Cu5N+tJ+y/XSxs7Y/kHRL7t9pefPBwAO235W0LnBV3rY88HvgMhbPCXYEsA/wT/4zJ1f5tfxY7vfXSCtCdyddf0hpMA4H/gFsB/xF0t9sF1dn3h3oB2wG3J/3jyk7zO7AS8X8Yvk7/AnwUD7mrcAZwEm2D5W0G3B4qS2lkac7gSHAfNL3cA3QoYtTxohNCHUy8rmRNE1qqnc3Qudwm6QZwCzgM8D/le3/k+1HbC+y/b7tMbafy5+fJQUQe+Sy84F1gX62F9oeZ3tW3rcI2FbSqrbftv18hb7sB0yyfUE+1mzbTwDkth63vcD2eOA3heOW+wpwl+17bZd+5FYlBRvljgVOs/2G7Q+AYcBBWvy20zDbc23Py335fe5bqfz2ktbKP65fBb5r+818DR7N5Q4H7rZ9d75295KCtH1zm3fZfsXJA8BoUrbzSkayeOA1KG/D9lTbN9t+z/ZsUjqJ8us03Pbz+VrOr3KMksOBP9u+MZefYrspH+t+23/P5/MMcH2FY52Z+/IMcHVZv6uy/S/bf7H9b9vvAhdVaLtYfpHt4fl7eZ/0veyQg+oOEyM2IdRR/579GTN4TL27EZYyDWn2D3KAA23fpzR/5gDggTyyMSnvn7hYe9JOwPmk0ZqVgJWBm/Lua0ijNddL6g6MIAUNcyV9BTgJ+F2+DXKi7RfL+tIbeKXieaS/4C8k/QW+Guk3ZFyVc+oFTCh9sL1I0kRgowpl+wC3SlpU2LYQ2KDw+cNrkK/TOcCXgfX4aBSnB+larFLlHPoAX5a0f2HbisBfc7v7AD8GNif94b8a8FyV87sfWDV/F5NIo0W35nZWIwUBnwfWzuW7SVre9sLy86lBc9/Jp0mZ2bfho38L15UVKx5rAlDT02iSepJGz/4fKUnqcqQRtWrll899OYj0XRS/l7m1HLM9xIhNCCF0Enl04RbSj/quxV1lRUeSEk72tr0WKbO2chvzbZ9pe2vS6Mh+pNtW2L7H9mdIk5RfJGUELzcR2LRKF3+V621me03gVKrfRnmLFEgAaZIs6Qf6zSrH3Md298JrFdvFssVrMIgUAO5Nuu3Wt3QYYArwfpVzmAhcU3ac1W2fL2ll4GbSyNIGtrsDd1c7v5xR/UbS6Mcg4M48OgMpSekWwE75Ou1e6F+l82lJc9/J9bnfpX8LV1boc/Epu01I3025Sv35Kel25nb5PAbT/DkcSRr9+m/S99Ivb28xsm9PEdiEEEInoeQA0l/5LzRTtBswzfb7kj5F+mEttbGnpO3yX8+zSLemFkraIE80XZ30YzWHFECVuxPoKen4PEm2Wx6VKB13FjBH0pbAN5rp443AFyTtJWlF0o/9B8CjFcr+GjhHUp98Duvl69Dc+X8ATCWNqpxb2pEDjt8DF0rqJWl5SZ/OgcsIYH9Jn8vbV5E0ME/sLY12TAYW5NGblh6BH0m65XZYfl/s3zxghtIk6B+30E5LRgCfl/S/klaQ1ENSKat78d/CzsAhFeqfoTShejvSHKwbKpR5F7Ck4kTmbqSRlplKSxCcVFbnHRaf+Fz+vZzTutNsHxHYhBBC/d0haQ4paDgHOKrK/JeS44CzJM0GfkQKIkp6kp6smkUKjh4g/TAuRwou3gKmkeZKHFfecB51+Axp4vAk4CVgz7z7JFIQNZs02lPpB7LUzj9Jc0N+QRpF2Z/0WHv5xGGAS0gjUKPzOT1OmvBbzR9It1TeJE2aLZ9sfRLpFtJT+Vx/CixneyJppOdUUgAzETg575tNmsR8IzA9n+ftzfSBPPdoLvnpoMKui0nziabkvv25uXYAJF0p6bIqx3mNdP1+kM/nadJEYUjB5Xn5up3K4v8WSh4mTYgeDZxn+/4Kx5hNuo30RH6iaQApIPsUMJN0LW4uq3YucGYufzxpwvRb+fU8lYPYpU527aNhAwYM8NixY9t0oIED03/HjGlT9RAazsDhAwFijs0yQNI42x36ZEgIkvqRnnTq0FtB9RaTh0Mo6MhHsJsmNdG/Z/+WC4YQQqhZ3IoKoSAewQ4hhK4tRmxCKNNRj2CXbkWFEMLSYPtlOviJpM4gRmxCCCGE0DAisAkhhBBCw4jAJoQQQqckaTdJ/6x3P0LXEoFNCCHUkaTxSlmh5yhl0r5K0hot12x8th+yvUW9+xG6lghsQgih/vbPWZU/CewInF5eIK9K3HD/n51XSA6h3cRTUSHUUdOkpng6KnzI9puSRpGSWyJpDPAIMJAU9Gwn6T7gaNv35TLDSJm8D5fUF3iNlNPnbNKy9hfZPieXHQ68Yfv0/HkgMML2xvnzD0ir765JWj32ONt/Ke9nbqeUj2ln0kq4R9qekPdvSVpxeAfSCr9n2L6xUHceKY/UHsABklYi5WjqTVox+SLbP6/Qv96kVYp3I/1hfp3tb0nalLQS8vak/EX3AN+0PaM11z80hoaL/kPoKgZtNygW6AuLyT/c+wJ/K2w+AhhKysMzoVK9CnYlJWHcC/iRpK1qOPYWwLeAHW13Az4HjG+mymGk4KkH0ARcm9tZHbiXlDtpfVKSyF9K2qZQdxApdUQ30nL/vwOOzcfdlpQ5u7x/y5PyWE0gJb3ciJQAEtIjzeeRUhtsRQqQhrV0zqExxYhNCHUydIehDN1haL27ETqAhrS4lMhtkhaQcvLcRSGpIzC8mDcqJclu0Zm25wHPSHqGNJLRXFJNSAkxVwa2ljTZ9vgWyt9l+8Hcp9P4KFHiLsB421flck9Luhk4iJQ/COBPth/J79+XND8f9xnb00m5msp9ihS4nGx7Qd72MHy4XsvLedtkSRey5IknQxcVIzYhhFB/B9rubruP7eNyUFIysQ3tTSq8fw9ocTJyDg6OJ410vCvpekm9mqnyYb9szyElZ+xFusW0U06MOEPSDNLoTs9KdbP/JY1UTZD0gKRPVzheb2BCIaj5kKT1c3/flDSLlPSzRwunHBpUBDYhhNC5lWcqnkuaO1PSk9o1W9f2SNu7koITk7JiV9O79CY/xbUOaV7OROCBHKiVXmvY/kbxUGXHfcr2AaRbV7dROUP1RGATSZXuNJyX2/yE7TVJWcWXuRV3QxKBTQghdC1NwCGSVpQ0gHSLpzV195W0jqSepBEaIM2xkfTfklYmTQyeR7o9Vc2+knbNE3/PBp6wPZE0D2ZzSUfkPq4oacdq83wkrSTpMElr2Z5Pmjxc6bhPAm8D50taXdIqkv5f3tcNmAPMkLQRcHIrrkloMBHYhBBC13IG6Wmk6cCZpEm6tboGeIY0KXg0cENh38rA+cAU0q2s9YFTm2lrJGkeyzTS00+HAdieDXwWOIQ0gjOJNPKzcjNtHQGMz7eRvk4acVmM7YXA/kA/4HXgDeArefeZpKfGSnOUbmnmWKHByS4f5axuwIABHjt2bJsONHBg+u+YMW2qHkKHKD163RFJMMOyQ9I42wPq3Y/2Uv7YeAidSYzYhBBCCKFhRGATQgghhIYR69iEEEJoFduD692HEKqJEZsQQgghNIwIbEIIIYTQMCKwCSGETkDSGEnT8zoyIYQ2isAmhBDqLGfl3o20eu4X69qZELq4CGxCCKH+jgQeB4YDRxV3SBou6XJJd0maLekJSZsW9l8iaaKkWZLGSdqtsG+YpBsl/SHXfT6vVlza/0NJr+R9/5D0P4V9/XLeppmSpkgqLua3mLwC8aM5N9RESYPz9rXysSdLmiDpdEnL5X2DJT0i6aJc71VJu+TtEyW9K+mowjGaa6tiXyX1leRiGoY8MnZ0a88xdB3xVFQIIRRcMe4KRj7XmsV828WRwIXAE8Djkjaw/U5h/6HA54GngauBc0gr+wI8BZxFWnX3u8BNkvrafj/v/yLwJWAI8BPgMmDnvO8V0kjRJODLwAhJ/Wy/TUqTMBrYE1gJqLjAoKRNgFHAUOCPwJp8lEfqF8BawMeBdXN7bwO/y/t3Aq7M+84ErgfuIK0uvAdws6Sbc5LN5tqqqa8VtLVe6MRixCaEEApGPjeSpklNHXY8SaWkkzfaHkcKNgaVFbvF9pM5s/W1QP/SDtsjbE+1vcD2BaTUBVsU6j5s++6ckuAaYPtC3Ztsv2V7ke0bgJeAT+Xd83O/etl+3/bDVU7hMOA+29fZnp/70iRpeVLKg1Nsz7Y9HriAlD6h5DXbV+W+3UAKiM6y/YHt0cC/gX41tFVrX8u1tV7oxGLEJoQQyvTv2b9d02poSLOJpo8CRtuekj+PzNsuKpSZVHj/HrDGh21LJwJHA71Ic3TWBHo0U3cVSSvYXiDpSOAEoG/ev0ah7vdJIxpPSpoOXGD79xX635sUjJXrQRoFmVDYNgHYqPC5OCo1D6BspGpeoU/NtVVrX8u1tV7oxCKwCSGEOpG0KnAwsLykUgCyMtBd0va2n2mh/m7AD4C9gOdtL8o/0M1GUrluH+C3ue5jthdKairVtT0JOCaX3RW4T9KDtl8ua2oiH43yFE3hoxGRf+RtmwBvttS31rZVra+k23MAq5GyhgP0LDXainMMXUjcigohhPo5EFgIbE26vdQf2Ap4iDTvpiXdgAXAZGAFST8ijdjUYnXSCM9kAElDgG1LOyV9WdLG+eP0XHZhhXauBfaWdLCkFSStK6l/vr10I3COpG45kDoBGFFj/z7UUlvV+mp7Min4OVzS8pK+SsqM3tpzDF1IBDYhhFA/RwFX2X7d9qTSizTB97Di0zxV3EOauPsv0q2Z90kjKC2y/Q/SPJXHSLeEtgMeKRTZEXhC0hzgduC7tl+r0M7rwL7AicA0oImP5vF8G5gLvAo8TLrN1tZbPc211VxfjwFOBqYC2wCPtvYcQ9ci2zUXHjBggMeOHdumAw0cmP47ZkybqofQIQYOHwjQrvMrQteyNP4NSBpnO564CaEDxIhNCCGEEBpGTB4OITScJVmLpmlSE/179m+5YAihU4oRmxBCw+notWhCCJ1HjNiEEBpSW9eiKc2xCSF0TTFiE0IIIYSGEYFNCCGEDiXpVElXFj7/T058OUfSf9Wzb6Hri8AmhBDqTNIgSWPzD/vbkkbllXAbku1zbR9d2PRz4Fu217D9t3r1KzSGCGxCCKGOJJ0AXAycC2xAShXwS+CAevarg/UBnm9LxRoWMQzLmPgHEUIIZZomNXXIJGJJawFnAUNs31LYdUd+IWll4KeknFKQUgv8wPYHef8BwJnAx0npEb5p+8+SegG/BnYlrQj8U9u/zXU+RQqeNiclmrzW9gkV+jeQlLbgIlJOqoXAqbavyvvHACNsX5k/DwaOtr1r/rwNKWjbgZTr6RLb50oaBvQDvkZaEXh54Bl4AYBUAAAgAElEQVRJk2xvmvv+C2B3YA5wke1Lc5vDSKkf3ge+CJwg6VngElI6innAzcAJtv9dy/cQGkuM2IQQQsGg7QZ15Do2nwZWAW5tpsxpwM6kPFLbkxJOng4fBih/IKUM6E4KBMbnetcBb5Cyfh8EnCtpr7zvElKQsSYpd9KNzRy/J7AWKZP214DLJa3d0olJ6gbcB/w596Ef8JdiGdsf2C5lKt8+BzXLkYK6Z/Ix9wKOl/S5QtUDgD/mc76WFHB9j5QF/NO5znEt9TE0phixCSGEgqE7DGXoDkPbtU0NqZpse11giu0FzVQ/DPi27XcBJJ0J/AY4gxRo/N72vbnsm7lMb9JIzX623wea8mTdI0jBxXygn6QetqcAjzdz/PnAWbmPd+e8Slu0UAdgP2CS7Qvy5/eBJ1qoAyl/03q2z8qfX5X0W+AQUm4sSNnIb8vv5wHjCvXHS/oNsAdptCgsY2LEJoQQ6mcq0KOFeSK9SAkuSybkbQC9gVeq1Jlme3ZZvY3y+6+RbkO9KOkpSfs118eywOs9YI1qhQuq9a0lfYBekmaUXsCppPlHJYsl+pS0uaQ7JU2SNIs0X6lHG44dGkCM2IQQurRK6RO6UFqEx0gjGQeSbq1U8haLT67dJG+D9AO/aZU660jqVghuNiGP6Nh+CTg03/b5EvBHSevantvK/s8FVit87ll4PxE4tJXtleq9ZnuzZsqUZ2/+FfA34FDbsyUdT7r9FpZBMWITQnbFuCtiGf4uqCunT7A9E/gRad7KgZJWk7SipH0k/SwXuw44XdJ6knrk8iPyvt8BQyTtJWk5SRtJ2tL2ROBR4DxJq0j6BGmU5loASYdLWs/2ImBGbmthG06hCfhS7ndpMnDJnUBPScdLWllSN0k71dDmk8AsST+QtKqk5SVtK2nHZup0A2YBcyRtCXyjDecSGkQENiFkI58bycwPZta7G6ENSukTSq8uMloDgO0LgRNIE4Ink0YsvgWU5pD8BBgLPAs8Bzydt2H7SWAI6amlmcADpNEdSKMlfUmjN7cCPy7Mxfk88HyeL3MJcEiei9NaFwH/Bt4BriYHTrlvs4HPAPsDk4CXgD1batD2wlynP/AaMAW4kjSBuZqTgEHAbOC3wA2tP5XQKOJWVAgFa63c3P93hrB02L6WQlBQtu994Dv5VWn/rVR4qsr2G6QJvJXqHF5jv8YAG5dt61t4PwX4bFm1YYX9fyc9oVTe7rCyzyr7/BZVbmOV183bHgS2LNv8o0r1Q+OLEZsQQgghNIwIbEIIIYTQMCKwCSGEEELDiMAmhBBCCA0jApsQQgidiqTDJI2udz9C1xSBTQgh1JGk8ZL2bsf2Bkp6o73aW9ok9ZXk4urLtq+1Xf60VQg1icAmhBBCCA0jApsQQuiEJK2d8x9NljQ9v9+4sH8dSVdJeivvv03S6sAoUq6lOfnVS9JwST8p1F1sVEdSb0m35GNNlXRZ3r6ppPvztimSrpXUvUp/JekiSe9KminpWUnb5n1fkPQ3SbMkTZQ0rFD1wfzfGbm/n5Y0WNLDue6vJf287Fh/knRCfv9DSa9Imi3pH5L+Zwkue2gAsUBfCKEhNU1qYuDwgfXuxpJYDrgKOBhYHvg9cBkprxTANcAcYJv8311sz5W0DzDCdjEIqnoQScuT0h/cT8r+vRAYUNoNnEcKPtYEbiYtwHd8haY+C+xOSq45k7RgXildw1zgSFK+q22BeyU15Qzdu5NWGO5eSrYpaYtCuyOBayWdbNuS1s7HKqVNeAXYjbS68ZeBEZL62X676kmHhhaBTQih4QzablC9u7DEbE8lBRIASDoH+Gt+vyGwD7Cu7em5yANtPNSnSNnATy5k8X449+Fl4OW8bbKkC4EfV2lnPiln05bAk7ZfKJzLmEK5ZyVdB+zBR2kjmvMQKenlbqQA6yDgsbw6MbZvKpS9QdIp+Zz+VEPboQFFYBNC6JQqZe2upFIm76E7DGXoDkOXVtdaTUOqj5hUrSOtRsrF9Hlg7by5Wx5h6Q1MKwQ1S6I3MKEQ1BT7sD5wKSmo6EYaRap4TNv351tYlwObSLoVOMn2rJz88nzSaM1KwMrATZXaqdCuJV1PSrHwICknVCkJKJKOJOXa6ps3rQH0qKXt0Jhijk0IoVPqylm728mJwBbATrbXJN2ygXR7aCKwTpX5Lq6wbS6wWuFzz8L7iaRApNIfuufl9j6R+3B4Pn5Fti+1vQPp9tjmwMl510jgdqC37bWAXxfaqdTfctcBB0nqA+xEHsnKn39LShq6ru3uwN+b62NofDFiE0LotEpZu5vTxefRlKwoaZXC5wWkEZJ5pEm161C4BWT7bUmjgF9K+iZpjs2nczLId4B1Ja1lu5Suvgk4MU8gXonF58g8CbwNnC/px6Q5NjvYfiT3YWbuw0Z8FKj8B0k7kv5YfpoUSL2f2yK3M832+5I+RRp1Ka1TMxlYBHwc+Feltm3/TdJkUpbve2yX5u6sTgqMJuc+DCGNCoVlWIzYhBBC/d1NCmJKr2HAxcCqwBTgceDPZXWOIM1reRF4lxys2H6RNMLxqqQZknqRJho/A4wnBRQ3lBqxvRDYH+gHvA68AXwl7z4T+CQpuLkLuKWZc1iTNHoyHZgATAVKTzMdB5wlaTYp6/aNheO/B5wDPJL7u3OV9q8D9iaN/pTq/gO4AHiMFNBtBzzSTB/DMiBGbEIIoY5s921m98Cyz78p1JsGHFWlza9W2PyVss8XFcq/zkdPWxXbeR7YoWzzBVWO+RfgE1X2/RH4Y6V9ef+PSAFPyePA8LIyZwNnV6h7GnBatbbDsidGbEIIIYTQMCKwCSGEEELDiMAmhBBCCA0jApsQQgghNIwIbEIIIYTQMCKwCSGETkzSMEkjWi7ZeRXPQdImOdnl8vXuV3Mk7Sbpn/XuR2i9CGxCCKHOJA2SNDb/4L8taZSkXevdr6XB9uu218jr53Rath+yvUXLJUNnE4FNCCHUkaQTSIvxnQtsAGwC/BI4oJ79CqGrigX6QghdXtOkpi6ZWkHSWsBZwBDbxVV978ivkpUk/QH4H9LqwEfZHpvb6AX8gpRLag5wke1L875hwNak9AaV6m4F/AroD7wJnGL79rxveK63KbAzKVXCkbYn5P27AJeQckL9C/iu7Ufzvo+RFtj7JGmxvQ9v6UjqC7wGrGh7gaTxwNG27yv0uZ/twwtlv5qv0xrAKcA44HekIHCE7W9Vub6r5vM7gJQ24qrcz43zfgOb5SzmpXN+w/bpkgbmtktlxwOXAUcCfUgrQR9l+/28/xjgB8A6pOzoXy9lIJd0CfAlYC3gJeB42w9V6nNYcjFiE0Lo0gZtN+g/snt3IZ8GVgFubaHcF4Hrge6kZJKXAUhajhQAPQNsBOwFHC/pczXUXTHXHQ2sD3wbuFZS8fbLYaTVfnuQ8k1dm+uuQ0qxcCmwLnAhcJekdXO9kaTgo0euX3GF5FbYCdiMtHryxaSVhvcmJds8WNIeVer9mBSYbQp8rh36cTAp2/rHSKssDwaQ9N+khKEHAxuSUkpcX6j3FCl4XId0bW4qyw0W2lGM2IQQurShOwxl6A5D692NZmlI1WTT6wJTbC9ooYmHbd8NIOkaPkpiuSOwnu2z8udXJf0WOAS4p4W6O5NGQM63vQi4X9KdwKGkXFUAd+XEmkg6DZgpqTcp1cNLtq/J5a6T9B1gf0n3537tbfsD4EFJxdGntjg7j4yMljQXuM72u7lfDwH/BTxQod7BwHE5/cQ0SZeyeOqG1rq0MApzBylYgRQA/t7203nfKcB0SX1tj7ddnPx9gaTTSZnbn1mCvoQqYsQmhBDqZyrQQ1JLf2ROKrx/D1gl1+kD9MrJI2dImgGcSpqr01LdXsDEHNSUTCCN/JRMLL2xPQeYluv1ymWLSnV7AdNtzy3btyTeKbyfV+HzGlXq9aJwDu3Qj/JrWTruYtcjX6up5Gsp6URJL0iamb+jtUijWWEpiMAmhBDq5zHSPJb/SEBZo4nAa7a7F17dbO9bQ923gN75dlbJJqS5NiW9S28krUG6lfJWfvUpa69U921gbUmrl+2rZi6wWuFzzxr6Xqu3KZxDhX68107HXux65HNfF3hT0m6kuTcHA2vb7k7Kll51GC8smQhsQgihTmzPJN0auVzSgZJWk7SipH0k/ayGJp4EZkn6gaRVJS0vaVtJO9ZQ9wlSUPH9fMyBwP4sPjdkX0m7SlqJNFfmCdsTgbuBzfNj6itI+gppkvKdeXLxWOBMSSvlx9b3b6YfTcAhuQ8DgINq6HutbgROkbS2pI1J84jKjz0oX7fPA9Xm6rRkJDBEUn9JK5OecHvC9nigG7AAmAysIOlHwJptPE6oQQQ2IYRQR7YvBE4ATif9+E0EvgXcVkPdhaSgoT/p6aEpwJWkWx0t1f03aWLxPrneL0lPPb1YKDaSNAF3GrADaS4JtqcC+wEnkm65fB/Yz/aUXG8QacLvtFz/D8105QzS5N7pwJn5mO3lTNItotdIk6SvKdv/XdL1m0E6txaveSW2/0I6j5tJo0SbkuY5QZrrNIr05NgE0gjdxArNhHYi2zUXHjBggMeOHdumAw0cmP47Zkybqoew1A0cPpCmSU3079mfMYPH1Ls7y7zS49uN8F1IGmd7QL370RrFR5/r3Zf2Uv4Id2hMMWITQgghhIYRgU0IIYQQGkasYxNCCOE/2B5c7z60N9tjgLgN1eBixCaEEEIIDSMCmxBC6IIk9ZXkGhb3qytJgyU9XO9+tIecdX1J0zKEpSwCmxBCqKO8TsyjeVXaaZIeqXEdmtDBbO9j++p69yM0r1NH+iGE0MgkrQncCXyDtJjcSsBuwAf17NeyStIKNeTtCp1cjNiEEEL9bA5g+zrbC23Psz3a9rOQsndLOl3SBEnvSvqDpIqL70kakvMRzZb0qqRjC/sGSnoj5yx6V9LbkoYU9o+RdHTh82K3j/Itr+MkvZTbP1vSppIekzRL0o15deJmSbpc0gVl2+6QdHx+31vSLZImS5oq6cMs5tWuQ+GW3FBJb+VzO7HQ/nKSfijpldzmjUrZyYt1vybpdVIi0FUkjchlZ0h6StIGVa7TMYVr/g9Jnyxcr36FcsMl/aTsu/h+4bs4UNK+kv6VR+1ObelahupixCaEENrRFeOuYORzNS+e+y9goaSrSakMHrc9vbB/cH7tCbxLWsH3MuCICm29S1oN+FVgd2CUpKdKGadJeZDWIiVm/AzwR0m3lR2vOZ8nrT7cG3ga2IW0Wu9UUs6rQ4GWbtNcDdwm6WTbiyT1APYCjpa0PGn06v58fguB0qKGtVyHPYHNgI+TApRnbN8HfIeUi2sP0srOlwKX5/6W7AFsBSwCjiJdp96kkbP+pESbi5H0ZVIW9ANJKSQ2Bea3cP4lPYFVSN/FYOC3wL2k67sJME7S9bZfrbG9UBAjNiGE0I5GPjeSpklNNZW1PQvYFTDpx22ypNtLIwSkwOFC26/mjNGnkPIq/ccfpbbvsv2KkwdIKQR2KxSZD5xle77tu4E5wBatOLWf2p5l+3ng78Do3K+ZpJQB/1XD+T5JSgC5V950CDDG9jvAp0hZsk+2Pdf2+7ZLo0a1XIczc73ngKv4KHA5FjjN9hu2PyAFIweV1R2W687L12ldoF8eRRuXv6dyRwM/s/1UvuYv5zxZtZgPnGN7Pimg7QFcYnt2vr7PA5+osa1QJkZsQijTNKnpw+X8Q/2U0lt0ReVpOTSkeiJn2y+Q/mpH0pbACOBi0g9zL1J+oZIJpP/f3oAykvYh5WXanPRH62rAc4UiU8vmj7wHrFHzScE7hffzKnyuNTP21cDhpBGKw4FL8vbewIQqc1xquQ4Ty/Zvl9/3AW6VtKiwf2Ezda/JfbleUnfS93FaDkKKegOvVDrBGkzNeb7go9Gg8uvZmu8mFMSITQgFG6y+QZf9MQ1dX05AORzYNm96i/TDXLIJKVN08UcQpYzSNwM/Bzaw3Z2Ugbt6RLW4uaRAqKTWIKUtRgAHSNqedPunlHhyIrBJlcfXa7kOvcv2v1Vodx/b3QuvVWy/WSj/YdLEPKJ1pu2tSbfb9gOOrNCniaTbT5W8R8ddz1AmRmxCKNiw24YNkXSxESwLo2Z5hOYLwA2235DUmzRS83guch3wA0mjSPNDzs1lF0iLxSwrASvnMgvy6M1nSbeMatEEfEnSlaTRka9RFjy1l3yeT5FGRm7Ot38AniRlxj5f0o9Joyo72H6E2q7DGZKOAT4GDCGNBgH8GjhH0lG2J0haD9jF9p8q9U/SnqRs5/8AZpFuGy2sUPRK4MI8yfpp8hybfDuqCRgk6XnSfKY9SPNwQgeIEZsQQqif2cBOwBOS5pICmr8Dpad6fk8KAB4EXgPeB75d3ojt2aRJsjcC04FBwO2t6MdFwL9JwczVwLVtOJfWuJp0q+ia0oZ8a2Z/oB/wOvAG8JW8u5br8ADwMvAX4Oe2R+ftl5CuxWhJs0nXeKdm+tYT+CMpqHkhtzuivJDtm4BzgJGk7/E2YJ28+7v5XGaQ5gfdVl4/LD2y3XKpbMCAAR47tm1B58CB6b9jxrSpeghLXWmEIEZsOoeu+n1U6rekcbYHVK6x7JG0OylY6Gt7UUvlW2irLynYWTHWoAkQIzYhhBA6kKQVSSMaVy5pUBNCJTHHJoQQmtHKdWm69NNcS5ukrUhzTZ4hzYMJod3FiE0IITSjNevShObZfsH26rZ3qbI2TFvaHG9bcRsqlHToiE1T00dzbULobJomXUz/Hx5f726ETqh8XZrmLAtPc4XQmXVYYDNoUEcdKYQQQgjLqg4LbIYOTa8QOquBw2O0JnQdjfw0kCQDm9l+ud59CV1PzLEJIYQ6krSrpEclzcyZnR+RtGO9+9VeyjOFh7C0xVNRIYRQJ5LWJGW0/gZpcb2VSIkrP6hnvxqFpBUabTQrtCxGbEIIoZ2VEqmWXs3YHMD2dTmT9Dzbo20/CyBpOUmnS5og6V1Jf5C0VqWGJA2R9IKk2ZJelXRsYd9ASW9IOjG387akIYX9YyQdXfi82CiLJEs6TtJLuf2zJW0q6TFJsyTdKGmlCn3aipTS4NOS5kiaUcvxsn3zeUyR9H+SliuU/2o+1+mS7pHUp7DPkr4p6SXgpbxtF0lP5VGxpyTtUnbsV/N5vSbpsLx9mKQRhXJ9c9srFM7hJ3m0bY6kOyStK+nafE2eyrcLQweLwCaEENrRoO0GtWYdm38BCyVdLWkfSWuX7R+cX3sCHydlfL6sSlvvkhI2rklaI+YiSZ8s7O8JrAVsRMoFdXmF4zXn88AOwM7A94ErSOkCepOSdh5aXiFnLv868JjtNXJyzlr9DzAA+CRwAPBVAEkHAqcCXwLWAx4i5ZIqOpCUNmFrSesAdwGXAusCFwJ35SBk9bx9H9vdSEkvW/Ns/yHAEaRruinwGHAVKbXCC6Rs66GDRWATQgjtaOgOQxkzeMxir2ryWi67krJL/xaYLOl2SRvkIocBF9p+1fYc4BTgkEoZsG3fZfsVJw8Ao0m3tUrmA2fl7NV3A3OALVpxaj+1Pcv286R8VqNzv2YCo4D/akVbtR5vmu3XgYv5KHA6Fjgvr4mzgJQQs39x1Cbvn5YTbH4BeMn2NbYX2L4OeJGUywlgEbCtpFVtv53Pr1ZX5Wteugav2L4v9+sm2v+ahBpEYBNCCHWUf6AH296YNPLRi/RDTn4/oVB8Amlu5AaUySM+j+cJyDOAfYEehSJTy+abvEcaAapVMdv3vAqfW9NWLSYW3k8gXQuAPsAlkmbk85wGiDRqUqlu+TUstbeR7bmkRJtfB96WdJdSxvVadfQ1CTWIycMhhIbW2pQI5ToyRYLtFyUNJ41KALxF+iEv2QRYQPoB3bi0UdLKwM3AkcCfbM+XdBvpB78Wc4HVCp97tukEKquUabmW4/UGSqMnm5CuBaSg5RzbzWUgLx6z/BqW2vszgO17gHskrQr8hDRytluNfQydUIzYhBAaWmdOiSBpyzyhd+P8uTfplsvjuch1wPckfUzSGqTbLjdUeNJnJWBlYDKwQNI+wGdb0ZUm4EuSVpPUjzQHp728A2xcNrm4luOdLGntfE2+C9yQt/8aOEXSNgCS1pL05WaOfzewuaRBklaQ9BVga+BOSRtI+mKea/MB6fbcwkIfd5e0SZ6wfUqbzj50uBixCSE0vNakRCi3lFMkzCZNcj1BUndgBunx75Pz/t+TbqU8CKwC3AN8u7wR27MlfYf0yPjKwB3A7a3ox0XAjqQg5FngWmDvNpxPJfeTRl4mSVpku0eNx/sTMI404Xk48DsA27fmIO/6PK9mJnAvaU7Lf7A9VdJ+wCXAr4CXgf1sT5G0IXAicA1plKcJOC7Xu1fSDbl/U4CfAl9c4qsRljrZlUYJKxswYIDHjh27FLsTQv2UfsDa+gMY2ld7fR9L2k579EPSONsD2txACKFmcSsqhBBCCA0jApsQQgghNIwIbEIIIYTQMCKwCSGEEELDiMAmhBA6ufK8RaHjlfJt1eG4zo/E111+9H2OpOXboa1Rko5qj36Vi8AmhBDqLCdifE7Se5ImSfpVfvw7hE7D9us559fClku32NY+tq9uj36Vi8AmhBDqSNKJpDVSTiat2bIzaaXceytlzO6sKuWvCl1Tpe+yK32/EdiEEEKdSFoTOBP4tu0/5wSV44GDScHN4VXq7Szp0Zwv6RlJAwv7hkh6QdJsSa9KOrZSG4XyX83lp0u6p5RMUlLffBtkhULZMZKOzu8HS3pE0kWSpgHDJC0n6XRJEyS9K+kPedXeYntDJb0l6e0c1JXa/pSkx/I5vS3psmJgl1dpvjfnwvqnpIML+74g6W+SZkmaKGlYM+fbQ9Kd+TjTJD0kabm8b7HbPpKGS/pJWf1TJU2RNF7SYYXt+0r6R77ub0o6KW8fKOmNavVao7nzrOH6Lifph5JekTRV0o1Kmc+Ldb8m6XXg/ha2rZDrjZF0dv53MFvSaEk9Csds7t/ph/+W2luXicBCCKFemiY1La0ViHchrSh8S3Gj7TmSRgGfIa0+/CFJGwF3AUeQ8h3tBdwsaUvbk4F3gf2AV4HdgVGSnrL9dPnBJR0InErKdP0S8ENSGoddauz/TsD1wPrAisDg/Noz9+MPwGW5ryV7ApsBHyf9WD5j+z5SKoPvAWNJebBGkVYBvlgp5cG9wI+AfYBPAKMlPZ+zcc8l5cl6npRI9F5JTbZvq9DnE4E3gPXy552pnM+qkp6kxKIb5Xp3Sxpr+5+klZEPtv2QpLWBj9VYrzVqOc9q1/c7wIHAHqTUG5cCl/NR1nTyvq1IGc83aGZb0SDSdzKR9J2dBPywhn+nS02M2IQQQjMGbTdoaSbB7AFMqZD7CeBtFs/OXXI4cLftu20vsn0vKRjYF8D2XbZfcfIAMJqU1LGSY4HzcobxBaRcVP1LozY1eMv2L2wvsD0POAy40PartueQ8isdosVvY5xpe67t54CryD+stsfZfjy3NR74DelHFVKgNt72VXn/06SknwflumNsP5evx7Ok4GwPKpsPbAj0ySNkD7k1S/DDGbY/yNf2LtLoWqndrSWtaXt6hUCyWr2a1XieFa8v6bs+zfYbtj8AhgEHlX03w3LdeS1sK7rK9r/y/huB0v9Ymv13ujTFiE0IodNqj5GSJc3OPXSHoQzdYegS9UFDqibZngL0kLRCheBmw7y/XB/gy5L2L2xbEfgrgFICzB8Dm5P+eF0NeK7K8fsAl0i6oNhd0sjCW5WrLGZi2edewITC5wmk35niX/oTy/Zvl/u9OXAhMCD3eQVSrqhSP3eSNKNQdwVSjick7QT/v70zj5Krqtb47wMhDIGEB8iUIKMDiKJEhMhTnygIMj2fA4hAUCbx+RQQBXEAVIyKgoqIggMS5hkHBBWCqEwJIFOYISQkgTCEDIAE+N4f+xRcKlXd1Z2urqZ7/9a6a1Xdc+65e597u++ufc49H+OJLEZNELShdhTwfeKhfrkkgF/YHt+1my/xpO0FdfavWT7/D/BVYLykW4DDbF/TwnEt06KfDfuX6MMLJb1YKX+B5temq31VZlU+Pw0Mr5yv6X3aTjJjkyTJgKTNmZKBwjWEqvSHqzvL0Mt2wF8bHDMNOM32yMq2vO3xkoYRmYxjgdVsjyTUrZtFVtOA/evaWtb2P4lhD4ggo8bqdcfXZzpmEA+0GmsDzxNilzVG15XXAqifAXcCG9pekRgiq9k9Dbiqzs7htj9Tys8gRD9H2x5BKIA39Nn2PNuH2F6PGII7WNLWpfjpbvxdqVybRey3fYPtnYlhuYuI7EW3x/WQVvxs1r/TgO3q+nAZ2w9X6jfKXPUkm1Wl6X3ay/ZaJjM2SZIMSPoiUwJtV+deLGw/Jeko4CeS5hKBzFrAicQ8kNMaHDYBuEHStsBfiF/BWxCq1U8Rv+JnA8+X7M02wG1NTDgJ+GaZp3G7YqLvNrbPtT1b0sPAJyX9HNgLWL8bl84EvlzmB80mhrbOtv18yY4AfE3SvsQclL15eYL0CsBcYL6kNwKfKW1AKJ6Pl7QHMacHYshjvu0p5dgnbD8raXNi3sfljQxUKH3fCdxXzvdC2SDUvT8h6XZiftN7iOGTKkdJ+goxv2gH4BuKSc4fBX5fruncSptNj2vWiV3Qip/N+vck4NuS9rI9VdKqwFjbF/fCjlZoep/abut6QJmxSZIk6SC2v0dkJ44lHrTXEb92ty5zIerrTwN2LsfMLnUPBZawPY+YJHoO8CTx4Luki3NfSLxqflZ5GN9GZIpq7FvafhzYGPhnN+78igjG/gY8ADwLfK6uzlVEEPZX4FjbtQfzF4u984CTgbMrds4jArRdiQzErGL3sFLlQOBoSfOICcbVbEk9GxIP2vlExuxE2xNL2eeJLM4cYr5Q/eTjWUS/zgBOBw6wfWcp2wN4sPTjAbzyjbamx5W3pS7twt4qrfjZrOENQPcAABFFSURBVH9/RNwLl5fjryWCrLbQ1X3arnPWUE/mTI0ZM8aTJtUHr0kyOKj9sp84bmJH7Uj6loFwXSVNtj2mYwYMACStQwQ7SzWZLD0oKa84T7A9qs3nWYch2L+NyIxNkiRJkiSDhgxskiRJkiQZNOTk4SRJkqTtlLVpmr73Plgp83faOgxVzvMgQ7B/G5EZmyRJkiRJBg0Z2CRJkiRJMmjIwCZJkuRVSL0g4VCgJijZaTterRSRyl4JcL6ayMAmSZKkg0jaqiggP6VQm/6HpHd02q5k8GF7G9und9qOdjNkIv0kSZKBhqQViVV1P0MstrY0IVi5yMJ8SdIVTfTGhiSZsUmSJOkcrwewfabtF2w/Y/vyotyMpCUkfVXSVEmPSvptkT1YBEl7S5oiaZ6k+yXtXyl7r6Tpkg4p7cyUtHelfKKkfSrfx0n6e+W7JR0o6Z7S/jclrS/pGklzJZ1TZAUa2TWuZKGOkzSn2Da27J9W7NmrUn+YpGMlPSTpEUknSVq2SdtvKrbPkXS7pJ0qZb+R9FNJfyg2Xydp/Ur5NpLuKpmyEyVdVdcHnyr9+aSky1RRPC/9cUDpjyfLeVTK1pd0haTHJT0m6XRJI5vYf5KkY+v2XSzp4PL5MEn3FfvvkPTfTfr1CULYs0sk/V3SuPJ5n+Jz9bpsU6n7aUkPVu6nXcv+DSVdWfHvtOo9We6zgyXdWvr2TIWGWe2cEyt1X1P6cp3yfYKkE0p/z5f0N0mrSfpJsXGKpLd252dmbJIkGfT0hUp4m7gbeEHSqYQG0rW2n6yUjyvbfwGPAr8FTiCW76/nUUKD6H7g3cClkm6wfWMpXx0YQWhRfQA4T9JFdefrig8CmxEiizcCYwnZgccJaYLdgFObHPtO4BRgZeCo4uvvgA0IPabzJZ1vez4hlbAeoQW1kBB+/DpweLVBSUuVNn5FyC1sBVwsaYztu0q13YrdNxbbvg3sKmkV4Dyiby8BPkvIR9TUwnchpAB2BO4BDiN0sMZWTNgBeAewIqFC/jvgT8Qr198hZCVWJERJjwS+0KBfzgBOl3SobUtaqfhSE/e8j8jgzSK0qCZI2sD2zEq/nkUIby7VoP3uGFv6ZWVCruGXwGhFJvGHwBjb90haA1ipHCPgW8DVwEjgQuBrhCRGjY8R99hC4t7Yg7j+rfAxYFtgCnAZIf1wBNF/3yakRz7QVQOZsUmSZFAzkFXCbc8lHsgm9JFmS7pE0mqlyu7AD23fXx76hxMP5kV+lNr+g+37HFxFiCP+Z6XKQuBo2wtt/5HQSnpDD8z9ru25tm8nNKUuL3Y9BVwKvK2LYx+w/WvbLxAaUKOLLf8uWkbPARuUrMe+wEG2nygaUccQGlH1bAEMB8bbfs72FcSw3m6VOhfYvr4M0ZxOBEsA2wO3276glP2YCB5q7A98x/aUUn4MsGk1a1POO8f2Q8CVtbZt32v7z8W32USA8J4m/XI1ce1r1+kjwDW2a4rh59qeYftF22cTQdbmleNn2P6J7edtP9PkHF1xn+1fletyKjCqBH0Uu94saRnbM23fUWy62/ZfS58/ChzXwL/jbc+y/ThxTXryB3i+7ZtsP0todc23fUbl3unqPgMyY5MkySCnr1TCFwft3XzdtKJOPQ5AoWo9ATieeECvCUytVJ9K/N9ejToUSt7fIIa3lgCWA26tVHm8bg7G00Rg0CqPVD4/0+D76j04Ftv1+4YDqxJ2T9bLauAClmzQ5prANNsvVvZNJTJSNarBStXfNQlRRoot1ivftnod8CNJP6jsU2m7dj0ati3ptUSg9J+EGvcShADmIpTznkVc678RIqATXjqhtCdwMLBO2TUcWKXSxDQWj3ofAIbbflDSbsAhwK8Vw5IH275b0uqEf+/iZf9m80rq2/2PHtjU3X3W7T2bGZskSZIBQlF8/g3w5rJrBvGQrbE28Dyv/GdPmcNwPpGmX832SOCPtL4S7QIioKjRVZDSTh4jHl4b2x5ZthG2Gz3MZhDDJtXn2NrAwy2cZyaV1YBLpqi6OvA0YP+KDSNtL2u7O3VziGEoA2+xvSKh8t3VdTgT+EjJBr2TuI6U7ycD/wusXK7pbXVtta5i3UNsX2r7/cAahFr4z0vRd4nJ7ZsU/8YxwO6zDGySJEk6hKQ3Kib0jirfRxO/3q8tVc4EDpK0rqThxJDI2Q3eflkaGEb8cn6+ZG+2oXVuBj4saTlJGwCf7r1XvadkX04GjiuZDyStJWnbBtWvIx6UX5K0lEJFe0dizkl3/AHYRNIuZVjvs7zyIXsScLikjYsNIyR9tEU3ViCG+eZIWgs4tKvKtm8irtspwGW255Si5YnAZXaxYW9eDnjbiqQ1JO0oaTlimHAB8EIpXqF8f6rcr19s0kwj/gW8RdImignh3+hLu2tkYJMkSdI55hG/0q+TtIAIaG4jhgAgJsaeRgxTPAA8C3yuvpEyF+X/iFfGnySGNC7pgR3HEQ+wR4i5Fp1c6+TLRIbgWklzgb/QYC6Q7eeAnYDtiEzPicCeJevVJbYfIybjfo+Y/LwRMInymr3tC4nMxFnFhtvKeVrhKODtwFNEAHVBC8ecCbyfmExcs/EO4AfE5NtHgE2Af3TViOLttzld1WmRJYmAbCbRP2OJzBFEMLI54d8llAxTKxSfjgEmAncR93WfI7v1TNaYMWM8adKkdtiRJB2n9tbMxHETO2pHMviQNNn2mE7bkTSmDGdNB3a3fWWn7UkWj8zYJEmSJEMOSdtKGlnmJ32FmCdybTeHJa8CMrBJkiRJhiJbEuvEPEbMzdmll69MJwOMfN07SZIkGXLYPpIWVutNXn1kxiZJkiRJkkFDBjZJkiSDHElfkXRK+bxO0ecZFBl71elcJUkGNkmSJB2kCA0+V1nKvrb/5qpA4OJg+xjbi/3wL68TT+++5qsP1Ql/Jq9eMrBJkiTpPA9Q0TiStAnQUNE6SZKuycAmSZKk85wG7Fn5vheh5P0Skj4k6SZJcyVNk3Rkpaw2vLSXpIckPSbpiEr5kZIm8Ep2b1J3mKTjJc0o2/Fl3/KE2OWakuaXbc16R0rdY0vbj0g6qawyi6RVJP1e0hxJT0i6uiaJIGm0pAskzZb0uKQTGtne3VCapE9JmiLpSUmXVYUry3EHSLqnlP9UwZuI1Ya3LH7N6a7Pk4HLoBhjTZK+4uZZN7+0UF+S9CPXAnuUB+zdwMcJ1e9vVeosIIKf24ml9f8s6WbbF1XqbEWs0vt64HpJFxSRzUY0q3sEoZy9KbGk/8XAV21/rUg1TLA9qkmbECv2rleOX0ispvt1Qpn8EGIhvFVL3S0AS1qSUIG+AtiDWL6/xwsaStqFWJNmR0IJ+zBiVd+xlWo7AO8AVgQmA7+z/SdJBwD72N6qUreVPk8GGJmxSZLCJzb5BJuuvmmnzUiGLrWszQeAO6kTc7Q90fattl+0fQvxwH5PXRtH2X7G9r8IXZ63dnG+ZnV3B462/ajt2YREwB6tOCBJwL7AQbafKFIPxwC7lioLCVHF19leaPtqx/L3mxOK24faXmD7Wdu9me+yP/Ad21OKntYxwKbVrA0w3vYc2w8BVxIBWENa7PNkgJEZmyQp7LfZfuy32X6dNiMZhGjvlsSPa5pQ61I3DAUg6Z3AeCJzUBO9PLeu2qzK56eBRqrY3dVdE5haKZta9rXCqoR68+SIccJ0QnsI4PvE2jGXl/Jf2B4PjAamNhD37CmvA34k6QeVfQLW4mWfWu6jFvs8GWBkxiZJkmQAYHsqMYl4exoLJ55BiA6Otj2CmBPSUsTUQ2YQAUKNtcs+iKGprngMeAbY2PbIso2wPRxCrNP2IbbXI4aLDpa0NTANWLvJvJkFRLBUY/UGdWpMA/avnHuk7WVt/7Mbu5v51l99nvQhGdgkSZIMHD4NvM/2ggZlKwBP2H5W0uaEgnc7OBP4qqRVyyvoXwdqk3cfAVaWNKLRgbZfBE4GjpP0WgBJa0natnzeQdIGZchqLjGX5gXgekJJeryk5SUtI+ldpdmbgXdLWruc9/AubD8JOFzSxuV8IyR9tEW/HwFGSVq6sq+/+jzpQzKwSZIkGSDYvs/2pCbFBwJHS5pHBBvntMmMbwGTgFuAW4Ebyz5s30kEPveXN5saDVF9GbgXuFbSXOAvxCRlgA3L9/nANcCJZR7LC0QGZwPgIWKC8cfLOf8MnF3smUxMMm6I7QuJyctnlXPfBmzXot9XEJOEZ0l6rOzrrz5P+hDFvK3WGDNmjCdNavY3lyRJkjRC0mTbPX7LJ0mSnpMZmyRJkiRJBg0Z2CRJkiRJMmjo0VCUpNm88jXAdrMKMct+qDEU/U6fhwZD0WeAN9heodNGJMlQoEfr2NhetftafYekSUNxXHoo+p0+Dw2Gos8QfnfahiQZKuRQVJIkSZIkg4YMbJIkSZIkGTQM9MDmF502oEMMRb/T56HBUPQZhq7fSdLv9GjycJIkSZIkyUBmoGdskiRJkiRJWqZfAxtJH5R0l6R7JR3WoHyYpLNL+XWS1in7V5Z0paT5kk6oO2YzSbeWY36siqTsQKC3Ppeyw8v+u2paK2X/QZJul3SbpDMlLdM/3rRGm3weKek8SXdKmiJpy/7xpnXa4XcpW1LSTZKaLiXfKfraZ0mjy9/6lHKPf77/vGmNNt3fXbaZJEkPsN0vGyFbfx+wHiH//i9go7o6BwInlc+7AmeXz8sDWwEHACfUHXM9sCWhuHopsF1/+dRmnzcq9YcB65Z2lgTWIhSAly31zgHGddrXdvpcyk4F9imflwZGdtrX/vC7lB9MqAz/vtN+9sP9vQbw9lJnBeDu+jYHoc/dtplbbrm1vvVnxmZz4F7b99t+DjgL2Lmuzs7EAwzgPGBrSbK9wPbfgWerlSWtAaxo+xrbBn4L7NJWL3pGr30u+8+y/W/bDxCicpuXeq8BlpX0GmA5YEab/egJfe6zpBWBdwO/BLD9nO05/eBLT2jLtZY0CvgQcEo/+NBT+txn2zNt3whgex4whQjmBwrtuM6ttJkkSYv0Z2CzFjCt8n06i/7DeqmO7eeBp4CVu2lzejdtdpLF8bnhsbYfBo4lFHBnAk/Zvrwt1veOPveZ+CU7G/h1GZI5RdLy7TG/17TDb4DjgS8BL/a9yYtNu3wGoAzhvA24rg9tXlza4XMrbSZJ0iL9Gdg0mvtS/0pWK3UWp35/szg+N9wvaSXi19y6wJrA8pI+uVhW9i197jORoXo78DPbbwMWAANtHkI7rvUOwKO2Jy+ucW2iHdc6DpKGA+cDX7A9t9cW9j3t8Hmg/x9LklcV/RnYTAdGV76PYtEhlJfqlGGWEcAT3bQ5qps2O8ni+Nzs2PcDD9iebXshcAEwti3W9452+DwdmG679sv9PCLQGUi0w+93ATtJepAYnnifpAntML6XtMNnJC1FBDWn276gLZb3nnbd3921mSRJi/RnYHMDsKGkdSUtTUyqu6SuziXAXuXzR4ArytyZhtieCcyTtEUZw94TuLjvTe81i+PzJcCu5Q2LdYENiYnSDwFbSFqu+Lw1MQ9hoNDnPtueBUyT9IZyzNbAHe12pIe0w+/DbY+yvU5p7wrbAyk71+c+l3v6l8AU2z/sFy96Rjv+pltpM0mSVunPmcrA9sRbDvcBR5R9RwM7lc/LAOcSk+quB9arHPsg8atnPvELZ6OyfwxwW2nzBMqigwNlW0yfjyjH3UXlbS/gKODO4vdpwLBO+9kPPm8KTAJuAS4CVuq0n/3hd6X8vQywt6La4TPx9qPLdb65bNt32s9+uL8XaTO33HLr3ZYrDydJkiRJMmjIlYeTJEmSJBk0ZGCTJEmSJMmgIQObJEmSJEkGDRnYJEmSJEkyaMjAJkmSJEmSQUMGNkmSJEmSDBoysEmSJEmSZNCQgU2SJEmSJIOG/wd0pSeIxiHLoQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gc_array = np.array(list(gcs.values()))\n",
    "gc_array = np.expand_dims(gc_array,-1)\n",
    "ytdist = pdist(aux2)\n",
    "Z = hierarchy.linkage(ytdist, 'single')\n",
    "plt.figure()\n",
    "dn = hierarchy.dendrogram(Z, orientation=\"left\", labels=list(gcs.keys()))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_genes(record):\n",
    "    genes_set = set()\n",
    "    genes_dict = {}\n",
    "    for ft in record.features[1:]:\n",
    "        gene = ft.qualifiers.get('gene')\n",
    "        if gene is None:\n",
    "            continue\n",
    "        else:\n",
    "            gene = gene[0]\n",
    "        if gene not in genes_set:\n",
    "            genes_dict[gene] = ft.extract(record)\n",
    "            genes_set.add(gene)\n",
    "    return genes_dict"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 78,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_codons(seq,frame=0):\n",
    "    codons = []\n",
    "    n = len(seq)\n",
    "    for i in range(frame,n,3):\n",
    "        if i+3 <= n:\n",
    "            codons.append(seq[i:i+3])\n",
    "    return codons\n",
    "\n",
    "def get_frecs(codons):\n",
    "    #cod_set = set(codons)\n",
    "    frecs = {}\n",
    "    for cod in codons:\n",
    "        if cod not in frecs:\n",
    "            frecs[str(cod)] = 1\n",
    "        else:\n",
    "            frecs[str(cod)] += 1\n",
    "    \n",
    "    # En caso de que la secuencia tenga codones ambigüos\n",
    "    for cod in list(frecs.keys()):\n",
    "        if \"N\" in cod:\n",
    "            del frecs[cod]\n",
    "    \n",
    "    # Normalizamos las frecuencias para que se pueda comparar entre especies\n",
    "    norm = sum(frecs.values())\n",
    "    for cod in frecs:\n",
    "        frecs[cod] = frecs[cod]/norm\n",
    "    \n",
    "    return frecs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [],
   "source": [
    "genes = get_genes(mercado[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'rps12': SeqRecord(seq=Seq('ATGCCAACTAGAAAACAACTTATTAGAAACACAAGACAGCCAATCCGGACTGTA...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'trnH-GUG': SeqRecord(seq=Seq('GGCGAACGACGGGAATTGAACCCGCGCATGGTGGATTCACAATCCACTGCCTTG...GCC', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'psbA': SeqRecord(seq=Seq('ATGACTGCAATTTTAGAGAGACGCGATAGCGAAAACCTATGGGGTCGCTTCTGT...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnK-UUU': SeqRecord(seq=Seq('GGGTTGCTAACTCAATGGTAGAGTACTCGGCTTTTAAGTGCGACTCCATTTTTT...CCA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'matK': SeqRecord(seq=Seq('ATGAAGGAATATCAAGTATATTTAGAACGATCTAGATCTCGCCAACAGGACTTC...TGA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rbcL': SeqRecord(seq=Seq('ATGTCACCACAAACAGAAACGAAAGCAAAGGTTGGGTTCAAAGCTGGTGTTAAA...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'atpB': SeqRecord(seq=Seq('ATGACAATAACTCCTCCCCCTTCTGATACTGAGGTTTCTGTACTTGAAAACAAA...TAG', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'atpE': SeqRecord(seq=Seq('ATGACTTTTAATCTTTGTGTACTGACCCCCAATCGAATTGTTTGGGATTCCGAA...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'trnM-CAU': SeqRecord(seq=Seq('ACCTACTTAACTCAGTGGTTAGAGTATTGCTTTCATACGGCAGGAGTCATTGGT...GTA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnV-UAC': SeqRecord(seq=Seq('AGGGCTATAGCTCAGTTAGGTAGAGCACCTCGTTTACACATGCGCCAATGTTTT...CTA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'ndhC': SeqRecord(seq=Seq('ATGTTTTTGCTTTATGAATATGATATTTTTTGGACATTTCTAATAATATCGATC...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'ndhK': SeqRecord(seq=Seq('ATGAATTCCATTGAGTTTCCCTTACTTGATCAAAAAACCAAAAATTCAGTTATT...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'ndhJ': SeqRecord(seq=Seq('ATGCAAGGTAATTTGTCTGCTTCGCTAGTCAAACACGGAATAGTTCATAGATCT...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'trnF-GAA': SeqRecord(seq=Seq('GTCGGGATAGCTCAGCTGGTAGAGCAGAGGACTGAAAATCCTCGTGTCACCAGT...GCG', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnL-UAA': SeqRecord(seq=Seq('GGGGATATGGCGAAATTGGTAGACGCTACGGACTTAATTGGATTGAGCCTTGGT...CCA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnT-UGU': SeqRecord(seq=Seq('GCCCGCTTAGCTCAGAGGTTAGAGCATCGCATTTGTAATGCGATGGTCATCGGT...GCT', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'rps4': SeqRecord(seq=Seq('ATGTCACGTTACAGGGGGCCTCGTTTCAAAAAAATACGGCGTCTGGGGGCTTTA...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'trnS-GGA': SeqRecord(seq=Seq('GGAGAGATGGCCGAGTGGTTCAAGGCGTAGCATTGGAACTGCTATGTAGGCTTT...CCG', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'ycf3': SeqRecord(seq=Seq('ATGACTAGATCTAGGATAAATGGAAATTTTATTGATAAGACCTTTACCATTGTA...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'ndhE': SeqRecord(seq=Seq('ATGATGCTCGAACATGTACTTGTTTTGAGTGCCTATTTATTTTCTATTGGTATC...TAG', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'psaC': SeqRecord(seq=Seq('ATGTCACATTCAGTAAAGATTTATGATACATGTATAGGATGTACTCAATGTGTC...TGA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'ndhD': SeqRecord(seq=Seq('ACGAATTATTTTCCTTGGTTAACAATCGTTGTTGTTTTGCCAATAGCCGCAGGT...TGA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'trnI-CAU': SeqRecord(seq=Seq('GCATCCATGGCTGAATGGTTAAAGCGCCCAACTCATAATTGGAGAAGTCGTAGG...GCA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'ycf2': SeqRecord(seq=Seq('ATGAAGAAGAATGAAGTAGACATTTGGATTTTAGAATTGAGAGAAATTTTTAGA...TAG', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'trnG-GCC': SeqRecord(seq=Seq('GCGGATATAGTTGAATGGTATAATTTCTCTTTGCCAAGGAGAAGACGCGGGTTC...GCC', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'trnfM-CAU': SeqRecord(seq=Seq('CGCGGAGTAGAGCAACTTGGTAGCTCGCAAGGCTCATAACCTTGAGGTTACGGG...GCA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rps14': SeqRecord(seq=Seq('ATGGCAAAGAAAAGTTTGATTGCCAGGGAGAAGAAGAGGAAAAAATTGGAAGAA...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'psaB': SeqRecord(seq=Seq('ATGGCATTAAGATTTCCAAGGTTTAGCCAAGGCTTAGCTCAGGACCCCACTACT...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'psaA': SeqRecord(seq=Seq('ATGATTATTCGTTCGCCGGAACCAAAAGTCCAAATTTTGGCAGATCCAGAAGTC...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'psbM': SeqRecord(seq=Seq('ATGGAAGTAAATATTCTCGCATTTATAGCTACTGCACTCTTCATTCTAGTTCCT...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnD-GUC': SeqRecord(seq=Seq('GGGATTGTAGTTCAATTGGTCAGAGCACCGCCCTGTCAAGGCGGAAGCTGCGGG...CCG', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnY-GUA': SeqRecord(seq=Seq('GGGTCGATGCCCGAGCGGTTAATGGGGACGGACTGTAAATTCGTTGGCAATATG...CCA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnE-UUC': SeqRecord(seq=Seq('GCCCCCATCGTCTAGCGGTTTAGGACATCTCTCTTTCAAGGAGGCAACGGGGAT...GTA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnT-GGU': SeqRecord(seq=Seq('GCCCTTTTAACTCAGTGGTAGAGTAACGCCATGGTAAGGCGTAAGTCGTCGGTT...GCT', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'psbD': SeqRecord(seq=Seq('ATGACTATAGCTCTTGGTAAATTTACCAAAGATCAAAATGATTTATTTGATATT...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'psbC': SeqRecord(seq=Seq('ATGAAAACCTTATATTCCCTGAGGAGGTTCTACCACGTGGAAACGCTCTTTAAT...TGA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'trnS-UGA': SeqRecord(seq=Seq('GGAGAGATGGCTGAGTGGTTGATAGCTCCGGTCTTGAAAACCGGTATAGTTCGA...CCT', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'psbZ': SeqRecord(seq=Seq('ATGACTATTGCGTTTCAATTAGCTGTTTTTGCATTAATTGTTACTTCATCCATT...TGA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'trnL-CAA': SeqRecord(seq=Seq('TGCCTTGGTGGTGAAATGGTAGACACGCGAGACTCAAAATCTCGTGCTAAACAG...GCA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnN-GUU': SeqRecord(seq=Seq('TCCCCAGTAGCTCAGTGGTAGAGCGGTCGGCTGTTAACCGATTGGTCGTAGGTT...GAG', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'trnR-ACG': SeqRecord(seq=Seq('GGGCCTGTAGCTCAGAGGATTAGAGCACGTGGCTACGAACCACGGTGTCGGGGG...CCA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rrn5': SeqRecord(seq=Seq('TATTCTGGTGTCCTTGGCGTAGAGGAACCACACCAATCCATCCCGAACTTGGTG...GAT', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rrn4.5': SeqRecord(seq=Seq('GAAGGTCACGGCGAGACGAGCCGTTTTTCATTAACGATAGGTGTCAAGTGGAAG...AAC', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rrn23': SeqRecord(seq=Seq('TTCAAACGAGGAAAGGCTTATGGTGGATACCTAGGCACCCAGAGACGAGGAAGG...CCT', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnA-UGC': SeqRecord(seq=Seq('GGGGATATAGCTCAGTTGGTAGAGCTCCGCTCTTGCAATTGGGTCGTTGCGATT...CCA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnI-GAU': SeqRecord(seq=Seq('GGGCTATTAGCTCAGTGGTAGAGCGCGCCCCTGATAATTGCGTCGTTGTGCCTG...CCA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rrn16': SeqRecord(seq=Seq('GCTCATGGAGAGTTTGATCCTGGCTCAGGATGAACGCTGGCGGCATGCTTTACA...TTT', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnV-GAC': SeqRecord(seq=Seq('AGGGGTATAACTCAGCGGTAGAGTGTCACCTTGACGTGGTGGAAGTCATCAGTT...CTA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rps7': SeqRecord(seq=Seq('ATGTCACGACGAGGTACTGCAGAAAAAAAAAAAATAGCAAAATCCGATCCAATT...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'ndhB': SeqRecord(seq=Seq('ATGAAAGCCTTTCATTTGCTTATTTTCGATGGAAGTTTTCTTTTTCCGGAATTG...TAG', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'ycf1': SeqRecord(seq=Seq('ATGATTTGCAAGTATTTTGTACGAGATTCTTTTAGAGGTTTCTGCCTTCAAATA...TGA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'rps15': SeqRecord(seq=Seq('ATGATAAAAAATTCTGTTATTTCAAAAAAAGAAGAAGAAAACCCGGGATCCGTT...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'ndhH': SeqRecord(seq=Seq('ATGAATATCCCAGCTACAAGAAAAGACCTCATGATAGTAAATATGGGGCCTCAA...TGA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'ndhA': SeqRecord(seq=Seq('ATGATAATTGATACAACGAAAGTCCAAGATATAAACTCCTTTTCCGGATTGGAA...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'ndhI': SeqRecord(seq=Seq('ATGTTCCTTATGGTAACTGGGTTCATGAATTATGGTCAACAAACACTACGAGCT...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'ndhG': SeqRecord(seq=Seq('ATGGATCTGCCTGAAACGATGCACGATTTTCTTTTAGTCTTTCTGGGATCGGGT...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'psbE': SeqRecord(seq=Seq('ATGTCTGGAAGCACAGGAGAACGTTCTTTTGCTGATATTATTACCAGTATTCGA...TAG', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'psbF': SeqRecord(seq=Seq('ATGACCATAGATAGGACCTATCCCATTTTTACAGTACGATGGTTGGCTGTTCAC...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'psbL': SeqRecord(seq=Seq('ATGACACAATCAAACCCAAATGAACAAAATGTTGAATTGAATCGTACCAGTCTC...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'psbJ': SeqRecord(seq=Seq('ATGGCCAATACTACTGGAAGAATTCCTCTTTGGATAATAGGTACTGTAGCTGGT...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'petA': SeqRecord(seq=Seq('ATGCAAACTAGAAATGCTTTTTCTTGGATAAAGAAAGAGATTACTCGATCTATT...TAG', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'cemA': SeqRecord(seq=Seq('ATGGCAAAAAAGAAAGCATTCATTCCCCTTCTATGTCTTACATCTATAGTCTTT...TGA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'psaI': SeqRecord(seq=Seq('ATGATAAACTTACCATCCCTTTTTGTTCCTTTAGTAGGCCTACTATTTCCGGCA...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'accD': SeqRecord(seq=Seq('ATGATAAATGAATACCCATCTAGTTTAACTGATATGGATAATAACATTGATAGC...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnQ-UUG': SeqRecord(seq=Seq('TGGGGCGTGGCCAAGCGGTAAGGCGCCGGGTTTTGATCCCGTGATTCGGAGGTT...CAG', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'psbK': SeqRecord(seq=Seq('ATGCTTAATATCTTTAGCTTGGTCTGTATTTGTATTAACTCTGCCCTTTATTCA...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'psbI': SeqRecord(seq=Seq('ATGCTTACTCTAAAACTATTTGTTTACACGATAGTGATATTCTTTGTTTCTCTT...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'trnS-GCU': SeqRecord(seq=Seq('GGAGAGATGGCTGAGTGGACTAAAGCGTCGGATTGCTAATCCGTTGTACGCATT...CCG', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnR-UCU': SeqRecord(seq=Seq('GCGTCCATTGTCTAATGGATAGGACAGAGGTCTTCTAAACCTTTGGTATAGGTT...GCA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'atpA': SeqRecord(seq=Seq('ATGGTAACCAGTCGTGCAGATGAAATTAGTCAAATTATCCGTAAACGTATTGAG...TGA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'atpF': SeqRecord(seq=Seq('ATGCAAAATATAACCGACGGAAGTTTCGGGTTTGATACCGATATTTTAGCAACA...TGA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'atpH': SeqRecord(seq=Seq('ATGAATCCACTGATTGCTGCTGCTTCTGTTATCGCTGCTGGGTTGGCCGTCGGG...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'atpI': SeqRecord(seq=Seq('ATGAATGTTCTATTATGTTACATCAACACATTAAACAGATTCTATGATATATCG...TGA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rps2': SeqRecord(seq=Seq('ATGACTAAAAGATATTGGAACATAACTTTTGAAGAGATGATGGAAGCTGGGGTT...TGA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rpoC2': SeqRecord(seq=Seq('ATGGCAGAACGGGCCAATCTAGTCTTTCACAATAACGTGATAGGTGGAACTGCC...TGA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rpoC1': SeqRecord(seq=Seq('ATGATTGATCAGTATAAACATCAACACCTTCGAATTGGATCAGTTTCTCCTGAA...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rpoB': SeqRecord(seq=Seq('ATGCTTGGGGTTGGGAATGAAGGAATGTCTACACTACCCGGATTGAATCAAATA...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnC-GCA': SeqRecord(seq=Seq('GGCGGCATGGCCGAGTGGTAAGGCGGGGGACTGCAAATCCTTTTTCCCCAGTTC...CCT', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'petN': SeqRecord(seq=Seq('ATGGATATAGTAAGTCTTGCTTGGGCTGCTTTAATGGTAGTTTTTACATTTTCC...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'clpP': SeqRecord(seq=Seq('ATGCCAGTAGGTGTTCCAAAGGTAGGCTTTGAAGTTCCTGGCGATGATGAAGCT...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'rpl20': SeqRecord(seq=Seq('ATGACCAGAACTAAACGAGGATCCATCGCTCGAAAACGGAGGACAAAAAGTACA...TAG', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'rps18': SeqRecord(seq=Seq('ATGGATAAATCTAAGCGACTGTTTGTTAAATCGAAGCAATCAATTCGTAGGAGT...TAG', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rpl33': SeqRecord(seq=Seq('ATGGCTAAGCGTAAAGATATCCGAATAACTGTTATTTTAGAATGTACTAGTTGT...TAG', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'psaJ': SeqRecord(seq=Seq('ATGCGAGATCTAAAAACATATCTCTCTGTGGCACCAGTGGCAAGTACTTTATGG...TGA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnP-UGG': SeqRecord(seq=Seq('AGGGATGTAGCGCAGCTTGGTAGCGCGTTTGTTTTGGGTACAAAATGTCACAGG...CTA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'trnW-CCA': SeqRecord(seq=Seq('GCGCTCTTAGTTCAGTTCGGTAGAACGCGGGTCTCCAAAACCCGATGTCGTAGG...GCG', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'petG': SeqRecord(seq=Seq('ATGATTGAAGTTTTTCTATTTGGAATCGTGTTAGGTCTAATTCCTATTACTTTG...TGA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'petL': SeqRecord(seq=Seq('ATGCCCACTCTAACCAGTTATTTCGGTTTTCTACTAGCAGTTTTAACTATAACC...TGA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'psbB': SeqRecord(seq=Seq('ATGGGTTTGCCTTGGTATCGTGTTCATACCGTTGTATTAAATGATCCCGGCCGT...TGA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'psbT': SeqRecord(seq=Seq('ATGGAAGCATTGGTTTATACATTTCTCTTAGTCTCGACTTTAGGAATTATTTTT...TGA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'psbN': SeqRecord(seq=Seq('ATGGAAACAGCAACCCTAATTGCCATCTCTATATCCGGTTTAATTGTAAGTTTT...TAG', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'psbH': SeqRecord(seq=Seq('ATGGCTACACAAACCGTTGAGGGTAGTTCTAGATCTGGTCCAAGACGAACTGTT...TAG', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'petB': SeqRecord(seq=Seq('ATGAGTGTGTGACTTGTTATAATGGATCCTATTGATAGTACAGAACATAAAATG...TAG', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'petD': SeqRecord(seq=Seq('ATGGGAGTGTGTGACTTGAACTATTGATTGGGCCATGCAGATATATGACTTAAT...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[]),\n",
       " 'rpoA': SeqRecord(seq=Seq('ATGATTCGAGAAAAACTAAAAGTATCGACTCAGACACTACAGTGGAAGTGTGTT...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rps11': SeqRecord(seq=Seq('ATGGCAAAATCTATACCAAAAATTGGTTCACGGAAAACTGGACGTATTGGTTCG...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rpl36': SeqRecord(seq=Seq('ATGAAAGTGGCGGCCTCTGTTCGTAAAATTTGTGAAAAATGTCGATTGATCCGT...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rps8': SeqRecord(seq=Seq('GTGGGTAAAGATACTATCGCTGATATACTAACCTTGATACGCAATGCTGACATG...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rpl14': SeqRecord(seq=Seq('ATGATTCAACCTCAAACCTATTTAAATGTAGCCGATAATAGTGGAGCGCGCGAA...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rpl16': SeqRecord(seq=Seq('ATGCTTAGTGTGTGACTCGTTGGTTTTTTTTTTATAGTGGTTAAATTTCTACAA...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rps3': SeqRecord(seq=Seq('ATGGGACAAAAAATACATCCGCTTGGTTTCAGACTTGGTGCAACTCAAAGTCAT...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rps19': SeqRecord(seq=Seq('ATGACACGTTCACGAAAAAAAAACCCTTTTGTAGCGAATCATTTATTAAGAAAA...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rpl2': SeqRecord(seq=Seq('ATGGCGATACATTTATCCAAAACTTCTAGCCCGAGCACACGCAATGGAGCCGTA...TAG', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rpl23': SeqRecord(seq=Seq('ATGAATGGAATCAAATATGCAGTATTTACAGAAAAAAGTATTCGGTTATGGGGA...CAT', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'ccsA': SeqRecord(seq=Seq('ATGATATTTTCAACTTTAGAGCATATATTAACTCATATATCGTTTTCCCTGATA...TGA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'trnL-UAG': SeqRecord(seq=Seq('GCCGCTATGGTGAAATTGGTAGACACGCTGCTCTTAGGAAGCAGTGCTAGAGCA...GCA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'rpl32': SeqRecord(seq=Seq('ATGCCAGTTCCAAAAAAACGCACTTCTATATCAAAAAAAAAAATTCGTAAGAAT...TAA', IUPACAmbiguousDNA()), id='<unknown id>', name='<unknown name>', description='<unknown description>', dbxrefs=[]),\n",
       " 'ndhF': SeqRecord(seq=Seq('ATGGAATATACACATCCATCTTCATGGATCATACCTTTCATCCCTCTTCCAGTT...TAA', IUPACAmbiguousDNA()), id='MG859922.1', name='MG859922', description='Pisum sativum subsp. elatius isolate CE1 chloroplast, complete genome', dbxrefs=[])}"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "genes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[Seq('ATG', IUPACAmbiguousDNA()),\n",
       " Seq('CCA', IUPACAmbiguousDNA()),\n",
       " Seq('ACT', IUPACAmbiguousDNA()),\n",
       " Seq('AGA', IUPACAmbiguousDNA()),\n",
       " Seq('AAA', IUPACAmbiguousDNA()),\n",
       " Seq('CAA', IUPACAmbiguousDNA()),\n",
       " Seq('CTT', IUPACAmbiguousDNA()),\n",
       " Seq('ATT', IUPACAmbiguousDNA()),\n",
       " Seq('AGA', IUPACAmbiguousDNA()),\n",
       " Seq('AAC', IUPACAmbiguousDNA()),\n",
       " Seq('ACA', IUPACAmbiguousDNA()),\n",
       " Seq('AGA', IUPACAmbiguousDNA()),\n",
       " Seq('CAG', IUPACAmbiguousDNA()),\n",
       " Seq('CCA', IUPACAmbiguousDNA()),\n",
       " Seq('ATC', IUPACAmbiguousDNA()),\n",
       " Seq('CGG', IUPACAmbiguousDNA()),\n",
       " Seq('ACT', IUPACAmbiguousDNA()),\n",
       " Seq('GTA', IUPACAmbiguousDNA()),\n",
       " Seq('ACA', IUPACAmbiguousDNA()),\n",
       " Seq('AAA', IUPACAmbiguousDNA()),\n",
       " Seq('TCT', IUPACAmbiguousDNA()),\n",
       " Seq('CCC', IUPACAmbiguousDNA()),\n",
       " Seq('GCT', IUPACAmbiguousDNA()),\n",
       " Seq('CTT', IUPACAmbiguousDNA()),\n",
       " Seq('CGG', IUPACAmbiguousDNA()),\n",
       " Seq('GGA', IUPACAmbiguousDNA()),\n",
       " Seq('TGC', IUPACAmbiguousDNA()),\n",
       " Seq('CCT', IUPACAmbiguousDNA()),\n",
       " Seq('CAA', IUPACAmbiguousDNA()),\n",
       " Seq('CGC', IUPACAmbiguousDNA()),\n",
       " Seq('CGA', IUPACAmbiguousDNA()),\n",
       " Seq('GGA', IUPACAmbiguousDNA()),\n",
       " Seq('ACA', IUPACAmbiguousDNA()),\n",
       " Seq('TGT', IUPACAmbiguousDNA()),\n",
       " Seq('ACT', IUPACAmbiguousDNA()),\n",
       " Seq('AGG', IUPACAmbiguousDNA()),\n",
       " Seq('GTG', IUPACAmbiguousDNA()),\n",
       " Seq('TAT', IUPACAmbiguousDNA()),\n",
       " Seq('GTG', IUPACAmbiguousDNA()),\n",
       " Seq('CGA', IUPACAmbiguousDNA()),\n",
       " Seq('CTC', IUPACAmbiguousDNA()),\n",
       " Seq('GTT', IUPACAmbiguousDNA()),\n",
       " Seq('ACT', IUPACAmbiguousDNA()),\n",
       " Seq('ATT', IUPACAmbiguousDNA()),\n",
       " Seq('AAC', IUPACAmbiguousDNA()),\n",
       " Seq('CCT', IUPACAmbiguousDNA()),\n",
       " Seq('AAA', IUPACAmbiguousDNA()),\n",
       " Seq('AAA', IUPACAmbiguousDNA()),\n",
       " Seq('CCA', IUPACAmbiguousDNA()),\n",
       " Seq('AAC', IUPACAmbiguousDNA()),\n",
       " Seq('TCT', IUPACAmbiguousDNA()),\n",
       " Seq('GCG', IUPACAmbiguousDNA()),\n",
       " Seq('TTA', IUPACAmbiguousDNA()),\n",
       " Seq('CGT', IUPACAmbiguousDNA()),\n",
       " Seq('AAA', IUPACAmbiguousDNA()),\n",
       " Seq('GTT', IUPACAmbiguousDNA()),\n",
       " Seq('GCC', IUPACAmbiguousDNA()),\n",
       " Seq('AGA', IUPACAmbiguousDNA()),\n",
       " Seq('GTA', IUPACAmbiguousDNA()),\n",
       " Seq('CGC', IUPACAmbiguousDNA()),\n",
       " Seq('TTA', IUPACAmbiguousDNA()),\n",
       " Seq('ACC', IUPACAmbiguousDNA()),\n",
       " Seq('TCT', IUPACAmbiguousDNA()),\n",
       " Seq('GGA', IUPACAmbiguousDNA()),\n",
       " Seq('TTT', IUPACAmbiguousDNA()),\n",
       " Seq('GAA', IUPACAmbiguousDNA()),\n",
       " Seq('ATC', IUPACAmbiguousDNA()),\n",
       " Seq('ACT', IUPACAmbiguousDNA()),\n",
       " Seq('GCT', IUPACAmbiguousDNA()),\n",
       " Seq('TAT', IUPACAmbiguousDNA()),\n",
       " Seq('ATA', IUPACAmbiguousDNA()),\n",
       " Seq('CCT', IUPACAmbiguousDNA()),\n",
       " Seq('GGT', IUPACAmbiguousDNA()),\n",
       " Seq('ATT', IUPACAmbiguousDNA()),\n",
       " Seq('GGC', IUPACAmbiguousDNA()),\n",
       " Seq('CAT', IUPACAmbiguousDNA()),\n",
       " Seq('AAT', IUPACAmbiguousDNA()),\n",
       " Seq('TTA', IUPACAmbiguousDNA()),\n",
       " Seq('CAA', IUPACAmbiguousDNA()),\n",
       " Seq('GAA', IUPACAmbiguousDNA()),\n",
       " Seq('CAT', IUPACAmbiguousDNA()),\n",
       " Seq('TCT', IUPACAmbiguousDNA()),\n",
       " Seq('GTA', IUPACAmbiguousDNA()),\n",
       " Seq('GTC', IUPACAmbiguousDNA()),\n",
       " Seq('TTA', IUPACAmbiguousDNA()),\n",
       " Seq('GTA', IUPACAmbiguousDNA()),\n",
       " Seq('AGA', IUPACAmbiguousDNA()),\n",
       " Seq('GGG', IUPACAmbiguousDNA()),\n",
       " Seq('GGA', IUPACAmbiguousDNA()),\n",
       " Seq('AGG', IUPACAmbiguousDNA()),\n",
       " Seq('GTT', IUPACAmbiguousDNA()),\n",
       " Seq('AAG', IUPACAmbiguousDNA()),\n",
       " Seq('GAT', IUPACAmbiguousDNA()),\n",
       " Seq('TTA', IUPACAmbiguousDNA()),\n",
       " Seq('CCT', IUPACAmbiguousDNA()),\n",
       " Seq('GGT', IUPACAmbiguousDNA()),\n",
       " Seq('GTG', IUPACAmbiguousDNA()),\n",
       " Seq('AGA', IUPACAmbiguousDNA()),\n",
       " Seq('TAT', IUPACAmbiguousDNA()),\n",
       " Seq('CAC', IUPACAmbiguousDNA()),\n",
       " Seq('ATT', IUPACAmbiguousDNA()),\n",
       " Seq('GTT', IUPACAmbiguousDNA()),\n",
       " Seq('CGA', IUPACAmbiguousDNA()),\n",
       " Seq('GGA', IUPACAmbiguousDNA()),\n",
       " Seq('ACC', IUPACAmbiguousDNA()),\n",
       " Seq('CTG', IUPACAmbiguousDNA()),\n",
       " Seq('GAT', IUPACAmbiguousDNA()),\n",
       " Seq('GCT', IUPACAmbiguousDNA()),\n",
       " Seq('GTC', IUPACAmbiguousDNA()),\n",
       " Seq('GCA', IUPACAmbiguousDNA()),\n",
       " Seq('GTA', IUPACAmbiguousDNA()),\n",
       " Seq('AAG', IUPACAmbiguousDNA()),\n",
       " Seq('GAT', IUPACAmbiguousDNA()),\n",
       " Seq('CGT', IUPACAmbiguousDNA()),\n",
       " Seq('CAA', IUPACAmbiguousDNA()),\n",
       " Seq('CAA', IUPACAmbiguousDNA()),\n",
       " Seq('GGG', IUPACAmbiguousDNA()),\n",
       " Seq('CGT', IUPACAmbiguousDNA()),\n",
       " Seq('TCA', IUPACAmbiguousDNA()),\n",
       " Seq('AAA', IUPACAmbiguousDNA()),\n",
       " Seq('TAT', IUPACAmbiguousDNA()),\n",
       " Seq('GGA', IUPACAmbiguousDNA()),\n",
       " Seq('GCC', IUPACAmbiguousDNA()),\n",
       " Seq('AAA', IUPACAmbiguousDNA()),\n",
       " Seq('AAG', IUPACAmbiguousDNA()),\n",
       " Seq('CCA', IUPACAmbiguousDNA()),\n",
       " Seq('AAA', IUPACAmbiguousDNA()),\n",
       " Seq('TAA', IUPACAmbiguousDNA())]"
      ]
     },
     "execution_count": 35,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "get_codons(genes[\"rps12\"].seq)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 260,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "128"
      ]
     },
     "execution_count": 260,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(get_codons(genes[\"rps12\"].seq))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 79,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/sivico26/miniconda3/lib/python3.6/site-packages/Bio/Seq.py:163: BiopythonWarning: Biopython Seq objects now use string comparison. Older versions of Biopython used object comparison. During this transition, please use hash(id(my_seq)) or my_dict[id(my_seq)] if you want the old behaviour, or use hash(str(my_seq)) or my_dict[str(my_seq)] for the new string hashing behaviour.\n",
      "  \"the new string hashing behaviour.\", BiopythonWarning)\n"
     ]
    }
   ],
   "source": [
    "frec_ex = get_frecs(get_codons(mercado[0].seq))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 81,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'GGC': 0.005033886651605933,\n",
       " 'GAA': 0.02615165504370887,\n",
       " 'CGA': 0.009674884588940183,\n",
       " 'CGG': 0.006580885964050683,\n",
       " 'TTG': 0.02146154601709066,\n",
       " 'AAC': 0.014512326883410274,\n",
       " 'CCG': 0.006458108240840782,\n",
       " 'CGC': 0.004469109124840389,\n",
       " 'ATG': 0.01699243689225027,\n",
       " 'GTG': 0.006777330321186524,\n",
       " 'GAT': 0.01991454670464591,\n",
       " 'TCA': 0.019742657892152047,\n",
       " 'CAA': 0.021289657204596798,\n",
       " 'TCC': 0.01532265985659562,\n",
       " 'ACT': 0.013431882919163148,\n",
       " 'GCC': 0.005500441999803556,\n",
       " 'ATC': 0.02006187997249779,\n",
       " 'CAC': 0.008054218642569493,\n",
       " 'GCT': 0.008471662901483155,\n",
       " 'ACA': 0.013259994106669285,\n",
       " 'CTC': 0.011197328356742952,\n",
       " 'TAA': 0.025193988802671645,\n",
       " 'CTA': 0.016525881544052646,\n",
       " 'TGT': 0.013014438660249485,\n",
       " 'AAT': 0.03661231706119242,\n",
       " 'TAC': 0.013751104999508889,\n",
       " 'ATA': 0.032634318829191634,\n",
       " 'TTT': 0.05082997740889893,\n",
       " 'TAT': 0.03290442982025341,\n",
       " 'ACC': 0.010386995383557607,\n",
       " 'TTA': 0.024923877811609863,\n",
       " 'TTC': 0.02799332089185738,\n",
       " 'TCT': 0.022983989784893428,\n",
       " 'CTT': 0.02016010215106571,\n",
       " 'AAA': 0.0461644239269227,\n",
       " 'AGA': 0.01895688046360868,\n",
       " 'CAT': 0.015862881838719183,\n",
       " 'CTG': 0.008766329437186917,\n",
       " 'GTA': 0.014831548963756016,\n",
       " 'AGG': 0.011295550535310873,\n",
       " 'GAG': 0.010509773106767508,\n",
       " 'CCA': 0.01404577153521265,\n",
       " 'TGA': 0.01988999116000393,\n",
       " 'ATT': 0.03644042824869856,\n",
       " 'CGT': 0.006580885964050683,\n",
       " 'GGA': 0.012989883115607503,\n",
       " 'GCA': 0.008643551713977016,\n",
       " 'GGG': 0.008152440821137414,\n",
       " 'TGC': 0.008348885178273254,\n",
       " 'CAG': 0.007440330026519988,\n",
       " 'CCT': 0.011663883704940576,\n",
       " 'TGG': 0.013382771829879186,\n",
       " 'GCG': 0.003658776151655044,\n",
       " 'AAG': 0.017532658874373833,\n",
       " 'TAG': 0.014954326686965917,\n",
       " 'GTT': 0.016722325901188488,\n",
       " 'ACG': 0.0062862194283469204,\n",
       " 'CCC': 0.00947844023180434,\n",
       " 'AGT': 0.012817994303113643,\n",
       " 'GTC': 0.008127885276495432,\n",
       " 'TCG': 0.009895884490718003,\n",
       " 'GAC': 0.007538552205087909,\n",
       " 'GGT': 0.008938218249680777,\n",
       " 'AGC': 0.007808663196149691}"
      ]
     },
     "execution_count": 81,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "frec_ex"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1.0"
      ]
     },
     "execution_count": 80,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sum(frec_ex.values())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 85,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Seq('E', ExtendedIUPACProtein())"
      ]
     },
     "execution_count": 85,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "get_codons(mercado[0].seq)[1].translate(table=11)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 86,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1.5625"
      ]
     },
     "execution_count": 86,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 82,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/sivico26/miniconda3/lib/python3.6/site-packages/Bio/Seq.py:163: BiopythonWarning: Biopython Seq objects now use string comparison. Older versions of Biopython used object comparison. During this transition, please use hash(id(my_seq)) or my_dict[id(my_seq)] if you want the old behaviour, or use hash(str(my_seq)) or my_dict[str(my_seq)] for the new string hashing behaviour.\n",
      "  \"the new string hashing behaviour.\", BiopythonWarning)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "64\n",
      "64\n",
      "64\n",
      "64\n",
      "64\n",
      "64\n",
      "64\n",
      "64\n",
      "64\n",
      "64\n",
      "64\n",
      "64\n",
      "64\n",
      "64\n"
     ]
    }
   ],
   "source": [
    "for plant in mercado:\n",
    "    print(len(get_frecs(get_codons(plant.seq))))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "122198\n",
      "155570\n",
      "154740\n",
      "157790\n",
      "125319\n",
      "152772\n",
      "122967\n",
      "155886\n",
      "155312\n",
      "159636\n",
      "151169\n",
      "153366\n",
      "155452\n",
      "161453\n"
     ]
    }
   ],
   "source": [
    "for plant in mercado:\n",
    "    print(len(plant.seq))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "CGG\n",
      "ATG\n",
      "TTA\n",
      "CTT\n",
      "ATC\n",
      "CCT\n",
      "GCC\n",
      "CCA\n",
      "GTG\n",
      "CGT\n",
      "AAA\n",
      "CAT\n",
      "ACA\n",
      "ATT\n",
      "TCT\n",
      "TCA\n",
      "GGC\n",
      "GAT\n",
      "TGT\n",
      "GCT\n",
      "GGA\n",
      "ATA\n",
      "AAT\n",
      "CGC\n",
      "TAA\n",
      "TTT\n",
      "CAA\n",
      "ACT\n",
      "GAA\n",
      "CAC\n",
      "AAG\n",
      "AGA\n",
      "TAT\n",
      "CTC\n",
      "GCG\n",
      "GTC\n",
      "GGG\n",
      "CAG\n",
      "AGG\n",
      "GGT\n",
      "TGC\n",
      "AAC\n",
      "GCA\n",
      "GTA\n",
      "CCC\n",
      "ACC\n",
      "GTT\n",
      "CTG\n",
      "CGA\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/sivico26/miniconda3/lib/python3.6/site-packages/Bio/Seq.py:163: BiopythonWarning: Biopython Seq objects now use string comparison. Older versions of Biopython used object comparison. During this transition, please use hash(id(my_seq)) or my_dict[id(my_seq)] if you want the old behaviour, or use hash(str(my_seq)) or my_dict[str(my_seq)] for the new string hashing behaviour.\n",
      "  \"the new string hashing behaviour.\", BiopythonWarning)\n"
     ]
    }
   ],
   "source": [
    "for i in set(get_codons(genes[\"rps12\"].seq)):\n",
    "    print(i)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 167,
   "metadata": {},
   "outputs": [],
   "source": [
    "rec = mercado[1]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 177,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "338"
      ]
     },
     "execution_count": 177,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(rec.features)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 225,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "OrderedDict([('gene', ['rps12']), ('trans_splicing', [''])])\n",
      "OrderedDict([('gene', ['rps12']), ('trans_splicing', ['']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein S12']), ('protein_id', ['AXQ38216.1']), ('translation', ['MPTIKQLIRNTRQPIRNVTKSPALRGCPQRRGTCTRVYTITPKKPNSALRKVARVRLTSGFEITAYIPGIGHNLQEHSVVLVRGGRVKDLPGVRYHIVRGTLDAVGVKDRQQGRSKYGVKKPK'])])\n",
      "OrderedDict([('gene', ['trnH'])])\n",
      "OrderedDict([('gene', ['trnH']), ('product', ['tRNA-His'])])\n",
      "OrderedDict([('gene', ['psbA'])])\n",
      "OrderedDict([('gene', ['psbA']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II protein D1']), ('protein_id', ['AXQ38247.1']), ('translation', ['MTAILERRESESLWGRFCNWITSTENRLYIGWFGVLMIPTLLTATSVFIIAFIAAPPVDIDGIREPVSGSLLYGNNIISGAIIPTSAAIGLHFYPIWEAASVDEWLYNGGPYELIVLHFLLGVACYMGREWELSFRLGMRPWIAVAYSAPVAAATAVFLIYPIGQGSFSDGMPLGISGTFNFMIVFQAEHNILMHPFHMLGVAGVFGGSLFSAMHGSLVTSSLIRETTENESANEGYRFGQEEETYNIVAAHGYFGRLIFQYASFNNSRSLHFFLAAWPVVGIWFTALGISTMAFNLNGFNFNQSVVDSQGRVINTWADIINRANLGMEVMHERNAHNFPLDLAAIEAPSTNG'])])\n",
      "OrderedDict([('gene', ['trnK'])])\n",
      "OrderedDict([('gene', ['trnK']), ('product', ['tRNA-Lys'])])\n",
      "OrderedDict([('gene', ['trnK'])])\n",
      "OrderedDict([('gene', ['trnK'])])\n",
      "OrderedDict([('gene', ['matK'])])\n",
      "OrderedDict([('gene', ['matK']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['maturase K']), ('protein_id', ['AXQ38231.1']), ('translation', ['MEEIHRYLQPDSSQPHNFLYPLIFQEYIYALAQDHGLNRNRSILLENSGYNNKFSFLIVKRLITRMYQQNHFLISTNDSNKNPFLGCNKSLYSQMISEGFSCIVEIPFSLRLISSLSSFEEKKIFKSHNLRSIHSTFPFLEDNLSHLNYVLDILIPYPVHLEILVQTLRYWVKDASSLHLLRFFLHEYCNLNSLITSKKPGYSFSKKNPKFFFFLYNSYVYECESTFVFLRKQSFHLRSTSFGALLERIYFYGKIERLVKVFAKDFQVTLWLFKDPLMHYVRYEGKSILASKGTFLLMNKWKFYLVNFWQCHFSMYFHTGRIHINQLSNHSRDFMGYRSSVRLNHSMVRSQMFENSFLINNPIKKFDTLVPIIPLIGSLAKAHFCTVLGHPISKPVWSDLSDSDIIDRFGRICRNLFHYYSGSSKKKTLYRIKYILRLSCARTLARKHKSTVRTFLKRSGSELLEEFLTSEEQVLSLTFPRASSSLWGVYRSRIWYLDIFCINDLANYQ'])])\n",
      "OrderedDict([('gene', ['trnK'])])\n",
      "OrderedDict([('gene', ['rps16'])])\n",
      "OrderedDict([('gene', ['rps16']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein S16']), ('protein_id', ['AXQ38246.1']), ('translation', ['MVKLRLKRCGRKQRAVYRIVAIDVRSRREGKDLQKVGFYDPIRNQTYLNVPAILYFLEKGAQPTETVQDILKKAEVFKELRPTKFN'])])\n",
      "OrderedDict([('gene', ['rps16'])])\n",
      "OrderedDict([('gene', ['rps16'])])\n",
      "OrderedDict([('gene', ['rps16'])])\n",
      "OrderedDict([('gene', ['trnQ'])])\n",
      "OrderedDict([('gene', ['trnQ']), ('product', ['tRNA-Gln'])])\n",
      "OrderedDict([('gene', ['psbK'])])\n",
      "OrderedDict([('gene', ['psbK']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II protein K']), ('protein_id', ['AXQ38283.1']), ('translation', ['MLNTFSLIGICLNSTLYSSSFFFGKLPEAYAFLNPIVDIMPVIPLFFFLLAFVWQAAVSFR'])])\n",
      "OrderedDict([('gene', ['psbI'])])\n",
      "OrderedDict([('gene', ['psbI']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II protein I']), ('protein_id', ['AXQ38292.1']), ('translation', ['MLTLKLFVYTVVIFFVSLFIFGFLSNDPGRNPGREE'])])\n",
      "OrderedDict([('gene', ['trnS'])])\n",
      "OrderedDict([('gene', ['trnS']), ('product', ['tRNA-Ser'])])\n",
      "OrderedDict([('gene', ['trnG'])])\n",
      "OrderedDict([('gene', ['trnG']), ('product', ['tRNA-Gly'])])\n",
      "OrderedDict([('gene', ['trnG'])])\n",
      "OrderedDict([('gene', ['trnG'])])\n",
      "OrderedDict([('gene', ['trnG'])])\n",
      "OrderedDict([('gene', ['trnR'])])\n",
      "OrderedDict([('gene', ['trnR']), ('product', ['tRNA-Arg'])])\n",
      "OrderedDict([('gene', ['atpA'])])\n",
      "OrderedDict([('gene', ['atpA']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ATP synthase CF1 alpha subunit']), ('protein_id', ['AXQ38233.1']), ('translation', ['MVTIRADEISNIIRERIEQYNREVKIVNTGTVLQVGDGIARIHGLDEVMAGELVEFEEGTIGIALNLESNNVGVVLMGDGLLIQEGSSVKATGRIAQIPVSEAYLGRVVNALAKPIDGRGEISASEFRLIESAAPGIISRRSVYEPLQTGLIAIDSMIPIGRGQRELIIGDRQTGKTAVATDTILNQQGQNVICVYVAIGQKASSVAQVVTTLQERGAMEYTIVVAETADSPATLQYLAPYTGAALAEYFMYRERHTLIIYDDLSKQAQAYRQMSLLLRRPPGREAYPGDVFYLHSRLLERAAKLSSSLGEGSMTALPIVETQSGDVSAYIPTNVISITDGQIFLSADLFNSGIRPAINVGISVSRVGSAAQIKAMKQVAGKLKLELAQFAELEAFAQFASDLDKATQNQLARGQRLRELLKQSQSAPLTVEEQIMTIYTGTNGYLDSLEVGQVRKFLVELRTYLKTTKPQFQEIISSTKTFTEEAEDLLKEAIQEQMDRFILQEQA'])])\n",
      "OrderedDict([('gene', ['atpF'])])\n",
      "OrderedDict([('gene', ['atpF']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ATP synthase CF0 subunit I']), ('protein_id', ['AXQ38243.1']), ('translation', ['MKNVTDSFVSLGHWPSAGSFGFNTDILATNPINLSVVLGVLIFFGKGVLSDLLDNRKQRILNTIRNSEELRGGAIEQLEKARSRLRKVETEAEQFRVNGYSEIEREKLNLINSTYKTLEQLENYKNETIQFEQQRAINQVRQRVFQQALRGALGTLNSCLNNELHLRTISANIGMLGTMKEITD'])])\n",
      "OrderedDict([('gene', ['atpF'])])\n",
      "OrderedDict([('gene', ['atpF'])])\n",
      "OrderedDict([('gene', ['atpF'])])\n",
      "OrderedDict([('gene', ['atpH'])])\n",
      "OrderedDict([('gene', ['atpH']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ATP synthase CF0 subunit III']), ('protein_id', ['AXQ38278.1']), ('translation', ['MNPLISAASVIAAGLAVGLASIGPGVGQGTAAGQAVEGIARQPEAEGKIRGTLLLSLAFMEALTIYGLVVALALLFANPFV'])])\n",
      "OrderedDict([('gene', ['atpI'])])\n",
      "OrderedDict([('gene', ['atpI']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ATP synthase CF0 subunit IV']), ('protein_id', ['AXQ38252.1']), ('translation', ['MNVLSCSINTLKGLYDISGVEVGQHFYWQIGGFQVHGQVLITSWVVIAILLGSATIAVRNPQTIPTGGQNFFEYVLEFIRDVSKTQIGEEYGPWVPFIGTMFLFIFVSNWSGALLPWKIIQLPHGELAAPTNDINTTVALALLTSVAYFYAGLTKKGLGYFGKYIQPTPILLPINILEDFTKPLSLSFRLFGNILADELVVVVLVSLVPLVVPIPVMLLGLFTSGIQALIFATLAAAYIGESMEGHH'])])\n",
      "OrderedDict([('gene', ['rps2'])])\n",
      "OrderedDict([('gene', ['rps2']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein S2']), ('protein_id', ['AXQ38254.1']), ('translation', ['MTRRYWNINLEEMMEAGVHFGHGTRKWNPKMAPYISAKRKGIHITNLTRTARFLSEACDLVFDAASRGKQFLIVGTKNKAADSVEWAAIRARCHYVNKKWLGGMLTNWSTTETRLHKFRDLRMEQKTGRLNRLPKRDAAMLKRQLSRLQTYLGGIKYMTGVPDIVIIVDQHEEYTALRECITLGIPTICLTDTNCDPYLADISIPANDDAISSIRLILNKLVFAICEGRSSYIRNP'])])\n",
      "OrderedDict([('gene', ['rpoC2'])])\n",
      "OrderedDict([('gene', ['rpoC2']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', [\"RNA polymerase beta' subunit\"]), ('protein_id', ['AXQ38220.1']), ('translation', ['MAERANLVFHNKAIDGTAMKRLISRLIEHFGMAYTSHILDQVKTLGFQQATATSISLGIDDLLTIPSKGWLVQDAEQQSLILEKHHHYGNVHAVEKLRQSIEIWYATSEYLRQEMNPNFRMTDPFNPVHIMSFSGARGNASQVHQLVGMRGLMSDPQGQMIDLPIQSNLREGLSLTEYIISCYGARKGVVDTAVRTSDAGYLTRRLVEVVQHIVVRRTDCGTARGISVSPRNGIMPERIFSQTLIGRVLADDIYMGSRCIATRNQAIGIGLVNRFITFRAQPISIRTPFTCRSTSWICRLCYGRSPTHGDLVELGEAVGIIAGQSIGEPGTQLTLRTFHTGGVFTGGTAEHVRAPSNGKIKFNEDLVHPTRTRHGHPAFLCSMDLYVTIESEDILHTVNIPPKSLLLVQNDQYVESEQVIAEIRAGISTLNFKETVRKHIYSDSDGEMHWSTDVYHAPQFTYGNVHLLPKTSHLWILLGRPCRSSLVYLSIHKDQDQMNAHSPSGKRRYTSNLSVTNDQARQKLFSSDFYGQKDDRIPDYSDLNRIICTGQYNLVYSPILHENSDLLSKRRRNKFIIPLHSIQELENELMPRSGISIEIPVNGSFRRNSILAYFDDPRYRRKSSGIIKYGTIETHSVIKKEDLIEYRGVKEFRPKYQMKVDRFFFIPEEVHILPGSSSIMVRNNSIVGVDTQITLNLRSRVGGLVRVEKKKKRIELKIFSGDIHFPGETDKISRHTGVLIPPGTGKRNSKEYKKVKNWIYVQRITPSKKRFFVLVRPVVTYEITDGINLATLFPPDPLQERDNVQLRIVNYILYENGKPIRGISDTSIQLVRTCLVLNWNQDKKSSSCEEARASFVEIRTNGLIRHFLRINLVKSPISYIGKRNDPSGSGLLSDNGSDCTNINPFSAIYSYSKAKMHQSLNQPQGTIHTLLNRNKECQSLIILSAANCSRMGPFKDVKYHSVIKESIKKDPLIPIRNSLGPLGTCLPIENFYSSYHLITHNQILVTKYLQLDNLKQTFQVIQLKYYLMDENGKIFNPDPCRNIILNPFNLNWFFLHHNYCAETSKIISLGQFICENVCIAKKGSPLKSGQVILVQVDSIVIRSAKPYLATPGATVHGHYGETLYEGDTLVTFIYEKSRSGDITQGLPKVEQVLEVRSIDSISMNLEKRVEVWNKCIPRILGIPWGFLIGAELTIAQSRISLVNKIQQVYRSQGVQIHNRHIEIIVRQITSKVLISEDGMSNVFSPGELIGLLRAERMGRALEEAICYRVVLLGITRASLNTQSFISEASFQETARVLAKAALRGRIDWLKGLKENVVLGGVIPVGTGFKGLVHPSKQHNNIPLETKKTNLFEGEMRDILFHHRKLFDSCLSKKFHDIPEQSFKSFIGFNHS'])])\n",
      "OrderedDict([('gene', ['rpoC1'])])\n",
      "OrderedDict([('gene', ['rpoC1']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['RNA polymerase beta']), ('protein_id', ['AXQ38222.1']), ('translation', ['MIDRYKHQQLRIGSVSPQQISAWATKILPNGEIVGEVTKPYTFHYKTNKPEKDGLFCERIFGPLKSGICACGNYRVIGDEKEDPKFCEQCGVEFVDSRIRRYQMGYIKLACPVTHVWYLKRLPSYIANLLDKPLKELEGLVYCDVYPNFSFARPITKKPTFLRLRGLFEYEIQSWKYSIPLFFTTQGFDTFRNREISTGAGAIREQLADLDLRIIIENSLVEWEELGEEGHTGNEWEDRKVGRRKDFLVRRVELAKHFIRTTIEPEWMVLCLLPVLPPELRPIIQIDGGKLMSSDINELYRRVIYRNNTLTDLLTTSRSTPGELVMCQEKLVQEAVDTLLDNGIRGQPMRDGHNKVYKSFSDVIEGKEGRFRETLLGKRVDYSGRSVIVVGPSLSLHRCGLPREIAIELFQTFVIRGLIRQHLASNIGVAKSKIREKEPIVWEILQEVMQGHPVLLNRAPTLHRLGIQAFQPVLVEGRAICLHPLVCKGFNADFDGDQMAVHVPLSLEAQVEARLLMFSHMNLLSPAIGDPISVPTQDMLIGLYVLTSGNHRGICVNRYNPCNRRNYQKRSDNSYYKYTKEPFFSNSYDAIGAYRQKRINLDSPLWLRWRLDQRVIASRETPVEVHYESLGTFYEIYGHYLIVRSLKKKILFVYIRTTVGHIALYREIEEAIQGFSRAYSYGT'])])\n",
      "OrderedDict([('gene', ['rpoC1'])])\n",
      "OrderedDict([('gene', ['rpoC1'])])\n",
      "OrderedDict([('gene', ['rpoC1'])])\n",
      "OrderedDict([('gene', ['rpoB'])])\n",
      "OrderedDict([('gene', ['rpoB']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['RNA polymerase beta subunit']), ('protein_id', ['AXQ38221.1']), ('translation', ['MLGDGNEGISTIPGFNQIQFEGFCRFIDQGLTEELYKFPKIEDTDQEIEFQLFVETYQLVEPLIKERDAVYESLTYSSELYVSAGLIWKNSRDMQEQTIFIGNIPLMNSLGTSIVNGIYRIVINQILQSPGIYYRSELDHNGISVYTGTIISDWGGRSELEIDRKARIWARVSRKQKISILVLSSAMGLNLREILENVCYPEIFLSFLNDKERKKIGSKENSILEFYQQFACVGGDPVFSESLCKELQKKFFQQRCELGRIGRRNMNRKLNLDIPQNNTFLLPRDILAAADHLIGLKFGMGALDDMNHLKNKRIRSVADLLQDQFGLALVRLENVVRGTICGAIRHKLIPTPQNLVTSPPLTTTYESFFGLHPLSQVLDRTNPLTQIVHGRKLSYLGPGGLTGRTASFRIRDIHPSHYGRICPIDTSEGINVGLIGSLAIHARIGHWGSLESPFYEISERSTGVRMLYLSPGSDEYYMVAAGNSLALNRDIQEEQVVPARYRQEFLTIAWEQVHLRSIFPFQYFSIGASLIPFIEHNDANRALMSSNMQRQAVPLSRSEKCIVGTGLERQAALDSGALAIAEREGRVVYTNTDKILLAGNGDILSIPLVIYQRSNKNTCMHQKLRVPRGKCIKKGQILADGAATVGGELALGKNVLVAYMPWEGYNFEDAVLISERLVYEDIYTSFHIRKYEIHTHVTSQGPEKVTNEIPHLEAHLLRNLDKKGIVMLGSWVETGDILVGKLTPQVVKESSYAPEDRLLRAILGIQVSTSKETCLKLPIGGRGRVIDVRWIQKRGGSSYNPETIRVYISQKREIKVGDKVAGRHGNKGIISKILPRQDMPYLQDGRSIDMVFNPLGVPSRMNVGQIFECSLGLAGSLLDRHYRIAPFDERYEQEASRKLVFSELYEASKQTANPWVFEPEYPGKSRIFDGRTGNPFEQPVIIGKPYILKLIHQVDDKIHGRSSGHYALVTQQPLRGRAKQGGQRVGEMEVWALEGFGVAHILQEMLTYKSDHIRARQEILGTTIIGGTIPNPEDAPESFRLLVRELRSLALELNHFLVSEKNFQINRKEA'])])\n",
      "OrderedDict([('gene', ['trnC-GCA'])])\n",
      "OrderedDict([('gene', ['trnC-GCA']), ('product', ['tRNA-Cys'])])\n",
      "OrderedDict([('gene', ['petN'])])\n",
      "OrderedDict([('gene', ['petN']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['cytochrome b6/f complex subunit VIII']), ('protein_id', ['AXQ38297.1']), ('translation', ['MDIVSLAWAALMVVFTFSLSLVVWGRSGL'])])\n",
      "OrderedDict([('gene', ['psbM'])])\n",
      "OrderedDict([('gene', ['psbM']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II protein M']), ('protein_id', ['AXQ38294.1']), ('translation', ['MEVNILAFIATALFILVPTAFLLIIYVKTVSQND'])])\n",
      "OrderedDict([('gene', ['trnD-GUC'])])\n",
      "OrderedDict([('gene', ['trnD-GUC']), ('product', ['tRNA-Asp'])])\n",
      "OrderedDict([('gene', ['trnY-GUA'])])\n",
      "OrderedDict([('gene', ['trnY-GUA']), ('product', ['tRNA-Tyr'])])\n",
      "OrderedDict([('gene', ['trnE-UUC'])])\n",
      "OrderedDict([('gene', ['trnE-UUC']), ('product', ['tRNA-Glu'])])\n",
      "OrderedDict([('gene', ['trnT-GGU'])])\n",
      "OrderedDict([('gene', ['trnT-GGU']), ('product', ['tRNA-Thr'])])\n",
      "OrderedDict([('gene', ['psbD'])])\n",
      "OrderedDict([('gene', ['psbD']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II protein D2']), ('protein_id', ['AXQ38248.1']), ('translation', ['MTIAIGKFTKEENDLFDIMDDWLRRDRFVFVGWSGLLLFPCAYFAVGGWFTGTTFVTSWYTHGLASSYLEGCNFLTAAVSTPANSLAHSLLLLWGPEAQGDFTRWCQLGGLWTFVALHGAFGLIGFMLRQFELARSVQLRPYNAIAFSGPIAVFVSVFLIYPLGQSGWFFAPSFGVAAIFRFILFFQGFHNWTLNPFHMMGVAGVLGAALLCAIHGATVENTLFEDGDGANTFRAFNPTQAEETYSMVTANRFWSQIFGVAFSNKRWLHFFMLFVPVTGLWMSALGVVGLALNLRAYDFVSQEIRAAEDPEFETFYTKNILLNEGIRAWMAAQDQPHENLIFPEEVLPRGNAL'])])\n",
      "OrderedDict([('gene', ['psbC'])])\n",
      "OrderedDict([('gene', ['psbC']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II CP43 chlorophyll apoprotein']), ('protein_id', ['AXQ38241.1']), ('translation', ['MKTLYSLRRFYHVETLFNGTLALAGRDQETTGFAWWAGNARLINLSGKLLGAHVAHAGLIVFWAGAMNLFEVAHFVPEKPMYEQGLILLPHLATLGWGVGPGGEVIDTFPYFVSGVLHLISSAVLGFGGIYHALLGPETLEESFPFFGYVWKDRNKMTTILGIHLILLGIGAFLLVFKALYFGGVYDTWAPGGGDVRKITNLTLSPSIIFGYLLKSPFGGEGWIVSVDDLEDIIGGHVWLGSICILGGIWHILTKPFAWARRALVWSGEAYLSYSLGALAVFGFIACCFVWFNNTAYPSEFYGPTGPEASQAQAFTFLVRDQRLGANVGSAQGPTGLGKYLMRSPTGEVIFGGETMRFWDLRAPWLEPLRGPNGLDLSRLKKDIQPWQERRSAEYMTHAPLGSLNSVGGVATEINAVNYVSPRSWLATSHFVLGFFFFVGHLWHAGRARAAAAGFEKGIDRDLEPVLFMTPLN'])])\n",
      "OrderedDict([('gene', ['trnS-UGA'])])\n",
      "OrderedDict([('gene', ['trnS-UGA']), ('product', ['tRNA-Ser'])])\n",
      "OrderedDict([('gene', ['psbZ'])])\n",
      "OrderedDict([('gene', ['psbZ']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II protein Z']), ('protein_id', ['AXQ38282.1']), ('translation', ['MTLAFQLAVFALIATSLILLISVPVVFASPDGWSSNKNVVFSGTSLWIGLVFLVGILNSLIS'])])\n",
      "OrderedDict([('gene', ['trnG-UCC'])])\n",
      "OrderedDict([('gene', ['trnG-UCC']), ('product', ['tRNA-Gly'])])\n",
      "OrderedDict([('gene', ['trnfM-CAU'])])\n",
      "OrderedDict([('gene', ['trnfM-CAU']), ('product', ['tRNA-Met'])])\n",
      "OrderedDict([('gene', ['rps14'])])\n",
      "OrderedDict([('gene', ['rps14']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein S14']), ('protein_id', ['AXQ38273.1']), ('translation', ['MARKSLIQREKKRQKLEQKYHSIRRSSKKEISKVPSLSDKWEIYGKLQSLPRNSAPTRLHRRCFLTGRPRANYRDFGLSGHILREMVHACLLPGATRSSW'])])\n",
      "OrderedDict([('gene', ['psaB'])])\n",
      "OrderedDict([('gene', ['psaB']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem I P700 apoprotein A2']), ('protein_id', ['AXQ38228.1']), ('translation', ['MALRFPRFSQGLAQDPTTRRIWFGIATAHDFESHDDITEERLYQNIFASHFGQLAIIFLWTSGNLFHVAWQGNFESWVQDPLHVRPIAHAIWDPHFGQPAVEAFTRGGALGPVNIAYSGVYQWWYTIGLRTNEDLYTGALFLLFLSAISLLAGWLHLQPKWKPSVSWFKNAESRLNHHLSGLFGVSSLAWTGHLVHVAIPASRGEYVRWNNFLDVLPHPQGLGPLFTGQWNLYAQNPDSSSHLFGTAQGAGTAILTLLGGFHPQTQSLWLTDIAHHHLAIAFIFLVAGHMYRTNFGIGHSMKDLLDAHIPPGGRLGRGHKGLYDTINNSLHFQLGLALASLGVITSLVAQHMYSLPAYAFIAQDFTTQAALYTHHQYIAGFIMTGAFAHGAIFFIRDYNPEQNEDNVLARMLDHKEAIISHLSWASLFLGFHTLGLYVHNDVMLAFGTPEKQILIEPIFAQWIQSAHGKTSYGFDVLLSSTSGPAFNAGRSIWLPGWLNAVNENSNSLFLTIGPGDFLVHHAIALGLHTTTLILVKGALDARGSKLMPDKKDFGYSFPCDGPGRGGTCDISAWDAFYLAVFWMLNTIGWVTFYWHWKHITLWQGNVSQFNESSTYLMGWLRDYLWLNSSQLINGYNPFGMNSLSVWAWMFLFGHLVWATGFMFLISWRGYWQELIETLAWAHERTPLANLIRWRDKPVALSIVQARLVGLAHFSVGYIFTYAAFLIASTSGKFG'])])\n",
      "OrderedDict([('gene', ['psaA'])])\n",
      "OrderedDict([('gene', ['psaA']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem I p700 apoprotein A1']), ('protein_id', ['AXQ38224.1']), ('translation', ['MIIRSPEPEVKILVDRDPVKTSFEEWARPGHFSRTIAKGPDTTTWIWNLHADAHDFDSHTSDLEEISRKVFSAHFGQLSIIFLWLSGMYFHGARFSNYEAWLSDPTHIGPSAQVVWPIVGQEILNGDVGGGFRGIQITSGFFQLWRASGITSELQLYCTAIGALIFAALMLFAGWFHYHKAAPKLAWFQDVESMLNHHLAGLLGLGSLSWAGHQVHVSLPINQFLNAGVDPKEIPLPHEFILNRDLLAQLYPSFAEGATPFFTLNWSKYADFLTFRGGLDPVTGGLWLTDIAHHHLAIAILFLIAGHMYRTNWGIGHGLKDILEAHKGPFTGQGHKGLYEILTTSWHAQLSLNLAMLGSLTIVVAHHMYSMPPYPYLATDYGTQLSLFTHHMWIGGFLIVGAAAHAAIFMVRDYDPTTRYNDLLDRVLRHRDAIISHLNWACIFLGFHSFGLYIHNDTMSALGRPQDMFSDTAIQLQPVFAQWIQNTHALAPSATAPGATASTSLTWGGGDLVAVGGKVALLPIPLGTADFLVHHIHAFTIHVTVLILLKGVLFARSSRLIPDKANLGFRFPCDGPGRGGTCQVSAWDHVFLGLFWMYNSISVVIFHFSWKMQSDVWGSVSDQGVVTHITGGNFAQSSITINGWLRDFLWAQASQVIQSYGSSLSAYGLFFLGAHFVWAFSLMFLFSGRGYWQELIESIVWAHNKLKVAPATQPRALSIIQGRAVGVTHYLLGGIATTWAFFLARIIAVG'])])\n",
      "OrderedDict([('gene', ['ycf3'])])\n",
      "OrderedDict([('gene', ['ycf3']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['hypothetical chloroplast RF3']), ('protein_id', ['AXQ38230.1']), ('translation', ['MPRSRINGNFIDKTFSIVANILLRVIPTTSGEKEAFTYYRDGAIMSAQSEGNYAEALQNYYEAMRLEIDPYDRSYILYNIGLIHTSNGEHTKALEYYFRALERNPFLPQAFNNMAVICHYRGEQAIQQGDSEIAEAWFDQAAEYWKQAIALTPGNYIEAHNWLKITRRFE'])])\n",
      "OrderedDict([('gene', ['ycf3'])])\n",
      "OrderedDict([('gene', ['ycf3'])])\n",
      "OrderedDict([('gene', ['ycf3'])])\n",
      "OrderedDict([('gene', ['ycf3'])])\n",
      "OrderedDict([('gene', ['ycf3'])])\n",
      "OrderedDict([('gene', ['trnS-GGA'])])\n",
      "OrderedDict([('gene', ['trnS-GGA']), ('product', ['tRNA-Ser'])])\n",
      "OrderedDict([('gene', ['rps4'])])\n",
      "OrderedDict([('gene', ['rps4']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein S4']), ('protein_id', ['AXQ38257.1']), ('translation', ['MSRYRGPRFKKIRRLGALPGLTNKKPRTGSDLRNQSRSGKKSQYRIRLEEKQKLRFHYGLTERQLLKYVRIARKAKGSTGQVLLQLLEMRLDNILFRLGMASTIPAARQLVNHRHILVNGRIVDIPSYRCKPRDIITAKDEQKSRALIQISLDSSPHEELPNHLTLQPFQYKGLVNQIIDSKWVGLKINELLVVEYYSRQT'])])\n",
      "OrderedDict([('gene', ['trnT-UGU'])])\n",
      "OrderedDict([('gene', ['trnT-UGU']), ('product', ['tRNA-Thr'])])\n",
      "OrderedDict([('gene', ['trnL-UAA'])])\n",
      "OrderedDict([('gene', ['trnL-UAA']), ('product', ['tRNA-Leu'])])\n",
      "OrderedDict([('gene', ['trnL-UAA'])])\n",
      "OrderedDict([('gene', ['trnL-UAA'])])\n",
      "OrderedDict([('gene', ['trnL-UAA'])])\n",
      "OrderedDict([('gene', ['trnF-GAA'])])\n",
      "OrderedDict([('gene', ['trnF-GAA']), ('product', ['tRNA-Phe'])])\n",
      "OrderedDict([('gene', ['ndhJ'])])\n",
      "OrderedDict([('gene', ['ndhJ']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['NADH-plastoquinone oxidoreductase subunit J']), ('protein_id', ['AXQ38261.1']), ('translation', ['MQGRLSAWLVKHGLIHRSLGFDYQGIETLQIKPEDWHSIAVIFYVYGYNYLRSQCAYDVAPGGLLASVYHLTRIEDGVDQPEEVCIKVFASRRNPRIPSVFWVWKSVDFQERESYDMLGISYDNHPRLKRILMPESWIGWPLRKDYIAPNFYEIQDAH'])])\n",
      "OrderedDict([('gene', ['ndhK'])])\n",
      "OrderedDict([('gene', ['ndhK']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['NADH-plastoquinone oxidoreductase subunit K']), ('protein_id', ['AXQ38253.1']), ('translation', ['MVLAPEYSDNKKKNGKNKIETVMNSIQFPLLDRTAPNSVISTTLNDLSNWSRLSSLWPLLYGTSCCFIEFASLIGSRFDFDRYGLVPRSSPRQADLILTAGTVTMKMAPSLVRLYEQMPEPKYVIAMGACTITGGMFSTDSYSTVRGVDKLIPVDVYLPGCPPKPEAVIDAITKLRKKISRELYEDRIRSQRANRCFTTNHKFHVRRSIHTGNYDQRVLYQPPSTSEIPTEIFFKYKNSVSSAELVN'])])\n",
      "OrderedDict([('gene', ['ndhC'])])\n",
      "OrderedDict([('gene', ['ndhC']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['NADH-plastoquinone oxidoreductase subunit 3']), ('protein_id', ['AXQ38270.1']), ('translation', ['MFLLYEYDFFWAFLIISILVPILAFFISGVLAPISKGPEKLSTYESGIEPMGDAWLQFRIRYYMFALVFVVFDVETVFLYPWAMSFDVLGVSVFIEALIFVLILIIGLVYAWRKGALEWS'])])\n",
      "OrderedDict([('gene', ['trnV-UAC'])])\n",
      "OrderedDict([('gene', ['trnV-UAC']), ('product', ['tRNA-Val'])])\n",
      "OrderedDict([('gene', ['trnV-UAC'])])\n",
      "OrderedDict([('gene', ['trnV-UAC'])])\n",
      "OrderedDict([('gene', ['trnV-UAC'])])\n",
      "OrderedDict([('gene', ['trnM-CAU'])])\n",
      "OrderedDict([('gene', ['trnM-CAU']), ('product', ['tRNA-Met'])])\n",
      "OrderedDict([('gene', ['atpE'])])\n",
      "OrderedDict([('gene', ['atpE']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ATP synthase CF1 epsilon subunit']), ('protein_id', ['AXQ38267.1']), ('translation', ['MTLNLSVLTPNRIVWDSEVEEIVLSTNSGQIGILPNHAPIATAVDIGILRIRLNDQWLTMALMGGFARIGNNEITVLVNDAERGSDIDPQEAQQTLEIAEANVKKAEGRRQKIEANLALRRARTRVEASNPIS'])])\n",
      "OrderedDict([('gene', ['atpB'])])\n",
      "OrderedDict([('gene', ['atpB']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ATP synthase CF1 beta subunit']), ('protein_id', ['AXQ38234.1']), ('translation', ['MRINPTTSGSGVSTLEKKNLGRVVQIIGPVLDVAFPPGKMPNIYNALVVQGRDSVGQPINVACEVQQLLGNNRVRAVAMSATDGLTRGMEVIDKGAPISVPVGGATLGRIFNVLGEPVDNLGPVDTSTTSPIHRSAPAFIQLDTKLSIFETGIKVVDLLAPYRRGGKIGLFGGAGVGKTVLIMELINNIAKAHGGVSVFGGVGERTREGNDLYMEMKESGVINKENIAESKVALVYGQMNEPPGARMRVGLTALTMAEYFRDVNEQDVLLFIDNIFRFVQAGSEVSALLGRMPSAVGYQPTLSTEMGSLQERITSTKEGSITSIQAVYVPADDLTDPAPATTFAHLDATTVLSRGLAAKGIYPAVDPLDSTSTMLQPRIVGEEHYETAQRVKQTLQRYKELQDIIAILGLDELSEEDRLLVARARKIERFLSQPFFVAEVFTGSPGKYVGLAETIRGFQLILSGELDGLPEQAFYLVGTIDEATAKAMNLEMESNLKK'])])\n",
      "OrderedDict([('gene', ['rbcL'])])\n",
      "OrderedDict([('gene', ['rbcL']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit']), ('protein_id', ['AXQ38239.1']), ('translation', ['MSPQTETKASVGFKAGVKEYKLTYYTPEYQTKDTDILAAFRVTPQPGVPPEEAGAAVAAESSTGTWTTVWTDGLTSLDRYKGRCYRIERVVGEKDQYIAYVAYPLDLFEEGSVTNMFTSIVGNVFGFKALRALRLEDLRIPPAYVKTFQGPPHGIQVERDKLNKYGRPLLGCTIKPKLGLSAKNYGRAVYECLRGGLDFTKDDENVNSQPFMRWRDRFCFCAEALFKAQAETGEIKGHYLNATAGTCEEMIKRAVFARELGVPIVMHDYLTGGFTANTTLAHYCRDNGLLLHIHRAMHAVIDRQKNHGIHFRVLAKALRMSGGDHIHAGTVVGKLEGERDITLGFVDLLRDDFVEQDRSRGIYFTQDWVSLPGVIPVASGGIHVWHMPALTEIFGDDSVLQFGGGTLGHPWGNAPGAVANRVALEACVKARNEGRDLAQEGNEIIREASKWSPELAAACEVWKEIVFNFAAVDVLDK'])])\n",
      "OrderedDict([('gene', ['accD'])])\n",
      "OrderedDict([('gene', ['accD']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['acetyl-CoA carboxyltransferase beta subunit']), ('protein_id', ['AXQ38238.1']), ('translation', ['MERWWFNSMLFKKEFERRCGLNKSMGSLGPSENTSEDLNLKVKNIHSCSNVDYLFSVKDIRNFISDDTFLVSDRNGDSYSIYFDIENQIFEIDNDYSFLRELESSFYSYRNSSYLNNGFRGEDPYYNSYMSYMYDTQYSWNNHINSCIDNYLQSQICIDTSIISGNENYGDSYIYRAICSGESRNSSEKQGSSRRTRTKGSDLTIRESSNDLEVTQKYKHLWVQCENCYGLNYKKFLKSKMNICEQCGYHLKMSSSDRIELLIDPGTWDPMDEDMVSLDPIEFHSEEEPYKDRIDSYQRKTGLTEAVQTGIGQLNGIPVAIGVMDFQFMGGSMGSVVGEKITRLIEHAANQKLPLIIVCASGGARMQEGSLSLMQMAKISSALYDYQFNKKLFYVSILTSPTTGGVTASFGMLGDIIIAEPNAYIAFAGKRVIEQTLNKTIPEGSQAAEYLFQKGLFDLIVPRNLLKSVLSELFKLHAFFPLNQKSSKIK'])])\n",
      "OrderedDict([('gene', ['psaI'])])\n",
      "OrderedDict([('gene', ['psaI']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem I subunit VIII']), ('protein_id', ['AXQ38293.1']), ('translation', ['MTNLNLPSIFVPLVGLVFPAIAMASLFLHVQKNKIV'])])\n",
      "OrderedDict([('gene', ['ycf4'])])\n",
      "OrderedDict([('gene', ['ycf4']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['hypothetical chloroplast RF4']), ('protein_id', ['AXQ38258.1']), ('translation', ['MTWRSDDIWIELITGSRKISNFCWAFILFLGSLGFLLVGTSSYLGRNLLSFFPPQQIIFFPQGIVMSFYGIAGLFISSYLWCTISWNVGSGYDRFDRKEGVVCIFRWGFPGKNRRIFLRFLIKDIQSVRIEVKEGIYARRVLYMDIRGQGSIPLTRTDENLTPREIEQKAAELAYFLRVPIEVF'])])\n",
      "OrderedDict([('gene', ['cemA'])])\n",
      "OrderedDict([('gene', ['cemA']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['chloroplast envelope membrane protein']), ('protein_id', ['AXQ38255.1']), ('translation', ['MAKKKAFTPLFYLASIVFLPWWISFSVNKCLESWVTNWWNTGQSEIFLNNIQEKSLLEKFIELEELLFLDEMIKEYSETHLEEFGIGIHKETIQLIKIQNENRIHTILHFSTNTICFIILSGYSILGNEKLVILNSWAQEFLYNLSDTVKAFSILLLTDLCIGFHSPQGWELIMCSIYKDFGFVHNDQIIAGLVSTFPVILDTIFKYWIFRYLNRLSPSLVVIYHSMND'])])\n",
      "OrderedDict([('gene', ['petA'])])\n",
      "OrderedDict([('gene', ['petA']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['cytochrome f']), ('protein_id', ['AXQ38250.1']), ('translation', ['MQTRNAFSWLKKQITRSISVSLMIYILTRTSISSAYPIFAQQGYENPREATGRIVCANCHLANKPVEIEVPQAVLPDTVFEAVVRIPYDMQLKQVLANGKKGGLNVGAVLILPEGFELAPPDRISPEMKEKIGNLSFQSYRPNKKNILVVGPVPGKKYSEITFPILSPDPATKKDVHFLKYPIYVGGNRGRGQIYPDGNKSNNTVYNATAAGIVSKIIRKEKGGYEITITDASEGRQVVDIIPPGPELLVSEGESIKFDQPLTSNPNVGGFGQGDAEIVLQDPLRVQGLLFFLASVILAQIFLVLKKKQFEKVQLAEMNF'])])\n",
      "OrderedDict([('gene', ['psbJ'])])\n",
      "OrderedDict([('gene', ['psbJ']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II protein J']), ('protein_id', ['AXQ38287.1']), ('translation', ['MADTTGRIPLWIIGTVAGILVIGLIGVFFYGSYSGLGSSL'])])\n",
      "OrderedDict([('gene', ['psbL'])])\n",
      "OrderedDict([('gene', ['psbL']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II protein L']), ('protein_id', ['AXQ38289.1']), ('translation', ['MTQSNPNEQNVELNRTSLYWGLLLIFVLAVLFSNYFFN'])])\n",
      "OrderedDict([('gene', ['psbF'])])\n",
      "OrderedDict([('gene', ['psbF']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II cytochrome b559 beta subunit']), ('protein_id', ['AXQ38288.1']), ('translation', ['MTIDRTYPIFTVRWLAVHGLAVPTVFFLGSISAMQFIQR'])])\n",
      "OrderedDict([('gene', ['psbE'])])\n",
      "OrderedDict([('gene', ['psbE']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II cytochrome b559 alpha subunit']), ('protein_id', ['AXQ38277.1']), ('translation', ['MSGSTGERSFADIITSIRYWVIHSITIPSLFIAGWLFVSTGLAYDVFGSPRPNEYFTESRQGIPLITGRFDPLEQLDEFSKSF'])])\n",
      "OrderedDict([('gene', ['petL'])])\n",
      "OrderedDict([('gene', ['petL']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['cytochrome b6/f complex subunit VI']), ('protein_id', ['AXQ38296.1']), ('translation', ['MLTITSYFGFLLAALTITSALFIGLSKIRLI'])])\n",
      "OrderedDict([('gene', ['petG'])])\n",
      "OrderedDict([('gene', ['petG']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['cytochrome b6/f complex subunit V']), ('protein_id', ['AXQ38290.1']), ('translation', ['MIEVFLFGIVLGLIPITLAGLFVTAYLQYRRGDQLDL'])])\n",
      "OrderedDict([('gene', ['trnW-CCA'])])\n",
      "OrderedDict([('gene', ['trnW-CCA']), ('product', ['tRNA-Trp'])])\n",
      "OrderedDict([('gene', ['trnP-UGG'])])\n",
      "OrderedDict([('gene', ['trnP-UGG']), ('product', ['tRNA-Pro'])])\n",
      "OrderedDict([('gene', ['psaJ'])])\n",
      "OrderedDict([('gene', ['psaJ']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem I subunit IX']), ('protein_id', ['AXQ38285.1']), ('translation', ['MRDLKTYLSVAPVLSTLWFGALAGLLIEINRFFPDALTFPFFSF'])])\n",
      "OrderedDict([('gene', ['rpl33'])])\n",
      "OrderedDict([('gene', ['rpl33']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein L33']), ('protein_id', ['AXQ38281.1']), ('translation', ['MAKGKDVRVTVILECTSCVRNSVDKVSRGISRYITQKNRHNTPNRFELKKFCPYCYKHTIHGEIKK'])])\n",
      "OrderedDict([('gene', ['rps18'])])\n",
      "OrderedDict([('gene', ['rps18']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein S18']), ('protein_id', ['AXQ38271.1']), ('translation', ['MDKSKRPFLKFKRSFRRRLPPIQSGDRIDYRNMSLISRFISEQGKILSRRVNRLTLKQQRLITLAIKQARILSLLPFLNNEKQFERTESTARTTGFKARNK'])])\n",
      "OrderedDict([('gene', ['rpl20'])])\n",
      "OrderedDict([('gene', ['rpl20']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein L20']), ('protein_id', ['AXQ38268.1']), ('translation', ['MTRIKRGYIARRRRTKIRLFASSFRAAHSRLTRTITQQKIRALVSAHRDRDRKKRDFRRLWITRINAVIRERGVSYSYSRLIHDLYKRQLLLNRKILAQIAILNRNCLYMISNEIINHKGSRLERIQQNNLN'])])\n",
      "OrderedDict([('gene', ['rps12']), ('trans_splicing', [''])])\n",
      "OrderedDict([('gene', ['rps12']), ('trans_splicing', ['']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein S12']), ('protein_id', ['AXQ38215.1']), ('translation', ['MPTIKQLIRNTRQPIRNVTKSPALRGCPQRRGTCTRVYTITPKKPNSALRKVARVRLTSGFEITAYIPGIGHNLQEHSVVLVRGGRVKDLPGVRYHIVRGTLDAVGVKDRQQGRSKYGVKKPK'])])\n",
      "OrderedDict([('gene', ['clpP'])])\n",
      "OrderedDict([('gene', ['clpP']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['clp protease proteolytic subunit']), ('protein_id', ['AXQ38229.1']), ('translation', ['MPIGIPKVKIKRRGIRRRLWTDIYNRLYRERLLFLGQGIGTGLGNQLICLLIYLSMEDEDKDLYMFINSPGGWVVPGLALYDTMQYVRPDIHTVCIGLAASMSSVVLVGGEVNKRIAFPRARVMMHEPFSGFYMAQVGEFVVEAGEMLKLRGILARIYEKKTGKPFWVISEDMERDVFMSATEALNYGIVDYVAVDGRKKANKMYANYAKLL'])])\n",
      "OrderedDict([('gene', ['clpP'])])\n",
      "OrderedDict([('gene', ['clpP'])])\n",
      "OrderedDict([('gene', ['clpP'])])\n",
      "OrderedDict([('gene', ['clpP'])])\n",
      "OrderedDict([('gene', ['clpP'])])\n",
      "OrderedDict([('gene', ['psbB'])])\n",
      "OrderedDict([('gene', ['psbB']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II CP47 chlorophyll apoprotein']), ('protein_id', ['AXQ38232.1']), ('translation', ['MGLPWYRVHTVVLNDPGRLLSVHIMHTALVAGWAGSMALYELAVFDPSDPVLDPMWRQGMFVIPFMTRLGITNSWGGWSITGGTVTNPGIWSYEGVAGAHIVFSGLCFLAAIWHWVYWDLEIFCDERTGKPSLDLPKIFGIHLFLSGVACFGFGAFHVTGLYGPGIWVSDPYGLTGKVQPVNPAWGVEGFDPFVPGGIASHHIAAGTLGILAGLFHLSVRPPQRLYKGLRMGNIETVLSSSIAAVFFAAFVVAGTMWYGSATTPIELFGPTRYQWDQGYFQQEIYRRVSAGLAENQSLSEAWSKIPEKLAFYDYIGNNPAKGGLFRAGSMDNGDGIAVGWLGHPIFRDKEGRELFVRRMPTFFETFPVVLVDGDGIVRADVPFRRAESKYSVEQVGVTVEFYGGELNGVSYSDPATVKKYARRAQLGEIFELDRATLKSDGVFRSSPRGWFTFGHASFALLFFFGHIWHGARTLFRDVFAGIDPDLDAQVEFGAFQKLGDPTTKRQAA'])])\n",
      "OrderedDict([('gene', ['psbT'])])\n",
      "OrderedDict([('gene', ['psbT']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II protein T']), ('protein_id', ['AXQ38295.1']), ('translation', ['MEALVYTFLLVSTLGIIFFAIFFREPPKVPTKKN'])])\n",
      "OrderedDict([('gene', ['psbN'])])\n",
      "OrderedDict([('gene', ['psbN']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II protein N']), ('protein_id', ['AXQ38286.1']), ('translation', ['METATLVAIFISGLLVSFTGYALYTAFGQPSQQLRDPFEEHGD'])])\n",
      "OrderedDict([('gene', ['psbH'])])\n",
      "OrderedDict([('gene', ['psbH']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem II phosphoprotein']), ('protein_id', ['AXQ38280.1']), ('translation', ['MATQTVENSSRSGPRRTAVGDLLKPLNSEYGKVAPGWGTTPLMGVAMALFAVFLSIILEIYNSSVLLDGISMN'])])\n",
      "OrderedDict([('gene', ['petB'])])\n",
      "OrderedDict([('gene', ['petB']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['cytochrome b6']), ('protein_id', ['AXQ38242.1']), ('translation', ['MSKVYDWFEERLEIQAIADDITSKYVPPHVNIFYCLGGITLTCFLVQVATGFAMTFYYRPTVTEAFASVQYIMTEANFGWLIRSVHRWSASMMVLMMILHVFRVYLTGGFKKPRELTWVTGVVLAVLTASFGVTGYSLPWDQIGYWAVKIVTGVPDAIPVIGSPLVELLRGSASVGQSTLTRFYSLHTFVLPLLTAVFMLMHFPMIRKQGISGPL'])])\n",
      "OrderedDict([('gene', ['petB'])])\n",
      "OrderedDict([('gene', ['petB'])])\n",
      "OrderedDict([('gene', ['petB'])])\n",
      "OrderedDict([('gene', ['petD'])])\n",
      "OrderedDict([('gene', ['petD']), ('note', ['cytochrome b6/f complex subunit 4']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['cytochrome b6/f complex subunit IV']), ('protein_id', ['AXQ38244.1']), ('translation', ['MGVTKKPDLNDPVLRAKLAKGMGHNYYGEPAWPNDLLYIFPVVILGTIACNVGLAVLEPSMIGEPADPFATPLEILPEWYFFPVFQILRTVPNKLLGVLLMVSVPVGLLTVPFLENVNKFQNPFRRPVATTVFLIGTAVALWLGIGATLPIDKSLTLGLF'])])\n",
      "OrderedDict([('gene', ['petD'])])\n",
      "OrderedDict([('gene', ['petD'])])\n",
      "OrderedDict([('gene', ['petD'])])\n",
      "OrderedDict([('gene', ['rpoA'])])\n",
      "OrderedDict([('gene', ['rpoA']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['RNA polymerase alpha subunit']), ('protein_id', ['AXQ38249.1']), ('translation', ['MVREKVTVSTRTLQWKCVESRTDSKRLYYGRFILSPLMKGQADTIGIAMRRALLGEIEGTCITRVKSEKVPHEYSTITGIHESVHEILMNLKEIVLRSKLYGTSDASICVKGPGYVTAQDIILPPYVEIVDNTQHIASLTEPIDFCIGLQIERNRGYLIKTPHNFQDGSYPIDAVFMPVRNANHSIHSYGNGNEKQEILFIEIWTNGSLTPKEALHEASRNLIDLFIPFLHMEEDNLYLQDNQHTVPLSPFTFHDKLAKLLKNKNKIALKSIFIDQSELSSRIYNCLKRSNIYTLVDLLNNSQEDLMKIEHFRSEDIKQILGILEKYFVIDLAKNKF'])])\n",
      "OrderedDict([('gene', ['rps11'])])\n",
      "OrderedDict([('gene', ['rps11']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein S11']), ('protein_id', ['AXQ38265.1']), ('translation', ['MAKAIPKISSRRNGRIGSRKGARRIPKGVIHVQASFNNTIVTVTDVRGRVVSWSSAGTSGFKGTRRGTPFAAQTAAANAIRTVVDQGMQRAEVMIKGPGLGRDAALRAIRRSGILLTFVRDVTPMPHNGCRPPKKRRV'])])\n",
      "OrderedDict([('gene', ['rpl36'])])\n",
      "OrderedDict([('gene', ['rpl36']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein L36']), ('protein_id', ['AXQ38291.1']), ('translation', ['MKIRASVRKICEKCRLIRRRGRIIVICSNPRHKQRQG'])])\n",
      "OrderedDict([('gene', ['infA']), ('pseudo', [''])])\n",
      "OrderedDict([('gene', ['rps8'])])\n",
      "OrderedDict([('gene', ['rps8']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein S8']), ('protein_id', ['AXQ38266.1']), ('translation', ['MGRDTIAEIITSIRNADMDRKRVVRIVSTNITENIVQLLLREGFIENVRKHRENNKYFLVLTLRHRRNRKRPYRNIFNLKRISRPGLRIYSNYQRIPRILGGIGIVILSTSRGIMTDREARLEGIGGEILCYIW'])])\n",
      "OrderedDict([('gene', ['rpl14'])])\n",
      "OrderedDict([('gene', ['rpl14']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein L14']), ('protein_id', ['AXQ38269.1']), ('translation', ['MIQPQTHLNVADNSGARELMCIRIIGASNRRYAHIGDVIVAVIKEAVPNMPLQRSEVVRAVIVRTCKELKRDNGMIIRYDDNAAVVIDQEGNPKGTRIFGAIARELRDLNFTKIVSLAPEVL'])])\n",
      "OrderedDict([('gene', ['rpl16'])])\n",
      "OrderedDict([('gene', ['rpl16']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein L16']), ('protein_id', ['AXQ38240.1']), ('translation', ['MLSPKRTRFRKQHRGRMKGISYRGNRISFGKYALQALEPAWITSRQIEAGRRAMTRNARRGGKIWVRIFPDKPVTLRPAETRMGSGKGSPEYWVAVVKPGRILYEMGGVTENIARRAISLAASKMPIRTQFIIS'])])\n",
      "OrderedDict([('gene', ['rpl16'])])\n",
      "OrderedDict([('gene', ['rpl16'])])\n",
      "OrderedDict([('gene', ['rpl16'])])\n",
      "OrderedDict([('gene', ['rps3'])])\n",
      "OrderedDict([('gene', ['rps3']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein S3']), ('protein_id', ['AXQ38256.1']), ('translation', ['MGQKINPLGFRLGTTQGHHSLWFSQPKNYSEGLQEDKKIRDCIKNYVQKNMRTSSGIEGIARIEIQKRIDLIQVIIFMGFPKLLIESRPRGIEELQMTLQKEFNCVNRKLNIAVTRIGKPYGNPNILAEFIAGQLKNRVSFRKAMKKAIELTEQADTKGIQIQIAGRIDGKEIARVEWIREGRVPLQTIQAKIDYCSYTVRTIYGILGIKIWIFLDEE'])])\n",
      "OrderedDict([('gene', ['rpl22'])])\n",
      "OrderedDict([('gene', ['rpl22']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein L22']), ('protein_id', ['AXQ38262.1']), ('translation', ['MLKKKKTEVYALGEHISMSADKARRVIDQIRGRSYEETLMILELMPYRACYAILKLVYSAAANASYNLGSSETNLVISKAEVNEGTIVKKLKPRARGRSFPIKRSTCHIIIVMKDISLDDEYGEMYSLKKTRWKKKSTAMTYRDMYNSGGLWDKK'])])\n",
      "OrderedDict([('gene', ['rps19']), ('pseudo', [''])])\n",
      "OrderedDict([('note', ['inverted repeat B (IRb)']), ('rpt_type', ['inverted'])])\n",
      "OrderedDict([('gene', ['rpl2'])])\n",
      "OrderedDict([('gene', ['rpl2']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein L2']), ('protein_id', ['AXQ38236.1']), ('translation', ['MAIHLYKTSTPSTRNETVDSQVKSNPRNNLIYGQRRCGKGRNARGIITARHRGGGHKRLYRKIDFRRNEKDIYGRIVTIEYDPNRNAYICLIHYGDGEKRYILHPRGAIIGDTIVSGTEVPIKMGNALPLSTDMPLGTAIHNIEITLGKGGQLARAAGAVAKLIAKEGKSATLKLPSGEVRLISKNCSATVGQVGNVGVNQKSLGRAGSKRWLGKRPVVRGVVMNPVDHPHGGGEGRAPIGRKKPTTPWGYPALGRRSRKRNKYSDNLILRRRSK'])])\n",
      "OrderedDict([('gene', ['rpl2'])])\n",
      "OrderedDict([('gene', ['rpl2'])])\n",
      "OrderedDict([('gene', ['rpl2'])])\n",
      "OrderedDict([('gene', ['rpl23'])])\n",
      "OrderedDict([('gene', ['rpl23']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein L23']), ('protein_id', ['AXQ38274.1']), ('translation', ['MDGIKYAVFTDKSIRLLGKNQYTSNVESGSTRTEIKHWVELFFGVKVIAMNSHRLPGKSRRMGPIMGHTMHYRRMIITLQPGYSIPPLRKKRT'])])\n",
      "OrderedDict([('gene', ['trnI-CAU'])])\n",
      "OrderedDict([('gene', ['trnI-CAU']), ('product', ['tRNA-Ile'])])\n",
      "OrderedDict([('gene', ['ycf2'])])\n",
      "OrderedDict([('gene', ['ycf2']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['hypothetical chloroplast RF2']), ('protein_id', ['AXQ38217.1']), ('translation', ['MRGHQFKSWIFELREILREIKNSHHFLDSWTQFNSVGSFIHIFFHQERFLKLFDPRIWSILLSRNSQGSTSNRYFTIKGVILFVVAVLIYRINNRNMVERKNLYLIGLLPIPMHSIGPRNDTLEESVGSSNINRLIVSLLYLPKGKKISESCFLNPKESTWVLPITKKCSMPESNWGSRWWRNWIGKKRDFSCKISNETVAGIEILFKEKDLKYLEFLFVYYMDDPIRKDHDWELFDRLSLRKSRNRINLNSGPLFEILVKHWISYLMSAFREKIPIEVEGFFKQQGAGSTIQSNDIEHVSHLFSRNKWAITLQNCAQFHMWQFRQDLFVSWGKNPPESDFLRNVSRENWIWLDNVWLVNKDRFFSKVQNVSSNIQYDSTRSSFVQVTDSSQLKGSSDQSRDHLDSISNEDSEYHTLINQREIQQRKERSILWDPSFLQTERKEIESGRFPKCLSGYSSMSRLFTEREKQMINHLFPEEIEEFLGNPTRSVRSFFSDRWSELHLGSNPTERSTRDQKLLKKQQDLSFVPSRRSENKEMVNIFKIITYLQNTVSIHPISSDPGCDMVPKDEPDMDSSNKISFLNKNPFFDLFHLFHDRNRGGYTLHYDFASEERFQEMADLFTLSITEPDLVYHKGFAFSIDSCGLDQKQFLNEARDESKKKSLLVLPPIFYEENESFSRRIRKKWVRISCGNDLEDPKPKIVVFASNNIMEAVTQYRLIRNLIQIQYSTYGYIRNVLNRFFLMNRSDRNFEYGIQRDQIGKDTLNHRTIMKYTINQYLSNLKKSQKKWFEPLILISRTERSMNRDPDAYRYKWSNGSKNFQEHLEQSVSEQKSRFQVVFDRLRINQYSIDWSEVIDKKDLSKSLRFFLSKSLRFFLSKSLLFLSKLLLFLSNSLPFFCVSFGNIPIHRSEIYIYELKGPNDQLCNQLLESIGLQIVHLKKLKPFLLDDHDTSQKSKFLINGGTISPFLFNKIPKWMIDSFHTRNNRRKSFDSPDSYFSMIFHDQDNWLNPVKPFHRSSLISSFYKANRLRFLNNPHHFCFYWNTRFPFSVEKARINNSDFTYGQFLNILFIRNKIFSLCVGKKKHAFWGRDTISPIESQVSHIFIPNDFPQSGDETYNLYKSFHFPSRSDPSVRRAIYSIADISGTPLTEGQIVNFERTYCQPLSDMNLSDSEGKNLHQYLNFNSNMGLIHTPCSEKDLSSEKRKKWSLCLKKCVEKGQMYRTFQRDSAFSTLSKWNLFQTYMPWFLTSTGYKYLNLIFLDTFSDLLPILSSSQKFVSIFPDIMHGSGISWRILQKKLCLPQWNLISEISSKCLHNLLLSEEMIHRNNESPLISTHLRSPNAREFLYSILFLLLVAGYLVRTHLLFVSRASSELQTEFEKVKSLMIPSSMIELRKLLDRYPTSEPNSFWLKNLFLVALEQLGDSLEEIRGSASGGNMLGPAYGVKSIRSKKKDWNINLIEIIDLIPNPINRITFSRNTRHLSHTSKEIYSLIRKRKNVNGDWIDEKIESWVANSDSIDDEEKEFLVHFSTLTTENRIDQILLSLTHSDHLSKNDSGYQMIEQPGAIYLRYLVDIHKKHLMNYEFNPSCLAERRIFLAHYQTITYSQTSCGENSFHFPSHGKPFSLRLALSPSRGILVIGSIGTGRSYLVKYLATNSYVPFITVFLNKFLDNKSKGFLLDEIDIDDSDDIDDSDNLDASDDIDRDLDTELELLTRMNGLTVDMMPEIDRFYITLQFELAKAMSPCIIWIPNIHDLDVNESNDLSLGLLVNHLSRDCERCSTRNILVIASTHIPQKVDPALIAPNKLNTCIKIRRLLIPQQRKHFFTLSYTRGFHLEKKMFHTNGFGSITMGSNARDLVALTNEVLSISITQKKSIIDTNTIRSALHRQTWDLRSQVRSVQDHGILFYQIGRAVAQNVLLSNCPIEPISIYMKKKSCNGGDSYLYKWYFELGTSMKRLTILLYLLSCSAGSVAQDLWSLSVPDEKNGITSYGLVENDSDLVHGLLEVEGALVGSSRTEKDCSQFDNDRVTLLLRPEPRNPLDMMQKGSWSILDQRFLYEKYESEFEEGEGEGALDPQEDLFNHIVWAPRIWRPWGFLFDCIERPNELGFPYWSRSFRGKRIIYDEEDELQENDSGFLQSGTMQYQTRDRSSKEQGLFRISQFIWDPADPLFFLFKDQPPGSVFSHRELFADEEMSKGLLTSQTDPPTSIYKRWFIKNTQEKHFELLINRQRWLRTNSSLSNGSFRSNTLSESYQYLSNLFLSNGTLLDQMTKTLLRKRWLFPDEMKIGFM'])])\n",
      "OrderedDict([('gene', ['ycf15']), ('pseudo', [''])])\n",
      "OrderedDict([('gene', ['ndhB'])])\n",
      "OrderedDict([('gene', ['ndhB']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['NADH-plastoquinone oxidoreductase subunit 2']), ('protein_id', ['AXQ38226.1']), ('translation', ['MIWHVQNENFILDSTRIFMKAFHLLLFDGSLIFPECILIFGLILLLMIDSTSDQKDIPWLYFISSTSLVMSITALLFRWREEPMISFSGNFQTNNFNEIFQFLILLCSTLCIPLSVEYIECTEMAITEFLLFVLTATLGGMFLCGANDLITIFVAPECFSLCSYLLSGYTKKDVRSNEATMKYLLMGGASSSILVHGFSWLYGSSGGEIELQEIVNGLINTQMYNSPGISIALIFITVGIGFKLSPAPSHQWTPDVYEGSPTPVVAFLSVTSKVAASASATRIFDIPFYFSSNEWHLLLEILAILSMILGNLIAITQTSMKRMLAYSSIGQIGYVIIGIIVGDSNDGYASMITYMLFYISMNLGTFACIVLFGLRTGTDNIRDYAGLYTKDPFLALSLALCLLSLGGLPPLAGFFGKLYLFWCGWQAGLYFLVLIGLLTSVVSIYYYLKIIKLLMTGRNQEITPHVRNYRRSPLRSNNSIELSMIVCVIASTIPGISMNPIIAIAQDSLF'])])\n",
      "OrderedDict([('gene', ['ndhB'])])\n",
      "OrderedDict([('gene', ['ndhB'])])\n",
      "OrderedDict([('gene', ['ndhB'])])\n",
      "OrderedDict([('gene', ['rps7'])])\n",
      "OrderedDict([('gene', ['rps7']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein S7']), ('protein_id', ['AXQ38263.1']), ('translation', ['MSRRGTAEKKTAKSDPIYRNRLVNMLVNRILKHGKKSLAYQIIYRAVKKIQQKTETNPLSVLRQAIRGVTPDITVKARRVGGSTHQVPIEIGSTQGKALAIRWLLAASRKRPGRNMAFKLSSELVDAAKGSGDAIRKKEETHRMAEANRAFAHFR'])])\n",
      "OrderedDict([('gene', ['rps12'])])\n",
      "OrderedDict([('gene', ['rps12'])])\n",
      "OrderedDict([('gene', ['rps12'])])\n",
      "OrderedDict([('gene', ['trnV-GAC'])])\n",
      "OrderedDict([('gene', ['trnV-GAC']), ('product', ['tRNA-Val'])])\n",
      "OrderedDict([('gene', ['rrn16'])])\n",
      "OrderedDict([('gene', ['rrn16']), ('product', ['16S ribosomal RNA'])])\n",
      "OrderedDict([('gene', ['trnI-GAU'])])\n",
      "OrderedDict([('gene', ['trnI-GAU']), ('product', ['tRNA-Ile'])])\n",
      "OrderedDict([('gene', ['trnI-GAU'])])\n",
      "OrderedDict([('gene', ['trnI-GAU'])])\n",
      "OrderedDict([('gene', ['trnI-GAU'])])\n",
      "OrderedDict([('gene', ['trnA-UGC'])])\n",
      "OrderedDict([('gene', ['trnA-UGC']), ('product', ['tRNA-Ala'])])\n",
      "OrderedDict([('gene', ['trnA-UGC'])])\n",
      "OrderedDict([('gene', ['trnA-UGC'])])\n",
      "OrderedDict([('gene', ['trnA-UGC'])])\n",
      "OrderedDict([('gene', ['rrn23'])])\n",
      "OrderedDict([('gene', ['rrn23']), ('product', ['23S ribosomal RNA'])])\n",
      "OrderedDict([('gene', ['rrn4.5'])])\n",
      "OrderedDict([('gene', ['rrn4.5']), ('product', ['4.5S ribosomal RNA'])])\n",
      "OrderedDict([('gene', ['rrn5'])])\n",
      "OrderedDict([('gene', ['rrn5']), ('product', ['5S ribosomal RNA'])])\n",
      "OrderedDict([('gene', ['trnR-ACG'])])\n",
      "OrderedDict([('gene', ['trnR-ACG']), ('product', ['tRNA-Arg'])])\n",
      "OrderedDict([('gene', ['trnN-GUU'])])\n",
      "OrderedDict([('gene', ['trnN-GUU']), ('product', ['tRNA-Asn'])])\n",
      "OrderedDict([('gene', ['ycf1']), ('pseudo', [''])])\n",
      "OrderedDict([('gene', ['ndhF'])])\n",
      "OrderedDict([('gene', ['ndhF']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['NADH-plastoquinone oxidoreductase subunit 5']), ('protein_id', ['AXQ38225.1']), ('translation', ['MEFMEQTYEYAWIIPFIPLPVPMLIGAGLILFPTATKSFRRMWAFQSVLLLSIVMIFSISLSIQQINSSSFYQYGWSWIINNDFSLDFGYLIDPLTSIMSILITTVGIMVLIYSDNYMAHDQGYLRFFVYMSFFSTSMLGLVTSSNLIQIYIFWELVGLCSYLLIGFWFTRPVAANACQKAFVTNRVGDFGLLLGILGFYWITGSFEFRDLFEIFTNLIYNNEVNFLFITLCAVLLFAGAVAKSAQFPLHIWLPDAMEGPTPISALIHAATMVAAGIFLVARLLPLFRVIPYIMYLISVIGIITVLLGATLALAQQDIKRGLAYSTMSQLGYMMLALGMGSYRSALFHLITHAYSKALLFLGSGSIIHSMETIVGYSPAKSQNMGLMGGLRKHVPITKITFLLGTLSLCGIPPLACFWSKDEILNDSWLYSPIFAIIAWATAGLTAFYMFRIYLLTFEGHLNVHFQNYRGKQKTPFYSLSLWGKNGIKKNSCLLPMNNNESTYFFSKTKYPIDKNERKMTRPFMTIAHFEHKAVFSYPYESDNTMLFPIFVLGLFTLFVGAIGIPFNQEGVNLDILSKWLASSINLLHPKSNNSMDWNEFFKDAVLSVSIAYFGIFIASFLYKPIYSSLKNLEFINSFVKKGPKRILWDKILNGIYDWSYNRAYIDAFYTRFFVGGIRGLAEFTHFFDRRVIDGMTNGVGVISFILGEGIKYIGGGRISSYLFLYLAYVSVFLLVYYLLFSTF'])])\n",
      "OrderedDict([('gene', ['rpl32'])])\n",
      "OrderedDict([('gene', ['rpl32']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein L32']), ('protein_id', ['AXQ38284.1']), ('translation', ['MAVPKKRTSTSKKRIRKNIWKGKGYWIALKAFSLAKSLSTGNSKSFFVRQTQINK'])])\n",
      "OrderedDict([('gene', ['sprA'])])\n",
      "OrderedDict([('ncRNA_class', ['other']), ('gene', ['sprA']), ('product', ['chloroplast small RNA']), ('note', ['possible role in the maturation of chloroplast 16S rRNA'])])\n",
      "OrderedDict([('gene', ['trnL-UAG'])])\n",
      "OrderedDict([('gene', ['trnL-UAG']), ('product', ['tRNA-Leu'])])\n",
      "OrderedDict([('gene', ['ccsA'])])\n",
      "OrderedDict([('gene', ['ccsA']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['cytochrome c heme attachment protein']), ('protein_id', ['AXQ38251.1']), ('translation', ['MIFSTLEHILTHISFSIVSIVITIHLITLLVEEIVKLYDSSEKGIIVTFFCITGLLVTRWISSGHFPLSDLYESLIFLSWSFSLIHIIPYFKKNVFVLSKITGPSAILTQGFATSGILTEIHQSGILVPALQSEWLIMHVSMMILGYAALLCGSLLSVALLVITFRKNRKLFYIYKSNGFLNASFFVGENVFENTSFFSAKNYYRFQLIQQLDYWSYRVISLGFTFLTIGILSGAVWANEAWGSYWNWDPKETWAFITWIVFAIYLHTRTNRNLRGANSAIVASIGFLIIWICYFGVNLLGIGLHSYGSFPSTFN'])])\n",
      "OrderedDict([('gene', ['ndhD'])])\n",
      "OrderedDict([('gene', ['ndhD']), ('exception', ['RNA editing']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['NADH-plastoquinone oxidoreductase subunit 4']), ('protein_id', ['AXQ38235.1']), ('translation', ['TNYFPWLTIIVVFPIFAGSLIFFLPHKGNRVIRWYTICICILELLLTTYAFCYHFQLDDPLIQLVEDYKWVDFFDFHWRLGIDGLSIGPILLTGFITTLATLAAWPVTRDSRLFHFLMLAMYSGQIGSFSSRDLLLFFIMWELELIPVYLLLAMWGGKKRLYSATKFILYTAGGSIFLLMGVLGVALYGSNEPTLNFETSVNQSYPVILEIIFYIGFFIAFAVKSPIIPLHTWLPDTHGEAHYSTCMLLAGILLKMGAYGLIRINMELLPHAHSIFSPWLMIIGTIQIIYAASTSLGQRNLKKRIAYSSVSHMGFIIIGISSLTDTGLNGALLQIISHGFIGAALFFLAGTTYDRIRLVYLDEMGGIAIPMPKMFTMFSSFSMASLALPGMSGFVAELLVFFGIITGQKYLLIPKLLITFVMAIGIILTPIYSLSMPRQMFYGYKLFNAPKDSFFDSGPRELFLSISIFLPVIGIGIYPDFVLSLAVDKVEVILSNFF'])])\n",
      "OrderedDict([('gene', ['psaC'])])\n",
      "OrderedDict([('gene', ['psaC']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['photosystem I iron-sulfur center']), ('protein_id', ['AXQ38279.1']), ('translation', ['MSHSVKIYDTCIGCTQCVRACPTDVLEMIPWDGCKAKQIASAPRTEDCVGCKRCESACPTDFLSVRVYLWHETTRSMGLAY'])])\n",
      "OrderedDict([('gene', ['ndhE'])])\n",
      "OrderedDict([('gene', ['ndhE']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['NADH-plastoquinone oxidoreductase subunit 4L']), ('protein_id', ['AXQ38272.1']), ('translation', ['MILEHVLVLSAYLFSIGIYGLITSRNMVRALMCLELLLNAVNINFVTFSDFFDNRQLKGDIFSIFVIAIAAAEAAIGLAIVSSIYRNRKSTRINQSNLLNN'])])\n",
      "OrderedDict([('gene', ['ndhG'])])\n",
      "OrderedDict([('gene', ['ndhG']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['NADH-plastoquinone oxidoreductase subunit 6']), ('protein_id', ['AXQ38259.1']), ('translation', ['MDLSEPIHDFLLVFLGSGLILGGLGVVLLPNPIYSAFSLGLVLVCTSLVYILSNSYFVAAAQLLIYVGAINVLIIFAVMFMNGSEYYKDFHLWTVGDGITSMVCISLFISLITIISDTSWYGIIWTTRSNQIIEQDFLSNSQQIGIHLSTDFFLPFELISIILLVALIGAIAVARQ'])])\n",
      "OrderedDict([('gene', ['ndhI'])])\n",
      "OrderedDict([('gene', ['ndhI']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['NADH-plastoquinone oxidoreductase subunit I']), ('protein_id', ['AXQ38260.1']), ('translation', ['MLPMITEFLNYGQQTIRAARYIGQGFIITLSHANRLPVTIQYPYEKLITSERFRGRIHFEFDKCIACEVCVRVCPIDLPVVDWKLETDIRKKRLLNYSIDFGICIFCGNCVEYCPTNCLSMTEEYELSTYDRHELNYNQIALGRLPMSVIDDYTIRTILNLPQIKNE'])])\n",
      "OrderedDict([('gene', ['ndhA'])])\n",
      "OrderedDict([('gene', ['ndhA']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['NADH-plastoquinone oxidoreductase subunit 1']), ('protein_id', ['AXQ38223.1']), ('translation', ['MIIDTTEMETINSFSKLESLKEVYGIIWILVPILTLVLGITIGVLVIVWLEREISAGIQQRIGPEYAGPLGILQALADGTKLLLKENLIPSTGDTRLFSIGPSIAVISIFLSYSVIPFGDHLVLADLSIGVFFWIAISSIAPVGLLMSGYGSNNKYSFLGGLRAAAQSISYEIPLALCVLSISLLSNSSSTVDIVEAQSKYGFWGWNLWRQPIGFIVFLISSLAECERLPFDLPEAEEELVAGYQTEYSGIKFGLFYIASYLNLLVSSLFVTVLYLGGWNLSIPYILVPELFGINKRGKVFGRLIGIFITLAKTYLFLFIPIATRWTLPRLRMDQLLNLGWKFLLPISLGNLLLTTSFQLLSL'])])\n",
      "OrderedDict([('gene', ['ndhA'])])\n",
      "OrderedDict([('gene', ['ndhA'])])\n",
      "OrderedDict([('gene', ['ndhA'])])\n",
      "OrderedDict([('gene', ['ndhH'])])\n",
      "OrderedDict([('gene', ['ndhH']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['NADH-plastoquinone oxidoreductase subunit 7']), ('protein_id', ['AXQ38245.1']), ('translation', ['MTAPTTRKDLMIVNMGPQHPSMHGVLRLIVTLDGEDVVDCEPILGYLHRGMEKIAENRTIIQYLPYVTRWDYLATMFTEAITINGPEQLGNIQVPKRASYIRVIMLELSRIASHLLWLGPFMADIGAQTPFFYIFRERELIYDLFEAATGMRMMHNYFRIGGVAADLPYGWIDKCLDFCDYFLTGVAEYQKLITRNPIFLERVEGVGIIGRDEAVNWGLSGPMLRASGIEWDLRKVDHYESYDEFDWQVQWQREGDSLARYLVRISEMTESIKIIQQALEGIPGGPYENLEIRRFDRLKDPEWNDFEYRFISKKPSPTFELSKQELYVRVEAPKGELGIFLIGDQSVFPWRWKIRPPGFINLQILPQLVKRMKLADIMTILGSIDIIMGEVDR'])])\n",
      "OrderedDict([('gene', ['rps15'])])\n",
      "OrderedDict([('gene', ['rps15']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein S15']), ('protein_id', ['AXQ38276.1']), ('translation', ['MVKNYVISVTSQEEKKGSVEFQVFNFTNKIQRLTSHLELHKKDYLSQRGLKKILGKRQRLLAYLAKKNRVRYKELINKLAIRETKTR'])])\n",
      "OrderedDict([('gene', ['ycf1'])])\n",
      "OrderedDict([('gene', ['ycf1']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['hypothetical chloroplast RF1']), ('protein_id', ['AXQ38219.1']), ('translation', ['MIFQSFLLGNLVSLCMKIINSVVVVGLYYGFLTTFSIGPSYLFLLRALVMEEGTEKKVSATTGFITGQLMMFISIYYAPLHLALGRPHTITVLALPYLLFHFFWNNHKHFFDYGSTTRNSMRNLSIQCVFLNNLIFQLFNHFILPSSMLARLVNIYLFRCNNKILFVTSGFVGWLIGHILFMKWLGLVLVWIRQNHSIRSNKYIRSNKYLVLELRNSMARIFSILLFITCVYYLGRIPSPILTKKLNKKLKEASKTEEGVESTGYNNINTSNSPIYDYQDSYLNNNNTGNLENFQLQLLDKKNENQEFDLLGFQKPFVSLLFDYNRWNRPFRYIKNDRFEQAVRTEMSQYFFDTCKSDGKQRISFTYPPSLSTFWKMIKRKIPLLSLQKTFPNELDNQWVSTNKEKSNTLNKEFVNRLEILDKESLSLDILETQTRLCNDDTKKEYMPRVYDPLLNGPYRGTIKKGVSSSIIKNTLLENRKEKVRLNRIHTIFLPNTDYQEFEQKAYTVDKKPLSTENDEFLTLINEFGNEGNEPKFSLNLKGLSLFSDQEQERVNSEKRTRFVKFVFKAIDPNEIKSGKQSIGIKEIGKKVSRWSHKLITELDQQMGEFQDRISMDHQLRLRKSKRVVVFTDKTNNNKNNATKDPDEEELVLISYSQLSDFRRGIIKGSMRAQRRKTSISKLFQANVHSPLFVDRITPFRLFSFDISELIKPIFRNWTGKGGKFKILASREEQTKREETKEKDEKEANKKKEETRIMIAEAWDSIILAQMVRGYMLITQSILRKYIILPSLIIAKNIGRMLFLQLPEWSEDFQEWNREMQIKCTYNGVQLSETEFPKNWLRDGIQIKILFPFCLKPWHISKLYPSRRELMKKKKQKDDFCFLTVWGMETELPFGSPRKRPSFFEPIFKELEKKTAKFKKKYFRTLKILKRKTKLFRRVFKEAKKWLIKSIGFLKKIKKELSKINLFILFRFKEISESNETEKDKDYLISNQIINESFSQIESGNRLNSSLIETKMKDLTDRTSTIKNQIERITKEKKKVTPEIDIKPNTNKTNNIKKFESPKNFFQKLKNRNTRVIRKFHYFRKLFIPRLYIDLFLSIINIPRITTQLFLESTNKLIEKSISNNEINQEKINNNKMHFIFISTIKKSVYNISKQNSHIFCDLSYLSQAYVFYKLSQTQVINLSKFKSVLQYNTMSCFLTTKIKDYFRTLGIFHSELKHKKLQSYRINQWKNWLRRHYQYDLSQIRWSRLMPKKWRTKVNQSCMTQNKNRNLNKWNSYEKDQLSHYKKENDLKLYSLSNEKDNFKKCYRYDLLASKSINYENKNDSFISRLPFQVKKNFEISYNFNPSKHNFVDMSGNLHINNYLRKGNILDIERNLDRKYFDWKILHFSLRQKGDIEAWVKIDTNSNPNTQIGINNYKIIDNYQIIEKRGVFYLTTHQNPEKTQKNSKKVFFDWMGMNEKIFNCPILNLEFWFFPEFLLLSNVYKIKPWIIPSKFLLLNLNTNEKVSQNQKQNFFIPSNKKKQNRIQEAKEPASQREHGSSIENKENRGPVFSKHQNDLEKDYVESDTKKGKTKKQYKSNTEAELDLFLQRYLLFQLKWKDALNQSLLENIKVYCLLLRLTNPTKITISSIQRREMCLDTMLIQANLPLTDLLKKGVLIIEPIRFSVKNNGEFLMYQTIGISLVHKSKHQTNQRYREQRYVAKKNFDEFILQPQTQRINTEQNHFDLLVPENILWSRRRRELRIRSFFNSLNWDGVDRNSVFCNETNVKNWSQFLDERKPFYKEKNELIKLKFFLWPNYRLEDLACMNRYWFDTNNGSRFNILRIHMYPRLKIN'])])\n",
      "OrderedDict([('note', ['inverted repeat A (IRa)']), ('rpt_type', ['inverted'])])\n",
      "OrderedDict([('gene', ['trnN-GUU'])])\n",
      "OrderedDict([('gene', ['trnN-GUU']), ('product', ['tRNA-Asn'])])\n",
      "OrderedDict([('gene', ['trnR-ACG'])])\n",
      "OrderedDict([('gene', ['trnR-ACG']), ('product', ['tRNA-Arg'])])\n",
      "OrderedDict([('gene', ['rrn5'])])\n",
      "OrderedDict([('gene', ['rrn5']), ('product', ['5S ribosomal RNA'])])\n",
      "OrderedDict([('gene', ['rrn4.5'])])\n",
      "OrderedDict([('gene', ['rrn4.5']), ('product', ['4.5S ribosomal RNA'])])\n",
      "OrderedDict([('gene', ['rrn23'])])\n",
      "OrderedDict([('gene', ['rrn23']), ('product', ['23S ribosomal RNA'])])\n",
      "OrderedDict([('gene', ['trnA-UGC'])])\n",
      "OrderedDict([('gene', ['trnA-UGC']), ('product', ['tRNA-Ala'])])\n",
      "OrderedDict([('gene', ['trnA-UGC'])])\n",
      "OrderedDict([('gene', ['trnA-UGC'])])\n",
      "OrderedDict([('gene', ['trnA-UGC'])])\n",
      "OrderedDict([('gene', ['trnI-GAU'])])\n",
      "OrderedDict([('gene', ['trnI-GAU']), ('product', ['tRNA-Ile'])])\n",
      "OrderedDict([('gene', ['trnI-GAU'])])\n",
      "OrderedDict([('gene', ['trnI-GAU'])])\n",
      "OrderedDict([('gene', ['trnI-GAU'])])\n",
      "OrderedDict([('gene', ['rrn16'])])\n",
      "OrderedDict([('gene', ['rrn16']), ('product', ['16S ribosomal RNA'])])\n",
      "OrderedDict([('gene', ['trnV-GAC'])])\n",
      "OrderedDict([('gene', ['trnV-GAC']), ('product', ['tRNA-Val'])])\n",
      "OrderedDict([('gene', ['rps12'])])\n",
      "OrderedDict([('gene', ['rps12'])])\n",
      "OrderedDict([('gene', ['rps12'])])\n",
      "OrderedDict([('gene', ['rps7'])])\n",
      "OrderedDict([('gene', ['rps7']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein S7']), ('protein_id', ['AXQ38264.1']), ('translation', ['MSRRGTAEKKTAKSDPIYRNRLVNMLVNRILKHGKKSLAYQIIYRAVKKIQQKTETNPLSVLRQAIRGVTPDITVKARRVGGSTHQVPIEIGSTQGKALAIRWLLAASRKRPGRNMAFKLSSELVDAAKGSGDAIRKKEETHRMAEANRAFAHFR'])])\n",
      "OrderedDict([('gene', ['ndhB'])])\n",
      "OrderedDict([('gene', ['ndhB']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['NADH-plastoquinone oxidoreductase subunit 2']), ('protein_id', ['AXQ38227.1']), ('translation', ['MIWHVQNENFILDSTRIFMKAFHLLLFDGSLIFPECILIFGLILLLMIDSTSDQKDIPWLYFISSTSLVMSITALLFRWREEPMISFSGNFQTNNFNEIFQFLILLCSTLCIPLSVEYIECTEMAITEFLLFVLTATLGGMFLCGANDLITIFVAPECFSLCSYLLSGYTKKDVRSNEATMKYLLMGGASSSILVHGFSWLYGSSGGEIELQEIVNGLINTQMYNSPGISIALIFITVGIGFKLSPAPSHQWTPDVYEGSPTPVVAFLSVTSKVAASASATRIFDIPFYFSSNEWHLLLEILAILSMILGNLIAITQTSMKRMLAYSSIGQIGYVIIGIIVGDSNDGYASMITYMLFYISMNLGTFACIVLFGLRTGTDNIRDYAGLYTKDPFLALSLALCLLSLGGLPPLAGFFGKLYLFWCGWQAGLYFLVLIGLLTSVVSIYYYLKIIKLLMTGRNQEITPHVRNYRRSPLRSNNSIELSMIVCVIASTIPGISMNPIIAIAQDSLF'])])\n",
      "OrderedDict([('gene', ['ndhB'])])\n",
      "OrderedDict([('gene', ['ndhB'])])\n",
      "OrderedDict([('gene', ['ndhB'])])\n",
      "OrderedDict([('gene', ['trnL-CAA'])])\n",
      "OrderedDict([('gene', ['trnL-CAA']), ('product', ['tRNA-Leu'])])\n",
      "OrderedDict([('gene', ['ycf15']), ('pseudo', [''])])\n",
      "OrderedDict([('gene', ['ycf2'])])\n",
      "OrderedDict([('gene', ['ycf2']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['hypothetical chloroplast RF2']), ('protein_id', ['AXQ38218.1']), ('translation', ['MRGHQFKSWIFELREILREIKNSHHFLDSWTQFNSVGSFIHIFFHQERFLKLFDPRIWSILLSRNSQGSTSNRYFTIKGVILFVVAVLIYRINNRNMVERKNLYLIGLLPIPMHSIGPRNDTLEESVGSSNINRLIVSLLYLPKGKKISESCFLNPKESTWVLPITKKCSMPESNWGSRWWRNWIGKKRDFSCKISNETVAGIEILFKEKDLKYLEFLFVYYMDDPIRKDHDWELFDRLSLRKSRNRINLNSGPLFEILVKHWISYLMSAFREKIPIEVEGFFKQQGAGSTIQSNDIEHVSHLFSRNKWAITLQNCAQFHMWQFRQDLFVSWGKNPPESDFLRNVSRENWIWLDNVWLVNKDRFFSKVQNVSSNIQYDSTRSSFVQVTDSSQLKGSSDQSRDHLDSISNEDSEYHTLINQREIQQRKERSILWDPSFLQTERKEIESGRFPKCLSGYSSMSRLFTEREKQMINHLFPEEIEEFLGNPTRSVRSFFSDRWSELHLGSNPTERSTRDQKLLKKQQDLSFVPSRRSENKEMVNIFKIITYLQNTVSIHPISSDPGCDMVPKDEPDMDSSNKISFLNKNPFFDLFHLFHDRNRGGYTLHYDFASEERFQEMADLFTLSITEPDLVYHKGFAFSIDSCGLDQKQFLNEARDESKKKSLLVLPPIFYEENESFSRRIRKKWVRISCGNDLEDPKPKIVVFASNNIMEAVTQYRLIRNLIQIQYSTYGYIRNVLNRFFLMNRSDRNFEYGIQRDQIGKDTLNHRTIMKYTINQYLSNLKKSQKKWFEPLILISRTERSMNRDPDAYRYKWSNGSKNFQEHLEQSVSEQKSRFQVVFDRLRINQYSIDWSEVIDKKDLSKSLRFFLSKSLRFFLSKSLLFLSKLLLFLSNSLPFFCVSFGNIPIHRSEIYIYELKGPNDQLCNQLLESIGLQIVHLKKLKPFLLDDHDTSQKSKFLINGGTISPFLFNKIPKWMIDSFHTRNNRRKSFDSPDSYFSMIFHDQDNWLNPVKPFHRSSLISSFYKANRLRFLNNPHHFCFYWNTRFPFSVEKARINNSDFTYGQFLNILFIRNKIFSLCVGKKKHAFWGRDTISPIESQVSHIFIPNDFPQSGDETYNLYKSFHFPSRSDPSVRRAIYSIADISGTPLTEGQIVNFERTYCQPLSDMNLSDSEGKNLHQYLNFNSNMGLIHTPCSEKDLSSEKRKKWSLCLKKCVEKGQMYRTFQRDSAFSTLSKWNLFQTYMPWFLTSTGYKYLNLIFLDTFSDLLPILSSSQKFVSIFPDIMHGSGISWRILQKKLCLPQWNLISEISSKCLHNLLLSEEMIHRNNESPLISTHLRSPNAREFLYSILFLLLVAGYLVRTHLLFVSRASSELQTEFEKVKSLMIPSSMIELRKLLDRYPTSEPNSFWLKNLFLVALEQLGDSLEEIRGSASGGNMLGPAYGVKSIRSKKKDWNINLIEIIDLIPNPINRITFSRNTRHLSHTSKEIYSLIRKRKNVNGDWIDEKIESWVANSDSIDDEEKEFLVHFSTLTTENRIDQILLSLTHSDHLSKNDSGYQMIEQPGAIYLRYLVDIHKKHLMNYEFNPSCLAERRIFLAHYQTITYSQTSCGENSFHFPSHGKPFSLRLALSPSRGILVIGSIGTGRSYLVKYLATNSYVPFITVFLNKFLDNKSKGFLLDEIDIDDSDDIDDSDNLDASDDIDRDLDTELELLTRMNGLTVDMMPEIDRFYITLQFELAKAMSPCIIWIPNIHDLDVNESNDLSLGLLVNHLSRDCERCSTRNILVIASTHIPQKVDPALIAPNKLNTCIKIRRLLIPQQRKHFFTLSYTRGFHLEKKMFHTNGFGSITMGSNARDLVALTNEVLSISITQKKSIIDTNTIRSALHRQTWDLRSQVRSVQDHGILFYQIGRAVAQNVLLSNCPIEPISIYMKKKSCNGGDSYLYKWYFELGTSMKRLTILLYLLSCSAGSVAQDLWSLSVPDEKNGITSYGLVENDSDLVHGLLEVEGALVGSSRTEKDCSQFDNDRVTLLLRPEPRNPLDMMQKGSWSILDQRFLYEKYESEFEEGEGEGALDPQEDLFNHIVWAPRIWRPWGFLFDCIERPNELGFPYWSRSFRGKRIIYDEEDELQENDSGFLQSGTMQYQTRDRSSKEQGLFRISQFIWDPADPLFFLFKDQPPGSVFSHRELFADEEMSKGLLTSQTDPPTSIYKRWFIKNTQEKHFELLINRQRWLRTNSSLSNGSFRSNTLSESYQYLSNLFLSNGTLLDQMTKTLLRKRWLFPDEMKIGFM'])])\n",
      "OrderedDict([('gene', ['trnI-CAU'])])\n",
      "OrderedDict([('gene', ['trnI-CAU']), ('product', ['tRNA-Ile'])])\n",
      "OrderedDict([('gene', ['rpl23'])])\n",
      "OrderedDict([('gene', ['rpl23']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein L23']), ('protein_id', ['AXQ38275.1']), ('translation', ['MDGIKYAVFTDKSIRLLGKNQYTSNVESGSTRTEIKHWVELFFGVKVIAMNSHRLPGKSRRMGPIMGHTMHYRRMIITLQPGYSIPPLRKKRT'])])\n",
      "OrderedDict([('gene', ['rpl2'])])\n",
      "OrderedDict([('gene', ['rpl2']), ('codon_start', ['1']), ('transl_table', ['11']), ('product', ['ribosomal protein L2']), ('protein_id', ['AXQ38237.1']), ('translation', ['MAIHLYKTSTPSTRNETVDSQVKSNPRNNLIYGQRRCGKGRNARGIITARHRGGGHKRLYRKIDFRRNEKDIYGRIVTIEYDPNRNAYICLIHYGDGEKRYILHPRGAIIGDTIVSGTEVPIKMGNALPLSTDMPLGTAIHNIEITLGKGGQLARAAGAVAKLIAKEGKSATLKLPSGEVRLISKNCSATVGQVGNVGVNQKSLGRAGSKRWLGKRPVVRGVVMNPVDHPHGGGEGRAPIGRKKPTTPWGYPALGRRSRKRNKYSDNLILRRRSK'])])\n",
      "OrderedDict([('gene', ['rpl2'])])\n",
      "OrderedDict([('gene', ['rpl2'])])\n",
      "OrderedDict([('gene', ['rpl2'])])\n",
      "OrderedDict([('gene', ['rps19']), ('pseudo', [''])])\n"
     ]
    }
   ],
   "source": [
    "for ft in rec.features[1:]:\n",
    "    print(ft.qualifiers)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 87,
   "metadata": {},
   "outputs": [],
   "source": [
    "import seaborn as sns\n",
    "sns.set()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 88,
   "metadata": {},
   "outputs": [],
   "source": [
    "exercise = sns.load_dataset(\"exercise\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 89,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "pandas.core.frame.DataFrame"
      ]
     },
     "execution_count": 89,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "type(exercise)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 90,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Unnamed: 0</th>\n",
       "      <th>id</th>\n",
       "      <th>diet</th>\n",
       "      <th>pulse</th>\n",
       "      <th>time</th>\n",
       "      <th>kind</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>low fat</td>\n",
       "      <td>85</td>\n",
       "      <td>1 min</td>\n",
       "      <td>rest</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>low fat</td>\n",
       "      <td>85</td>\n",
       "      <td>15 min</td>\n",
       "      <td>rest</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>low fat</td>\n",
       "      <td>88</td>\n",
       "      <td>30 min</td>\n",
       "      <td>rest</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>low fat</td>\n",
       "      <td>90</td>\n",
       "      <td>1 min</td>\n",
       "      <td>rest</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>2</td>\n",
       "      <td>low fat</td>\n",
       "      <td>92</td>\n",
       "      <td>15 min</td>\n",
       "      <td>rest</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Unnamed: 0  id     diet  pulse    time  kind\n",
       "0           0   1  low fat     85   1 min  rest\n",
       "1           1   1  low fat     85  15 min  rest\n",
       "2           2   1  low fat     88  30 min  rest\n",
       "3           3   2  low fat     90   1 min  rest\n",
       "4           4   2  low fat     92  15 min  rest"
      ]
     },
     "execution_count": 90,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "exercise.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [],
   "source": [
    "import itertools as it"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Con Biopython\n",
    "def count_bp(archivo,formato):\n",
    "    from Bio import SeqIO\n",
    "    records = SeqIO.parse(archivo,formato)\n",
    "    return sum(len(read) for read in records)\n",
    "def count_reads(archivo,formato):\n",
    "    from Bio import SeqIO\n",
    "    records = SeqIO.parse(archivo,formato)\n",
    "    return sum(1 for read in records)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Cantidad de información: 98574 bp\n",
      "Número de reads: 1000\n",
      "Longitud promedio: 98.574 bp\n"
     ]
    }
   ],
   "source": [
    "bases = count_bp(\"../../Biopython/reads.fastq\", \"fastq\")\n",
    "nreads = count_reads(\"../../Biopython/reads.fastq\", \"fastq\")\n",
    "print(f'Cantidad de información: {bases} bp')\n",
    "print(f'Número de reads: {nreads}')\n",
    "print(f'Longitud promedio: {bases/nreads} bp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Sin Biopython\n",
    "def count_bp2(archivo,step=4):\n",
    "    import itertools as it\n",
    "    f = open(archivo,\"r\")\n",
    "    bases = sum(len(line.strip(\"\\n\")) for line in it.islice(f,1,None,step))\n",
    "    f.close()\n",
    "    return bases\n",
    "def count_reads2(archivo,step=4):\n",
    "    import itertools as it\n",
    "    f = open(archivo,\"r\")\n",
    "    nreads = sum(1 for line in it.islice(f,1,None,step))\n",
    "    f.close()\n",
    "    return nreads\n",
    "\n",
    "def len_prom(archivo,step=4):\n",
    "    import itertools as it\n",
    "    with open(archivo,\"r\") as f:\n",
    "        bases = ((len(line.strip(\"\\n\")),1) for line in it.islice(f,1,None,step))\n",
    "    return bases"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 74,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Cantidad de información: 98574 bp\n",
      "Número de reads: 1000\n",
      "Longitud promedio: 98.574 bp\n"
     ]
    }
   ],
   "source": [
    "bases = count_bp2(\"../../Biopython/reads.fastq\")\n",
    "nreads = count_reads2(\"../../Biopython/reads.fastq\")\n",
    "print(f'Cantidad de información: {bases} bp')\n",
    "print(f'Número de reads: {nreads}')\n",
    "print(f'Longitud promedio: {bases/nreads} bp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 146,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Cantidad de información: 0.098574 Mbp\n",
      "Número de reads: 1000\n",
      "Longitud promedio: 98.574 bp\n",
      "CPU times: user 4.21 ms, sys: 127 µs, total: 4.33 ms\n",
      "Wall time: 3.58 ms\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "import itertools as it\n",
    "import os\n",
    "\n",
    "def generate_filenames(folder):\n",
    "    \"\"\"\n",
    "    generates a sequence of opened files\n",
    "    matching a specific extension\n",
    "    \"\"\"\n",
    "    for dir_path, dir_names, file_names in os.walk(folder):\n",
    "        for file_name in file_names:\n",
    "            if file_name.endswith('.fastq'):\n",
    "                yield open(os.path.join(dir_path, file_name))\n",
    "\n",
    "def cat_files(files):\n",
    "    \"\"\"\n",
    "    takes in an iterable of filenames\n",
    "    \"\"\"\n",
    "    for fname in files:\n",
    "        for line in fname:\n",
    "            yield line\n",
    "            \n",
    "def count_bp3(lines,step=4):\n",
    "    bases = sum(len(line.strip(\"\\n\")) for line in it.islice(lines,1,None,step))\n",
    "    return bases\n",
    "\n",
    "def count_reads3(lines,step=4):\n",
    "    nreads = sum(1 for line in it.islice(lines,1,None,step))\n",
    "    return nreads\n",
    "\n",
    "folder = \"../../Biopython\"\n",
    "\n",
    "bases = count_bp3(cat_files(generate_filenames(folder)))\n",
    "nreads = count_reads3(cat_files(generate_filenames(folder)))\n",
    "print(f'Cantidad de información: {bases/1000000} Mbp')\n",
    "print(f'Número de reads: {nreads}')\n",
    "print(f'Longitud promedio: {bases/nreads} bp')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 164,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Cantidad de información: 0.098574 Mbp\n",
      "Número de reads: 1000\n",
      "Longitud promedio: 98.574 bp\n",
      "CPU times: user 4.22 ms, sys: 76 µs, total: 4.29 ms\n",
      "Wall time: 3.31 ms\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "import itertools as it\n",
    "import os\n",
    "\n",
    "def generate_filenames(folder):\n",
    "    \"\"\"\n",
    "    generates a sequence of opened files\n",
    "    matching a specific extension\n",
    "    \"\"\"\n",
    "    for dir_path, dir_names, file_names in os.walk(folder):\n",
    "        for file_name in file_names:\n",
    "            if file_name.endswith('.fastq'):\n",
    "                yield open(os.path.join(dir_path, file_name))\n",
    "\n",
    "def cat_files(files):\n",
    "    \"\"\"\n",
    "    takes in an iterable of filenames\n",
    "    \"\"\"\n",
    "    for fname in files:\n",
    "        for line in fname:\n",
    "            yield line\n",
    "\n",
    "def reads_stats(lines,step=4):\n",
    "    len_read = lambda read: len(read.rstrip(\"\\n\"))\n",
    "    road = ((len_read(line),1) for line in it.islice(lines,1,None,step))\n",
    "    bases, nreads = map(sum, zip(*road))\n",
    "    return bases, nreads\n",
    "\n",
    "folder = \"../../Biopython\"\n",
    "\n",
    "bases,nreads = reads_stats(cat_files(generate_filenames(folder)))\n",
    "print(f'Cantidad de información: {bases/1000000} Mbp')\n",
    "print(f'Número de reads: {nreads}')\n",
    "print(f'Longitud promedio: {bases/nreads} bp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 104,
   "metadata": {},
   "outputs": [],
   "source": [
    "aux = [(i,1) for i in \"ABCDEF\"]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 120,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[('A', 'B', 'C', 'D', 'E', 'F'), (1, 1, 1, 1, 1, 1)]"
      ]
     },
     "execution_count": 120,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "list(zip(*((i,1) for i in \"ABCDEF\")))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 123,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[('B', 1), ('D', 1), ('F', 1)]"
      ]
     },
     "execution_count": 123,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "list(it.islice(((i,1) for i in \"ABCDEF\"),1,None,2))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 71,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "98574\n"
     ]
    }
   ],
   "source": [
    "bases = count_bp2(\"../../Biopython/reads.fastq\")\n",
    "print(bases)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 69,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "58\n",
      "101\n",
      "101\n",
      "91\n",
      "101\n",
      "101\n",
      "101\n",
      "59\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "44\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "48\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "81\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "98\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "43\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "99\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "63\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "59\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "52\n",
      "101\n",
      "101\n",
      "101\n",
      "44\n",
      "101\n",
      "57\n",
      "101\n",
      "101\n",
      "49\n",
      "101\n",
      "70\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "75\n",
      "66\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "93\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "57\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "69\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "86\n",
      "62\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "41\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "99\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "61\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "60\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "63\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "47\n",
      "68\n",
      "101\n",
      "64\n",
      "101\n",
      "101\n",
      "85\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "53\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "64\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "84\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "49\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "69\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "58\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "47\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "53\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "87\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "63\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "50\n",
      "101\n",
      "101\n",
      "100\n",
      "46\n",
      "101\n",
      "62\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "93\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "95\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "87\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "50\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "71\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "94\n",
      "101\n",
      "101\n",
      "101\n",
      "62\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "99\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "40\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "97\n",
      "101\n",
      "77\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "46\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "92\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "71\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "40\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "53\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "36\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "78\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "42\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "57\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "39\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "100\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n",
      "50\n",
      "101\n",
      "101\n",
      "101\n",
      "101\n"
     ]
    }
   ],
   "source": [
    "import itertools as it\n",
    "f = open(\"../../Biopython/reads.fastq\",\"r\")\n",
    "handle =(len(line.strip(\"\\n\")) for line in it.islice(f,1,None,4))\n",
    "for i in handle:\n",
    "    print(i)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 89,
   "metadata": {},
   "outputs": [
    {
     "ename": "TypeError",
     "evalue": "unsupported operand type(s) for +: 'int' and 'tuple'",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mTypeError\u001b[0m                                 Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-89-68154ce4fa5d>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m      2\u001b[0m \u001b[0mtest\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0menumerate\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mstring\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      3\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 4\u001b[0;31m \u001b[0msum\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mit\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mislice\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mtest\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;36m1\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;32mNone\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;36m2\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;31mTypeError\u001b[0m: unsupported operand type(s) for +: 'int' and 'tuple'"
     ]
    }
   ],
   "source": [
    "string = \"ABCDEFGHIJ\"\n",
    "test = enumerate(string)\n",
    "\n",
    "sum(it.islice(test,1,None,2))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 94,
   "metadata": {},
   "outputs": [
    {
     "ename": "TypeError",
     "evalue": "islice expected at most 4 arguments, got 12",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mTypeError\u001b[0m                                 Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-94-38e37e562f01>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m      3\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      4\u001b[0m \u001b[0;31m#print(test)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 5\u001b[0;31m \u001b[0msum\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mit\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mislice\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mtest\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;32mNone\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;36m2\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;31mTypeError\u001b[0m: islice expected at most 4 arguments, got 12"
     ]
    }
   ],
   "source": [
    "string = \"ABCDEFGHIJ\"\n",
    "test = list(enumerate(string))\n",
    "\n",
    "#print(test)\n",
    "sum(it.islice(*test,None,2))"
   ]
  },
  {
   "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.6.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}