8_generate_cluster_consensus_sequences.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 34,
   "id": "7f27ccce",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Imports\n",
    "import os\n",
    "import subprocess\n",
    "from Bio.Seq import Seq\n",
    "from Bio.motifs import Motif\n",
    "from Bio import SeqIO, AlignIO\n",
    "from Bio.SeqRecord import SeqRecord\n",
    "from collections import defaultdict\n",
    "from Bio.Align import MultipleSeqAlignment"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "id": "e8e94df9",
   "metadata": {},
   "outputs": [],
   "source": [
    "def write_temp_fasta(records, filepath):\n",
    "    SeqIO.write(records, filepath, \"fasta\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "id": "0641362a",
   "metadata": {},
   "outputs": [],
   "source": [
    "def process_sequences_by_class(input_fasta):\n",
    "    sequences_by_class = defaultdict(list)\n",
    "    for record in SeqIO.parse(input_fasta, \"fasta\"):\n",
    "        sequence_class = record.description.split('|')[-1]\n",
    "        sequences_by_class[sequence_class].append(record)\n",
    "    return sequences_by_class"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "id": "18343f2a",
   "metadata": {},
   "outputs": [],
   "source": [
    "def align_sequences(muscle_path, input_fasta, output_fasta):\n",
    "    # Construct the MUSCLE command\n",
    "    command = f'\"{muscle_path}\" -in \"{input_fasta}\" -out \"{output_fasta}\" -maxiters 100'\n",
    "    print(f\"Executing command: {command}\")\n",
    "    \n",
    "    try:\n",
    "        # Execute the command using subprocess.run\n",
    "        subprocess.run(command, check=True, shell=True, text=True, stderr=subprocess.PIPE)\n",
    "        print(f\"MUSCLE successfully executed for {input_fasta}\")\n",
    "    \n",
    "    except subprocess.CalledProcessError as e:\n",
    "        # Handle the error and display the stderr output\n",
    "        print(f\"Error executing MUSCLE for {input_fasta}: {e.stderr}\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "id": "d1d1c9cb",
   "metadata": {},
   "outputs": [],
   "source": [
    "def read_alignment(aligned_fasta):\n",
    "    return AlignIO.read(aligned_fasta, \"fasta\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "id": "95e3e3dd",
   "metadata": {},
   "outputs": [],
   "source": [
    "def generate_consensus(align):\n",
    "    alignment = align.alignment\n",
    "    # Include all possible characters (including '-' for gaps) in the alphabet\n",
    "    alphabet = set(''.join(str(record.seq) for record in align))\n",
    "    motif = Motif(''.join(sorted(alphabet)), alignment)\n",
    "    consensus = motif.counts.calculate_consensus(identity=0.7)\n",
    "    return consensus"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "id": "c1370bef",
   "metadata": {},
   "outputs": [],
   "source": [
    "def generate_consensus_per_class(input_fasta, output_fasta, muscle_path):\n",
    "    sequences_by_class = process_sequences_by_class(input_fasta)\n",
    "    print(sequences_by_class)\n",
    "    consensus_records = []\n",
    "    for cls, records in sequences_by_class.items():\n",
    "        temp_fasta = f'temp_{cls}.fasta'\n",
    "        aligned_fasta = f'{cls}_aligned.fasta'\n",
    "        \n",
    "        if not os.path.exists(aligned_fasta):\n",
    "            write_temp_fasta(records, temp_fasta)\n",
    "            align_sequences(muscle_path, temp_fasta, aligned_fasta)\n",
    "            if os.path.exists(temp_fasta):\n",
    "                os.remove(temp_fasta)\n",
    "\n",
    "        if os.path.exists(aligned_fasta):\n",
    "            align = read_alignment(aligned_fasta)\n",
    "            os.remove(aligned_fasta)\n",
    "        else:\n",
    "            print(f\"Failed to generate alignment for class {cls}\")\n",
    "            continue\n",
    "\n",
    "        consensus = generate_consensus(align)\n",
    "        # Create a Seq object from the consensus string\n",
    "        consensus_seq = Seq(str(consensus))\n",
    "        consensus_records.append(SeqRecord(consensus_seq, id=cls, description=\"\"))\n",
    "\n",
    "    SeqIO.write(consensus_records, output_fasta, \"fasta\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "id": "d3bb2818",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "defaultdict(<class 'list'>, {'C5': [SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGATCCGCAATACCAGCACCCCGCACGTC...TGA'), id='39.B1.W11.12|C5', name='39.B1.W11.12|C5', description='39.B1.W11.12|C5', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGATCCGCAATACCAGCACCCCGCACGTC...TGA'), id='GQ396662.1|C5', name='GQ396662.1|C5', description='GQ396662.1|C5', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGATCCGCAATACCAGCACCCCGCACGTC...TGA'), id='02-520-S1a|C5', name='02-520-S1a|C5', description='02-520-S1a|C5', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGATCCGCAATACCAGCACCCCGCACGTC...TGA'), id='KP745648.1|C5', name='KP745648.1|C5', description='KP745648.1|C5', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGATCCGCAATACCAGCACCCCGCACGTC...TGA'), id='KP745677.1|C5', name='KP745677.1|C5', description='KP745677.1|C5', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGATCCGCAATACCAGCACCCCGCACGTC...TGA'), id='KY490070.1|C5', name='KY490070.1|C5', description='KY490070.1|C5', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGATCCGCAATACCAGCACCCCGCACGTC...TGA'), id='KP745721.1|C5', name='KP745721.1|C5', description='KP745721.1|C5', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGATCCGCAATACCAGCACCCCGCACGTC...TGA'), id='GU180094.1|C5', name='GU180094.1|C5', description='GU180094.1|C5', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGATCCGCAATACCAGCACCCCGCACGTC...TGA'), id='KP745685.1|C5', name='KP745685.1|C5', description='KP745685.1|C5', dbxrefs=[])], 'TOWNE': [SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='30.B1.W11.08|TOWNE', name='30.B1.W11.08|TOWNE', description='30.B1.W11.08|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='30.M1.W04.09|TOWNE', name='30.M1.W04.09|TOWNE', description='30.M1.W04.09|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KY490072.1|TOWNE', name='KY490072.1|TOWNE', description='KY490072.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KX544836.1|TOWNE', name='KX544836.1|TOWNE', description='KX544836.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KX544831.1|TOWNE', name='KX544831.1|TOWNE', description='KX544831.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KX101024.1|TOWNE', name='KX101024.1|TOWNE', description='KX101024.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KX101023.1|TOWNE', name='KX101023.1|TOWNE', description='KX101023.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KX101022.1|TOWNE', name='KX101022.1|TOWNE', description='KX101022.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KJ361971.1|TOWNE', name='KJ361971.1|TOWNE', description='KJ361971.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KF493877.1|TOWNE', name='KF493877.1|TOWNE', description='KF493877.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KF493876.1|TOWNE', name='KF493876.1|TOWNE', description='KF493876.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='GU980198.1|TOWNE', name='GU980198.1|TOWNE', description='GU980198.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='FJ616285.1|TOWNE', name='FJ616285.1|TOWNE', description='FJ616285.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='GQ121041.1|TOWNE', name='GQ121041.1|TOWNE', description='GQ121041.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KP745697.1|TOWNE', name='KP745697.1|TOWNE', description='KP745697.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACGCCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KT726941.2|TOWNE', name='KT726941.2|TOWNE', description='KT726941.2|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KT726940.2|TOWNE', name='KT726940.2|TOWNE', description='KT726940.2|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='OU342917.1|TOWNE', name='OU342917.1|TOWNE', description='OU342917.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='OU342914.1|TOWNE', name='OU342914.1|TOWNE', description='OU342914.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='OU342913.1|TOWNE', name='OU342913.1|TOWNE', description='OU342913.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='MW528459.1|TOWNE', name='MW528459.1|TOWNE', description='MW528459.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='MT649474.1|TOWNE', name='MT649474.1|TOWNE', description='MT649474.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='MT070139.1|TOWNE', name='MT070139.1|TOWNE', description='MT070139.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='MT044481.1|TOWNE', name='MT044481.1|TOWNE', description='MT044481.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='MT044477.1|TOWNE', name='MT044477.1|TOWNE', description='MT044477.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='MK290743.1|TOWNE', name='MK290743.1|TOWNE', description='MK290743.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KY490080.1|TOWNE', name='KY490080.1|TOWNE', description='KY490080.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KY490063.1|TOWNE', name='KY490063.1|TOWNE', description='KY490063.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KU550090.1|TOWNE', name='KU550090.1|TOWNE', description='KU550090.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KR534203.1|TOWNE', name='KR534203.1|TOWNE', description='KR534203.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KR534196.1|TOWNE', name='KR534196.1|TOWNE', description='KR534196.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KP745718.1|TOWNE', name='KP745718.1|TOWNE', description='KP745718.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KP745717.1|TOWNE', name='KP745717.1|TOWNE', description='KP745717.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KP745710.1|TOWNE', name='KP745710.1|TOWNE', description='KP745710.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KP745701.1|TOWNE', name='KP745701.1|TOWNE', description='KP745701.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KP745694.1|TOWNE', name='KP745694.1|TOWNE', description='KP745694.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KP745676.1|TOWNE', name='KP745676.1|TOWNE', description='KP745676.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KP745673.1|TOWNE', name='KP745673.1|TOWNE', description='KP745673.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KP745671.1|TOWNE', name='KP745671.1|TOWNE', description='KP745671.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KP745667.1|TOWNE', name='KP745667.1|TOWNE', description='KP745667.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KP745663.1|TOWNE', name='KP745663.1|TOWNE', description='KP745663.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KP745641.1|TOWNE', name='KP745641.1|TOWNE', description='KP745641.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KP745633.1|TOWNE', name='KP745633.1|TOWNE', description='KP745633.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KJ361949.1|TOWNE', name='KJ361949.1|TOWNE', description='KJ361949.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='KJ361947.1|TOWNE', name='KJ361947.1|TOWNE', description='KJ361947.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='JX512197.1|TOWNE', name='JX512197.1|TOWNE', description='JX512197.1|TOWNE', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCAGTACTCCGCATGTC...TGA'), id='JX512201.1|TOWNE', name='JX512201.1|TOWNE', description='JX512201.1|TOWNE', dbxrefs=[])], 'MERLIN': [SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACGAACCGCAGTACCGATACTCCGCACATCAAT...TGA'), id='OK000912.1|MERLIN', name='OK000912.1|MERLIN', description='OK000912.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='02-519-Pc|MERLIN', name='02-519-Pc|MERLIN', description='02-519-Pc|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='02-520-Pa|MERLIN', name='02-520-Pa|MERLIN', description='02-520-Pa|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='02-520-Pc|MERLIN', name='02-520-Pc|MERLIN', description='02-520-Pc|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='02-519-Pd|MERLIN', name='02-519-Pd|MERLIN', description='02-519-Pd|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='OV100760.1|MERLIN', name='OV100760.1|MERLIN', description='OV100760.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='OU342918.1|MERLIN', name='OU342918.1|MERLIN', description='OU342918.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='ON119191.1|MERLIN', name='ON119191.1|MERLIN', description='ON119191.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KT726945.2|MERLIN', name='KT726945.2|MERLIN', description='KT726945.2|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745702.1|MERLIN', name='KP745702.1|MERLIN', description='KP745702.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745659.1|MERLIN', name='KP745659.1|MERLIN', description='KP745659.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KJ872540.1|MERLIN', name='KJ872540.1|MERLIN', description='KJ872540.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745634.1|MERLIN', name='KP745634.1|MERLIN', description='KP745634.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCGTATGCTGTTCCGCCCAGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='OK000909.1|MERLIN', name='OK000909.1|MERLIN', description='OK000909.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='NC_006273.2|MERLIN', name='NC_006273.2|MERLIN', description='NC_006273.2|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KU221100.1|MERLIN', name='KU221100.1|MERLIN', description='KU221100.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KU221099.1|MERLIN', name='KU221099.1|MERLIN', description='KU221099.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KU221098.1|MERLIN', name='KU221098.1|MERLIN', description='KU221098.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KU221097.1|MERLIN', name='KU221097.1|MERLIN', description='KU221097.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KU221096.1|MERLIN', name='KU221096.1|MERLIN', description='KU221096.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KU221095.1|MERLIN', name='KU221095.1|MERLIN', description='KU221095.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KU221094.1|MERLIN', name='KU221094.1|MERLIN', description='KU221094.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KU221093.1|MERLIN', name='KU221093.1|MERLIN', description='KU221093.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KU221092.1|MERLIN', name='KU221092.1|MERLIN', description='KU221092.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KU221091.1|MERLIN', name='KU221091.1|MERLIN', description='KU221091.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KU221090.1|MERLIN', name='KU221090.1|MERLIN', description='KU221090.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KT726952.2|MERLIN', name='KT726952.2|MERLIN', description='KT726952.2|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973642.1|MERLIN', name='KP973642.1|MERLIN', description='KP973642.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973641.1|MERLIN', name='KP973641.1|MERLIN', description='KP973641.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973640.1|MERLIN', name='KP973640.1|MERLIN', description='KP973640.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973639.1|MERLIN', name='KP973639.1|MERLIN', description='KP973639.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973638.1|MERLIN', name='KP973638.1|MERLIN', description='KP973638.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973637.1|MERLIN', name='KP973637.1|MERLIN', description='KP973637.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973636.1|MERLIN', name='KP973636.1|MERLIN', description='KP973636.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973635.1|MERLIN', name='KP973635.1|MERLIN', description='KP973635.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973634.1|MERLIN', name='KP973634.1|MERLIN', description='KP973634.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973633.1|MERLIN', name='KP973633.1|MERLIN', description='KP973633.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973632.1|MERLIN', name='KP973632.1|MERLIN', description='KP973632.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973631.1|MERLIN', name='KP973631.1|MERLIN', description='KP973631.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973630.1|MERLIN', name='KP973630.1|MERLIN', description='KP973630.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973629.1|MERLIN', name='KP973629.1|MERLIN', description='KP973629.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973628.1|MERLIN', name='KP973628.1|MERLIN', description='KP973628.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973627.1|MERLIN', name='KP973627.1|MERLIN', description='KP973627.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973626.1|MERLIN', name='KP973626.1|MERLIN', description='KP973626.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973625.1|MERLIN', name='KP973625.1|MERLIN', description='KP973625.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973624.1|MERLIN', name='KP973624.1|MERLIN', description='KP973624.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP973623.1|MERLIN', name='KP973623.1|MERLIN', description='KP973623.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KM192302.1|MERLIN', name='KM192302.1|MERLIN', description='KM192302.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KM192301.1|MERLIN', name='KM192301.1|MERLIN', description='KM192301.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KM192300.1|MERLIN', name='KM192300.1|MERLIN', description='KM192300.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KM192299.1|MERLIN', name='KM192299.1|MERLIN', description='KM192299.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KM192298.1|MERLIN', name='KM192298.1|MERLIN', description='KM192298.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='AY446894.2|MERLIN', name='AY446894.2|MERLIN', description='AY446894.2|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='GU179001.1|MERLIN', name='GU179001.1|MERLIN', description='GU179001.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='02-579-Pa|MERLIN', name='02-579-Pa|MERLIN', description='02-579-Pa|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='02-505-Pa|MERLIN', name='02-505-Pa|MERLIN', description='02-505-Pa|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='63.M1.W02.30|MERLIN', name='63.M1.W02.30|MERLIN', description='63.M1.W02.30|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='63.B1.W09.29|MERLIN', name='63.B1.W09.29|MERLIN', description='63.B1.W09.29|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='05.M1.W08.03|MERLIN', name='05.M1.W08.03|MERLIN', description='05.M1.W08.03|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='05.B1.W14.04|MERLIN', name='05.B1.W14.04|MERLIN', description='05.B1.W14.04|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='OV100762.1|MERLIN', name='OV100762.1|MERLIN', description='OV100762.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='OU342912.1|MERLIN', name='OU342912.1|MERLIN', description='OU342912.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='ON119199.1|MERLIN', name='ON119199.1|MERLIN', description='ON119199.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='ON119197.1|MERLIN', name='ON119197.1|MERLIN', description='ON119197.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='ON119195.1|MERLIN', name='ON119195.1|MERLIN', description='ON119195.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='ON119190.1|MERLIN', name='ON119190.1|MERLIN', description='ON119190.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='MW980585.1|MERLIN', name='MW980585.1|MERLIN', description='MW980585.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='MW197157.1|MERLIN', name='MW197157.1|MERLIN', description='MW197157.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='MT649472.1|MERLIN', name='MT649472.1|MERLIN', description='MT649472.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='MT649471.1|MERLIN', name='MT649471.1|MERLIN', description='MT649471.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='MT649469.1|MERLIN', name='MT649469.1|MERLIN', description='MT649469.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='MT649468.1|MERLIN', name='MT649468.1|MERLIN', description='MT649468.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='MT044483.1|MERLIN', name='MT044483.1|MERLIN', description='MT044483.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='MF084223.1|MERLIN', name='MF084223.1|MERLIN', description='MF084223.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KY490086.1|MERLIN', name='KY490086.1|MERLIN', description='KY490086.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KY490084.1|MERLIN', name='KY490084.1|MERLIN', description='KY490084.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATTATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KY490083.1|MERLIN', name='KY490083.1|MERLIN', description='KY490083.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KY490081.1|MERLIN', name='KY490081.1|MERLIN', description='KY490081.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KY490068.1|MERLIN', name='KY490068.1|MERLIN', description='KY490068.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KY490065.1|MERLIN', name='KY490065.1|MERLIN', description='KY490065.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KY490064.1|MERLIN', name='KY490064.1|MERLIN', description='KY490064.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KY490062.1|MERLIN', name='KY490062.1|MERLIN', description='KY490062.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KY123649.1|MERLIN', name='KY123649.1|MERLIN', description='KY123649.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KX544838.1|MERLIN', name='KX544838.1|MERLIN', description='KX544838.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KX544837.1|MERLIN', name='KX544837.1|MERLIN', description='KX544837.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KU550089.1|MERLIN', name='KU550089.1|MERLIN', description='KU550089.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KU550087.1|MERLIN', name='KU550087.1|MERLIN', description='KU550087.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KT959235.1|MERLIN', name='KT959235.1|MERLIN', description='KT959235.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KT726954.2|MERLIN', name='KT726954.2|MERLIN', description='KT726954.2|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KT726949.2|MERLIN', name='KT726949.2|MERLIN', description='KT726949.2|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KT726946.2|MERLIN', name='KT726946.2|MERLIN', description='KT726946.2|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KT634296.1|MERLIN', name='KT634296.1|MERLIN', description='KT634296.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KR534213.1|MERLIN', name='KR534213.1|MERLIN', description='KR534213.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KR534212.1|MERLIN', name='KR534212.1|MERLIN', description='KR534212.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KR534207.1|MERLIN', name='KR534207.1|MERLIN', description='KR534207.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KR534205.1|MERLIN', name='KR534205.1|MERLIN', description='KR534205.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745724.1|MERLIN', name='KP745724.1|MERLIN', description='KP745724.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745712.1|MERLIN', name='KP745712.1|MERLIN', description='KP745712.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745711.1|MERLIN', name='KP745711.1|MERLIN', description='KP745711.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745706.1|MERLIN', name='KP745706.1|MERLIN', description='KP745706.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745695.1|MERLIN', name='KP745695.1|MERLIN', description='KP745695.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745691.1|MERLIN', name='KP745691.1|MERLIN', description='KP745691.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745684.1|MERLIN', name='KP745684.1|MERLIN', description='KP745684.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAATCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745672.1|MERLIN', name='KP745672.1|MERLIN', description='KP745672.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745664.1|MERLIN', name='KP745664.1|MERLIN', description='KP745664.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745662.1|MERLIN', name='KP745662.1|MERLIN', description='KP745662.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745661.1|MERLIN', name='KP745661.1|MERLIN', description='KP745661.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745653.1|MERLIN', name='KP745653.1|MERLIN', description='KP745653.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745647.1|MERLIN', name='KP745647.1|MERLIN', description='KP745647.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745646.1|MERLIN', name='KP745646.1|MERLIN', description='KP745646.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745643.1|MERLIN', name='KP745643.1|MERLIN', description='KP745643.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745638.1|MERLIN', name='KP745638.1|MERLIN', description='KP745638.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745637.1|MERLIN', name='KP745637.1|MERLIN', description='KP745637.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KP745635.1|MERLIN', name='KP745635.1|MERLIN', description='KP745635.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KJ872539.1|MERLIN', name='KJ872539.1|MERLIN', description='KJ872539.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KJ361958.1|MERLIN', name='KJ361958.1|MERLIN', description='KJ361958.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KJ361954.1|MERLIN', name='KJ361954.1|MERLIN', description='KJ361954.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KJ361953.1|MERLIN', name='KJ361953.1|MERLIN', description='KJ361953.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KJ361948.1|MERLIN', name='KJ361948.1|MERLIN', description='KJ361948.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KC519322.1|MERLIN', name='KC519322.1|MERLIN', description='KC519322.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='KC519320.1|MERLIN', name='KC519320.1|MERLIN', description='KC519320.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='JX512202.1|MERLIN', name='JX512202.1|MERLIN', description='JX512202.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='GU179291.1|MERLIN', name='GU179291.1|MERLIN', description='GU179291.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='GQ466044.1|MERLIN', name='GQ466044.1|MERLIN', description='GQ466044.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCGATACTCCGCATGTC...TGA'), id='GU179288.1|MERLIN', name='GU179288.1|MERLIN', description='GU179288.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACGCCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='MT044485.1|MERLIN', name='MT044485.1|MERLIN', description='MT044485.1|MERLIN', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGAACCGCAGTACCAGTACTCCGCATGTC...TGA'), id='39.M1.W05.14|MERLIN', name='39.M1.W05.14|MERLIN', description='39.M1.W05.14|MERLIN', dbxrefs=[])], 'AD169': [SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCTGCAACAACACTCCTCCGCACATCAAT...TGA'), id='OV100763.1|AD169', name='OV100763.1|AD169', description='OV100763.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='PP803491.1|AD169', name='PP803491.1|AD169', description='PP803491.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='OQ466312.1|AD169', name='OQ466312.1|AD169', description='OQ466312.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='OK000910.1|AD169', name='OK000910.1|AD169', description='OK000910.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MW528461.1|AD169', name='MW528461.1|AD169', description='MW528461.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MW528460.1|AD169', name='MW528460.1|AD169', description='MW528460.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MW528458.1|AD169', name='MW528458.1|AD169', description='MW528458.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MW197154.1|AD169', name='MW197154.1|AD169', description='MW197154.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MT070142.1|AD169', name='MT070142.1|AD169', description='MT070142.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MT070141.1|AD169', name='MT070141.1|AD169', description='MT070141.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MT070138.1|AD169', name='MT070138.1|AD169', description='MT070138.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MT044482.1|AD169', name='MT044482.1|AD169', description='MT044482.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MT044479.1|AD169', name='MT044479.1|AD169', description='MT044479.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MK425187.1|AD169', name='MK425187.1|AD169', description='MK425187.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MK422176.1|AD169', name='MK422176.1|AD169', description='MK422176.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MF084224.1|AD169', name='MF084224.1|AD169', description='MF084224.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY490088.1|AD169', name='KY490088.1|AD169', description='KY490088.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY490087.1|AD169', name='KY490087.1|AD169', description='KY490087.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY490085.1|AD169', name='KY490085.1|AD169', description='KY490085.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY490082.1|AD169', name='KY490082.1|AD169', description='KY490082.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY490079.1|AD169', name='KY490079.1|AD169', description='KY490079.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY490078.1|AD169', name='KY490078.1|AD169', description='KY490078.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY490075.1|AD169', name='KY490075.1|AD169', description='KY490075.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY490074.1|AD169', name='KY490074.1|AD169', description='KY490074.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY490073.1|AD169', name='KY490073.1|AD169', description='KY490073.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY490069.1|AD169', name='KY490069.1|AD169', description='KY490069.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY490061.1|AD169', name='KY490061.1|AD169', description='KY490061.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY123650.1|AD169', name='KY123650.1|AD169', description='KY123650.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KX544841.1|AD169', name='KX544841.1|AD169', description='KX544841.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KT726955.2|AD169', name='KT726955.2|AD169', description='KT726955.2|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KT726947.2|AD169', name='KT726947.2|AD169', description='KT726947.2|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KT726943.2|AD169', name='KT726943.2|AD169', description='KT726943.2|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATTATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KR534210.1|AD169', name='KR534210.1|AD169', description='KR534210.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KR534200.1|AD169', name='KR534200.1|AD169', description='KR534200.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KR534199.1|AD169', name='KR534199.1|AD169', description='KR534199.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745728.1|AD169', name='KP745728.1|AD169', description='KP745728.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745727.1|AD169', name='KP745727.1|AD169', description='KP745727.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745719.1|AD169', name='KP745719.1|AD169', description='KP745719.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745709.1|AD169', name='KP745709.1|AD169', description='KP745709.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745705.1|AD169', name='KP745705.1|AD169', description='KP745705.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745696.1|AD169', name='KP745696.1|AD169', description='KP745696.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745693.1|AD169', name='KP745693.1|AD169', description='KP745693.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745692.1|AD169', name='KP745692.1|AD169', description='KP745692.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745690.1|AD169', name='KP745690.1|AD169', description='KP745690.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745688.1|AD169', name='KP745688.1|AD169', description='KP745688.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745683.1|AD169', name='KP745683.1|AD169', description='KP745683.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745682.1|AD169', name='KP745682.1|AD169', description='KP745682.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745680.1|AD169', name='KP745680.1|AD169', description='KP745680.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745675.1|AD169', name='KP745675.1|AD169', description='KP745675.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745670.1|AD169', name='KP745670.1|AD169', description='KP745670.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745668.1|AD169', name='KP745668.1|AD169', description='KP745668.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745666.1|AD169', name='KP745666.1|AD169', description='KP745666.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745665.1|AD169', name='KP745665.1|AD169', description='KP745665.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745654.1|AD169', name='KP745654.1|AD169', description='KP745654.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745652.1|AD169', name='KP745652.1|AD169', description='KP745652.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745645.1|AD169', name='KP745645.1|AD169', description='KP745645.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745644.1|AD169', name='KP745644.1|AD169', description='KP745644.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745640.1|AD169', name='KP745640.1|AD169', description='KP745640.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745636.1|AD169', name='KP745636.1|AD169', description='KP745636.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KJ361969.1|AD169', name='KJ361969.1|AD169', description='KJ361969.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KJ361965.1|AD169', name='KJ361965.1|AD169', description='KJ361965.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KJ361962.1|AD169', name='KJ361962.1|AD169', description='KJ361962.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KJ361961.1|AD169', name='KJ361961.1|AD169', description='KJ361961.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KJ361955.1|AD169', name='KJ361955.1|AD169', description='KJ361955.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KF021605.1|AD169', name='KF021605.1|AD169', description='KF021605.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='JX512208.1|AD169', name='JX512208.1|AD169', description='JX512208.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='JX512207.1|AD169', name='JX512207.1|AD169', description='JX512207.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='JX512205.1|AD169', name='JX512205.1|AD169', description='JX512205.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='JX512203.1|AD169', name='JX512203.1|AD169', description='JX512203.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='HQ380895.1|AD169', name='HQ380895.1|AD169', description='HQ380895.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='GQ221975.1|AD169', name='GQ221975.1|AD169', description='GQ221975.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='GU179289.1|AD169', name='GU179289.1|AD169', description='GU179289.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KR534204.1|AD169', name='KR534204.1|AD169', description='KR534204.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='JX512206.1|AD169', name='JX512206.1|AD169', description='JX512206.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745660.1|AD169', name='KP745660.1|AD169', description='KP745660.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MK290742.1|AD169', name='MK290742.1|AD169', description='MK290742.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745726.1|AD169', name='KP745726.1|AD169', description='KP745726.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KC519321.1|AD169', name='KC519321.1|AD169', description='KC519321.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KJ872541.1|AD169', name='KJ872541.1|AD169', description='KJ872541.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745722.1|AD169', name='KP745722.1|AD169', description='KP745722.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='OU342916.1|AD169', name='OU342916.1|AD169', description='OU342916.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='OU342911.1|AD169', name='OU342911.1|AD169', description='OU342911.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MZ327301.1|AD169', name='MZ327301.1|AD169', description='MZ327301.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MW439039.1|AD169', name='MW439039.1|AD169', description='MW439039.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MW439038.1|AD169', name='MW439038.1|AD169', description='MW439038.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MT649470.1|AD169', name='MT649470.1|AD169', description='MT649470.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MT070140.1|AD169', name='MT070140.1|AD169', description='MT070140.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MT044484.1|AD169', name='MT044484.1|AD169', description='MT044484.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MT044478.1|AD169', name='MT044478.1|AD169', description='MT044478.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='MT044476.1|AD169', name='MT044476.1|AD169', description='MT044476.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY490071.1|AD169', name='KY490071.1|AD169', description='KY490071.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KY123653.1|AD169', name='KY123653.1|AD169', description='KY123653.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KX544835.1|AD169', name='KX544835.1|AD169', description='KX544835.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KX544834.1|AD169', name='KX544834.1|AD169', description='KX544834.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KX544832.1|AD169', name='KX544832.1|AD169', description='KX544832.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KR534198.1|AD169', name='KR534198.1|AD169', description='KR534198.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745723.1|AD169', name='KP745723.1|AD169', description='KP745723.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745720.1|AD169', name='KP745720.1|AD169', description='KP745720.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745714.1|AD169', name='KP745714.1|AD169', description='KP745714.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745699.1|AD169', name='KP745699.1|AD169', description='KP745699.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745681.1|AD169', name='KP745681.1|AD169', description='KP745681.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KP745657.1|AD169', name='KP745657.1|AD169', description='KP745657.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KJ872542.1|AD169', name='KJ872542.1|AD169', description='KJ872542.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KJ361959.1|AD169', name='KJ361959.1|AD169', description='KJ361959.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KJ361957.1|AD169', name='KJ361957.1|AD169', description='KJ361957.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KJ361956.1|AD169', name='KJ361956.1|AD169', description='KJ361956.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KF297339.1|AD169', name='KF297339.1|AD169', description='KF297339.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KC519323.1|AD169', name='KC519323.1|AD169', description='KC519323.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='KC519319.1|AD169', name='KC519319.1|AD169', description='KC519319.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='JX512204.1|AD169', name='JX512204.1|AD169', description='JX512204.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='JX512198.1|AD169', name='JX512198.1|AD169', description='JX512198.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='GQ221974.1|AD169', name='GQ221974.1|AD169', description='GQ221974.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='GQ221973.1|AD169', name='GQ221973.1|AD169', description='GQ221973.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='BK000394.5|AD169', name='BK000394.5|AD169', description='BK000394.5|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACAACACTCCTCCGCACATCAAT...TGA'), id='FJ527563.1|AD169', name='FJ527563.1|AD169', description='FJ527563.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGACCCGCAACAACACTCCTCCGCACATC...TGA'), id='KJ361950.1|AD169', name='KJ361950.1|AD169', description='KJ361950.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCACGACCCGCAACAACACTCCTCCGCACATC...TGA'), id='KJ361951.1|AD169', name='KJ361951.1|AD169', description='KJ361951.1|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCTCGACCCGCAACAACACTCCTCCGCACATC...TGA'), id='04.B1.W14.01|AD169', name='04.B1.W14.01|AD169', description='04.B1.W14.01|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCTCGACCCGCAACAACACTCCTCCGCACATC...TGA'), id='56.B1.W09.22|AD169', name='56.B1.W09.22|AD169', description='56.B1.W09.22|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCTCGACCCGCAACAACACTCCTCCGCACATC...TGA'), id='66.M1.W02.24|AD169', name='66.M1.W02.24|AD169', description='66.M1.W02.24|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCTCGACCCGCAACACCACTCCTCCGCACATC...TGA'), id='04.M1.W09.02|AD169', name='04.M1.W09.02|AD169', description='04.M1.W09.02|AD169', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCTCGACCCGCAACACCACTCCTCCGCACATC...TGA'), id='56.M1.W03.21|AD169', name='56.M1.W03.21|AD169', description='56.M1.W03.21|AD169', dbxrefs=[])], 'TOLEDO': [SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745725.1|TOLEDO', name='KP745725.1|TOLEDO', description='KP745725.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='OQ466311.1|TOLEDO', name='OQ466311.1|TOLEDO', description='OQ466311.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='02-579-M|TOLEDO', name='02-579-M|TOLEDO', description='02-579-M|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='02-519-S1a|TOLEDO', name='02-519-S1a|TOLEDO', description='02-519-S1a|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='02-579-S1b|TOLEDO', name='02-579-S1b|TOLEDO', description='02-579-S1b|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='02-579-S1a|TOLEDO', name='02-579-S1a|TOLEDO', description='02-579-S1a|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='MT649473.1|TOLEDO', name='MT649473.1|TOLEDO', description='MT649473.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='MT044480.1|TOLEDO', name='MT044480.1|TOLEDO', description='MT044480.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KY490076.1|TOLEDO', name='KY490076.1|TOLEDO', description='KY490076.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KY002201.1|TOLEDO', name='KY002201.1|TOLEDO', description='KY002201.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KY002200.1|TOLEDO', name='KY002200.1|TOLEDO', description='KY002200.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KX101021.1|TOLEDO', name='KX101021.1|TOLEDO', description='KX101021.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KT726950.2|TOLEDO', name='KT726950.2|TOLEDO', description='KT726950.2|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KT726944.2|TOLEDO', name='KT726944.2|TOLEDO', description='KT726944.2|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745713.1|TOLEDO', name='KP745713.1|TOLEDO', description='KP745713.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745669.1|TOLEDO', name='KP745669.1|TOLEDO', description='KP745669.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745656.1|TOLEDO', name='KP745656.1|TOLEDO', description='KP745656.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='GU937742.2|TOLEDO', name='GU937742.2|TOLEDO', description='GU937742.2|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KJ426589.1|TOLEDO', name='KJ426589.1|TOLEDO', description='KJ426589.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='MW197156.1|TOLEDO', name='MW197156.1|TOLEDO', description='MW197156.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745715.1|TOLEDO', name='KP745715.1|TOLEDO', description='KP745715.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KR534209.1|TOLEDO', name='KR534209.1|TOLEDO', description='KR534209.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='MW528463.1|TOLEDO', name='MW528463.1|TOLEDO', description='MW528463.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='OU342915.1|TOLEDO', name='OU342915.1|TOLEDO', description='OU342915.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='OU342900.1|TOLEDO', name='OU342900.1|TOLEDO', description='OU342900.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAATACCACTACTCCGCACATCAAT...TGA'), id='OK000911.1|TOLEDO', name='OK000911.1|TOLEDO', description='OK000911.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='MW528464.1|TOLEDO', name='MW528464.1|TOLEDO', description='MW528464.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='MW528462.1|TOLEDO', name='MW528462.1|TOLEDO', description='MW528462.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='MW197155.1|TOLEDO', name='MW197155.1|TOLEDO', description='MW197155.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='MK290744.1|TOLEDO', name='MK290744.1|TOLEDO', description='MK290744.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KY490077.1|TOLEDO', name='KY490077.1|TOLEDO', description='KY490077.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KY490067.1|TOLEDO', name='KY490067.1|TOLEDO', description='KY490067.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KY490066.1|TOLEDO', name='KY490066.1|TOLEDO', description='KY490066.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KY123652.1|TOLEDO', name='KY123652.1|TOLEDO', description='KY123652.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KY123651.1|TOLEDO', name='KY123651.1|TOLEDO', description='KY123651.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KX544833.1|TOLEDO', name='KX544833.1|TOLEDO', description='KX544833.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KU550088.1|TOLEDO', name='KU550088.1|TOLEDO', description='KU550088.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KT726953.2|TOLEDO', name='KT726953.2|TOLEDO', description='KT726953.2|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KT726951.2|TOLEDO', name='KT726951.2|TOLEDO', description='KT726951.2|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KT726948.2|TOLEDO', name='KT726948.2|TOLEDO', description='KT726948.2|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KT726942.2|TOLEDO', name='KT726942.2|TOLEDO', description='KT726942.2|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KR534211.1|TOLEDO', name='KR534211.1|TOLEDO', description='KR534211.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KR534208.1|TOLEDO', name='KR534208.1|TOLEDO', description='KR534208.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KR534206.1|TOLEDO', name='KR534206.1|TOLEDO', description='KR534206.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KR534202.1|TOLEDO', name='KR534202.1|TOLEDO', description='KR534202.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KR534201.1|TOLEDO', name='KR534201.1|TOLEDO', description='KR534201.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KR534197.1|TOLEDO', name='KR534197.1|TOLEDO', description='KR534197.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745716.1|TOLEDO', name='KP745716.1|TOLEDO', description='KP745716.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745708.1|TOLEDO', name='KP745708.1|TOLEDO', description='KP745708.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745707.1|TOLEDO', name='KP745707.1|TOLEDO', description='KP745707.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745704.1|TOLEDO', name='KP745704.1|TOLEDO', description='KP745704.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745703.1|TOLEDO', name='KP745703.1|TOLEDO', description='KP745703.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745700.1|TOLEDO', name='KP745700.1|TOLEDO', description='KP745700.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745698.1|TOLEDO', name='KP745698.1|TOLEDO', description='KP745698.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745689.1|TOLEDO', name='KP745689.1|TOLEDO', description='KP745689.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745687.1|TOLEDO', name='KP745687.1|TOLEDO', description='KP745687.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745686.1|TOLEDO', name='KP745686.1|TOLEDO', description='KP745686.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745679.1|TOLEDO', name='KP745679.1|TOLEDO', description='KP745679.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745678.1|TOLEDO', name='KP745678.1|TOLEDO', description='KP745678.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745674.1|TOLEDO', name='KP745674.1|TOLEDO', description='KP745674.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745658.1|TOLEDO', name='KP745658.1|TOLEDO', description='KP745658.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745655.1|TOLEDO', name='KP745655.1|TOLEDO', description='KP745655.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745651.1|TOLEDO', name='KP745651.1|TOLEDO', description='KP745651.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745650.1|TOLEDO', name='KP745650.1|TOLEDO', description='KP745650.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745649.1|TOLEDO', name='KP745649.1|TOLEDO', description='KP745649.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745642.1|TOLEDO', name='KP745642.1|TOLEDO', description='KP745642.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KP745639.1|TOLEDO', name='KP745639.1|TOLEDO', description='KP745639.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KJ361970.1|TOLEDO', name='KJ361970.1|TOLEDO', description='KJ361970.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KJ361968.1|TOLEDO', name='KJ361968.1|TOLEDO', description='KJ361968.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KJ361967.1|TOLEDO', name='KJ361967.1|TOLEDO', description='KJ361967.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KJ361966.1|TOLEDO', name='KJ361966.1|TOLEDO', description='KJ361966.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KJ361964.1|TOLEDO', name='KJ361964.1|TOLEDO', description='KJ361964.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KJ361963.1|TOLEDO', name='KJ361963.1|TOLEDO', description='KJ361963.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KJ361960.1|TOLEDO', name='KJ361960.1|TOLEDO', description='KJ361960.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KJ361952.1|TOLEDO', name='KJ361952.1|TOLEDO', description='KJ361952.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='KJ361946.1|TOLEDO', name='KJ361946.1|TOLEDO', description='KJ361946.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='JX512200.1|TOLEDO', name='JX512200.1|TOLEDO', description='JX512200.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='JX512199.1|TOLEDO', name='JX512199.1|TOLEDO', description='JX512199.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='GQ396663.1|TOLEDO', name='GQ396663.1|TOLEDO', description='GQ396663.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACTCGACCCGCAACACCACTACTCCGCACATCAAT...TGA'), id='GU179290.1|TOLEDO', name='GU179290.1|TOLEDO', description='GU179290.1|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCTCGACCCGCAACACCACTACTCCGCACATC...TGA'), id='27.M1.W10.07|TOLEDO', name='27.M1.W10.07|TOLEDO', description='27.M1.W10.07|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACATCTCGACCCGCAACACCACTACTCCGCACATC...TGA'), id='53.B1.W14.17|TOLEDO', name='53.B1.W14.17|TOLEDO', description='53.B1.W14.17|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCTCGACCCGCAACACCACTACTCCGCACATC...TGA'), id='27.B1.W13.06|TOLEDO', name='27.B1.W13.06|TOLEDO', description='27.B1.W13.06|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACATCTCGACCCGCAATACCACTACTCCGCACATC...TGA'), id='53.M1.W07.16|TOLEDO', name='53.M1.W07.16|TOLEDO', description='53.M1.W07.16|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACATCTCGACCCGCAACACCACTACTCCGCACATC...TGA'), id='66.B1.W09.25|TOLEDO', name='66.B1.W09.25|TOLEDO', description='66.B1.W09.25|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACACCTCGACCCGCAACACCACTACTCCGCACATC...TGA'), id='38.B1.W10.11|TOLEDO', name='38.B1.W10.11|TOLEDO', description='38.B1.W10.11|TOLEDO', dbxrefs=[]), SeqRecord(seq=Seq('ATGGACACCATCATCCACAACATCTCGACCCGCAACACCACTACTCCGCACATC...TGA'), id='38.M1.W03.10|TOLEDO', name='38.M1.W03.10|TOLEDO', description='38.M1.W03.10|TOLEDO', dbxrefs=[])]})\n",
      "Executing command: \"./lib/muscle.exe\" -in \"temp_C5.fasta\" -out \"C5_aligned.fasta\" -maxiters 100\n",
      "MUSCLE successfully executed for temp_C5.fasta\n",
      "Executing command: \"./lib/muscle.exe\" -in \"temp_TOWNE.fasta\" -out \"TOWNE_aligned.fasta\" -maxiters 100\n",
      "MUSCLE successfully executed for temp_TOWNE.fasta\n",
      "Executing command: \"./lib/muscle.exe\" -in \"temp_MERLIN.fasta\" -out \"MERLIN_aligned.fasta\" -maxiters 100\n",
      "MUSCLE successfully executed for temp_MERLIN.fasta\n",
      "Executing command: \"./lib/muscle.exe\" -in \"temp_AD169.fasta\" -out \"AD169_aligned.fasta\" -maxiters 100\n",
      "MUSCLE successfully executed for temp_AD169.fasta\n",
      "Executing command: \"./lib/muscle.exe\" -in \"temp_TOLEDO.fasta\" -out \"TOLEDO_aligned.fasta\" -maxiters 100\n",
      "MUSCLE successfully executed for temp_TOLEDO.fasta\n"
     ]
    }
   ],
   "source": [
    "# Usage\n",
    "input_fasta = 'data/fasta/UL33_all_nucleotide_sequences_annotated.fasta'\n",
    "output_fasta = 'data/fasta/UL33_cluster_consensus.fasta'\n",
    "muscle_path = './lib/muscle.exe'\n",
    "generate_consensus_per_class(input_fasta, output_fasta, muscle_path)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.12.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}