Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,530,252 C→A intergenic (‑9/+3) sgcC ← / ← sgcB putative PTS enzyme IIC component SgcC/putative PTS enzyme IIB component SgcB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,530,2520CA88.6% 169.6 / 6.6 70intergenic (‑9/+3)sgcC/sgcBputative PTS enzyme IIC component SgcC/putative PTS enzyme IIB component SgcB
Reads supporting (aligned to +/- strand):  ref base C (6/2);  new base A (41/21);  total (47/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGAC  >  NC_000913/4530005‑4530498
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
ccTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCAttt                                                                                                                                                                                                                                                      >  1:607342/1‑251 (MQ=255)
  tGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCAtttt                                                                                                                                                                                                                                                     >  1:189702/1‑250 (MQ=255)
        aaCGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGCTCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGAACGGATGTAATCAAACATAATGCCTCATTTTATTGAGt                                                                                                                                                                                                                                              >  1:403656/1‑251 (MQ=255)
        aaCGTAAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CACAACCGTTACCCATGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGCTAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGt                                                                                                                                                                                                                                              >  1:194302/1‑251 (MQ=255)
         aCGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGtt                                                                                                                                                                                                                                             >  1:188481/1‑251 (MQ=255)
                  cctggaagggttgggcagAAACACTCAATCCTCATTTTGATGGGGGGACGGAGGCTGCCGATCCACATAACGA‑GACCCGTCTCCATACCGACAAAACCAAT‑CCCAACGGTTACCCCGTCATTAAGCGCCGGAAGCCAGGGGAGGCG‑AAGAATAAGGCCAAT‑CAGGAGTATGTTGATAGGGACGAAAACGGTATCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAATAGAGc‑c                                                                                                                                                                                                                                   <  2:514754/234‑1 (MQ=255)
                   cATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑Ct                                                                                                                                                                                                                                  >  1:91302/1‑251 (MQ=255)
                              gggCCAAAACGCTCAATCATCCTTTTGAGAGGGGAACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATCCCGACAAAACCAAC‑CCCACCCGTTCCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTtgttgt                                                                                                                                                                                                                       <  1:8603/250‑1 (MQ=255)
                                         cTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGCCAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTCTTGAGTTAACAGAGC‑CTTGGTTTGTTGTTTTAACTcgtc                                                                                                                                                                                                            >  2:131965/1‑251 (MQ=255)
                                          tCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAACGCCAAT‑CACGATCATGATGATAGGGACGAAAACAGTACCGCCGAGAGACAGAATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTCAcgctgtg                                                                                                                                                                                                           >  2:120122/1‑247 (MQ=255)
                                              tcatcaCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGcgtcgtcgt                                                                                                                                                                                                       >  1:529630/1‑251 (MQ=255)
                                              tcatcaCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGcgtcgtcgt                                                                                                                                                                                                       >  1:30721/1‑251 (MQ=255)
                                               caacatttgtagggggggccagAGGCTGTCAACCCCCATAACAATCGCCAGTC‑CCATACCG‑CAAAACCAATACCCAACCTTTCCCCCGGCTTTAATCGCCGGAAGCCAGGGGATGCGGAAGA‑TAGGGCCAAT‑CACGATCATGATGAGAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGtt                                                                                                                                                                                                      <  2:191839/230‑1 (MQ=255)
                                                 tcaCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCTTCGTCGTTg                                                                                                                                                                                                     >  1:414105/1‑250 (MQ=255)
                                                     tttGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTATTTGAGTTAACAGAGC‑CTTGATTTGTGGTTTTAAAGCGTCGTCGTTGATccc                                                                                                                                                                                                >  2:310611/1‑251 (MQ=255)
                                                       tgaaaGGGGGACTGAGGCTGTCGATCGCCATGCCAATCACCAGTC‑CCATACCGACAAACCCAAT‑CCCAACCGTTACCCCTGCCTTAACCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCAAAC‑CGCTATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCt                                                                                                                                                                                               <  1:598758/246‑1 (MQ=255)
                                                          taagaggACTGAGGCTGGCGTTCGCCATGACGATCACCAGTC‑ACATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCAGGAGGATAGGGACGACAACGGTACCGCCGATGGACAGGATGTAAACGAACAAAATGCATCATTATATTGAGTTAACAG‑GCACTTGGTTGGTGGTTTTAAAGCGTAGGCGTTGACCCCagtgt                                                                                                                                                                                           >  2:21247/6‑250 (MQ=255)
                                                                  cTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAAAATAATACATCATTTTATTGATTTAACTGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGGGGGCAGggcg                                                                                                                                                                                   >  2:141740/1‑251 (MQ=255)
                                                                    gggacTGTCGATCGACGTGCCGATCACCAGCT‑CTATCCGGACAAACCAAAT‑CTCACCTGTTCCCACTGCCTTAAACGCTTGCAGCCGTGGGATGCGGAAGA‑TAAGGCAAAT‑CACGATCATGATGATAGGGACGAAAACGTAACCGCCGAGAGACAGGATGTAAGCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATGTGTTGTTTTAACGCTTCGTCGTTGATCCCTGTGAGCAGggcggc                                                                                                                                                                                 <  2:518334/247‑1 (MQ=255)
                                                                    gAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGggcggc                                                                                                                                                                                 >  1:484360/1‑251 (MQ=255)
                                                                          gTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCAt                                                                                                                                                                            <  1:216475/250‑1 (MQ=255)
                                                                            cGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTa                                                                                                                                                                          >  1:157130/1‑250 (MQ=255)
                                                                                   cATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAgg                                                                                                                                                                  >  1:62154/1‑251 (MQ=255)
                                                                                         gATCACCAGTC‑CCATACCGAAAAAAACAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAACAATAATGCCTCATTTTATTGAGATAACAGGGC‑ATTGATTTGTTGTTTTAACGCGTCGTCGTTGATCACTGTGAACAGGGCGGCGCCATTAAggggtgggaagcc                                                                                                                                                            >  2:380386/1‑239 (MQ=255)
                                                                                           tCACCAGTC‑CCATCCCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAAGCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGC‑CTTGATTTgttgtt                                                                                                                                                                                                                      <  1:429414/191‑1 (MQ=255)
                                                                                           tCACCAGTA‑CAATAACGACAAAACCAAT‑CCCAAACGTTACACCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGC‑CTTGATTTgttgtt                                                                                                                                                                                                                      >  2:429414/1‑191 (MQ=255)
                                                                                               catcc‑ccATACCGACAAAACCAAT‑CCCAACCGTTCCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAC‑CCCGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCTCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATc                                                                                                                                                      <  1:179897/247‑1 (MQ=255)
                                                                                                   c‑ccATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTg                                                                                                                                                  >  1:494341/1‑251 (MQ=255)
                                                                                                      cATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGtt                                                                                                                                                >  1:545170/1‑251 (MQ=255)
                                                                                                        tACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGAAGGGATGCCGTAATCACTGTTgg                                                                                                                                              >  2:623286/1‑251 (MQ=255)
                                                                                                              cAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGc                                                                                                                                         >  1:558659/1‑250 (MQ=255)
                                                                                                                     aaT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACCGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGC‑CTTGTTTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTACGCGTAGGGTTGCCGTAATCACTGTTGGTTCgactgacgat                                                                                                                                 >  2:482524/1‑242 (MQ=255)
                                                                                                                      aaccccAACCGTGACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑AAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATgg                                                                                                                                     <  1:482524/244‑1 (MQ=255)
                                                                                                                      aT‑CCCAACCGTTACCCCTGCCTTAATCGCCTCCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAACCGGTACCGCCGAGAGACAGGATGTAATCAAACATACTGCCTCATTTTATTGTGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTACGCGTAGGGATGCCGTACTCACTGTTGGTTCGCATGGAGGtt                                                                                                                                >  2:348164/1‑251 (MQ=255)
                                                                                                                         cccAACCGTTCCCCCTGCCTTAATCGCCAGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTAc                                                                                                                              <  1:387913/251‑1 (MQ=255)
                                                                                                                          ccAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCACACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTAc                                                                                                                              >  1:381915/1‑250 (MQ=255)
                                                                                                                            aacctttcccccTGCTTAAATCGCCTGCAGCCAGGGGATGCGGAAGA‑AAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTTCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAAt                                                                                                                           <  1:680142/243‑1 (MQ=255)
                                                                                                                              ccGTTACCCCTGCCATAACCGCCTGCAGCCAGGGGATGCGGAGGA‑TAAGGCCAAT‑CACGAGCAGGAGGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAAGCACTGTTGGTTCGCATGGAGGTTACAATaa                                                                                                                         <  1:302554/251‑1 (MQ=255)
                                                                                                                                           cTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGtt                                                                                                            >  1:568564/1‑251 (MQ=255)
                                                                                                                                             tAATCGCCTGCAGCCAGGGGCTGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGAACCGCCGAGAGAACGGATGTAATCAAACATAATGCCTCATATTATTGAGTGAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCCGGGCGGCGCCATTAAGCGTAGGGATGACGTAATCACTGTTGGGTCGCATGGAGGATACAATAAGGCCCATACCGTTAc                                                                                                          >  2:355785/1‑251 (MQ=255)
                                                                                                                                                   ccTGCAGCAAGGGGATGAGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGCGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACACAGC‑CTTGATTTGTTGTTTTAACGCGTCGGCGTTGATCCCTGTGCGCAGGTCGGCGCCATTAAGCGTAGGGCTCCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATACGGTCCATACCGTTACAGTTa                                                                                                     >  2:676107/1‑250 (MQ=255)
                                                                                                                                                     ccccccccgacgctcttccgatctAAGGCAAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTATGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAag                                                                                                   <  2:490193/228‑1 (MQ=255)
                                                                                                                                                        aGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACAATCATGATGATAAGGTTGAAAACGGTGCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAggg                                                                                               <  1:635498/251‑1 (MQ=255)
                                                                                                                                                         gCCAGGGGAGGTAGAAGA‑TAAGGCAAATTCA‑GTTCAAGATGATAGGGACGAAATCGGTACCGCCGAGAGACAGGATTTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGa                                                                                              <  2:549440/251‑1 (MQ=255)
                                                                                                                                                              gagaggCGGAAGA‑AAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGGACCGCAGAGAGACAGGATGAAATAAAACATAATGACTCATTTTATTGAGTTAAACGAGC‑CTTGATTTGTTGATTTAAAGCGTGGTCGTTGAACGATGTGAGCAGGGCGGCGCCATTAAACGAAGGGATGCCGAAATCACTGTTGGTTGGCATGGAGGTGACAATAAGGGCAATACCGtaccagttcaggggggtcgc                                                                                          >  2:479984/6‑231 (MQ=255)
                                                                                                                                                                     ggAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTACTTTTATTGAGGTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATACCTGGGAGCAGGGCGGCGCCATTCAGCGTAGGGATGCCGTCATCACTGTTGGATCGCATGGGGGTTACAATAAGGTCCATACCGTAACAGTTAAGAGGGATCTCATTCAGAc                                                                                  >  2:500899/1‑251 (MQ=255)
                                                                                                                                                                      gAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGAc                                                                                  >  1:225224/1‑250 (MQ=255)
                                                                                                                                                                             aaGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTATGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAag                                                                                                   >  1:490193/1‑227 (MQ=255)
                                                                                                                                                                             aaGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTg                                                                           >  1:60209/1‑251 (MQ=255)
                                                                                                                                                                                  cAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTGACAGTTAAGAGGGCTCTCATTACGACAGACTTGGGCgg                                                                      >  2:600500/1‑251 (MQ=255)
                                                                                                                                                                                   aaT‑CACGAACATGCTGATAGGATCGAAAACGGTACCGCAGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGCTCCATGTGTGTAGGGCGGCGCCTTTAAGAGTAGGGAGGCCGTAATAACTGTTTGATCGAGTGGAGGTTACAATAAGGTTCATACCGTTACCGTTAAGAGGGGtctgatgacgcaagcattgagcggg                                                                     >  2:637980/1‑229 (MQ=255)
                                                                                                                                                                                    aaccacgataaggagtatggggACGAAACCGGTGCCGCCGAGAGACAGGATGTAGTCAAACATAATGCCTCATTTTATTTAGTTAACAGTGC‑CTTGATTTGTTTTTTTAACGCGTCGTCGTTGGTCCCTGTTAGCAGGGCGGCTCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCTTACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCgg                                                                      <  1:637980/231‑1 (MQ=255)
                                                                                                                                                                                          gaacgtaatgatAGGGACGAAACCGTTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGc                                                               <  2:62154/244‑1 (MQ=255)
                                                                                                                                                                                          gATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGc                                                               >  1:231245/1‑251 (MQ=255)
                                                                                                                                                                                            tCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGAAATCAAACATAATGCCTCATTTTCTTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTg                                                             >  2:120549/1‑251 (MQ=255)
                                                                                                                                                                                                      gggaaaaaaacatgaactgccgAGAGTCAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGAAAGCATTGGGCGGTAGTGGCTGAAATGccgtg                                                   <  2:115541/233‑4 (MQ=255)
                                                                                                                                                                                                      aGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACCTAATGCCACATTTTGTGGAGTTAAACGAGC‑CGTGATTTGTTGTTTTAACGCGTAGTCGTTGATCACTGTGAGCAGGGCGGAGACATTAAGAGTAGGGATGCCGTCATAACAGTTGTTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGATAAGCGGGGTCTCATTCAGACAGCATTGGGCgggaggagctgcaataccctg                                                   >  2:154291/1‑232 (MQ=255)
                                                                                                                                                                                                        ggACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCCTTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACCATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGAAAGCATTGGGCGGTAGTGGCTGAAATGCCCTGtt                                                 >  2:528165/1‑251 (MQ=255)
                                                                                                                                                                                                                       cgccgAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCtt                                  >  1:545470/1‑251 (MQ=255)
                                                                                                                                                                                                                        tccgAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGCGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATAAAGCGTAGGGAGGCCATAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACATTTAAGAGGGATCTACTTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTg                                 >  2:313472/2‑251 (MQ=255)
                                                                                                                                                                                                                              agaCAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTAGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACCGTTAAGAGGGATCTCATTCAGACAGCATTGGGAGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGCACTCTTGCAGCtt                           >  1:197555/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATa                                                                                                                          >  2:118405/1‑150 (MQ=255)
                                                                                                                                                                                                                                    gATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATa                                                                                                                          <  1:118405/150‑1 (MQ=255)
                                                                                                                                                                                                                                       gTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCACCTTGTGTGCAAt                  <  1:122259/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aaCATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGtt           <  1:78475/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aCATAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATATCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGACCTGTTCAGTCAGGAACTCTTTAAGCTTGTGTGCAAACAGGGTTg          >  2:164496/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cATAATGCCTCATTTTATTGGGTTAACAGAGA‑CTTGCTTTGTTGTTTTAACGCGTCGTCTTTGATCCCTGTGAGCAGAGCGGCGCCATTAAGCGTCGGGATCCCGTACTCACTGTTGGTTCGCATGGAGGTTACAATATGGTACCTACCGTTCCAGTTAAGGAGGACCTAATTCAGACACCATTGGACGGTAGTGGCTGTAATGCCCTGTTCAGTCAGGATATCTTGCAGATTGTGTGccatactgtatg          >  2:229796/1‑239 (MQ=255)
                                                                                                                                                                                                                                                 aTAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGTATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGTGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtcagtc                                            <  2:490944/215‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTAATGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGTGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtcagtc                                            >  1:490944/1‑215 (MQ=255)
                                                                                                                                                                                                                                                   aaTGCCTCATTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGt      <  1:533454/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       cctcatTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGAc  <  1:446500/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
CCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTC‑CCATACCGACAAAACCAAT‑CCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGA‑TAAGGCCAAT‑CACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGC‑CTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGAC  >  NC_000913/4530005‑4530498

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑