sourmash_flu_test.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "from Bio.SeqIO.FastaIO import SimpleFastaParser\n",
    "from Bio import SeqIO\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import sourmash\n",
    "import time"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\r",
      "\u001b[Ksourmash version 2.0.0\r\n",
      "\r",
      "\u001b[K- loaded from path: /home/brian/anaconda3/lib/python3.7/site-packages/sourmash\r\n",
      "\r",
      "\u001b[K\r\n"
     ]
    }
   ],
   "source": [
    "!sourmash info"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "with open('influenza.fna') as fasta_file: \n",
    "    cgi_titles = []\n",
    "    cgi_sequences = []\n",
    "    for title, sequence in SimpleFastaParser(fasta_file):\n",
    "        cgi_titles.append(title)\n",
    "        cgi_sequences.append(sequence)\n",
    "influenza_na = pd.DataFrame(list(zip(cgi_titles, cgi_sequences)), columns=['titles','sequences'])\n",
    "uniflna = influenza_na.drop_duplicates('sequences')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "uniflna_HA = uniflna.loc[uniflna.iloc[:,0].str.contains(r'\\(HA\\) | hemagglutinin | (segment 4)')]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "np.random.seed(12478)\n",
    "uniflna_HA_sub1 = uniflna_HA.sample(n=1000)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "scrolled": true
   },
   "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>titles</th>\n",
       "      <th>sequences</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>252112</th>\n",
       "      <td>gi|469936972|gb|KC738886|Influenza A virus (A/...</td>\n",
       "      <td>ACAGCTACATATGCAGACACAATATGTATAGGCTACCATGCCAACA...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>259962</th>\n",
       "      <td>gi|484849114|gb|KC865637|Influenza A virus (A/...</td>\n",
       "      <td>AAAGCAGGGGATAATTCTATTAACCATGAAGACTATCATTGCTTTG...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>428295</th>\n",
       "      <td>gi|984690389|gb|KU590388|Influenza A virus (A/...</td>\n",
       "      <td>ATGAAGACTATCATTGCTTTGAGCTACATTCTATGTCTGGTTTTCG...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>628621</th>\n",
       "      <td>gi|1409674136|gb|MH540925|Influenza A virus (A...</td>\n",
       "      <td>GGATAATTCTATTAACCATGAAGACTATCATTGCTTTGAGCTACAT...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>244546</th>\n",
       "      <td>gi|451789264|gb|KC535453|Influenza A virus (A/...</td>\n",
       "      <td>ATGAAGACTATCATTGCTTTGAGCTACATTCTATGTCTGGTTTTCG...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                   titles  \\\n",
       "252112  gi|469936972|gb|KC738886|Influenza A virus (A/...   \n",
       "259962  gi|484849114|gb|KC865637|Influenza A virus (A/...   \n",
       "428295  gi|984690389|gb|KU590388|Influenza A virus (A/...   \n",
       "628621  gi|1409674136|gb|MH540925|Influenza A virus (A...   \n",
       "244546  gi|451789264|gb|KC535453|Influenza A virus (A/...   \n",
       "\n",
       "                                                sequences  \n",
       "252112  ACAGCTACATATGCAGACACAATATGTATAGGCTACCATGCCAACA...  \n",
       "259962  AAAGCAGGGGATAATTCTATTAACCATGAAGACTATCATTGCTTTG...  \n",
       "428295  ATGAAGACTATCATTGCTTTGAGCTACATTCTATGTCTGGTTTTCG...  \n",
       "628621  GGATAATTCTATTAACCATGAAGACTATCATTGCTTTGAGCTACAT...  \n",
       "244546  ATGAAGACTATCATTGCTTTGAGCTACATTCTATGTCTGGTTTTCG...  "
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "uniflna_HA_sub1.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1000"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(uniflna_HA_sub1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [],
   "source": [
    "# make it back into text FASTA file\n",
    "with open (\"uniflna_HA_sub1.fasta\", 'w') as output:\n",
    "    for i in range(len(uniflna_HA_sub1)):\n",
    "        output.write(\">\" + list(uniflna_HA_sub1.titles.values)[i] + \"\\n\" + list(uniflna_HA_sub1.sequences.values)[i] + \"\\n\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "\u001b[Kcomputing signatures for files: uniflna_HA_sub1.fasta\n",
      "\u001b[KComputing signature for ksizes: [5]\n",
      "\u001b[KComputing only nucleotide (and not protein) signatures.\n",
      "\u001b[KComputing a total of 1 signature(s).\n",
      "\u001b[Kcalculated 1000 signatures for 1000 sequences in uniflna_HA_sub1.fasta\n",
      "\u001b[Ksaved 1000 signature(s). Note: signature license is CC0.\n",
      "finished in 1.0879552364349365seconds\n"
     ]
    }
   ],
   "source": [
    "start = time.time()\n",
    "!sourmash compute uniflna_HA_sub1.fasta  -k 5 --singleton\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "...sig loading 999_HA_sub1.fasta.sig\n",
      "\u001b[Kloaded 1000 signatures total.                                                  \n",
      "\u001b[K\n",
      "min similarity in matrix: 0.221\n",
      "\u001b[Ksaving labels to: uniflna_HA_sub1_cmp.labels.txt\n",
      "\u001b[Ksaving distance matrix to: uniflna_HA_sub1_cmp\n",
      "finished in 63.475311279296875seconds\n"
     ]
    }
   ],
   "source": [
    "start = time.time()\n",
    "!sourmash compare uniflna_HA_sub1.fasta.sig -o uniflna_HA_sub1_cmp\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "\u001b[Kloading comparison matrix from uniflna_HA_sub1_cmp...\n",
      "\u001b[K...got 1000 x 1000 matrix.\n",
      "\u001b[Kloading labels from uniflna_HA_sub1_cmp.labels.txt\n",
      "\u001b[Ksaving histogram of matrix values => uniflna_HA_sub1_cmp.hist.png\n",
      "\u001b[Kwrote dendrogram to: uniflna_HA_sub1_cmp.dendro.png\n",
      "\u001b[Kwrote numpy distance matrix to: uniflna_HA_sub1_cmp.matrix.png\n",
      "finished in 10.18021273612976seconds\n"
     ]
    }
   ],
   "source": [
    "start = time.time()\n",
    "!sourmash plot uniflna_HA_sub1_cmp\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "tree:  \n",
    "![dendogram](uniflna_HA_sub1_cmp.dendro.png)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "matrix and heatmap:  \n",
    "![matrix](uniflna_HA_sub1_cmp.matrix.png)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "...sig loading 999_HA_sub1.fasta.sig\n",
      "\u001b[Kloaded 1000 signatures total.                                                  \n",
      "\u001b[K\n",
      "min similarity in matrix: 0.221\n",
      "finished in 67.94533371925354seconds\n"
     ]
    }
   ],
   "source": [
    "# comparison to csv for plotting in R\n",
    "start = time.time()\n",
    "!sourmash compare uniflna_HA_sub1.fasta.sig --csv uniflna_HA_sub1_cmp.csv\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## not just HA, but all seqs..."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "np.random.seed(8903)\n",
    "uniflna_sub1 = uniflna.sample(n=1000)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# make it back into text FASTA file\n",
    "with open (\"uniflna_sub1.fasta\", 'w') as output:\n",
    "    for i in range(len(uniflna_sub1)):\n",
    "        output.write(\">\" + list(uniflna_sub1.titles.values)[i] + \"\\n\" + list(uniflna_sub1.sequences.values)[i] + \"\\n\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "\u001b[Kcomputing signatures for files: uniflna_sub1.fasta\n",
      "\u001b[KComputing signature for ksizes: [5]\n",
      "\u001b[KComputing only nucleotide (and not protein) signatures.\n",
      "\u001b[KComputing a total of 1 signature(s).\n",
      "\u001b[Kcalculated 1000 signatures for 1000 sequences in uniflna_sub1.fasta\n",
      "\u001b[Ksaved 1000 signature(s). Note: signature license is CC0.\n",
      "finished in 1.0115101337432861seconds\n"
     ]
    }
   ],
   "source": [
    "start = time.time()\n",
    "!sourmash compute uniflna_sub1.fasta  -k 5 --singleton\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "...sig loading 999_sub1.fasta.sig\n",
      "\u001b[Kloaded 1000 signatures total.                                                  \n",
      "\u001b[K\n",
      "min similarity in matrix: 0.079\n",
      "finished in 65.50607323646545seconds\n"
     ]
    }
   ],
   "source": [
    "# comparison to csv for plotting in R\n",
    "start = time.time()\n",
    "!sourmash compare uniflna_sub1.fasta.sig --csv uniflna_sub1_cmp.csv\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "...sig loading 999_sub1.fasta.sig\n",
      "\u001b[Kloaded 1000 signatures total.                                                  \n",
      "\u001b[K\n",
      "min similarity in matrix: 0.079\n",
      "\u001b[Ksaving labels to: uniflna_sub1_cmp.labels.txt\n",
      "\u001b[Ksaving distance matrix to: uniflna_sub1_cmp\n",
      "finished in 67.71590304374695seconds\n"
     ]
    }
   ],
   "source": [
    "# comparison to matrx\n",
    "start = time.time()\n",
    "!sourmash compare uniflna_sub1.fasta.sig -o uniflna_sub1_cmp\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "\u001b[Kloading comparison matrix from uniflna_sub1_cmp...\n",
      "\u001b[K...got 1000 x 1000 matrix.\n",
      "\u001b[Kloading labels from uniflna_sub1_cmp.labels.txt\n",
      "\u001b[Ksaving histogram of matrix values => uniflna_sub1_cmp.hist.png\n",
      "\u001b[Kwrote dendrogram to: uniflna_sub1_cmp.dendro.png\n",
      "\u001b[Kwrote numpy distance matrix to: uniflna_sub1_cmp.matrix.png\n"
     ]
    }
   ],
   "source": [
    "!sourmash plot uniflna_sub1_cmp"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "tree:  \n",
    "![dendogram](uniflna_sub1_cmp.dendro.png)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "\u001b[Kcomputing signatures for files: uniflna_sub1.fasta\n",
      "\u001b[KComputing signature for ksizes: [7]\n",
      "\u001b[KComputing only nucleotide (and not protein) signatures.\n",
      "\u001b[KComputing a total of 1 signature(s).\n",
      "\u001b[Kcalculated 1000 signatures for 1000 sequences in uniflna_sub1.fasta\n",
      "\u001b[Ksaved 1000 signature(s). Note: signature license is CC0.\n",
      "finished in 1.1013703346252441seconds\n"
     ]
    }
   ],
   "source": [
    "start = time.time()\n",
    "!sourmash compute uniflna_sub1.fasta  -k 7 --singleton\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "...sig loading 999_sub1.fasta.sig\n",
      "\u001b[Kloaded 1000 signatures total.                                                  \n",
      "\u001b[K\n",
      "min similarity in matrix: 0.002\n",
      "\u001b[Ksaving labels to: uniflna_sub1_k7_cmp.labels.txt\n",
      "\u001b[Ksaving distance matrix to: uniflna_sub1_k7_cmp\n",
      "finished in 57.62817907333374seconds\n"
     ]
    }
   ],
   "source": [
    "# comparison to matrx\n",
    "start = time.time()\n",
    "!sourmash compare uniflna_sub1.fasta.sig -o uniflna_sub1_k7_cmp\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "\u001b[Kloading comparison matrix from uniflna_sub1_k7_cmp...\n",
      "\u001b[K...got 1000 x 1000 matrix.\n",
      "\u001b[Kloading labels from uniflna_sub1_k7_cmp.labels.txt\n",
      "\u001b[Ksaving histogram of matrix values => uniflna_sub1_k7_cmp.hist.png\n",
      "\u001b[Kwrote dendrogram to: uniflna_sub1_k7_cmp.dendro.png\n",
      "\u001b[Kwrote numpy distance matrix to: uniflna_sub1_k7_cmp.matrix.png\n"
     ]
    }
   ],
   "source": [
    "!sourmash plot uniflna_sub1_k7_cmp"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "tree:  \n",
    "![dendogram](uniflna_sub1_k7_cmp.dendro.png)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "...sig loading 999_sub1.fasta.sig\n",
      "\u001b[Kloaded 1000 signatures total.                                                  \n",
      "\u001b[K\n",
      "min similarity in matrix: 0.002\n",
      "finished in 57.83396553993225seconds\n"
     ]
    }
   ],
   "source": [
    "# comparison to csv for plotting in R\n",
    "start = time.time()\n",
    "!sourmash compare uniflna_sub1.fasta.sig --csv uniflna_sub1_k7_cmp.csv\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### add back gene labels in new col"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "np.random.seed(78453)\n",
    "uniflna_test1 = uniflna.sample(n=1000)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['gi|1410306180|gb|MH546639|Influenza A virus (A/mallard/Minnesota/15-040386-2/2015(H11N1)) segment 8 nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|1370449901|gb|MG957554|Influenza A virus (A/duck/Bangladesh/32003/2017(H5N1)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|1048081328|gb|KX616653|Influenza B virus (B/Texas/62/2016) segment 3 polymerase PA (PA) gene, complete cds',\n",
       "       'gi|525553452|gb|KF314561|Influenza A virus (A/India/K0914850/2009(H3N2)) segment 6 neuraminidase (NA) gene, partial cds',\n",
       "       'gi|237874843|gb|FJ985245|Influenza A virus (A/turkey/France/05045/2005(H1N2)) segment 1 polymerase PB2 (PB2) gene, partial cds',\n",
       "       'gi|315440734|gb|CY079457|Influenza A virus (A/American black duck/Illinois/08OS2688/2008(H5N2)) segment 3, complete sequence',\n",
       "       'gi|937172275|gb|KT854258|Influenza B virus (B/Texas/09/2015) segment 7 matrix protein 1 (M1) and BM2 protein (BM2) genes, complete cds',\n",
       "       'gi|594215428|gb|CY174927|Influenza B virus (B/New York/1157/2008) polymerase PB1 (PB1) gene, complete cds',\n",
       "       'gi|198285780|gb|CY030920|Influenza A virus (A/goose/Yunnan/4389/2006(H5N1)) segment 3 sequence',\n",
       "       'gi|1120605462|gb|KY347748|Influenza A virus (A/wild bird feces/Anhui/S93/2014(H9N2)) segment 8 nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|660873930|gb|KM064239|Influenza A virus (A/Washington/10/2014(H3N2)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|1219455324|gb|CY239151|Influenza B virus (B/Michigan/78/2017) polymerase PA (PA) gene, complete cds',\n",
       "       'gi|1390209297|gb|MH363175|Influenza B virus (B/Louisiana/05/2018) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|512772000|gb|KF216939|Influenza B virus (B/Georgia/06/2013) segment 6 neuraminidase (NA) gene, complete cds; and NB protein (NB) gene, partial cds',\n",
       "       'gi|1421412631|gb|MH598206|Influenza A virus (A/ruddy turnstone/New Jersey/AI00-2296/2000(H12N5)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|1041633881|gb|LC123364|Influenza C virus (C/Yamagata/1/2012) PB2 gene for polymerase 2, complete cds',\n",
       "       'gi|324345671|gb|CY084441|Influenza A virus (A/New York/6064/2009(H1N1)) nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|336286441|gb|JN087059|Influenza A virus (A/environment/Korea/ESD3-3/2004(H3N8)) segment 7 matrix protein 2 (M2) gene, partial cds; and matrix protein 1 (M1) gene, complete cds',\n",
       "       'gi|755195545|gb|KP456628|Influenza A virus (A/Alaska/29/2009(H3N2)) segment 4 hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|974741015|gb|CY195307|Influenza A virus (A/swine/Manitoba/D0349/2014(H3N2)) nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|283551383|gb|CY047996|Influenza A virus (A/chicken/Nigeria/08VIR4337-377/2006(H5N1)) segment 5 sequence',\n",
       "       'gi|1237086641|gb|MF630455|Influenza A virus (A/duck/Hunan/S11682/2015(H7N9)) segment 6 neuraminidase (NA) gene, complete cds',\n",
       "       'gi|563400421|gb|CY158761|Influenza A virus (A/swine/Manitoba/02947/2010(H1N1)) hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|607259322|gb|KJ579959|Influenza A virus (A/pigeon/Anhui/08/2013(H3N8)) segment 5 nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|755202521|gb|KP460116|Influenza B virus (B/Kochi/1/2008) segment 4 hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|827499629|gb|KR701147|Influenza A virus (A/swine/Netherlands/Wanroij-CVI2117A/2012(H1N2)) segment 6 neuraminidase (NA) gene, complete cds',\n",
       "       'gi|1176711220|gb|KY855523|Influenza A virus (A/chicken/Guangdong/J2/2017(H7N9)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|1246295632|gb|MF801575|Influenza A virus (A/swine/Iowa/A01667088/2017(H1N1)) segment 2 polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|77747051|gb|CY003463|Influenza A virus (A/New York/439/2000(H3N2)) segment 1, complete sequence',\n",
       "       'gi|1408067743|gb|MH341870|Influenza A virus (A/mallard/Minnesota/UGAI15-6021/2015(H1N1)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|336286770|gb|JN087168|Influenza A virus (A/duck/Korea/U12/2007(H3N2)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|946524900|gb|KT889009|Influenza A virus (A/Philippines/FLU_AFPA_00013/2011(H3N2)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|1390210015|gb|MH363477|Influenza B virus (B/Texas/19/2018) segment 1 polymerase PB1 (PB1) gene, complete cds',\n",
       "       'gi|1435279360|gb|MH560176|Influenza A virus (A/chicken/Vietnam/H7F-LC4-371/2014(H9N2)) segment 8 nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|73665566|gb|CY000866|Influenza A virus (A/New York/193/2003(H3N2)) segment 7, complete sequence',\n",
       "       'gi|56785893|gb|AY738455|Influenza A virus (A/chicken/Talmei Elazar/1304/03(H9N2)) hemagglutinin precursor, gene, complete cds',\n",
       "       'gi|297632072|gb|AB537453|Influenza A virus (A/duck/Shiga/10413-24/2006(H3N8)) HA gene for hemagglutinin, complete cds',\n",
       "       'gi|336286387|gb|JN087276|Influenza A virus (A/duck/Korea/GJ87/2007(H3N8)) segment 8 nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|332079329|gb|JF789635|Influenza A virus (A/Astana/S2/2009(H1N1)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|393712449|gb|CY124530|Influenza A virus (A/Singapore/GP39/2011(H1N1)) polymerase PB1 (PB1) gene, complete cds; and PB1-F2 gene, complete sequence',\n",
       "       'gi|1073679892|gb|KX919563|Influenza A virus (A/Pennsylvania/12/2016(H1N1)) segment 3 polymerase PA (PA) and PA-X protein (PA-X) genes, complete cds',\n",
       "       'gi|660873760|gb|KM064199|Influenza A virus (A/Massachusetts/16/2013(H3N2)) segment 6 neuraminidase (NA) gene, complete cds',\n",
       "       'gi|37785075|gb|AY210036|Influenza A virus (A/Panama/1/61(H2N2)) membrane ion channel M2 and matrix protein M1 (M) gene, complete cds',\n",
       "       'gi|1060100525|gb|KX684034|Influenza A virus (A/swine/Iowa/A01777664/2016(H3N2)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|82654283|gb|CY005452|Influenza A virus (A/Chicken/Nanchang/7-010/2000(H3N6)) segment 2, complete sequence',\n",
       "       'gi|469935292|gb|KC738250|Influenza A virus (A/Taiwan/1151/1992(H3N2)) segment 2 polymerase PB1 (PB1) gene, partial cds; and PB1-F2 protein (PB1-F2) gene, complete cds',\n",
       "       'gi|383278558|gb|CY113883|Influenza A virus (A/Paris/614/1992(H3N2)) polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|1028349584|gb|KX005595|Influenza A virus (A/Texas/07/2016(H1N1)) clone 3000416252_285445_v1_5 segment 5 nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|1128979015|gb|CY212029|Influenza A virus (A/Wisconsin/97/2016(H3N2)) matrix protein 2 (M2) and matrix protein 1 (M1) genes, complete cds',\n",
       "       'gi|169822690|gb|CY030784|Influenza B virus (B/Tennessee/UR06-0079/2007) segment 7, complete sequence',\n",
       "       'gi|910300283|gb|CY191631|Influenza B virus (B/New York/WC-LVD-14-037/2014) nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|1009090830|gb|KU877383|Influenza A virus (A/swine/Illinois/A01731417/2016(H1N2)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|549434742|gb|CY156096|Influenza B virus (B/Wellington/9/2007) polymerase PB1 (PB1) gene, complete cds',\n",
       "       'gi|745630148|gb|KP285522|Influenza A virus (A/duck/Jiangxi/33641/2013(H10N3)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|93008624|gb|DQ493099|Influenza A virus (A/chicken/Kupang-1-NTT/BPPV6/2004(H5N1)) nucleoprotein (NP) gene, partial cds',\n",
       "       'gi|1181550534|gb|KY828763|Influenza A virus (A/northern shoveler/California/LDC185/2014(H8N4)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|452832157|gb|CY134861|Influenza A virus (A/Texas/JMM_19/2012(H3N2)) matrix protein 2 (M2) and matrix protein 1 (M1) genes, complete cds',\n",
       "       'gi|1168686406|gb|CY223564|Influenza A virus (A/Oklahoma/08/2017(H1N1)) polymerase PA (PA) and PA-X protein (PA-X) genes, complete cds',\n",
       "       'gi|1162229649|gb|KY678589|Influenza A virus (A/swine/Minnesota/A01781391/2016(H1N1)) segment 5 nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|1015396647|gb|KU982582|Influenza A virus (A/swine/Hong Kong/9/98-MA(H9N2)) escape mutant RAm8C4(1) segment 4 hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|310699813|gb|CY076324|Influenza A virus (A/northern shoveler/Washington/44249-783/2006(H7N3)) polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|1487102559|gb|MK025712|Influenza A virus (A/Florida/77/2018(H1N1)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|984688745|gb|KU589729|Influenza A virus (A/Ohio/28/2015(H3N2)) segment 2 polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|303303804|gb|CY070957|Influenza A virus (A/Managua/1181.03/2010(H3N2)) segment 2, complete sequence',\n",
       "       'gi|1215244063|gb|CY236525|Influenza B virus (B/Utah/08/2017) polymerase PA (PA) gene, complete cds',\n",
       "       'gi|27465228|gb|AY180404|Influenza A virus (A/Bantam/Nanchang/9-366/2000(H3N3)) hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|974988724|gb|CY196181|Influenza A virus (A/American green-winged teal/Interior Alaska/11BM03479/2011(H4N6)) nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|937145279|gb|KT836974|Influenza A virus (A/Michigan/13/2014(H3N2)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|937151569|gb|KT839367|Influenza A virus (A/Kentucky/14/2014(H3N2)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|727202210|gb|CY103829|Influenza A virus (A/duck/Viet Nam/TMU020/2009(H5N1)) matrix protein 2 (M2) and matrix protein 1 (M1) genes, complete cds',\n",
       "       'gi|411392031|gb|CY126159|Influenza A virus (A/chicken/Qalubia/CAI12/2008(H5N1)) polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|1276993959|gb|CY251499|Influenza A virus (A/Victoria/504/2013(H3N2)) nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|120434185|gb|CY019122|Influenza A virus (A/Memphis/3/1996(H1N1)) segment 3, complete sequence',\n",
       "       'gi|1474769382|gb|MH885281|Influenza A virus (A/Texas/6867/2017(H3N2)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|380836483|gb|CY108700|Influenza A virus (A/Jilin-Chaoyang/133/2007(H3N2)) hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|1248251053|gb|MF955655|Influenza B virus (B/Hong Kong/CUHK37079/2010) segment 4 hemagglutinin gene, complete cds',\n",
       "       'gi|1140316078|gb|KY463977|Influenza A virus (A/mallard/Ohio/14OS2695/2014(mixed)) segment 2 polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|1168684613|gb|CY222823|Influenza A virus (A/Massachusetts/04/2017(H3N2)) neuraminidase (NA) gene, complete cds',\n",
       "       'gi|474455949|gb|KC780219|Influenza A virus (A/Utah/33/2009(H1N1)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|589912749|gb|KJ439872|Influenza A virus (A/duck/Zhejiang/6D7/2013(H3N2)) segment 6 neuraminidase (NA) gene, complete cds',\n",
       "       'gi|293599381|gb|CY059418|Influenza A virus (A/Shenfang/262/1995(H3N2)) segment 7 sequence',\n",
       "       'gi|357382150|gb|CY101657|Influenza A virus (A/mallard/Alberta/220/1990(H3N8)) polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|161138879|gb|CY027957|Influenza A virus (A/Oregon/UR06-0185/2007(H1N1)) segment 6, complete sequence',\n",
       "       'gi|831182692|gb|KR821025|Influenza A virus (A/Mexico/INMEGEN-INER 09/2012(H1N1)) segment 5 nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|984688450|gb|KU589610|Influenza A virus (A/Florida/80/2015(H3N2)) segment 2 polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|83940846|gb|DQ299844|Influenza A virus (A/chicken/Korea/S24/04(H9N2)) polymerase PB2 (PB2) gene, partial cds',\n",
       "       'gi|354681921|gb|AB676809|Influenza A virus (A/goshawk/Tochigi/64/2011(H5N1)) PA gene for polymerase acidic protein, partial cds',\n",
       "       'gi|1237085517|gb|MF630005|Influenza A virus (A/chicken/Guangdong/SC086/2013(H7N9)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|1149046518|gb|CY219243|Influenza A virus (A/California/11/2017(H3N2)) polymerase PA (PA) and PA-X protein (PA-X) genes, complete cds',\n",
       "       'gi|1205299466|gb|MF147727|Influenza A virus (A/mallard duck/Netherlands/24/2011(H2N3)) segment 6 neuraminidase (NA) gene, complete cds',\n",
       "       'gi|552052680|gb|KF695222|Influenza A virus (A/mallard/Sweden/54/2003(H6N1)) segment 3 polymerase PA (PA) and PA-X protein (PA-X) genes, complete cds',\n",
       "       'gi|1350895264|gb|MG856252|Influenza A virus (A/Mexico/InDRE2362/2007(H3N2)) segment 5 nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|188504541|gb|CY031660|Influenza A virus (A/environment/New York/23164-1/2005(H7N2)) segment 7, complete sequence',\n",
       "       'gi|1464622050|gb|MH798603|Influenza A virus (A/New York/PV37/2017(H3N2)) segment 3 polymerase PA (PA) and PA-X protein (PA-X) genes, complete cds',\n",
       "       'gi|1421640142|gb|MH600510|Influenza A virus (A/West Virginia/25/2018(H1N1)) segment 3 polymerase PA (PA) and PA-X protein (PA-X) genes, complete cds',\n",
       "       'gi|662490809|gb|KM027976|Influenza A virus (A/swine/Guangxi/2841/2010(H1N1)) segment 5 nucleocapsid protein (NP) gene, partial cds',\n",
       "       'gi|1378290338|gb|MH216554|Influenza B virus (B/Alberta/006/2017) segment 4 hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|488464962|gb|KC891723|Influenza B virus (B/New Jersey/01/2012) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|574455772|gb|KF975499|Influenza A virus (A/chicken/Pakistan/288-102/2010(H9N2)) segment 8 nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds'],\n",
       "      dtype=object)"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# uniflna_test1.titles.values"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {
    "collapsed": true
   },
   "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>titles</th>\n",
       "      <th>sequences</th>\n",
       "      <th>gene</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>629174</th>\n",
       "      <td>gi|1410306180|gb|MH546639|Influenza A virus (A...</td>\n",
       "      <td>GTGACAAAAACATAATGGATTCCAACACGATAACCTCGTTTCAGGT...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>610132</th>\n",
       "      <td>gi|1370449901|gb|MG957554|Influenza A virus (A...</td>\n",
       "      <td>AATATGGAGAGAATAAAAGAACTAAGAGATTTGATGTCGCAGTCTC...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>470806</th>\n",
       "      <td>gi|1048081328|gb|KX616653|Influenza B virus (B...</td>\n",
       "      <td>GTTTGATTTATCATAATGGATACTTTTATTACAAGAAACTTCCAGA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>278994</th>\n",
       "      <td>gi|525553452|gb|KF314561|Influenza A virus (A/...</td>\n",
       "      <td>AATAACGATTGGCTCTGTTTCTCTCACCATTTCCACAATATGCTTC...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>86479</th>\n",
       "      <td>gi|237874843|gb|FJ985245|Influenza A virus (A/...</td>\n",
       "      <td>ATTCAATATGGAAAGAATAAAAGAATTAAGAGATCTAATGTCACAG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>151351</th>\n",
       "      <td>gi|315440734|gb|CY079457|Influenza A virus (A/...</td>\n",
       "      <td>TACTGATTCAAAATGGAAGACTTTGTGCGACAATGCTTCAATCCAA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>412958</th>\n",
       "      <td>gi|937172275|gb|KT854258|Influenza B virus (B/...</td>\n",
       "      <td>ATGTCGCTGTTTGGAGACACAATTGCCTACCTGCTTTCATTGACAG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>315125</th>\n",
       "      <td>gi|594215428|gb|CY174927|Influenza B virus (B/...</td>\n",
       "      <td>ATGAATATAAATCCTTATTTTCTCTTCATAGATGTACCCATACAGG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74069</th>\n",
       "      <td>gi|198285780|gb|CY030920|Influenza A virus (A/...</td>\n",
       "      <td>ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>505337</th>\n",
       "      <td>gi|1120605462|gb|KY347748|Influenza A virus (A...</td>\n",
       "      <td>ATGGATTCCAACACGATAACCTCGTTTCAGGTAGATTGTTATCTAT...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>340273</th>\n",
       "      <td>gi|660873930|gb|KM064239|Influenza A virus (A/...</td>\n",
       "      <td>ATGAAGACTATCATTGCTTTGAGCTACACTCTATGTCTGGTTTTCG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>556365</th>\n",
       "      <td>gi|1219455324|gb|CY239151|Influenza B virus (B...</td>\n",
       "      <td>GTTTGATTTGTCATAATGGATACTTTTATTACAAGAAACTTCCAGA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>623475</th>\n",
       "      <td>gi|1390209297|gb|MH363175|Influenza B virus (B...</td>\n",
       "      <td>ATTTTCTAATATCCACAGAATGAAGGCAATAATTGTACTACTCATG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>273912</th>\n",
       "      <td>gi|512772000|gb|KF216939|Influenza B virus (B/...</td>\n",
       "      <td>ATGCTACCTTCAACTATACAAACGTTAACCCTATTTCTCACATCAG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>636364</th>\n",
       "      <td>gi|1421412631|gb|MH598206|Influenza A virus (A...</td>\n",
       "      <td>GGTCACAATGGAAAAATTCATCATTTTGAGCACTGTCCTGGCAGTA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>467742</th>\n",
       "      <td>gi|1041633881|gb|LC123364|Influenza C virus (C...</td>\n",
       "      <td>ATGTCTCTTCTATTGACAATAGCAAAGGAATACAAAAGACTATGCC...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>156141</th>\n",
       "      <td>gi|324345671|gb|CY084441|Influenza A virus (A/...</td>\n",
       "      <td>ATGAGTGACATCGAAGCCATGGCGTCTCAAGGCACCAAACGATCAT...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>165106</th>\n",
       "      <td>gi|336286441|gb|JN087059|Influenza A virus (A/...</td>\n",
       "      <td>ATATTGAAAGATGAGTCTTCTAACCGAGGTCGAAACGTACGTTCTC...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>371494</th>\n",
       "      <td>gi|755195545|gb|KP456628|Influenza A virus (A/...</td>\n",
       "      <td>ATGAAGACTATCATTGCTTTGAGCTACATTCTATGTCTGGTTTTCG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>424904</th>\n",
       "      <td>gi|974741015|gb|CY195307|Influenza A virus (A/...</td>\n",
       "      <td>ATGGACTCCAATACTGTGTCAAGCTTTCAGGTAGACTGTTTCCTTT...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>111975</th>\n",
       "      <td>gi|283551383|gb|CY047996|Influenza A virus (A/...</td>\n",
       "      <td>GGAGAACGCCAGAATGCTACTGAGATCAGGGCATCTGTTGGAAGAA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>561122</th>\n",
       "      <td>gi|1237086641|gb|MF630455|Influenza A virus (A...</td>\n",
       "      <td>ATGAATCCAAATCAGAAGATTCTATGCACTTCAGCCACTGCTATCC...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>298789</th>\n",
       "      <td>gi|563400421|gb|CY158761|Influenza A virus (A/...</td>\n",
       "      <td>ATGAAGGCAATACTATTAGTCTTGTTATATACATTTGTACCCGCAA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>328210</th>\n",
       "      <td>gi|607259322|gb|KJ579959|Influenza A virus (A/...</td>\n",
       "      <td>TTATTCGTCTCAGGGAGCAAAAGCAGGGTAGATAATCACTCACTGA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>374982</th>\n",
       "      <td>gi|755202521|gb|KP460116|Influenza B virus (B/...</td>\n",
       "      <td>GATCGAATCTGCACTGGGATAACATCTTCAAACTCACCTCATGTGG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>389631</th>\n",
       "      <td>gi|827499629|gb|KR701147|Influenza A virus (A/...</td>\n",
       "      <td>ATGAATCCAAATCAGAAGATAATAACAATTGGCTCTGTTTCTCTAT...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>529407</th>\n",
       "      <td>gi|1176711220|gb|KY855523|Influenza A virus (A...</td>\n",
       "      <td>AGCGAAAGCAGGTCAAATATATTCAATATGGAAAGAATAAAAGAAC...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>567789</th>\n",
       "      <td>gi|1246295632|gb|MF801575|Influenza A virus (A...</td>\n",
       "      <td>ATGGATGTCAATCCGACCCTACTTTTCTTAAAGGTTCCAGCGCAAA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18713</th>\n",
       "      <td>gi|77747051|gb|CY003463|Influenza A virus (A/N...</td>\n",
       "      <td>AGCAAAAGCAGGTCAATTATATTCAGTATGGAAAGAATAAAAGAAC...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>627939</th>\n",
       "      <td>gi|1408067743|gb|MH341870|Influenza A virus (A...</td>\n",
       "      <td>TCAAATCAACCAAGATGGAAGCAAAACTGTTTGTACTATTCTGTAC...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>355132</th>\n",
       "      <td>gi|727202210|gb|CY103829|Influenza A virus (A/...</td>\n",
       "      <td>AGCAAAAGCAGGTAGATATTGAAAAATGAGTCTTCTAACTGAGGTC...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>226790</th>\n",
       "      <td>gi|411392031|gb|CY126159|Influenza A virus (A/...</td>\n",
       "      <td>TCAAATATATTCAATATGGAGAGAATAAAGGAATTAAGAGACCTAA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>577113</th>\n",
       "      <td>gi|1276993959|gb|CY251499|Influenza A virus (A...</td>\n",
       "      <td>GTGACAAAGACATAATGGATTCCAACACTGTGTCAAGTTTCCAGGT...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>40450</th>\n",
       "      <td>gi|120434185|gb|CY019122|Influenza A virus (A/...</td>\n",
       "      <td>ATTCGAAATGGAAGATTTTGTGCGACAGTGCTTCAATCCGATGATT...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>650629</th>\n",
       "      <td>gi|1474769382|gb|MH885281|Influenza A virus (A...</td>\n",
       "      <td>ATGAAGACTATCATTGCTTTGAGCTACATTCTATGTCTGGTTTTCG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>196481</th>\n",
       "      <td>gi|380836483|gb|CY108700|Influenza A virus (A/...</td>\n",
       "      <td>CAAAAACTTCCCGGAAATGACAACAGCACGGCAACGCTGTGCCTTG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>568332</th>\n",
       "      <td>gi|1248251053|gb|MF955655|Influenza B virus (B...</td>\n",
       "      <td>ATGAAGGCAATAATTGTACTACTCATGGTAGTAACATCCAATGCAG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>512342</th>\n",
       "      <td>gi|1140316078|gb|KY463977|Influenza A virus (A...</td>\n",
       "      <td>ATGGATGTCAATCCGACTTTACTTTTCTTGAAAGTTCCAGCGCAAA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>526752</th>\n",
       "      <td>gi|1168684613|gb|CY222823|Influenza A virus (A...</td>\n",
       "      <td>AGTAAAGATGAATCCAAATCAAAAGATAATAACGATTGGCTCTGTT...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>254049</th>\n",
       "      <td>gi|474455949|gb|KC780219|Influenza A virus (A/...</td>\n",
       "      <td>ATGAAGGCAATACTAGTAGTTCTGCTATATACATTTGCAACCGCAA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>308720</th>\n",
       "      <td>gi|589912749|gb|KJ439872|Influenza A virus (A/...</td>\n",
       "      <td>ATGAATCCAAATCAGAAGATAATAACAATTGGCTCTGTCTCTCTAA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>120107</th>\n",
       "      <td>gi|293599381|gb|CY059418|Influenza A virus (A/...</td>\n",
       "      <td>CCCCCTCAAAGCCGAAATCGCGCAGAGACTTGAAGATGTCTTTGCT...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>182482</th>\n",
       "      <td>gi|357382150|gb|CY101657|Influenza A virus (A/...</td>\n",
       "      <td>CAAACCATTTGAATGGATGTCAATCCGACTTTACTTTTCTTGAAAG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>58027</th>\n",
       "      <td>gi|161138879|gb|CY027957|Influenza A virus (A/...</td>\n",
       "      <td>TTTAAAATGAATCCAAATCAAAAAATAATAACCATTGGATCAATCA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>390631</th>\n",
       "      <td>gi|831182692|gb|KR821025|Influenza A virus (A/...</td>\n",
       "      <td>ATGGCGTCTCAAGGCACCAAACGATCATATGAACAAATGGAGACTG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>427517</th>\n",
       "      <td>gi|984688450|gb|KU589610|Influenza A virus (A/...</td>\n",
       "      <td>ATGGATGTCAATCCGACTCTACTGTTCTTAAAAGTTCCAGCGCAAA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23280</th>\n",
       "      <td>gi|83940846|gb|DQ299844|Influenza A virus (A/c...</td>\n",
       "      <td>AGCAAAAGCAGGTCAATTATATTCAATATGGAGAGAATAAAAGAAC...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>180865</th>\n",
       "      <td>gi|354681921|gb|AB676809|Influenza A virus (A/...</td>\n",
       "      <td>TATGGGGAAGATCCGAAAATCGAAACAAACAAATTTGCCGCAATAT...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>560672</th>\n",
       "      <td>gi|1237085517|gb|MF630005|Influenza A virus (A...</td>\n",
       "      <td>ATGAACACTCAAATCCTGGTATTCGCTCTGATTGCGATCATTCCAA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>520172</th>\n",
       "      <td>gi|1149046518|gb|CY219243|Influenza A virus (A...</td>\n",
       "      <td>TACTGATTCAAAATGGAAGATTTTGTGCGACAATGCTTCAACCCGA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>546639</th>\n",
       "      <td>gi|1205299466|gb|MF147727|Influenza A virus (A...</td>\n",
       "      <td>AGCAAAAGCAGGTGCGAAATGAATCCAAATCAGAAGATAATAACAA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>293165</th>\n",
       "      <td>gi|552052680|gb|KF695222|Influenza A virus (A/...</td>\n",
       "      <td>TACTGATCCAAAATGGAAGACTTTGTGCGACAATGCTTCAATCCAA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>604028</th>\n",
       "      <td>gi|1350895264|gb|MG856252|Influenza A virus (A...</td>\n",
       "      <td>AGCAAAAGCAGGGTTAATAATCACTCACTGAGTGACATCAAAATCA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>66867</th>\n",
       "      <td>gi|188504541|gb|CY031660|Influenza A virus (A/...</td>\n",
       "      <td>AAGATGAGTCTTCTAACCGAGGTCGAAACGTACGTTCTCTCTATCA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>649040</th>\n",
       "      <td>gi|1464622050|gb|MH798603|Influenza A virus (A...</td>\n",
       "      <td>CTGATTCAAAATGGAAGATTTTGTGCGACAATGCTTCAACCCGATG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>636950</th>\n",
       "      <td>gi|1421640142|gb|MH600510|Influenza A virus (A...</td>\n",
       "      <td>TACTGATCCAAAATGGAAGACTTTGTGCGACAATGCTTCAATCCAA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>341110</th>\n",
       "      <td>gi|662490809|gb|KM027976|Influenza A virus (A/...</td>\n",
       "      <td>ATGGCGTCTCAAGGCACCAAACGGTCTTATGAACAAATGGAAACTG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>616949</th>\n",
       "      <td>gi|1378290338|gb|MH216554|Influenza B virus (B...</td>\n",
       "      <td>ATGAAGGCAATAATTGTACTACTCATGGTAGTAACATCCAATGCAG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>263002</th>\n",
       "      <td>gi|488464962|gb|KC891723|Influenza B virus (B/...</td>\n",
       "      <td>ATGAAGGCAATAATTGTACTACTCATGGTAGTAACATCCAATGCAG...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>304678</th>\n",
       "      <td>gi|574455772|gb|KF975499|Influenza A virus (A/...</td>\n",
       "      <td>TGACAAAAACATAATGGATTCCAACACTGTGTCAAGCTTTCAGGTA...</td>\n",
       "      <td>not</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>99 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                   titles  \\\n",
       "629174  gi|1410306180|gb|MH546639|Influenza A virus (A...   \n",
       "610132  gi|1370449901|gb|MG957554|Influenza A virus (A...   \n",
       "470806  gi|1048081328|gb|KX616653|Influenza B virus (B...   \n",
       "278994  gi|525553452|gb|KF314561|Influenza A virus (A/...   \n",
       "86479   gi|237874843|gb|FJ985245|Influenza A virus (A/...   \n",
       "151351  gi|315440734|gb|CY079457|Influenza A virus (A/...   \n",
       "412958  gi|937172275|gb|KT854258|Influenza B virus (B/...   \n",
       "315125  gi|594215428|gb|CY174927|Influenza B virus (B/...   \n",
       "74069   gi|198285780|gb|CY030920|Influenza A virus (A/...   \n",
       "505337  gi|1120605462|gb|KY347748|Influenza A virus (A...   \n",
       "340273  gi|660873930|gb|KM064239|Influenza A virus (A/...   \n",
       "556365  gi|1219455324|gb|CY239151|Influenza B virus (B...   \n",
       "623475  gi|1390209297|gb|MH363175|Influenza B virus (B...   \n",
       "273912  gi|512772000|gb|KF216939|Influenza B virus (B/...   \n",
       "636364  gi|1421412631|gb|MH598206|Influenza A virus (A...   \n",
       "467742  gi|1041633881|gb|LC123364|Influenza C virus (C...   \n",
       "156141  gi|324345671|gb|CY084441|Influenza A virus (A/...   \n",
       "165106  gi|336286441|gb|JN087059|Influenza A virus (A/...   \n",
       "371494  gi|755195545|gb|KP456628|Influenza A virus (A/...   \n",
       "424904  gi|974741015|gb|CY195307|Influenza A virus (A/...   \n",
       "111975  gi|283551383|gb|CY047996|Influenza A virus (A/...   \n",
       "561122  gi|1237086641|gb|MF630455|Influenza A virus (A...   \n",
       "298789  gi|563400421|gb|CY158761|Influenza A virus (A/...   \n",
       "328210  gi|607259322|gb|KJ579959|Influenza A virus (A/...   \n",
       "374982  gi|755202521|gb|KP460116|Influenza B virus (B/...   \n",
       "389631  gi|827499629|gb|KR701147|Influenza A virus (A/...   \n",
       "529407  gi|1176711220|gb|KY855523|Influenza A virus (A...   \n",
       "567789  gi|1246295632|gb|MF801575|Influenza A virus (A...   \n",
       "18713   gi|77747051|gb|CY003463|Influenza A virus (A/N...   \n",
       "627939  gi|1408067743|gb|MH341870|Influenza A virus (A...   \n",
       "...                                                   ...   \n",
       "355132  gi|727202210|gb|CY103829|Influenza A virus (A/...   \n",
       "226790  gi|411392031|gb|CY126159|Influenza A virus (A/...   \n",
       "577113  gi|1276993959|gb|CY251499|Influenza A virus (A...   \n",
       "40450   gi|120434185|gb|CY019122|Influenza A virus (A/...   \n",
       "650629  gi|1474769382|gb|MH885281|Influenza A virus (A...   \n",
       "196481  gi|380836483|gb|CY108700|Influenza A virus (A/...   \n",
       "568332  gi|1248251053|gb|MF955655|Influenza B virus (B...   \n",
       "512342  gi|1140316078|gb|KY463977|Influenza A virus (A...   \n",
       "526752  gi|1168684613|gb|CY222823|Influenza A virus (A...   \n",
       "254049  gi|474455949|gb|KC780219|Influenza A virus (A/...   \n",
       "308720  gi|589912749|gb|KJ439872|Influenza A virus (A/...   \n",
       "120107  gi|293599381|gb|CY059418|Influenza A virus (A/...   \n",
       "182482  gi|357382150|gb|CY101657|Influenza A virus (A/...   \n",
       "58027   gi|161138879|gb|CY027957|Influenza A virus (A/...   \n",
       "390631  gi|831182692|gb|KR821025|Influenza A virus (A/...   \n",
       "427517  gi|984688450|gb|KU589610|Influenza A virus (A/...   \n",
       "23280   gi|83940846|gb|DQ299844|Influenza A virus (A/c...   \n",
       "180865  gi|354681921|gb|AB676809|Influenza A virus (A/...   \n",
       "560672  gi|1237085517|gb|MF630005|Influenza A virus (A...   \n",
       "520172  gi|1149046518|gb|CY219243|Influenza A virus (A...   \n",
       "546639  gi|1205299466|gb|MF147727|Influenza A virus (A...   \n",
       "293165  gi|552052680|gb|KF695222|Influenza A virus (A/...   \n",
       "604028  gi|1350895264|gb|MG856252|Influenza A virus (A...   \n",
       "66867   gi|188504541|gb|CY031660|Influenza A virus (A/...   \n",
       "649040  gi|1464622050|gb|MH798603|Influenza A virus (A...   \n",
       "636950  gi|1421640142|gb|MH600510|Influenza A virus (A...   \n",
       "341110  gi|662490809|gb|KM027976|Influenza A virus (A/...   \n",
       "616949  gi|1378290338|gb|MH216554|Influenza B virus (B...   \n",
       "263002  gi|488464962|gb|KC891723|Influenza B virus (B/...   \n",
       "304678  gi|574455772|gb|KF975499|Influenza A virus (A/...   \n",
       "\n",
       "                                                sequences gene  \n",
       "629174  GTGACAAAAACATAATGGATTCCAACACGATAACCTCGTTTCAGGT...  not  \n",
       "610132  AATATGGAGAGAATAAAAGAACTAAGAGATTTGATGTCGCAGTCTC...  not  \n",
       "470806  GTTTGATTTATCATAATGGATACTTTTATTACAAGAAACTTCCAGA...  not  \n",
       "278994  AATAACGATTGGCTCTGTTTCTCTCACCATTTCCACAATATGCTTC...  not  \n",
       "86479   ATTCAATATGGAAAGAATAAAAGAATTAAGAGATCTAATGTCACAG...  not  \n",
       "151351  TACTGATTCAAAATGGAAGACTTTGTGCGACAATGCTTCAATCCAA...  not  \n",
       "412958  ATGTCGCTGTTTGGAGACACAATTGCCTACCTGCTTTCATTGACAG...  not  \n",
       "315125  ATGAATATAAATCCTTATTTTCTCTTCATAGATGTACCCATACAGG...  not  \n",
       "74069   ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...  not  \n",
       "505337  ATGGATTCCAACACGATAACCTCGTTTCAGGTAGATTGTTATCTAT...  not  \n",
       "340273  ATGAAGACTATCATTGCTTTGAGCTACACTCTATGTCTGGTTTTCG...  not  \n",
       "556365  GTTTGATTTGTCATAATGGATACTTTTATTACAAGAAACTTCCAGA...  not  \n",
       "623475  ATTTTCTAATATCCACAGAATGAAGGCAATAATTGTACTACTCATG...  not  \n",
       "273912  ATGCTACCTTCAACTATACAAACGTTAACCCTATTTCTCACATCAG...  not  \n",
       "636364  GGTCACAATGGAAAAATTCATCATTTTGAGCACTGTCCTGGCAGTA...  not  \n",
       "467742  ATGTCTCTTCTATTGACAATAGCAAAGGAATACAAAAGACTATGCC...  not  \n",
       "156141  ATGAGTGACATCGAAGCCATGGCGTCTCAAGGCACCAAACGATCAT...  not  \n",
       "165106  ATATTGAAAGATGAGTCTTCTAACCGAGGTCGAAACGTACGTTCTC...  not  \n",
       "371494  ATGAAGACTATCATTGCTTTGAGCTACATTCTATGTCTGGTTTTCG...  not  \n",
       "424904  ATGGACTCCAATACTGTGTCAAGCTTTCAGGTAGACTGTTTCCTTT...  not  \n",
       "111975  GGAGAACGCCAGAATGCTACTGAGATCAGGGCATCTGTTGGAAGAA...  not  \n",
       "561122  ATGAATCCAAATCAGAAGATTCTATGCACTTCAGCCACTGCTATCC...  not  \n",
       "298789  ATGAAGGCAATACTATTAGTCTTGTTATATACATTTGTACCCGCAA...  not  \n",
       "328210  TTATTCGTCTCAGGGAGCAAAAGCAGGGTAGATAATCACTCACTGA...  not  \n",
       "374982  GATCGAATCTGCACTGGGATAACATCTTCAAACTCACCTCATGTGG...  not  \n",
       "389631  ATGAATCCAAATCAGAAGATAATAACAATTGGCTCTGTTTCTCTAT...  not  \n",
       "529407  AGCGAAAGCAGGTCAAATATATTCAATATGGAAAGAATAAAAGAAC...  not  \n",
       "567789  ATGGATGTCAATCCGACCCTACTTTTCTTAAAGGTTCCAGCGCAAA...  not  \n",
       "18713   AGCAAAAGCAGGTCAATTATATTCAGTATGGAAAGAATAAAAGAAC...  not  \n",
       "627939  TCAAATCAACCAAGATGGAAGCAAAACTGTTTGTACTATTCTGTAC...  not  \n",
       "...                                                   ...  ...  \n",
       "355132  AGCAAAAGCAGGTAGATATTGAAAAATGAGTCTTCTAACTGAGGTC...  not  \n",
       "226790  TCAAATATATTCAATATGGAGAGAATAAAGGAATTAAGAGACCTAA...  not  \n",
       "577113  GTGACAAAGACATAATGGATTCCAACACTGTGTCAAGTTTCCAGGT...  not  \n",
       "40450   ATTCGAAATGGAAGATTTTGTGCGACAGTGCTTCAATCCGATGATT...  not  \n",
       "650629  ATGAAGACTATCATTGCTTTGAGCTACATTCTATGTCTGGTTTTCG...  not  \n",
       "196481  CAAAAACTTCCCGGAAATGACAACAGCACGGCAACGCTGTGCCTTG...  not  \n",
       "568332  ATGAAGGCAATAATTGTACTACTCATGGTAGTAACATCCAATGCAG...  not  \n",
       "512342  ATGGATGTCAATCCGACTTTACTTTTCTTGAAAGTTCCAGCGCAAA...  not  \n",
       "526752  AGTAAAGATGAATCCAAATCAAAAGATAATAACGATTGGCTCTGTT...  not  \n",
       "254049  ATGAAGGCAATACTAGTAGTTCTGCTATATACATTTGCAACCGCAA...  not  \n",
       "308720  ATGAATCCAAATCAGAAGATAATAACAATTGGCTCTGTCTCTCTAA...  not  \n",
       "120107  CCCCCTCAAAGCCGAAATCGCGCAGAGACTTGAAGATGTCTTTGCT...  not  \n",
       "182482  CAAACCATTTGAATGGATGTCAATCCGACTTTACTTTTCTTGAAAG...  not  \n",
       "58027   TTTAAAATGAATCCAAATCAAAAAATAATAACCATTGGATCAATCA...  not  \n",
       "390631  ATGGCGTCTCAAGGCACCAAACGATCATATGAACAAATGGAGACTG...  not  \n",
       "427517  ATGGATGTCAATCCGACTCTACTGTTCTTAAAAGTTCCAGCGCAAA...  not  \n",
       "23280   AGCAAAAGCAGGTCAATTATATTCAATATGGAGAGAATAAAAGAAC...  not  \n",
       "180865  TATGGGGAAGATCCGAAAATCGAAACAAACAAATTTGCCGCAATAT...  not  \n",
       "560672  ATGAACACTCAAATCCTGGTATTCGCTCTGATTGCGATCATTCCAA...  not  \n",
       "520172  TACTGATTCAAAATGGAAGATTTTGTGCGACAATGCTTCAACCCGA...  not  \n",
       "546639  AGCAAAAGCAGGTGCGAAATGAATCCAAATCAGAAGATAATAACAA...  not  \n",
       "293165  TACTGATCCAAAATGGAAGACTTTGTGCGACAATGCTTCAATCCAA...  not  \n",
       "604028  AGCAAAAGCAGGGTTAATAATCACTCACTGAGTGACATCAAAATCA...  not  \n",
       "66867   AAGATGAGTCTTCTAACCGAGGTCGAAACGTACGTTCTCTCTATCA...  not  \n",
       "649040  CTGATTCAAAATGGAAGATTTTGTGCGACAATGCTTCAACCCGATG...  not  \n",
       "636950  TACTGATCCAAAATGGAAGACTTTGTGCGACAATGCTTCAATCCAA...  not  \n",
       "341110  ATGGCGTCTCAAGGCACCAAACGGTCTTATGAACAAATGGAAACTG...  not  \n",
       "616949  ATGAAGGCAATAATTGTACTACTCATGGTAGTAACATCCAATGCAG...  not  \n",
       "263002  ATGAAGGCAATAATTGTACTACTCATGGTAGTAACATCCAATGCAG...  not  \n",
       "304678  TGACAAAAACATAATGGATTCCAACACTGTGTCAAGCTTTCAGGTA...  not  \n",
       "\n",
       "[99 rows x 3 columns]"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "uniflna_test1\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "def row_test(row):\n",
    "    if 'HA' in row['titles']:\n",
    "        return('HA')\n",
    "    if 'NA' in row['titles']:\n",
    "        return('NA')\n",
    "    if 'PA' in row['titles']:\n",
    "        return('PA')\n",
    "    if 'PB2' in row['titles']:\n",
    "        return('PB2')\n",
    "    if 'NEP' in row['titles']:\n",
    "        return('NEP')\n",
    "    if 'M1' in row['titles']:\n",
    "        return('M1')\n",
    "    if 'PB1' in row['titles']:\n",
    "        return('PB1')\n",
    "    if 'NP' in row['titles']:\n",
    "        return('NP')\n",
    "    else:\n",
    "        return('uncl')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [],
   "source": [
    "uniflna_test1['gene'] = uniflna_test1.apply(row_test, axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {
    "scrolled": true
   },
   "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>titles</th>\n",
       "      <th>sequences</th>\n",
       "      <th>gene</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>629174</th>\n",
       "      <td>gi|1410306180|gb|MH546639|Influenza A virus (A...</td>\n",
       "      <td>GTGACAAAAACATAATGGATTCCAACACGATAACCTCGTTTCAGGT...</td>\n",
       "      <td>NEP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>610132</th>\n",
       "      <td>gi|1370449901|gb|MG957554|Influenza A virus (A...</td>\n",
       "      <td>AATATGGAGAGAATAAAAGAACTAAGAGATTTGATGTCGCAGTCTC...</td>\n",
       "      <td>PB2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>470806</th>\n",
       "      <td>gi|1048081328|gb|KX616653|Influenza B virus (B...</td>\n",
       "      <td>GTTTGATTTATCATAATGGATACTTTTATTACAAGAAACTTCCAGA...</td>\n",
       "      <td>PA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>278994</th>\n",
       "      <td>gi|525553452|gb|KF314561|Influenza A virus (A/...</td>\n",
       "      <td>AATAACGATTGGCTCTGTTTCTCTCACCATTTCCACAATATGCTTC...</td>\n",
       "      <td>NA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>86479</th>\n",
       "      <td>gi|237874843|gb|FJ985245|Influenza A virus (A/...</td>\n",
       "      <td>ATTCAATATGGAAAGAATAAAAGAATTAAGAGATCTAATGTCACAG...</td>\n",
       "      <td>PB2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>151351</th>\n",
       "      <td>gi|315440734|gb|CY079457|Influenza A virus (A/...</td>\n",
       "      <td>TACTGATTCAAAATGGAAGACTTTGTGCGACAATGCTTCAATCCAA...</td>\n",
       "      <td>uncl</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>412958</th>\n",
       "      <td>gi|937172275|gb|KT854258|Influenza B virus (B/...</td>\n",
       "      <td>ATGTCGCTGTTTGGAGACACAATTGCCTACCTGCTTTCATTGACAG...</td>\n",
       "      <td>M1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>315125</th>\n",
       "      <td>gi|594215428|gb|CY174927|Influenza B virus (B/...</td>\n",
       "      <td>ATGAATATAAATCCTTATTTTCTCTTCATAGATGTACCCATACAGG...</td>\n",
       "      <td>PB1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74069</th>\n",
       "      <td>gi|198285780|gb|CY030920|Influenza A virus (A/...</td>\n",
       "      <td>ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...</td>\n",
       "      <td>uncl</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>505337</th>\n",
       "      <td>gi|1120605462|gb|KY347748|Influenza A virus (A...</td>\n",
       "      <td>ATGGATTCCAACACGATAACCTCGTTTCAGGTAGATTGTTATCTAT...</td>\n",
       "      <td>NEP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>340273</th>\n",
       "      <td>gi|660873930|gb|KM064239|Influenza A virus (A/...</td>\n",
       "      <td>ATGAAGACTATCATTGCTTTGAGCTACACTCTATGTCTGGTTTTCG...</td>\n",
       "      <td>HA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>556365</th>\n",
       "      <td>gi|1219455324|gb|CY239151|Influenza B virus (B...</td>\n",
       "      <td>GTTTGATTTGTCATAATGGATACTTTTATTACAAGAAACTTCCAGA...</td>\n",
       "      <td>PA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>623475</th>\n",
       "      <td>gi|1390209297|gb|MH363175|Influenza B virus (B...</td>\n",
       "      <td>ATTTTCTAATATCCACAGAATGAAGGCAATAATTGTACTACTCATG...</td>\n",
       "      <td>HA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>273912</th>\n",
       "      <td>gi|512772000|gb|KF216939|Influenza B virus (B/...</td>\n",
       "      <td>ATGCTACCTTCAACTATACAAACGTTAACCCTATTTCTCACATCAG...</td>\n",
       "      <td>NA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>636364</th>\n",
       "      <td>gi|1421412631|gb|MH598206|Influenza A virus (A...</td>\n",
       "      <td>GGTCACAATGGAAAAATTCATCATTTTGAGCACTGTCCTGGCAGTA...</td>\n",
       "      <td>HA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>467742</th>\n",
       "      <td>gi|1041633881|gb|LC123364|Influenza C virus (C...</td>\n",
       "      <td>ATGTCTCTTCTATTGACAATAGCAAAGGAATACAAAAGACTATGCC...</td>\n",
       "      <td>PB2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>156141</th>\n",
       "      <td>gi|324345671|gb|CY084441|Influenza A virus (A/...</td>\n",
       "      <td>ATGAGTGACATCGAAGCCATGGCGTCTCAAGGCACCAAACGATCAT...</td>\n",
       "      <td>NP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>165106</th>\n",
       "      <td>gi|336286441|gb|JN087059|Influenza A virus (A/...</td>\n",
       "      <td>ATATTGAAAGATGAGTCTTCTAACCGAGGTCGAAACGTACGTTCTC...</td>\n",
       "      <td>M1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>371494</th>\n",
       "      <td>gi|755195545|gb|KP456628|Influenza A virus (A/...</td>\n",
       "      <td>ATGAAGACTATCATTGCTTTGAGCTACATTCTATGTCTGGTTTTCG...</td>\n",
       "      <td>HA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>424904</th>\n",
       "      <td>gi|974741015|gb|CY195307|Influenza A virus (A/...</td>\n",
       "      <td>ATGGACTCCAATACTGTGTCAAGCTTTCAGGTAGACTGTTTCCTTT...</td>\n",
       "      <td>NEP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>111975</th>\n",
       "      <td>gi|283551383|gb|CY047996|Influenza A virus (A/...</td>\n",
       "      <td>GGAGAACGCCAGAATGCTACTGAGATCAGGGCATCTGTTGGAAGAA...</td>\n",
       "      <td>uncl</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>561122</th>\n",
       "      <td>gi|1237086641|gb|MF630455|Influenza A virus (A...</td>\n",
       "      <td>ATGAATCCAAATCAGAAGATTCTATGCACTTCAGCCACTGCTATCC...</td>\n",
       "      <td>NA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>298789</th>\n",
       "      <td>gi|563400421|gb|CY158761|Influenza A virus (A/...</td>\n",
       "      <td>ATGAAGGCAATACTATTAGTCTTGTTATATACATTTGTACCCGCAA...</td>\n",
       "      <td>HA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>328210</th>\n",
       "      <td>gi|607259322|gb|KJ579959|Influenza A virus (A/...</td>\n",
       "      <td>TTATTCGTCTCAGGGAGCAAAAGCAGGGTAGATAATCACTCACTGA...</td>\n",
       "      <td>NP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>374982</th>\n",
       "      <td>gi|755202521|gb|KP460116|Influenza B virus (B/...</td>\n",
       "      <td>GATCGAATCTGCACTGGGATAACATCTTCAAACTCACCTCATGTGG...</td>\n",
       "      <td>HA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>389631</th>\n",
       "      <td>gi|827499629|gb|KR701147|Influenza A virus (A/...</td>\n",
       "      <td>ATGAATCCAAATCAGAAGATAATAACAATTGGCTCTGTTTCTCTAT...</td>\n",
       "      <td>NA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>529407</th>\n",
       "      <td>gi|1176711220|gb|KY855523|Influenza A virus (A...</td>\n",
       "      <td>AGCGAAAGCAGGTCAAATATATTCAATATGGAAAGAATAAAAGAAC...</td>\n",
       "      <td>PB2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>567789</th>\n",
       "      <td>gi|1246295632|gb|MF801575|Influenza A virus (A...</td>\n",
       "      <td>ATGGATGTCAATCCGACCCTACTTTTCTTAAAGGTTCCAGCGCAAA...</td>\n",
       "      <td>PB1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18713</th>\n",
       "      <td>gi|77747051|gb|CY003463|Influenza A virus (A/N...</td>\n",
       "      <td>AGCAAAAGCAGGTCAATTATATTCAGTATGGAAAGAATAAAAGAAC...</td>\n",
       "      <td>uncl</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>627939</th>\n",
       "      <td>gi|1408067743|gb|MH341870|Influenza A virus (A...</td>\n",
       "      <td>TCAAATCAACCAAGATGGAAGCAAAACTGTTTGTACTATTCTGTAC...</td>\n",
       "      <td>HA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4629</th>\n",
       "      <td>gi|15555102|gb|AX225029|Sequence 39 from Paten...</td>\n",
       "      <td>AGCAAAAGCAGGGTGACAAAAACATAATGGATTCCAACACTGTGTC...</td>\n",
       "      <td>uncl</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>388575</th>\n",
       "      <td>gi|827497043|gb|KR700091|Influenza A virus (A/...</td>\n",
       "      <td>ATGAATCCAAATCAAAAGATAATAACAATTGGCTCTGTCTCCCTAA...</td>\n",
       "      <td>NA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>400891</th>\n",
       "      <td>gi|923095095|gb|KT374233|Influenza A virus (A/...</td>\n",
       "      <td>ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...</td>\n",
       "      <td>PA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>495092</th>\n",
       "      <td>gi|1102154012|gb|KX979673|Influenza A virus (A...</td>\n",
       "      <td>AGCAAAAGCAGGCAAACCATTTGAATGGATGTCAATCCGACTTTAC...</td>\n",
       "      <td>PB1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>49282</th>\n",
       "      <td>gi|153972782|gb|EU050209|Influenza A virus (A/...</td>\n",
       "      <td>CACTCACTGAGTGACATTCACGCCATGGCGCTTCAAGGCACCAAAC...</td>\n",
       "      <td>NP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>76167</th>\n",
       "      <td>gi|209866820|gb|CY035609|Influenza A virus (A/...</td>\n",
       "      <td>ATGATCGTCGAGCTTGCGGAAAAGGCAATGAAAGAATATGGGGAGG...</td>\n",
       "      <td>uncl</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>164042</th>\n",
       "      <td>gi|334189508|gb|CY091572|Influenza A virus (A/...</td>\n",
       "      <td>ATTATATTCAGTATGGAAAGAATAAAAGAACTACGGAATCTGATGT...</td>\n",
       "      <td>PB2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>639802</th>\n",
       "      <td>gi|1423484309|gb|MH606905|Influenza B virus (B...</td>\n",
       "      <td>ATTTTCTTGTGAACTTCAAGCACCAGTGAAAGAACTGAAAATCAAA...</td>\n",
       "      <td>NP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>432587</th>\n",
       "      <td>gi|1000373120|gb|CY197474|Influenza B virus (B...</td>\n",
       "      <td>TTTCTTAAAATGTCGCTGTTTGGAGACACAATTGCCTACCTGCTTT...</td>\n",
       "      <td>M1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>396439</th>\n",
       "      <td>gi|908275253|gb|KT283502|Influenza B virus (B/...</td>\n",
       "      <td>ATGAAGGCAATAATTGTACTACTCATGGTAGTAACATCCAATGCAG...</td>\n",
       "      <td>HA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26965</th>\n",
       "      <td>gi|90571769|gb|CY010022|Influenza A virus (A/N...</td>\n",
       "      <td>AATCCAAATCAAAAGATAATAACGATTGGCTCTGTTTCTCTCACTA...</td>\n",
       "      <td>uncl</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>314234</th>\n",
       "      <td>gi|594211014|gb|CY174083|Influenza B virus (B/...</td>\n",
       "      <td>CTTCTCAAAACTGAGGCAAATAGGCCAAAAATGAACAATGCTACCT...</td>\n",
       "      <td>NA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63374</th>\n",
       "      <td>gi|183222959|gb|EU555227|Influenza B virus (B/...</td>\n",
       "      <td>ATTGTACTACTCATGGTAGTAACATCCAATGCAGATCGAATCTGCA...</td>\n",
       "      <td>uncl</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>60259</th>\n",
       "      <td>gi|167735711|gb|AB378686|Influenza A virus (A/...</td>\n",
       "      <td>CAAATTTCTCAATTCTATTTGATCAACTAGAGACGTTAGTCTCTCT...</td>\n",
       "      <td>uncl</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>170937</th>\n",
       "      <td>gi|341610267|gb|JN244211|Influenza A virus (A/...</td>\n",
       "      <td>TATTTGCATTATCTGGAGTGGCAGTAGCACTTAGTGTACTGAACTT...</td>\n",
       "      <td>NA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>122382</th>\n",
       "      <td>gi|294861120|gb|CY061602|Influenza A virus (A/...</td>\n",
       "      <td>ACAAAATGAATACTCAAATTTTAATACTCACTCTTGTGGCGGTCAT...</td>\n",
       "      <td>HA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4695</th>\n",
       "      <td>gi|16076700|gb|AF222811|Influenza A virus (A/S...</td>\n",
       "      <td>ATGGAAGCAGAATCACTAATAACTATACTACTAGTAGTAACAGCAA...</td>\n",
       "      <td>uncl</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>263090</th>\n",
       "      <td>gi|488465200|gb|KC891811|Influenza B virus (B/...</td>\n",
       "      <td>ATGAAGGCAATAATTGTACTACTAATGGTAGTAACATCCAATGCAG...</td>\n",
       "      <td>HA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>568420</th>\n",
       "      <td>gi|1248280161|gb|MF927815|Influenza A virus (A...</td>\n",
       "      <td>ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...</td>\n",
       "      <td>PA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>479411</th>\n",
       "      <td>gi|1072743785|gb|KX851867|Influenza A virus (A...</td>\n",
       "      <td>ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...</td>\n",
       "      <td>PA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>532585</th>\n",
       "      <td>gi|1181553196|gb|KY935553|Influenza A virus (A...</td>\n",
       "      <td>AAAGCAGGGTGACAAAAACATAATGGAATCCAACACCATGTCAAGC...</td>\n",
       "      <td>NEP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>493894</th>\n",
       "      <td>gi|1101014193|gb|KY089804|Influenza A virus (A...</td>\n",
       "      <td>ATGAATCCAAATCAAAAGATAATAACGATTGGCTCTGTTTCTCTCA...</td>\n",
       "      <td>NA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>286402</th>\n",
       "      <td>gi|549395916|gb|CY150729|Influenza B virus (B/...</td>\n",
       "      <td>TTTCTAATATCCACAAAATGAAGGCAATAATTGTACTACTCATGGT...</td>\n",
       "      <td>HA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>216098</th>\n",
       "      <td>gi|396578002|gb|AB735609|Influenza A virus (A/...</td>\n",
       "      <td>ATCCAAATCAAAAGATAACATGCATTTCAGCAACAGGAATGACACT...</td>\n",
       "      <td>NA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>628780</th>\n",
       "      <td>gi|1410305195|gb|MH546245|Influenza A virus (A...</td>\n",
       "      <td>AGTGAAAATGAATCCAAATCAGAAAATAATAACAATTGGCTCTGTC...</td>\n",
       "      <td>NA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18303</th>\n",
       "      <td>gi|76574387|gb|DQ145544|Influenza A virus (A/s...</td>\n",
       "      <td>ATGGATGTCAATCCGACCCTACTTTTCCTAAAAGTTCCAGCGCAAA...</td>\n",
       "      <td>PB1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>293124</th>\n",
       "      <td>gi|552052068|gb|KF695274|Influenza A virus (A/...</td>\n",
       "      <td>AGCAGGGTGACAAAGACATAATGGATTCCAACACGATAACCTCGTT...</td>\n",
       "      <td>NEP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>186510</th>\n",
       "      <td>gi|370321083|gb|EU026031|Influenza A virus (A/...</td>\n",
       "      <td>TCAAAATGAATCCAAATCAAAAGATAATAACTATTGGGTCAATCTG...</td>\n",
       "      <td>NA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>379609</th>\n",
       "      <td>gi|766559216|gb|KP767578|Influenza A virus (A/...</td>\n",
       "      <td>ACTGTGCCTAGGACACCATGCAGTATCGAATGGGACAATAGTGAAG...</td>\n",
       "      <td>HA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>580342</th>\n",
       "      <td>gi|1277263174|gb|CY254736|Influenza A virus (A...</td>\n",
       "      <td>ATTATATTCAGTATGGAAAGAATAAAAGAACTACGGAATCTGATGT...</td>\n",
       "      <td>PB2</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1000 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                   titles  \\\n",
       "629174  gi|1410306180|gb|MH546639|Influenza A virus (A...   \n",
       "610132  gi|1370449901|gb|MG957554|Influenza A virus (A...   \n",
       "470806  gi|1048081328|gb|KX616653|Influenza B virus (B...   \n",
       "278994  gi|525553452|gb|KF314561|Influenza A virus (A/...   \n",
       "86479   gi|237874843|gb|FJ985245|Influenza A virus (A/...   \n",
       "151351  gi|315440734|gb|CY079457|Influenza A virus (A/...   \n",
       "412958  gi|937172275|gb|KT854258|Influenza B virus (B/...   \n",
       "315125  gi|594215428|gb|CY174927|Influenza B virus (B/...   \n",
       "74069   gi|198285780|gb|CY030920|Influenza A virus (A/...   \n",
       "505337  gi|1120605462|gb|KY347748|Influenza A virus (A...   \n",
       "340273  gi|660873930|gb|KM064239|Influenza A virus (A/...   \n",
       "556365  gi|1219455324|gb|CY239151|Influenza B virus (B...   \n",
       "623475  gi|1390209297|gb|MH363175|Influenza B virus (B...   \n",
       "273912  gi|512772000|gb|KF216939|Influenza B virus (B/...   \n",
       "636364  gi|1421412631|gb|MH598206|Influenza A virus (A...   \n",
       "467742  gi|1041633881|gb|LC123364|Influenza C virus (C...   \n",
       "156141  gi|324345671|gb|CY084441|Influenza A virus (A/...   \n",
       "165106  gi|336286441|gb|JN087059|Influenza A virus (A/...   \n",
       "371494  gi|755195545|gb|KP456628|Influenza A virus (A/...   \n",
       "424904  gi|974741015|gb|CY195307|Influenza A virus (A/...   \n",
       "111975  gi|283551383|gb|CY047996|Influenza A virus (A/...   \n",
       "561122  gi|1237086641|gb|MF630455|Influenza A virus (A...   \n",
       "298789  gi|563400421|gb|CY158761|Influenza A virus (A/...   \n",
       "328210  gi|607259322|gb|KJ579959|Influenza A virus (A/...   \n",
       "374982  gi|755202521|gb|KP460116|Influenza B virus (B/...   \n",
       "389631  gi|827499629|gb|KR701147|Influenza A virus (A/...   \n",
       "529407  gi|1176711220|gb|KY855523|Influenza A virus (A...   \n",
       "567789  gi|1246295632|gb|MF801575|Influenza A virus (A...   \n",
       "18713   gi|77747051|gb|CY003463|Influenza A virus (A/N...   \n",
       "627939  gi|1408067743|gb|MH341870|Influenza A virus (A...   \n",
       "...                                                   ...   \n",
       "4629    gi|15555102|gb|AX225029|Sequence 39 from Paten...   \n",
       "388575  gi|827497043|gb|KR700091|Influenza A virus (A/...   \n",
       "400891  gi|923095095|gb|KT374233|Influenza A virus (A/...   \n",
       "495092  gi|1102154012|gb|KX979673|Influenza A virus (A...   \n",
       "49282   gi|153972782|gb|EU050209|Influenza A virus (A/...   \n",
       "76167   gi|209866820|gb|CY035609|Influenza A virus (A/...   \n",
       "164042  gi|334189508|gb|CY091572|Influenza A virus (A/...   \n",
       "639802  gi|1423484309|gb|MH606905|Influenza B virus (B...   \n",
       "432587  gi|1000373120|gb|CY197474|Influenza B virus (B...   \n",
       "396439  gi|908275253|gb|KT283502|Influenza B virus (B/...   \n",
       "26965   gi|90571769|gb|CY010022|Influenza A virus (A/N...   \n",
       "314234  gi|594211014|gb|CY174083|Influenza B virus (B/...   \n",
       "63374   gi|183222959|gb|EU555227|Influenza B virus (B/...   \n",
       "60259   gi|167735711|gb|AB378686|Influenza A virus (A/...   \n",
       "170937  gi|341610267|gb|JN244211|Influenza A virus (A/...   \n",
       "122382  gi|294861120|gb|CY061602|Influenza A virus (A/...   \n",
       "4695    gi|16076700|gb|AF222811|Influenza A virus (A/S...   \n",
       "263090  gi|488465200|gb|KC891811|Influenza B virus (B/...   \n",
       "568420  gi|1248280161|gb|MF927815|Influenza A virus (A...   \n",
       "479411  gi|1072743785|gb|KX851867|Influenza A virus (A...   \n",
       "532585  gi|1181553196|gb|KY935553|Influenza A virus (A...   \n",
       "493894  gi|1101014193|gb|KY089804|Influenza A virus (A...   \n",
       "286402  gi|549395916|gb|CY150729|Influenza B virus (B/...   \n",
       "216098  gi|396578002|gb|AB735609|Influenza A virus (A/...   \n",
       "628780  gi|1410305195|gb|MH546245|Influenza A virus (A...   \n",
       "18303   gi|76574387|gb|DQ145544|Influenza A virus (A/s...   \n",
       "293124  gi|552052068|gb|KF695274|Influenza A virus (A/...   \n",
       "186510  gi|370321083|gb|EU026031|Influenza A virus (A/...   \n",
       "379609  gi|766559216|gb|KP767578|Influenza A virus (A/...   \n",
       "580342  gi|1277263174|gb|CY254736|Influenza A virus (A...   \n",
       "\n",
       "                                                sequences  gene  \n",
       "629174  GTGACAAAAACATAATGGATTCCAACACGATAACCTCGTTTCAGGT...   NEP  \n",
       "610132  AATATGGAGAGAATAAAAGAACTAAGAGATTTGATGTCGCAGTCTC...   PB2  \n",
       "470806  GTTTGATTTATCATAATGGATACTTTTATTACAAGAAACTTCCAGA...    PA  \n",
       "278994  AATAACGATTGGCTCTGTTTCTCTCACCATTTCCACAATATGCTTC...    NA  \n",
       "86479   ATTCAATATGGAAAGAATAAAAGAATTAAGAGATCTAATGTCACAG...   PB2  \n",
       "151351  TACTGATTCAAAATGGAAGACTTTGTGCGACAATGCTTCAATCCAA...  uncl  \n",
       "412958  ATGTCGCTGTTTGGAGACACAATTGCCTACCTGCTTTCATTGACAG...    M1  \n",
       "315125  ATGAATATAAATCCTTATTTTCTCTTCATAGATGTACCCATACAGG...   PB1  \n",
       "74069   ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...  uncl  \n",
       "505337  ATGGATTCCAACACGATAACCTCGTTTCAGGTAGATTGTTATCTAT...   NEP  \n",
       "340273  ATGAAGACTATCATTGCTTTGAGCTACACTCTATGTCTGGTTTTCG...    HA  \n",
       "556365  GTTTGATTTGTCATAATGGATACTTTTATTACAAGAAACTTCCAGA...    PA  \n",
       "623475  ATTTTCTAATATCCACAGAATGAAGGCAATAATTGTACTACTCATG...    HA  \n",
       "273912  ATGCTACCTTCAACTATACAAACGTTAACCCTATTTCTCACATCAG...    NA  \n",
       "636364  GGTCACAATGGAAAAATTCATCATTTTGAGCACTGTCCTGGCAGTA...    HA  \n",
       "467742  ATGTCTCTTCTATTGACAATAGCAAAGGAATACAAAAGACTATGCC...   PB2  \n",
       "156141  ATGAGTGACATCGAAGCCATGGCGTCTCAAGGCACCAAACGATCAT...    NP  \n",
       "165106  ATATTGAAAGATGAGTCTTCTAACCGAGGTCGAAACGTACGTTCTC...    M1  \n",
       "371494  ATGAAGACTATCATTGCTTTGAGCTACATTCTATGTCTGGTTTTCG...    HA  \n",
       "424904  ATGGACTCCAATACTGTGTCAAGCTTTCAGGTAGACTGTTTCCTTT...   NEP  \n",
       "111975  GGAGAACGCCAGAATGCTACTGAGATCAGGGCATCTGTTGGAAGAA...  uncl  \n",
       "561122  ATGAATCCAAATCAGAAGATTCTATGCACTTCAGCCACTGCTATCC...    NA  \n",
       "298789  ATGAAGGCAATACTATTAGTCTTGTTATATACATTTGTACCCGCAA...    HA  \n",
       "328210  TTATTCGTCTCAGGGAGCAAAAGCAGGGTAGATAATCACTCACTGA...    NP  \n",
       "374982  GATCGAATCTGCACTGGGATAACATCTTCAAACTCACCTCATGTGG...    HA  \n",
       "389631  ATGAATCCAAATCAGAAGATAATAACAATTGGCTCTGTTTCTCTAT...    NA  \n",
       "529407  AGCGAAAGCAGGTCAAATATATTCAATATGGAAAGAATAAAAGAAC...   PB2  \n",
       "567789  ATGGATGTCAATCCGACCCTACTTTTCTTAAAGGTTCCAGCGCAAA...   PB1  \n",
       "18713   AGCAAAAGCAGGTCAATTATATTCAGTATGGAAAGAATAAAAGAAC...  uncl  \n",
       "627939  TCAAATCAACCAAGATGGAAGCAAAACTGTTTGTACTATTCTGTAC...    HA  \n",
       "...                                                   ...   ...  \n",
       "4629    AGCAAAAGCAGGGTGACAAAAACATAATGGATTCCAACACTGTGTC...  uncl  \n",
       "388575  ATGAATCCAAATCAAAAGATAATAACAATTGGCTCTGTCTCCCTAA...    NA  \n",
       "400891  ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...    PA  \n",
       "495092  AGCAAAAGCAGGCAAACCATTTGAATGGATGTCAATCCGACTTTAC...   PB1  \n",
       "49282   CACTCACTGAGTGACATTCACGCCATGGCGCTTCAAGGCACCAAAC...    NP  \n",
       "76167   ATGATCGTCGAGCTTGCGGAAAAGGCAATGAAAGAATATGGGGAGG...  uncl  \n",
       "164042  ATTATATTCAGTATGGAAAGAATAAAAGAACTACGGAATCTGATGT...   PB2  \n",
       "639802  ATTTTCTTGTGAACTTCAAGCACCAGTGAAAGAACTGAAAATCAAA...    NP  \n",
       "432587  TTTCTTAAAATGTCGCTGTTTGGAGACACAATTGCCTACCTGCTTT...    M1  \n",
       "396439  ATGAAGGCAATAATTGTACTACTCATGGTAGTAACATCCAATGCAG...    HA  \n",
       "26965   AATCCAAATCAAAAGATAATAACGATTGGCTCTGTTTCTCTCACTA...  uncl  \n",
       "314234  CTTCTCAAAACTGAGGCAAATAGGCCAAAAATGAACAATGCTACCT...    NA  \n",
       "63374   ATTGTACTACTCATGGTAGTAACATCCAATGCAGATCGAATCTGCA...  uncl  \n",
       "60259   CAAATTTCTCAATTCTATTTGATCAACTAGAGACGTTAGTCTCTCT...  uncl  \n",
       "170937  TATTTGCATTATCTGGAGTGGCAGTAGCACTTAGTGTACTGAACTT...    NA  \n",
       "122382  ACAAAATGAATACTCAAATTTTAATACTCACTCTTGTGGCGGTCAT...    HA  \n",
       "4695    ATGGAAGCAGAATCACTAATAACTATACTACTAGTAGTAACAGCAA...  uncl  \n",
       "263090  ATGAAGGCAATAATTGTACTACTAATGGTAGTAACATCCAATGCAG...    HA  \n",
       "568420  ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...    PA  \n",
       "479411  ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...    PA  \n",
       "532585  AAAGCAGGGTGACAAAAACATAATGGAATCCAACACCATGTCAAGC...   NEP  \n",
       "493894  ATGAATCCAAATCAAAAGATAATAACGATTGGCTCTGTTTCTCTCA...    NA  \n",
       "286402  TTTCTAATATCCACAAAATGAAGGCAATAATTGTACTACTCATGGT...    HA  \n",
       "216098  ATCCAAATCAAAAGATAACATGCATTTCAGCAACAGGAATGACACT...    NA  \n",
       "628780  AGTGAAAATGAATCCAAATCAGAAAATAATAACAATTGGCTCTGTC...    NA  \n",
       "18303   ATGGATGTCAATCCGACCCTACTTTTCCTAAAAGTTCCAGCGCAAA...   PB1  \n",
       "293124  AGCAGGGTGACAAAGACATAATGGATTCCAACACGATAACCTCGTT...   NEP  \n",
       "186510  TCAAAATGAATCCAAATCAAAAGATAATAACTATTGGGTCAATCTG...    NA  \n",
       "379609  ACTGTGCCTAGGACACCATGCAGTATCGAATGGGACAATAGTGAAG...    HA  \n",
       "580342  ATTATATTCAGTATGGAAAGAATAAAAGAACTACGGAATCTGATGT...   PB2  \n",
       "\n",
       "[1000 rows x 3 columns]"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "uniflna_test1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {
    "collapsed": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['gi|1410306180|gb|MH546639|Influenza A virus (A/mallard/Minnesota/15-040386-2/2015(H11N1)) segment 8 nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|1370449901|gb|MG957554|Influenza A virus (A/duck/Bangladesh/32003/2017(H5N1)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|1048081328|gb|KX616653|Influenza B virus (B/Texas/62/2016) segment 3 polymerase PA (PA) gene, complete cds',\n",
       "       'gi|525553452|gb|KF314561|Influenza A virus (A/India/K0914850/2009(H3N2)) segment 6 neuraminidase (NA) gene, partial cds',\n",
       "       'gi|237874843|gb|FJ985245|Influenza A virus (A/turkey/France/05045/2005(H1N2)) segment 1 polymerase PB2 (PB2) gene, partial cds',\n",
       "       'gi|315440734|gb|CY079457|Influenza A virus (A/American black duck/Illinois/08OS2688/2008(H5N2)) segment 3, complete sequence',\n",
       "       'gi|937172275|gb|KT854258|Influenza B virus (B/Texas/09/2015) segment 7 matrix protein 1 (M1) and BM2 protein (BM2) genes, complete cds',\n",
       "       'gi|594215428|gb|CY174927|Influenza B virus (B/New York/1157/2008) polymerase PB1 (PB1) gene, complete cds',\n",
       "       'gi|198285780|gb|CY030920|Influenza A virus (A/goose/Yunnan/4389/2006(H5N1)) segment 3 sequence',\n",
       "       'gi|1120605462|gb|KY347748|Influenza A virus (A/wild bird feces/Anhui/S93/2014(H9N2)) segment 8 nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|660873930|gb|KM064239|Influenza A virus (A/Washington/10/2014(H3N2)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|1219455324|gb|CY239151|Influenza B virus (B/Michigan/78/2017) polymerase PA (PA) gene, complete cds',\n",
       "       'gi|1390209297|gb|MH363175|Influenza B virus (B/Louisiana/05/2018) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|512772000|gb|KF216939|Influenza B virus (B/Georgia/06/2013) segment 6 neuraminidase (NA) gene, complete cds; and NB protein (NB) gene, partial cds',\n",
       "       'gi|1421412631|gb|MH598206|Influenza A virus (A/ruddy turnstone/New Jersey/AI00-2296/2000(H12N5)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|1041633881|gb|LC123364|Influenza C virus (C/Yamagata/1/2012) PB2 gene for polymerase 2, complete cds',\n",
       "       'gi|324345671|gb|CY084441|Influenza A virus (A/New York/6064/2009(H1N1)) nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|336286441|gb|JN087059|Influenza A virus (A/environment/Korea/ESD3-3/2004(H3N8)) segment 7 matrix protein 2 (M2) gene, partial cds; and matrix protein 1 (M1) gene, complete cds',\n",
       "       'gi|755195545|gb|KP456628|Influenza A virus (A/Alaska/29/2009(H3N2)) segment 4 hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|974741015|gb|CY195307|Influenza A virus (A/swine/Manitoba/D0349/2014(H3N2)) nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|283551383|gb|CY047996|Influenza A virus (A/chicken/Nigeria/08VIR4337-377/2006(H5N1)) segment 5 sequence',\n",
       "       'gi|1237086641|gb|MF630455|Influenza A virus (A/duck/Hunan/S11682/2015(H7N9)) segment 6 neuraminidase (NA) gene, complete cds',\n",
       "       'gi|563400421|gb|CY158761|Influenza A virus (A/swine/Manitoba/02947/2010(H1N1)) hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|607259322|gb|KJ579959|Influenza A virus (A/pigeon/Anhui/08/2013(H3N8)) segment 5 nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|755202521|gb|KP460116|Influenza B virus (B/Kochi/1/2008) segment 4 hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|827499629|gb|KR701147|Influenza A virus (A/swine/Netherlands/Wanroij-CVI2117A/2012(H1N2)) segment 6 neuraminidase (NA) gene, complete cds',\n",
       "       'gi|1176711220|gb|KY855523|Influenza A virus (A/chicken/Guangdong/J2/2017(H7N9)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|1246295632|gb|MF801575|Influenza A virus (A/swine/Iowa/A01667088/2017(H1N1)) segment 2 polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|77747051|gb|CY003463|Influenza A virus (A/New York/439/2000(H3N2)) segment 1, complete sequence',\n",
       "       'gi|1408067743|gb|MH341870|Influenza A virus (A/mallard/Minnesota/UGAI15-6021/2015(H1N1)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|336286770|gb|JN087168|Influenza A virus (A/duck/Korea/U12/2007(H3N2)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|946524900|gb|KT889009|Influenza A virus (A/Philippines/FLU_AFPA_00013/2011(H3N2)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|1390210015|gb|MH363477|Influenza B virus (B/Texas/19/2018) segment 1 polymerase PB1 (PB1) gene, complete cds',\n",
       "       'gi|1435279360|gb|MH560176|Influenza A virus (A/chicken/Vietnam/H7F-LC4-371/2014(H9N2)) segment 8 nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|73665566|gb|CY000866|Influenza A virus (A/New York/193/2003(H3N2)) segment 7, complete sequence',\n",
       "       'gi|56785893|gb|AY738455|Influenza A virus (A/chicken/Talmei Elazar/1304/03(H9N2)) hemagglutinin precursor, gene, complete cds',\n",
       "       'gi|297632072|gb|AB537453|Influenza A virus (A/duck/Shiga/10413-24/2006(H3N8)) HA gene for hemagglutinin, complete cds',\n",
       "       'gi|336286387|gb|JN087276|Influenza A virus (A/duck/Korea/GJ87/2007(H3N8)) segment 8 nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|332079329|gb|JF789635|Influenza A virus (A/Astana/S2/2009(H1N1)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|393712449|gb|CY124530|Influenza A virus (A/Singapore/GP39/2011(H1N1)) polymerase PB1 (PB1) gene, complete cds; and PB1-F2 gene, complete sequence',\n",
       "       'gi|1073679892|gb|KX919563|Influenza A virus (A/Pennsylvania/12/2016(H1N1)) segment 3 polymerase PA (PA) and PA-X protein (PA-X) genes, complete cds',\n",
       "       'gi|660873760|gb|KM064199|Influenza A virus (A/Massachusetts/16/2013(H3N2)) segment 6 neuraminidase (NA) gene, complete cds',\n",
       "       'gi|37785075|gb|AY210036|Influenza A virus (A/Panama/1/61(H2N2)) membrane ion channel M2 and matrix protein M1 (M) gene, complete cds',\n",
       "       'gi|1060100525|gb|KX684034|Influenza A virus (A/swine/Iowa/A01777664/2016(H3N2)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|82654283|gb|CY005452|Influenza A virus (A/Chicken/Nanchang/7-010/2000(H3N6)) segment 2, complete sequence',\n",
       "       'gi|469935292|gb|KC738250|Influenza A virus (A/Taiwan/1151/1992(H3N2)) segment 2 polymerase PB1 (PB1) gene, partial cds; and PB1-F2 protein (PB1-F2) gene, complete cds',\n",
       "       'gi|383278558|gb|CY113883|Influenza A virus (A/Paris/614/1992(H3N2)) polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|1028349584|gb|KX005595|Influenza A virus (A/Texas/07/2016(H1N1)) clone 3000416252_285445_v1_5 segment 5 nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|1128979015|gb|CY212029|Influenza A virus (A/Wisconsin/97/2016(H3N2)) matrix protein 2 (M2) and matrix protein 1 (M1) genes, complete cds',\n",
       "       'gi|169822690|gb|CY030784|Influenza B virus (B/Tennessee/UR06-0079/2007) segment 7, complete sequence',\n",
       "       'gi|910300283|gb|CY191631|Influenza B virus (B/New York/WC-LVD-14-037/2014) nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|1009090830|gb|KU877383|Influenza A virus (A/swine/Illinois/A01731417/2016(H1N2)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|549434742|gb|CY156096|Influenza B virus (B/Wellington/9/2007) polymerase PB1 (PB1) gene, complete cds',\n",
       "       'gi|745630148|gb|KP285522|Influenza A virus (A/duck/Jiangxi/33641/2013(H10N3)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|93008624|gb|DQ493099|Influenza A virus (A/chicken/Kupang-1-NTT/BPPV6/2004(H5N1)) nucleoprotein (NP) gene, partial cds',\n",
       "       'gi|1181550534|gb|KY828763|Influenza A virus (A/northern shoveler/California/LDC185/2014(H8N4)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|452832157|gb|CY134861|Influenza A virus (A/Texas/JMM_19/2012(H3N2)) matrix protein 2 (M2) and matrix protein 1 (M1) genes, complete cds',\n",
       "       'gi|1168686406|gb|CY223564|Influenza A virus (A/Oklahoma/08/2017(H1N1)) polymerase PA (PA) and PA-X protein (PA-X) genes, complete cds',\n",
       "       'gi|1162229649|gb|KY678589|Influenza A virus (A/swine/Minnesota/A01781391/2016(H1N1)) segment 5 nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|1015396647|gb|KU982582|Influenza A virus (A/swine/Hong Kong/9/98-MA(H9N2)) escape mutant RAm8C4(1) segment 4 hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|310699813|gb|CY076324|Influenza A virus (A/northern shoveler/Washington/44249-783/2006(H7N3)) polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|1487102559|gb|MK025712|Influenza A virus (A/Florida/77/2018(H1N1)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|984688745|gb|KU589729|Influenza A virus (A/Ohio/28/2015(H3N2)) segment 2 polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|303303804|gb|CY070957|Influenza A virus (A/Managua/1181.03/2010(H3N2)) segment 2, complete sequence',\n",
       "       'gi|1215244063|gb|CY236525|Influenza B virus (B/Utah/08/2017) polymerase PA (PA) gene, complete cds',\n",
       "       'gi|27465228|gb|AY180404|Influenza A virus (A/Bantam/Nanchang/9-366/2000(H3N3)) hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|974988724|gb|CY196181|Influenza A virus (A/American green-winged teal/Interior Alaska/11BM03479/2011(H4N6)) nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|937145279|gb|KT836974|Influenza A virus (A/Michigan/13/2014(H3N2)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|937151569|gb|KT839367|Influenza A virus (A/Kentucky/14/2014(H3N2)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|727202210|gb|CY103829|Influenza A virus (A/duck/Viet Nam/TMU020/2009(H5N1)) matrix protein 2 (M2) and matrix protein 1 (M1) genes, complete cds',\n",
       "       'gi|411392031|gb|CY126159|Influenza A virus (A/chicken/Qalubia/CAI12/2008(H5N1)) polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|1276993959|gb|CY251499|Influenza A virus (A/Victoria/504/2013(H3N2)) nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds',\n",
       "       'gi|120434185|gb|CY019122|Influenza A virus (A/Memphis/3/1996(H1N1)) segment 3, complete sequence',\n",
       "       'gi|1474769382|gb|MH885281|Influenza A virus (A/Texas/6867/2017(H3N2)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|380836483|gb|CY108700|Influenza A virus (A/Jilin-Chaoyang/133/2007(H3N2)) hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|1248251053|gb|MF955655|Influenza B virus (B/Hong Kong/CUHK37079/2010) segment 4 hemagglutinin gene, complete cds',\n",
       "       'gi|1140316078|gb|KY463977|Influenza A virus (A/mallard/Ohio/14OS2695/2014(mixed)) segment 2 polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|1168684613|gb|CY222823|Influenza A virus (A/Massachusetts/04/2017(H3N2)) neuraminidase (NA) gene, complete cds',\n",
       "       'gi|474455949|gb|KC780219|Influenza A virus (A/Utah/33/2009(H1N1)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|589912749|gb|KJ439872|Influenza A virus (A/duck/Zhejiang/6D7/2013(H3N2)) segment 6 neuraminidase (NA) gene, complete cds',\n",
       "       'gi|293599381|gb|CY059418|Influenza A virus (A/Shenfang/262/1995(H3N2)) segment 7 sequence',\n",
       "       'gi|357382150|gb|CY101657|Influenza A virus (A/mallard/Alberta/220/1990(H3N8)) polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|161138879|gb|CY027957|Influenza A virus (A/Oregon/UR06-0185/2007(H1N1)) segment 6, complete sequence',\n",
       "       'gi|831182692|gb|KR821025|Influenza A virus (A/Mexico/INMEGEN-INER 09/2012(H1N1)) segment 5 nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|984688450|gb|KU589610|Influenza A virus (A/Florida/80/2015(H3N2)) segment 2 polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|83940846|gb|DQ299844|Influenza A virus (A/chicken/Korea/S24/04(H9N2)) polymerase PB2 (PB2) gene, partial cds',\n",
       "       'gi|354681921|gb|AB676809|Influenza A virus (A/goshawk/Tochigi/64/2011(H5N1)) PA gene for polymerase acidic protein, partial cds',\n",
       "       'gi|1237085517|gb|MF630005|Influenza A virus (A/chicken/Guangdong/SC086/2013(H7N9)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|1149046518|gb|CY219243|Influenza A virus (A/California/11/2017(H3N2)) polymerase PA (PA) and PA-X protein (PA-X) genes, complete cds',\n",
       "       'gi|1205299466|gb|MF147727|Influenza A virus (A/mallard duck/Netherlands/24/2011(H2N3)) segment 6 neuraminidase (NA) gene, complete cds',\n",
       "       'gi|552052680|gb|KF695222|Influenza A virus (A/mallard/Sweden/54/2003(H6N1)) segment 3 polymerase PA (PA) and PA-X protein (PA-X) genes, complete cds',\n",
       "       'gi|1350895264|gb|MG856252|Influenza A virus (A/Mexico/InDRE2362/2007(H3N2)) segment 5 nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|188504541|gb|CY031660|Influenza A virus (A/environment/New York/23164-1/2005(H7N2)) segment 7, complete sequence',\n",
       "       'gi|1464622050|gb|MH798603|Influenza A virus (A/New York/PV37/2017(H3N2)) segment 3 polymerase PA (PA) and PA-X protein (PA-X) genes, complete cds',\n",
       "       'gi|1421640142|gb|MH600510|Influenza A virus (A/West Virginia/25/2018(H1N1)) segment 3 polymerase PA (PA) and PA-X protein (PA-X) genes, complete cds',\n",
       "       'gi|662490809|gb|KM027976|Influenza A virus (A/swine/Guangxi/2841/2010(H1N1)) segment 5 nucleocapsid protein (NP) gene, partial cds',\n",
       "       'gi|1378290338|gb|MH216554|Influenza B virus (B/Alberta/006/2017) segment 4 hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|488464962|gb|KC891723|Influenza B virus (B/New Jersey/01/2012) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|574455772|gb|KF975499|Influenza A virus (A/chicken/Pakistan/288-102/2010(H9N2)) segment 8 nuclear export protein (NEP) and nonstructural protein 1 (NS1) genes, complete cds'],\n",
       "      dtype=object)"
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#uniflna_test1.titles.values"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "def seg_test(row):\n",
    "    if 'segment 1' in row['titles']:\n",
    "        return(1)\n",
    "    if 'segment 2' in row['titles']:\n",
    "        return(2)\n",
    "    if 'segment 3' in row['titles']:\n",
    "        return(3)\n",
    "    if 'segment 4' in row['titles']:\n",
    "        return(4)\n",
    "    if 'segment 5' in row['titles']:\n",
    "        return(5)\n",
    "    if 'segment 6' in row['titles']:\n",
    "        return(6)\n",
    "    if 'segment 7' in row['titles']:\n",
    "        return(7)\n",
    "    if 'segment 8' in row['titles']:\n",
    "        return(8)\n",
    "    else:\n",
    "        return(int(0))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [],
   "source": [
    "uniflna_test1['segment'] = uniflna_test1.apply(seg_test, axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "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>titles</th>\n",
       "      <th>sequences</th>\n",
       "      <th>gene</th>\n",
       "      <th>segment</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>629174</th>\n",
       "      <td>gi|1410306180|gb|MH546639|Influenza A virus (A...</td>\n",
       "      <td>GTGACAAAAACATAATGGATTCCAACACGATAACCTCGTTTCAGGT...</td>\n",
       "      <td>NEP</td>\n",
       "      <td>8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>610132</th>\n",
       "      <td>gi|1370449901|gb|MG957554|Influenza A virus (A...</td>\n",
       "      <td>AATATGGAGAGAATAAAAGAACTAAGAGATTTGATGTCGCAGTCTC...</td>\n",
       "      <td>PB2</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>470806</th>\n",
       "      <td>gi|1048081328|gb|KX616653|Influenza B virus (B...</td>\n",
       "      <td>GTTTGATTTATCATAATGGATACTTTTATTACAAGAAACTTCCAGA...</td>\n",
       "      <td>PA</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>278994</th>\n",
       "      <td>gi|525553452|gb|KF314561|Influenza A virus (A/...</td>\n",
       "      <td>AATAACGATTGGCTCTGTTTCTCTCACCATTTCCACAATATGCTTC...</td>\n",
       "      <td>NA</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>86479</th>\n",
       "      <td>gi|237874843|gb|FJ985245|Influenza A virus (A/...</td>\n",
       "      <td>ATTCAATATGGAAAGAATAAAAGAATTAAGAGATCTAATGTCACAG...</td>\n",
       "      <td>PB2</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>151351</th>\n",
       "      <td>gi|315440734|gb|CY079457|Influenza A virus (A/...</td>\n",
       "      <td>TACTGATTCAAAATGGAAGACTTTGTGCGACAATGCTTCAATCCAA...</td>\n",
       "      <td>uncl</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>412958</th>\n",
       "      <td>gi|937172275|gb|KT854258|Influenza B virus (B/...</td>\n",
       "      <td>ATGTCGCTGTTTGGAGACACAATTGCCTACCTGCTTTCATTGACAG...</td>\n",
       "      <td>M1</td>\n",
       "      <td>7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>315125</th>\n",
       "      <td>gi|594215428|gb|CY174927|Influenza B virus (B/...</td>\n",
       "      <td>ATGAATATAAATCCTTATTTTCTCTTCATAGATGTACCCATACAGG...</td>\n",
       "      <td>PB1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74069</th>\n",
       "      <td>gi|198285780|gb|CY030920|Influenza A virus (A/...</td>\n",
       "      <td>ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...</td>\n",
       "      <td>uncl</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>505337</th>\n",
       "      <td>gi|1120605462|gb|KY347748|Influenza A virus (A...</td>\n",
       "      <td>ATGGATTCCAACACGATAACCTCGTTTCAGGTAGATTGTTATCTAT...</td>\n",
       "      <td>NEP</td>\n",
       "      <td>8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>340273</th>\n",
       "      <td>gi|660873930|gb|KM064239|Influenza A virus (A/...</td>\n",
       "      <td>ATGAAGACTATCATTGCTTTGAGCTACACTCTATGTCTGGTTTTCG...</td>\n",
       "      <td>HA</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>556365</th>\n",
       "      <td>gi|1219455324|gb|CY239151|Influenza B virus (B...</td>\n",
       "      <td>GTTTGATTTGTCATAATGGATACTTTTATTACAAGAAACTTCCAGA...</td>\n",
       "      <td>PA</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>623475</th>\n",
       "      <td>gi|1390209297|gb|MH363175|Influenza B virus (B...</td>\n",
       "      <td>ATTTTCTAATATCCACAGAATGAAGGCAATAATTGTACTACTCATG...</td>\n",
       "      <td>HA</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>273912</th>\n",
       "      <td>gi|512772000|gb|KF216939|Influenza B virus (B/...</td>\n",
       "      <td>ATGCTACCTTCAACTATACAAACGTTAACCCTATTTCTCACATCAG...</td>\n",
       "      <td>NA</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>636364</th>\n",
       "      <td>gi|1421412631|gb|MH598206|Influenza A virus (A...</td>\n",
       "      <td>GGTCACAATGGAAAAATTCATCATTTTGAGCACTGTCCTGGCAGTA...</td>\n",
       "      <td>HA</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>467742</th>\n",
       "      <td>gi|1041633881|gb|LC123364|Influenza C virus (C...</td>\n",
       "      <td>ATGTCTCTTCTATTGACAATAGCAAAGGAATACAAAAGACTATGCC...</td>\n",
       "      <td>PB2</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>156141</th>\n",
       "      <td>gi|324345671|gb|CY084441|Influenza A virus (A/...</td>\n",
       "      <td>ATGAGTGACATCGAAGCCATGGCGTCTCAAGGCACCAAACGATCAT...</td>\n",
       "      <td>NP</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>165106</th>\n",
       "      <td>gi|336286441|gb|JN087059|Influenza A virus (A/...</td>\n",
       "      <td>ATATTGAAAGATGAGTCTTCTAACCGAGGTCGAAACGTACGTTCTC...</td>\n",
       "      <td>M1</td>\n",
       "      <td>7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>371494</th>\n",
       "      <td>gi|755195545|gb|KP456628|Influenza A virus (A/...</td>\n",
       "      <td>ATGAAGACTATCATTGCTTTGAGCTACATTCTATGTCTGGTTTTCG...</td>\n",
       "      <td>HA</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>424904</th>\n",
       "      <td>gi|974741015|gb|CY195307|Influenza A virus (A/...</td>\n",
       "      <td>ATGGACTCCAATACTGTGTCAAGCTTTCAGGTAGACTGTTTCCTTT...</td>\n",
       "      <td>NEP</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>111975</th>\n",
       "      <td>gi|283551383|gb|CY047996|Influenza A virus (A/...</td>\n",
       "      <td>GGAGAACGCCAGAATGCTACTGAGATCAGGGCATCTGTTGGAAGAA...</td>\n",
       "      <td>uncl</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>561122</th>\n",
       "      <td>gi|1237086641|gb|MF630455|Influenza A virus (A...</td>\n",
       "      <td>ATGAATCCAAATCAGAAGATTCTATGCACTTCAGCCACTGCTATCC...</td>\n",
       "      <td>NA</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>298789</th>\n",
       "      <td>gi|563400421|gb|CY158761|Influenza A virus (A/...</td>\n",
       "      <td>ATGAAGGCAATACTATTAGTCTTGTTATATACATTTGTACCCGCAA...</td>\n",
       "      <td>HA</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>328210</th>\n",
       "      <td>gi|607259322|gb|KJ579959|Influenza A virus (A/...</td>\n",
       "      <td>TTATTCGTCTCAGGGAGCAAAAGCAGGGTAGATAATCACTCACTGA...</td>\n",
       "      <td>NP</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>374982</th>\n",
       "      <td>gi|755202521|gb|KP460116|Influenza B virus (B/...</td>\n",
       "      <td>GATCGAATCTGCACTGGGATAACATCTTCAAACTCACCTCATGTGG...</td>\n",
       "      <td>HA</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>389631</th>\n",
       "      <td>gi|827499629|gb|KR701147|Influenza A virus (A/...</td>\n",
       "      <td>ATGAATCCAAATCAGAAGATAATAACAATTGGCTCTGTTTCTCTAT...</td>\n",
       "      <td>NA</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>529407</th>\n",
       "      <td>gi|1176711220|gb|KY855523|Influenza A virus (A...</td>\n",
       "      <td>AGCGAAAGCAGGTCAAATATATTCAATATGGAAAGAATAAAAGAAC...</td>\n",
       "      <td>PB2</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>567789</th>\n",
       "      <td>gi|1246295632|gb|MF801575|Influenza A virus (A...</td>\n",
       "      <td>ATGGATGTCAATCCGACCCTACTTTTCTTAAAGGTTCCAGCGCAAA...</td>\n",
       "      <td>PB1</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18713</th>\n",
       "      <td>gi|77747051|gb|CY003463|Influenza A virus (A/N...</td>\n",
       "      <td>AGCAAAAGCAGGTCAATTATATTCAGTATGGAAAGAATAAAAGAAC...</td>\n",
       "      <td>uncl</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>627939</th>\n",
       "      <td>gi|1408067743|gb|MH341870|Influenza A virus (A...</td>\n",
       "      <td>TCAAATCAACCAAGATGGAAGCAAAACTGTTTGTACTATTCTGTAC...</td>\n",
       "      <td>HA</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4629</th>\n",
       "      <td>gi|15555102|gb|AX225029|Sequence 39 from Paten...</td>\n",
       "      <td>AGCAAAAGCAGGGTGACAAAAACATAATGGATTCCAACACTGTGTC...</td>\n",
       "      <td>uncl</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>388575</th>\n",
       "      <td>gi|827497043|gb|KR700091|Influenza A virus (A/...</td>\n",
       "      <td>ATGAATCCAAATCAAAAGATAATAACAATTGGCTCTGTCTCCCTAA...</td>\n",
       "      <td>NA</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>400891</th>\n",
       "      <td>gi|923095095|gb|KT374233|Influenza A virus (A/...</td>\n",
       "      <td>ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...</td>\n",
       "      <td>PA</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>495092</th>\n",
       "      <td>gi|1102154012|gb|KX979673|Influenza A virus (A...</td>\n",
       "      <td>AGCAAAAGCAGGCAAACCATTTGAATGGATGTCAATCCGACTTTAC...</td>\n",
       "      <td>PB1</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>49282</th>\n",
       "      <td>gi|153972782|gb|EU050209|Influenza A virus (A/...</td>\n",
       "      <td>CACTCACTGAGTGACATTCACGCCATGGCGCTTCAAGGCACCAAAC...</td>\n",
       "      <td>NP</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>76167</th>\n",
       "      <td>gi|209866820|gb|CY035609|Influenza A virus (A/...</td>\n",
       "      <td>ATGATCGTCGAGCTTGCGGAAAAGGCAATGAAAGAATATGGGGAGG...</td>\n",
       "      <td>uncl</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>164042</th>\n",
       "      <td>gi|334189508|gb|CY091572|Influenza A virus (A/...</td>\n",
       "      <td>ATTATATTCAGTATGGAAAGAATAAAAGAACTACGGAATCTGATGT...</td>\n",
       "      <td>PB2</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>639802</th>\n",
       "      <td>gi|1423484309|gb|MH606905|Influenza B virus (B...</td>\n",
       "      <td>ATTTTCTTGTGAACTTCAAGCACCAGTGAAAGAACTGAAAATCAAA...</td>\n",
       "      <td>NP</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>432587</th>\n",
       "      <td>gi|1000373120|gb|CY197474|Influenza B virus (B...</td>\n",
       "      <td>TTTCTTAAAATGTCGCTGTTTGGAGACACAATTGCCTACCTGCTTT...</td>\n",
       "      <td>M1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>396439</th>\n",
       "      <td>gi|908275253|gb|KT283502|Influenza B virus (B/...</td>\n",
       "      <td>ATGAAGGCAATAATTGTACTACTCATGGTAGTAACATCCAATGCAG...</td>\n",
       "      <td>HA</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26965</th>\n",
       "      <td>gi|90571769|gb|CY010022|Influenza A virus (A/N...</td>\n",
       "      <td>AATCCAAATCAAAAGATAATAACGATTGGCTCTGTTTCTCTCACTA...</td>\n",
       "      <td>uncl</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>314234</th>\n",
       "      <td>gi|594211014|gb|CY174083|Influenza B virus (B/...</td>\n",
       "      <td>CTTCTCAAAACTGAGGCAAATAGGCCAAAAATGAACAATGCTACCT...</td>\n",
       "      <td>NA</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63374</th>\n",
       "      <td>gi|183222959|gb|EU555227|Influenza B virus (B/...</td>\n",
       "      <td>ATTGTACTACTCATGGTAGTAACATCCAATGCAGATCGAATCTGCA...</td>\n",
       "      <td>uncl</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>60259</th>\n",
       "      <td>gi|167735711|gb|AB378686|Influenza A virus (A/...</td>\n",
       "      <td>CAAATTTCTCAATTCTATTTGATCAACTAGAGACGTTAGTCTCTCT...</td>\n",
       "      <td>uncl</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>170937</th>\n",
       "      <td>gi|341610267|gb|JN244211|Influenza A virus (A/...</td>\n",
       "      <td>TATTTGCATTATCTGGAGTGGCAGTAGCACTTAGTGTACTGAACTT...</td>\n",
       "      <td>NA</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>122382</th>\n",
       "      <td>gi|294861120|gb|CY061602|Influenza A virus (A/...</td>\n",
       "      <td>ACAAAATGAATACTCAAATTTTAATACTCACTCTTGTGGCGGTCAT...</td>\n",
       "      <td>HA</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4695</th>\n",
       "      <td>gi|16076700|gb|AF222811|Influenza A virus (A/S...</td>\n",
       "      <td>ATGGAAGCAGAATCACTAATAACTATACTACTAGTAGTAACAGCAA...</td>\n",
       "      <td>uncl</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>263090</th>\n",
       "      <td>gi|488465200|gb|KC891811|Influenza B virus (B/...</td>\n",
       "      <td>ATGAAGGCAATAATTGTACTACTAATGGTAGTAACATCCAATGCAG...</td>\n",
       "      <td>HA</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>568420</th>\n",
       "      <td>gi|1248280161|gb|MF927815|Influenza A virus (A...</td>\n",
       "      <td>ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...</td>\n",
       "      <td>PA</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>479411</th>\n",
       "      <td>gi|1072743785|gb|KX851867|Influenza A virus (A...</td>\n",
       "      <td>ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...</td>\n",
       "      <td>PA</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>532585</th>\n",
       "      <td>gi|1181553196|gb|KY935553|Influenza A virus (A...</td>\n",
       "      <td>AAAGCAGGGTGACAAAAACATAATGGAATCCAACACCATGTCAAGC...</td>\n",
       "      <td>NEP</td>\n",
       "      <td>8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>493894</th>\n",
       "      <td>gi|1101014193|gb|KY089804|Influenza A virus (A...</td>\n",
       "      <td>ATGAATCCAAATCAAAAGATAATAACGATTGGCTCTGTTTCTCTCA...</td>\n",
       "      <td>NA</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>286402</th>\n",
       "      <td>gi|549395916|gb|CY150729|Influenza B virus (B/...</td>\n",
       "      <td>TTTCTAATATCCACAAAATGAAGGCAATAATTGTACTACTCATGGT...</td>\n",
       "      <td>HA</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>216098</th>\n",
       "      <td>gi|396578002|gb|AB735609|Influenza A virus (A/...</td>\n",
       "      <td>ATCCAAATCAAAAGATAACATGCATTTCAGCAACAGGAATGACACT...</td>\n",
       "      <td>NA</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>628780</th>\n",
       "      <td>gi|1410305195|gb|MH546245|Influenza A virus (A...</td>\n",
       "      <td>AGTGAAAATGAATCCAAATCAGAAAATAATAACAATTGGCTCTGTC...</td>\n",
       "      <td>NA</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18303</th>\n",
       "      <td>gi|76574387|gb|DQ145544|Influenza A virus (A/s...</td>\n",
       "      <td>ATGGATGTCAATCCGACCCTACTTTTCCTAAAAGTTCCAGCGCAAA...</td>\n",
       "      <td>PB1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>293124</th>\n",
       "      <td>gi|552052068|gb|KF695274|Influenza A virus (A/...</td>\n",
       "      <td>AGCAGGGTGACAAAGACATAATGGATTCCAACACGATAACCTCGTT...</td>\n",
       "      <td>NEP</td>\n",
       "      <td>8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>186510</th>\n",
       "      <td>gi|370321083|gb|EU026031|Influenza A virus (A/...</td>\n",
       "      <td>TCAAAATGAATCCAAATCAAAAGATAATAACTATTGGGTCAATCTG...</td>\n",
       "      <td>NA</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>379609</th>\n",
       "      <td>gi|766559216|gb|KP767578|Influenza A virus (A/...</td>\n",
       "      <td>ACTGTGCCTAGGACACCATGCAGTATCGAATGGGACAATAGTGAAG...</td>\n",
       "      <td>HA</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>580342</th>\n",
       "      <td>gi|1277263174|gb|CY254736|Influenza A virus (A...</td>\n",
       "      <td>ATTATATTCAGTATGGAAAGAATAAAAGAACTACGGAATCTGATGT...</td>\n",
       "      <td>PB2</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1000 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                   titles  \\\n",
       "629174  gi|1410306180|gb|MH546639|Influenza A virus (A...   \n",
       "610132  gi|1370449901|gb|MG957554|Influenza A virus (A...   \n",
       "470806  gi|1048081328|gb|KX616653|Influenza B virus (B...   \n",
       "278994  gi|525553452|gb|KF314561|Influenza A virus (A/...   \n",
       "86479   gi|237874843|gb|FJ985245|Influenza A virus (A/...   \n",
       "151351  gi|315440734|gb|CY079457|Influenza A virus (A/...   \n",
       "412958  gi|937172275|gb|KT854258|Influenza B virus (B/...   \n",
       "315125  gi|594215428|gb|CY174927|Influenza B virus (B/...   \n",
       "74069   gi|198285780|gb|CY030920|Influenza A virus (A/...   \n",
       "505337  gi|1120605462|gb|KY347748|Influenza A virus (A...   \n",
       "340273  gi|660873930|gb|KM064239|Influenza A virus (A/...   \n",
       "556365  gi|1219455324|gb|CY239151|Influenza B virus (B...   \n",
       "623475  gi|1390209297|gb|MH363175|Influenza B virus (B...   \n",
       "273912  gi|512772000|gb|KF216939|Influenza B virus (B/...   \n",
       "636364  gi|1421412631|gb|MH598206|Influenza A virus (A...   \n",
       "467742  gi|1041633881|gb|LC123364|Influenza C virus (C...   \n",
       "156141  gi|324345671|gb|CY084441|Influenza A virus (A/...   \n",
       "165106  gi|336286441|gb|JN087059|Influenza A virus (A/...   \n",
       "371494  gi|755195545|gb|KP456628|Influenza A virus (A/...   \n",
       "424904  gi|974741015|gb|CY195307|Influenza A virus (A/...   \n",
       "111975  gi|283551383|gb|CY047996|Influenza A virus (A/...   \n",
       "561122  gi|1237086641|gb|MF630455|Influenza A virus (A...   \n",
       "298789  gi|563400421|gb|CY158761|Influenza A virus (A/...   \n",
       "328210  gi|607259322|gb|KJ579959|Influenza A virus (A/...   \n",
       "374982  gi|755202521|gb|KP460116|Influenza B virus (B/...   \n",
       "389631  gi|827499629|gb|KR701147|Influenza A virus (A/...   \n",
       "529407  gi|1176711220|gb|KY855523|Influenza A virus (A...   \n",
       "567789  gi|1246295632|gb|MF801575|Influenza A virus (A...   \n",
       "18713   gi|77747051|gb|CY003463|Influenza A virus (A/N...   \n",
       "627939  gi|1408067743|gb|MH341870|Influenza A virus (A...   \n",
       "...                                                   ...   \n",
       "4629    gi|15555102|gb|AX225029|Sequence 39 from Paten...   \n",
       "388575  gi|827497043|gb|KR700091|Influenza A virus (A/...   \n",
       "400891  gi|923095095|gb|KT374233|Influenza A virus (A/...   \n",
       "495092  gi|1102154012|gb|KX979673|Influenza A virus (A...   \n",
       "49282   gi|153972782|gb|EU050209|Influenza A virus (A/...   \n",
       "76167   gi|209866820|gb|CY035609|Influenza A virus (A/...   \n",
       "164042  gi|334189508|gb|CY091572|Influenza A virus (A/...   \n",
       "639802  gi|1423484309|gb|MH606905|Influenza B virus (B...   \n",
       "432587  gi|1000373120|gb|CY197474|Influenza B virus (B...   \n",
       "396439  gi|908275253|gb|KT283502|Influenza B virus (B/...   \n",
       "26965   gi|90571769|gb|CY010022|Influenza A virus (A/N...   \n",
       "314234  gi|594211014|gb|CY174083|Influenza B virus (B/...   \n",
       "63374   gi|183222959|gb|EU555227|Influenza B virus (B/...   \n",
       "60259   gi|167735711|gb|AB378686|Influenza A virus (A/...   \n",
       "170937  gi|341610267|gb|JN244211|Influenza A virus (A/...   \n",
       "122382  gi|294861120|gb|CY061602|Influenza A virus (A/...   \n",
       "4695    gi|16076700|gb|AF222811|Influenza A virus (A/S...   \n",
       "263090  gi|488465200|gb|KC891811|Influenza B virus (B/...   \n",
       "568420  gi|1248280161|gb|MF927815|Influenza A virus (A...   \n",
       "479411  gi|1072743785|gb|KX851867|Influenza A virus (A...   \n",
       "532585  gi|1181553196|gb|KY935553|Influenza A virus (A...   \n",
       "493894  gi|1101014193|gb|KY089804|Influenza A virus (A...   \n",
       "286402  gi|549395916|gb|CY150729|Influenza B virus (B/...   \n",
       "216098  gi|396578002|gb|AB735609|Influenza A virus (A/...   \n",
       "628780  gi|1410305195|gb|MH546245|Influenza A virus (A...   \n",
       "18303   gi|76574387|gb|DQ145544|Influenza A virus (A/s...   \n",
       "293124  gi|552052068|gb|KF695274|Influenza A virus (A/...   \n",
       "186510  gi|370321083|gb|EU026031|Influenza A virus (A/...   \n",
       "379609  gi|766559216|gb|KP767578|Influenza A virus (A/...   \n",
       "580342  gi|1277263174|gb|CY254736|Influenza A virus (A...   \n",
       "\n",
       "                                                sequences  gene  segment  \n",
       "629174  GTGACAAAAACATAATGGATTCCAACACGATAACCTCGTTTCAGGT...   NEP        8  \n",
       "610132  AATATGGAGAGAATAAAAGAACTAAGAGATTTGATGTCGCAGTCTC...   PB2        1  \n",
       "470806  GTTTGATTTATCATAATGGATACTTTTATTACAAGAAACTTCCAGA...    PA        3  \n",
       "278994  AATAACGATTGGCTCTGTTTCTCTCACCATTTCCACAATATGCTTC...    NA        6  \n",
       "86479   ATTCAATATGGAAAGAATAAAAGAATTAAGAGATCTAATGTCACAG...   PB2        1  \n",
       "151351  TACTGATTCAAAATGGAAGACTTTGTGCGACAATGCTTCAATCCAA...  uncl        3  \n",
       "412958  ATGTCGCTGTTTGGAGACACAATTGCCTACCTGCTTTCATTGACAG...    M1        7  \n",
       "315125  ATGAATATAAATCCTTATTTTCTCTTCATAGATGTACCCATACAGG...   PB1        0  \n",
       "74069   ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...  uncl        3  \n",
       "505337  ATGGATTCCAACACGATAACCTCGTTTCAGGTAGATTGTTATCTAT...   NEP        8  \n",
       "340273  ATGAAGACTATCATTGCTTTGAGCTACACTCTATGTCTGGTTTTCG...    HA        4  \n",
       "556365  GTTTGATTTGTCATAATGGATACTTTTATTACAAGAAACTTCCAGA...    PA        0  \n",
       "623475  ATTTTCTAATATCCACAGAATGAAGGCAATAATTGTACTACTCATG...    HA        4  \n",
       "273912  ATGCTACCTTCAACTATACAAACGTTAACCCTATTTCTCACATCAG...    NA        6  \n",
       "636364  GGTCACAATGGAAAAATTCATCATTTTGAGCACTGTCCTGGCAGTA...    HA        4  \n",
       "467742  ATGTCTCTTCTATTGACAATAGCAAAGGAATACAAAAGACTATGCC...   PB2        0  \n",
       "156141  ATGAGTGACATCGAAGCCATGGCGTCTCAAGGCACCAAACGATCAT...    NP        0  \n",
       "165106  ATATTGAAAGATGAGTCTTCTAACCGAGGTCGAAACGTACGTTCTC...    M1        7  \n",
       "371494  ATGAAGACTATCATTGCTTTGAGCTACATTCTATGTCTGGTTTTCG...    HA        4  \n",
       "424904  ATGGACTCCAATACTGTGTCAAGCTTTCAGGTAGACTGTTTCCTTT...   NEP        0  \n",
       "111975  GGAGAACGCCAGAATGCTACTGAGATCAGGGCATCTGTTGGAAGAA...  uncl        5  \n",
       "561122  ATGAATCCAAATCAGAAGATTCTATGCACTTCAGCCACTGCTATCC...    NA        6  \n",
       "298789  ATGAAGGCAATACTATTAGTCTTGTTATATACATTTGTACCCGCAA...    HA        0  \n",
       "328210  TTATTCGTCTCAGGGAGCAAAAGCAGGGTAGATAATCACTCACTGA...    NP        5  \n",
       "374982  GATCGAATCTGCACTGGGATAACATCTTCAAACTCACCTCATGTGG...    HA        4  \n",
       "389631  ATGAATCCAAATCAGAAGATAATAACAATTGGCTCTGTTTCTCTAT...    NA        6  \n",
       "529407  AGCGAAAGCAGGTCAAATATATTCAATATGGAAAGAATAAAAGAAC...   PB2        1  \n",
       "567789  ATGGATGTCAATCCGACCCTACTTTTCTTAAAGGTTCCAGCGCAAA...   PB1        2  \n",
       "18713   AGCAAAAGCAGGTCAATTATATTCAGTATGGAAAGAATAAAAGAAC...  uncl        1  \n",
       "627939  TCAAATCAACCAAGATGGAAGCAAAACTGTTTGTACTATTCTGTAC...    HA        4  \n",
       "...                                                   ...   ...      ...  \n",
       "4629    AGCAAAAGCAGGGTGACAAAAACATAATGGATTCCAACACTGTGTC...  uncl        0  \n",
       "388575  ATGAATCCAAATCAAAAGATAATAACAATTGGCTCTGTCTCCCTAA...    NA        6  \n",
       "400891  ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...    PA        3  \n",
       "495092  AGCAAAAGCAGGCAAACCATTTGAATGGATGTCAATCCGACTTTAC...   PB1        2  \n",
       "49282   CACTCACTGAGTGACATTCACGCCATGGCGCTTCAAGGCACCAAAC...    NP        0  \n",
       "76167   ATGATCGTCGAGCTTGCGGAAAAGGCAATGAAAGAATATGGGGAGG...  uncl        3  \n",
       "164042  ATTATATTCAGTATGGAAAGAATAAAAGAACTACGGAATCTGATGT...   PB2        0  \n",
       "639802  ATTTTCTTGTGAACTTCAAGCACCAGTGAAAGAACTGAAAATCAAA...    NP        5  \n",
       "432587  TTTCTTAAAATGTCGCTGTTTGGAGACACAATTGCCTACCTGCTTT...    M1        0  \n",
       "396439  ATGAAGGCAATAATTGTACTACTCATGGTAGTAACATCCAATGCAG...    HA        4  \n",
       "26965   AATCCAAATCAAAAGATAATAACGATTGGCTCTGTTTCTCTCACTA...  uncl        6  \n",
       "314234  CTTCTCAAAACTGAGGCAAATAGGCCAAAAATGAACAATGCTACCT...    NA        0  \n",
       "63374   ATTGTACTACTCATGGTAGTAACATCCAATGCAGATCGAATCTGCA...  uncl        0  \n",
       "60259   CAAATTTCTCAATTCTATTTGATCAACTAGAGACGTTAGTCTCTCT...  uncl        0  \n",
       "170937  TATTTGCATTATCTGGAGTGGCAGTAGCACTTAGTGTACTGAACTT...    NA        6  \n",
       "122382  ACAAAATGAATACTCAAATTTTAATACTCACTCTTGTGGCGGTCAT...    HA        0  \n",
       "4695    ATGGAAGCAGAATCACTAATAACTATACTACTAGTAGTAACAGCAA...  uncl        0  \n",
       "263090  ATGAAGGCAATAATTGTACTACTAATGGTAGTAACATCCAATGCAG...    HA        4  \n",
       "568420  ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...    PA        3  \n",
       "479411  ATGGAAGACTTTGTGCGACAATGCTTCAATCCAATGATCGTCGAGC...    PA        3  \n",
       "532585  AAAGCAGGGTGACAAAAACATAATGGAATCCAACACCATGTCAAGC...   NEP        8  \n",
       "493894  ATGAATCCAAATCAAAAGATAATAACGATTGGCTCTGTTTCTCTCA...    NA        6  \n",
       "286402  TTTCTAATATCCACAAAATGAAGGCAATAATTGTACTACTCATGGT...    HA        0  \n",
       "216098  ATCCAAATCAAAAGATAACATGCATTTCAGCAACAGGAATGACACT...    NA        0  \n",
       "628780  AGTGAAAATGAATCCAAATCAGAAAATAATAACAATTGGCTCTGTC...    NA        6  \n",
       "18303   ATGGATGTCAATCCGACCCTACTTTTCCTAAAAGTTCCAGCGCAAA...   PB1        0  \n",
       "293124  AGCAGGGTGACAAAGACATAATGGATTCCAACACGATAACCTCGTT...   NEP        8  \n",
       "186510  TCAAAATGAATCCAAATCAAAAGATAATAACTATTGGGTCAATCTG...    NA        6  \n",
       "379609  ACTGTGCCTAGGACACCATGCAGTATCGAATGGGACAATAGTGAAG...    HA        4  \n",
       "580342  ATTATATTCAGTATGGAAAGAATAAAAGAACTACGGAATCTGATGT...   PB2        0  \n",
       "\n",
       "[1000 rows x 4 columns]"
      ]
     },
     "execution_count": 30,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "uniflna_test1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "344"
      ]
     },
     "execution_count": 31,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sum(uniflna_test1['segment']==0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "153"
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sum(uniflna_test1['gene']=='uncl')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "22"
      ]
     },
     "execution_count": 33,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sum((uniflna_test1['segment']==0) & (uniflna_test1['gene']=='uncl'))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.022"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "22/1000"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### try sourmash, then export, then paint clusters"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 85,
   "metadata": {},
   "outputs": [],
   "source": [
    "# make it back into text FASTA file\n",
    "with open (\"uniflna_test1.fasta\", 'w') as output:\n",
    "    for i in range(len(uniflna_test1)):\n",
    "        output.write(\">\" + list(uniflna_test1.titles.values)[i] + \"\\n\" + list(uniflna_test1.sequences.values)[i] + \"\\n\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 86,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "\u001b[Kcomputing signatures for files: uniflna_test1.fasta\n",
      "\u001b[KComputing signature for ksizes: [7]\n",
      "\u001b[KComputing only nucleotide (and not protein) signatures.\n",
      "\u001b[KComputing a total of 1 signature(s).\n",
      "\u001b[Kcalculated 1000 signatures for 1000 sequences in uniflna_test1.fasta\n",
      "\u001b[Ksaved 1000 signature(s). Note: signature license is CC0.\n",
      "finished in 1.0520622730255127seconds\n"
     ]
    }
   ],
   "source": [
    "start = time.time()\n",
    "!sourmash compute uniflna_test1.fasta  -k 7 --singleton\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 87,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "...sig loading 999_test1.fasta.sig\n",
      "\u001b[Kloaded 1000 signatures total.                                                  \n",
      "\u001b[K\n",
      "min similarity in matrix: 0.007\n",
      "\u001b[Ksaving labels to: uniflna_test1_k7_cmp.labels.txt\n",
      "\u001b[Ksaving distance matrix to: uniflna_test1_k7_cmp\n",
      "finished in 56.853424310684204seconds\n"
     ]
    }
   ],
   "source": [
    "# comparison to matrx\n",
    "start = time.time()\n",
    "!sourmash compare uniflna_test1.fasta.sig -o uniflna_test1_k7_cmp\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "...sig loading 999_test1.fasta.sig\n",
      "\u001b[Kloaded 1000 signatures total.                                                  \n",
      "\u001b[K\n",
      "min similarity in matrix: 0.007\n",
      "finished in 61.660585165023804seconds\n"
     ]
    }
   ],
   "source": [
    "# comparison to csv for plotting in R\n",
    "start = time.time()\n",
    "!sourmash compare uniflna_test1.fasta.sig --csv uniflna_test1_cmp.csv\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [],
   "source": [
    "# export csv of dataframe\n",
    "uniflna_test1.to_csv(\"uniflna_test1.csv\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "larger 2nd sample"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "np.random.seed(23874)\n",
    "uniflna_test2 = uniflna.sample(n=2500)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "# export csv of dataframe\n",
    "uniflna_test2.to_csv(\"uniflna_test2.csv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "uniflna_test2['gene'] = uniflna_test2.apply(row_test, axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "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>titles</th>\n",
       "      <th>sequences</th>\n",
       "      <th>gene</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>136410</th>\n",
       "      <td>gi|302375385|gb|CY068989|Influenza A virus (A/...</td>\n",
       "      <td>TCACTCAATGAGTGACATCGAAGCCATGGCGTCTCAAGGCACCAAA...</td>\n",
       "      <td>uncl</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>219699</th>\n",
       "      <td>gi|402485249|gb|JX454740|Influenza A virus (A/...</td>\n",
       "      <td>AGCAAAAGCAGGGTAGATAATCACTCACTGAGTGACATCAACATCA...</td>\n",
       "      <td>NP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>76150</th>\n",
       "      <td>gi|209866781|gb|CY035592|Influenza A virus (A/...</td>\n",
       "      <td>TAAGAGATCTAATGTCACAGTCCCGCACTCGCGAGATACTCACTAA...</td>\n",
       "      <td>uncl</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>96229</th>\n",
       "      <td>gi|256595781|gb|CY045071|Influenza A virus (A/...</td>\n",
       "      <td>AAATGAATCCAAACCAAAAGATAATAACCATTGGTTCGGTCTGTAT...</td>\n",
       "      <td>uncl</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>611239</th>\n",
       "      <td>gi|1371434282|gb|MH126459|Influenza A virus (A...</td>\n",
       "      <td>CAAACCATTTGAATGGATGTCAATCCGACTCTACTGTTCTTAAAAG...</td>\n",
       "      <td>PB1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                   titles  \\\n",
       "136410  gi|302375385|gb|CY068989|Influenza A virus (A/...   \n",
       "219699  gi|402485249|gb|JX454740|Influenza A virus (A/...   \n",
       "76150   gi|209866781|gb|CY035592|Influenza A virus (A/...   \n",
       "96229   gi|256595781|gb|CY045071|Influenza A virus (A/...   \n",
       "611239  gi|1371434282|gb|MH126459|Influenza A virus (A...   \n",
       "\n",
       "                                                sequences  gene  \n",
       "136410  TCACTCAATGAGTGACATCGAAGCCATGGCGTCTCAAGGCACCAAA...  uncl  \n",
       "219699  AGCAAAAGCAGGGTAGATAATCACTCACTGAGTGACATCAACATCA...    NP  \n",
       "76150   TAAGAGATCTAATGTCACAGTCCCGCACTCGCGAGATACTCACTAA...  uncl  \n",
       "96229   AAATGAATCCAAACCAAAAGATAATAACCATTGGTTCGGTCTGTAT...  uncl  \n",
       "611239  CAAACCATTTGAATGGATGTCAATCCGACTCTACTGTTCTTAAAAG...   PB1  "
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "uniflna_test2.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "387"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sum(uniflna_test2['gene'] == 'uncl')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "# make it back into text FASTA file\n",
    "with open (\"uniflna_test2.fasta\", 'w') as output:\n",
    "    for i in range(len(uniflna_test2)):\n",
    "        output.write(\">\" + list(uniflna_test2.titles.values)[i] + \"\\n\" + list(uniflna_test2.sequences.values)[i] + \"\\n\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "\u001b[Kcomputing signatures for files: uniflna_test2.fasta\n",
      "\u001b[KComputing signature for ksizes: [7]\n",
      "\u001b[KComputing only nucleotide (and not protein) signatures.\n",
      "\u001b[KComputing a total of 1 signature(s).\n",
      "\u001b[Kcalculated 2500 signatures for 2500 sequences in uniflna_test2.fasta\n",
      "\u001b[Ksaved 2500 signature(s). Note: signature license is CC0.\n",
      "finished in 2.0141425132751465seconds\n"
     ]
    }
   ],
   "source": [
    "start = time.time()\n",
    "!sourmash compute uniflna_test2.fasta  -k 7 --singleton\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[K== This is sourmash version 2.0.0. ==\n",
      "\u001b[K== Please cite Brown and Irber (2016), doi:10.21105/joss.00027. ==\n",
      "\n",
      "...sig loading 2,499est2.fasta.sig\n",
      "\u001b[Kloaded 2500 signatures total.                                                  \n",
      "\u001b[K\n",
      "min similarity in matrix: 0.004\n",
      "finished in 352.2169532775879seconds\n"
     ]
    }
   ],
   "source": [
    "# comparison to csv for plotting in R\n",
    "start = time.time()\n",
    "!sourmash compare uniflna_test2.fasta.sig --csv uniflna_test2_cmp.csv\n",
    "print(\"finished in \" + str(time.time()-start) + \"seconds\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "more labels"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [],
   "source": [
    "# small test\n",
    "np.random.seed(4351)\n",
    "uniflna_smalltest = uniflna.sample(n=25)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {
    "collapsed": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['gi|301638037|gb|CY050131|Influenza A virus (A/Qingdao/1364/2009(H1N1)) segment 6 sequence',\n",
       "       'gi|496493392|gb|KF021599|Influenza A virus (A/Shanghai/02/2013(H7N9)) segment 6 neuraminidase (NA) gene, complete cds',\n",
       "       'gi|1039130038|gb|KX414423|Influenza A virus (A/Florida/08/2016(H3N2)) segment 2 polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|134048028|gb|CY021252|Influenza A virus (A/american black duck/Ohio/1823/2005(H11N9)) segment 1, complete sequence',\n",
       "       'gi|1162228474|gb|KY678645|Influenza A virus (A/green-winged teal/Ohio/15OS5497/2015(H4N6)) segment 6 neuraminidase (NA) gene, complete cds',\n",
       "       'gi|353251209|gb|JN804006|Influenza A virus (A/duck/Anhui/B27/2010(H9N2)) segment 4 hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|1148882880|gb|KY550827|Influenza A virus (A/American green-winged teal/California/AH0079678/2016(H7N3)) segment 6 neuraminidase (NA) gene, complete cds',\n",
       "       'gi|1017035011|gb|KU976528|Influenza A virus (A/swine/Mexico/AVX16/2012(H1N2)) segment 2 polymerase PB1 (PB1) and PB1-F2 protein (PB1-F2) genes, complete cds',\n",
       "       'gi|954172452|gb|KU160851|Influenza A virus (A/duck/Chongqing/S2086/2012(H4N8)) segment 4 hemagglutinin (HA) gene, complete cds',\n",
       "       'gi|1479548419|gb|MH932122|Influenza A virus (A/mallard/Alberta/201/2017(H4N6)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|158188049|gb|EU199324|Influenza A virus (A/New Jersey/04/2006(H1N1)) segment 4 hemagglutinin (HA) gene, partial cds',\n",
       "       'gi|1300791881|gb|CY259192|Influenza A virus (A/New York/A-WC-LVD-17-027/2017(H3N2)) matrix protein 2 (M2) and matrix protein 1 (M1) genes, complete cds',\n",
       "       'gi|1304402427|gb|CY260205|Influenza B virus (B/Iowa/30/2017) matrix protein 1 (M1) and BM2 protein (BM2) genes, complete cds',\n",
       "       'gi|1387697499|gb|MH306713|Influenza A virus (A/South Dakota/20/2018(H3N2)) segment 1 polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|973426422|gb|KU299845|Influenza A virus (A/Cambodia/W1023353/2012(H3N2)) segment 7 matrix protein 2 (M2) and matrix protein 1 (M1) genes, complete cds',\n",
       "       'gi|49356905|gb|AY633163|Influenza A virus (A/mallard/Alberta/136/00(H4N6)) polymerase basic protein 2 (PB2) gene, partial cds',\n",
       "       'gi|469937682|gb|KC739241|Influenza A virus (A/Taiwan/3/2007(H3N2)) segment 6 neuraminidase (NA) gene, partial cds',\n",
       "       'gi|1000381783|gb|CY197920|Influenza B virus (B/South Australia/12/2013) polymerase PB2 (PB2) gene, complete cds',\n",
       "       'gi|1017034988|gb|KU976946|Influenza A virus (A/swine/Mexico/AVX28/2012(H3N2)) segment 5 nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|115608201|gb|CY016947|Influenza A virus (A/chicken/Nigeria/1047-34/2006(H5N1)) segment 4, complete sequence',\n",
       "       'gi|662494622|gb|KM029544|Influenza A virus (A/swine/Hong Kong/2569/2010(H1N2)) segment 5 nucleocapsid protein (NP) gene, complete cds',\n",
       "       'gi|304441125|gb|HQ156772|Influenza A virus (A/chicken/Bangladesh/1151-10/2010(H5N1)) segment 1 polymerase PB2 (PB2) gene, partial cds',\n",
       "       'gi|61970814|gb|CY000699|Influenza A virus (A/New York/174/2000(H3N2)) segment 6, complete sequence',\n",
       "       'gi|82622954|gb|DQ280242|Influenza A virus (A/swine/Ontario/23866/04(H1N1)) neuraminidase (NA) gene, complete cds',\n",
       "       'gi|393395100|gb|CY114995|Influenza A virus (A/swine/Illinois/A01134451/2011(H1N1)) matrix protein 2 (M2) gene, partial cds; and matrix protein 1 (M1) gene, complete cds'],\n",
       "      dtype=object)"
      ]
     },
     "execution_count": 38,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "uniflna_smalltest.titles.values"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "' A virus (A/Qingdao/1364/2009(H1N1)) segment 6 sequence'"
      ]
     },
     "execution_count": 44,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "uniflna_smalltest.titles.values[0].split('Influenza',1)[1].split('')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [],
   "source": [
    "def virus_type(row):\n",
    "    if 'Influenza A virus' in row['titles']:\n",
    "        return('A')\n",
    "    if 'Influenza B virus' in row['titles']:\n",
    "        return('B')\n",
    "    if 'Influenza C virus' in row['titles']:\n",
    "        return('C')\n",
    "    else:\n",
    "        return('unk')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {},
   "outputs": [],
   "source": [
    "uniflna_smalltest['type'] = uniflna_smalltest.apply(virus_type, axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {},
   "outputs": [],
   "source": [
    "import re\n",
    "def HN_type(row):\n",
    "    #print(row['titles'])\n",
    "    gtm = re.search(r'(?:\\(H)(.*)(\\){2})', row['titles'])\n",
    "    if gtm==None:\n",
    "        return('none')\n",
    "    else:\n",
    "        return(re.sub(r'[\\(\\)]', '', gtm.group(0)))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "metadata": {},
   "outputs": [],
   "source": [
    "uniflna_smalltest['HNtype'] = uniflna_smalltest.apply(HN_type, axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "metadata": {
    "scrolled": true
   },
   "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>titles</th>\n",
       "      <th>sequences</th>\n",
       "      <th>type</th>\n",
       "      <th>HNtype</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>135018</th>\n",
       "      <td>gi|301638037|gb|CY050131|Influenza A virus (A/...</td>\n",
       "      <td>ATGAATCCAAACCAAAAGATAATAACCATTGGTTCGGTCTGTATGA...</td>\n",
       "      <td>A</td>\n",
       "      <td>H1N1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>264854</th>\n",
       "      <td>gi|496493392|gb|KF021599|Influenza A virus (A/...</td>\n",
       "      <td>ATGAATCCAAATCAGAAGATTCTATGCACTTCAGCCACTGCTATCA...</td>\n",
       "      <td>A</td>\n",
       "      <td>H7N9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>458734</th>\n",
       "      <td>gi|1039130038|gb|KX414423|Influenza A virus (A...</td>\n",
       "      <td>CAAACCATTTGAATGGATGTCAATCCGACTCTACTGTTCTTAAAAG...</td>\n",
       "      <td>A</td>\n",
       "      <td>H3N2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>44580</th>\n",
       "      <td>gi|134048028|gb|CY021252|Influenza A virus (A/...</td>\n",
       "      <td>AATATATTCAATATGGAGAGAATAAAAGAACTAAGAGATCTAATGT...</td>\n",
       "      <td>A</td>\n",
       "      <td>H11N9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>526161</th>\n",
       "      <td>gi|1162228474|gb|KY678645|Influenza A virus (A...</td>\n",
       "      <td>ATGAATCCAAATCAGAAGATAATATGCATCTCAGCAACAGGAATGA...</td>\n",
       "      <td>A</td>\n",
       "      <td>H4N6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>179856</th>\n",
       "      <td>gi|353251209|gb|JN804006|Influenza A virus (A/...</td>\n",
       "      <td>ATGGAGACAGTATCACTAATAACTATACTACTAGTAGTAACAGTAA...</td>\n",
       "      <td>A</td>\n",
       "      <td>H9N2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>517819</th>\n",
       "      <td>gi|1148882880|gb|KY550827|Influenza A virus (A...</td>\n",
       "      <td>TGCGAGATGAATCCAAATCAGAAGATAATAACAATCGGGGTAGTGA...</td>\n",
       "      <td>A</td>\n",
       "      <td>H7N3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>444053</th>\n",
       "      <td>gi|1017035011|gb|KU976528|Influenza A virus (A...</td>\n",
       "      <td>AGCGAAAGCAGGCAAACCATTTGAATGGATGTCAATCCGACTCTAC...</td>\n",
       "      <td>A</td>\n",
       "      <td>H1N2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>419418</th>\n",
       "      <td>gi|954172452|gb|KU160851|Influenza A virus (A/...</td>\n",
       "      <td>ATGCTATCAATTGTGATTATGTTTCTGCTTGTTGCAGAGAGCTCTT...</td>\n",
       "      <td>A</td>\n",
       "      <td>H4N8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>651420</th>\n",
       "      <td>gi|1479548419|gb|MH932122|Influenza A virus (A...</td>\n",
       "      <td>AAATATATTCAATATGGAGAGAATAAAAGAACTAAGAGATTTAATG...</td>\n",
       "      <td>A</td>\n",
       "      <td>H4N6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55128</th>\n",
       "      <td>gi|158188049|gb|EU199324|Influenza A virus (A/...</td>\n",
       "      <td>ATGAAAGTAAAACTACTGGTTCTGTTATGTACATTTACAGCTACAT...</td>\n",
       "      <td>A</td>\n",
       "      <td>H1N1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>587361</th>\n",
       "      <td>gi|1300791881|gb|CY259192|Influenza A virus (A...</td>\n",
       "      <td>ATATTGAAAGATGAGCCTTCTAACCGAGGTCGAAACGTATGTTCTC...</td>\n",
       "      <td>A</td>\n",
       "      <td>H3N2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>588190</th>\n",
       "      <td>gi|1304402427|gb|CY260205|Influenza B virus (B...</td>\n",
       "      <td>TCTTAAAATGTCGCTGTTTGGAGACACAATTGCCTACCTGCTTTCA...</td>\n",
       "      <td>B</td>\n",
       "      <td>none</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>621327</th>\n",
       "      <td>gi|1387697499|gb|MH306713|Influenza A virus (A...</td>\n",
       "      <td>TCAATTATATTCAGCATGGAAAGAATAAAAGAACTACGGAATCTAA...</td>\n",
       "      <td>A</td>\n",
       "      <td>H3N2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>421901</th>\n",
       "      <td>gi|973426422|gb|KU299845|Influenza A virus (A/...</td>\n",
       "      <td>ATGAGCCTTCTAACCGAGGTCGAAACGTATGTTCTCTCTATCGTTC...</td>\n",
       "      <td>A</td>\n",
       "      <td>H3N2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10392</th>\n",
       "      <td>gi|49356905|gb|AY633163|Influenza A virus (A/m...</td>\n",
       "      <td>ATGTCACAGTCTCGCACCCGCGAGATACTCACCAAAACCACCGTGG...</td>\n",
       "      <td>A</td>\n",
       "      <td>H4N6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>252467</th>\n",
       "      <td>gi|469937682|gb|KC739241|Influenza A virus (A/...</td>\n",
       "      <td>GGCTCTGTTTCTCTCACCATTTCCACAATATGCTTCTTCATGCAAA...</td>\n",
       "      <td>A</td>\n",
       "      <td>H3N2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>433033</th>\n",
       "      <td>gi|1000381783|gb|CY197920|Influenza B virus (B...</td>\n",
       "      <td>ATGACATTGGCTAAAATTGAATTGTTAAAACAACTGTTAAGGGACA...</td>\n",
       "      <td>B</td>\n",
       "      <td>none</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>444043</th>\n",
       "      <td>gi|1017034988|gb|KU976946|Influenza A virus (A...</td>\n",
       "      <td>AGCGAAAGCAGGGTAGATAATCACTCAATGAGTGACATCGAAGCCA...</td>\n",
       "      <td>A</td>\n",
       "      <td>H3N2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>36005</th>\n",
       "      <td>gi|115608201|gb|CY016947|Influenza A virus (A/...</td>\n",
       "      <td>AAATGGAAAAAATAGTGCTTCTTTTTGCAATAGTCAGTCTTGTTAA...</td>\n",
       "      <td>A</td>\n",
       "      <td>H5N1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>342677</th>\n",
       "      <td>gi|662494622|gb|KM029544|Influenza A virus (A/...</td>\n",
       "      <td>ATGGCATCTCAAGGTACCAAACGATCATATGAACAAATGGAGACTG...</td>\n",
       "      <td>A</td>\n",
       "      <td>H1N2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>140445</th>\n",
       "      <td>gi|304441125|gb|HQ156772|Influenza A virus (A/...</td>\n",
       "      <td>TCACAGTCTCGCACTCGCGAGATACTAACAAAGACCACCGTGGACC...</td>\n",
       "      <td>A</td>\n",
       "      <td>H5N1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13605</th>\n",
       "      <td>gi|61970814|gb|CY000699|Influenza A virus (A/N...</td>\n",
       "      <td>CAAAAGCAGGAGTAAAGATGAATCCAAATCAAAAGATAATAACGAT...</td>\n",
       "      <td>A</td>\n",
       "      <td>H3N2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21340</th>\n",
       "      <td>gi|82622954|gb|DQ280242|Influenza A virus (A/s...</td>\n",
       "      <td>ATGAATACAAATCAAAGAATAATAACCATTGGGACGGTTTGCCTGA...</td>\n",
       "      <td>A</td>\n",
       "      <td>H1N1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>213246</th>\n",
       "      <td>gi|393395100|gb|CY114995|Influenza A virus (A/...</td>\n",
       "      <td>ATGAGTCTTCTAACCGAGGTCGAAACGTACGTTCTTTCTATCATCC...</td>\n",
       "      <td>A</td>\n",
       "      <td>H1N1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                   titles  \\\n",
       "135018  gi|301638037|gb|CY050131|Influenza A virus (A/...   \n",
       "264854  gi|496493392|gb|KF021599|Influenza A virus (A/...   \n",
       "458734  gi|1039130038|gb|KX414423|Influenza A virus (A...   \n",
       "44580   gi|134048028|gb|CY021252|Influenza A virus (A/...   \n",
       "526161  gi|1162228474|gb|KY678645|Influenza A virus (A...   \n",
       "179856  gi|353251209|gb|JN804006|Influenza A virus (A/...   \n",
       "517819  gi|1148882880|gb|KY550827|Influenza A virus (A...   \n",
       "444053  gi|1017035011|gb|KU976528|Influenza A virus (A...   \n",
       "419418  gi|954172452|gb|KU160851|Influenza A virus (A/...   \n",
       "651420  gi|1479548419|gb|MH932122|Influenza A virus (A...   \n",
       "55128   gi|158188049|gb|EU199324|Influenza A virus (A/...   \n",
       "587361  gi|1300791881|gb|CY259192|Influenza A virus (A...   \n",
       "588190  gi|1304402427|gb|CY260205|Influenza B virus (B...   \n",
       "621327  gi|1387697499|gb|MH306713|Influenza A virus (A...   \n",
       "421901  gi|973426422|gb|KU299845|Influenza A virus (A/...   \n",
       "10392   gi|49356905|gb|AY633163|Influenza A virus (A/m...   \n",
       "252467  gi|469937682|gb|KC739241|Influenza A virus (A/...   \n",
       "433033  gi|1000381783|gb|CY197920|Influenza B virus (B...   \n",
       "444043  gi|1017034988|gb|KU976946|Influenza A virus (A...   \n",
       "36005   gi|115608201|gb|CY016947|Influenza A virus (A/...   \n",
       "342677  gi|662494622|gb|KM029544|Influenza A virus (A/...   \n",
       "140445  gi|304441125|gb|HQ156772|Influenza A virus (A/...   \n",
       "13605   gi|61970814|gb|CY000699|Influenza A virus (A/N...   \n",
       "21340   gi|82622954|gb|DQ280242|Influenza A virus (A/s...   \n",
       "213246  gi|393395100|gb|CY114995|Influenza A virus (A/...   \n",
       "\n",
       "                                                sequences type HNtype  \n",
       "135018  ATGAATCCAAACCAAAAGATAATAACCATTGGTTCGGTCTGTATGA...    A   H1N1  \n",
       "264854  ATGAATCCAAATCAGAAGATTCTATGCACTTCAGCCACTGCTATCA...    A   H7N9  \n",
       "458734  CAAACCATTTGAATGGATGTCAATCCGACTCTACTGTTCTTAAAAG...    A   H3N2  \n",
       "44580   AATATATTCAATATGGAGAGAATAAAAGAACTAAGAGATCTAATGT...    A  H11N9  \n",
       "526161  ATGAATCCAAATCAGAAGATAATATGCATCTCAGCAACAGGAATGA...    A   H4N6  \n",
       "179856  ATGGAGACAGTATCACTAATAACTATACTACTAGTAGTAACAGTAA...    A   H9N2  \n",
       "517819  TGCGAGATGAATCCAAATCAGAAGATAATAACAATCGGGGTAGTGA...    A   H7N3  \n",
       "444053  AGCGAAAGCAGGCAAACCATTTGAATGGATGTCAATCCGACTCTAC...    A   H1N2  \n",
       "419418  ATGCTATCAATTGTGATTATGTTTCTGCTTGTTGCAGAGAGCTCTT...    A   H4N8  \n",
       "651420  AAATATATTCAATATGGAGAGAATAAAAGAACTAAGAGATTTAATG...    A   H4N6  \n",
       "55128   ATGAAAGTAAAACTACTGGTTCTGTTATGTACATTTACAGCTACAT...    A   H1N1  \n",
       "587361  ATATTGAAAGATGAGCCTTCTAACCGAGGTCGAAACGTATGTTCTC...    A   H3N2  \n",
       "588190  TCTTAAAATGTCGCTGTTTGGAGACACAATTGCCTACCTGCTTTCA...    B   none  \n",
       "621327  TCAATTATATTCAGCATGGAAAGAATAAAAGAACTACGGAATCTAA...    A   H3N2  \n",
       "421901  ATGAGCCTTCTAACCGAGGTCGAAACGTATGTTCTCTCTATCGTTC...    A   H3N2  \n",
       "10392   ATGTCACAGTCTCGCACCCGCGAGATACTCACCAAAACCACCGTGG...    A   H4N6  \n",
       "252467  GGCTCTGTTTCTCTCACCATTTCCACAATATGCTTCTTCATGCAAA...    A   H3N2  \n",
       "433033  ATGACATTGGCTAAAATTGAATTGTTAAAACAACTGTTAAGGGACA...    B   none  \n",
       "444043  AGCGAAAGCAGGGTAGATAATCACTCAATGAGTGACATCGAAGCCA...    A   H3N2  \n",
       "36005   AAATGGAAAAAATAGTGCTTCTTTTTGCAATAGTCAGTCTTGTTAA...    A   H5N1  \n",
       "342677  ATGGCATCTCAAGGTACCAAACGATCATATGAACAAATGGAGACTG...    A   H1N2  \n",
       "140445  TCACAGTCTCGCACTCGCGAGATACTAACAAAGACCACCGTGGACC...    A   H5N1  \n",
       "13605   CAAAAGCAGGAGTAAAGATGAATCCAAATCAAAAGATAATAACGAT...    A   H3N2  \n",
       "21340   ATGAATACAAATCAAAGAATAATAACCATTGGGACGGTTTGCCTGA...    A   H1N1  \n",
       "213246  ATGAGTCTTCTAACCGAGGTCGAAACGTACGTTCTTTCTATCATCC...    A   H1N1  "
      ]
     },
     "execution_count": 60,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "uniflna_smalltest"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 61,
   "metadata": {},
   "outputs": [],
   "source": [
    "uniflna_test2['type'] = uniflna_test2.apply(virus_type, axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 62,
   "metadata": {},
   "outputs": [],
   "source": [
    "uniflna_test2['HNtype'] = uniflna_test2.apply(HN_type, axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "metadata": {},
   "outputs": [],
   "source": [
    "# export csv of dataframe\n",
    "uniflna_test2.to_csv(\"uniflna_test2.csv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}