Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,773,176 C→T S299F (TCT→TTT)  mtlA → mannitol‑specific PTS enzyme IICBA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,773,1760CT95.3% 225.7 / 3.1 65coding (896/1914 nt)mtlAmannitol‑specific PTS enzyme IICBA component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/28);  minor base . (1/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTAAAATC  >  NC_000913/3772932‑3773420
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
gAAGCTACCCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCCCTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCtt                                                                                                                                                                                                                                               <  2:378756/251‑1 (MQ=255)
gAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAACCAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCt                                                                                                                                                                                                                                                <  2:168578/250‑1 (MQ=255)
       aCCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGGCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTAc                                                                                                                                                                                                                                        >  2:192544/1‑251 (MQ=255)
        cccAGTTCCAGGGATGGGCGGGCGGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGGCTGCGGGCGGTGTGGCAATCATCCACGGCCTTGGGGGTATCCACGAAGTCTACTTCCCGTATGTGCTGATGAATGCTCGTCTGATCCTCGCAGTCAGCCGTGGCGGTAGGACTGGCGTGTTCACGCTGACTATCTTGGGGGGTGGTCTGGTTTCTCAGGCATCGCCGGGGT‑TATCCTTGCTGTACt                                                                                                                                                                                                                                       <  1:356396/251‑1 (MQ=255)
            ggtccaggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAACCAGTCTGCGGGCGGTGCGGCAATCATCCACTTCATGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGGGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCGCCGG‑TTTTATCCTTGCTGTCCTGGCg                                                                                                                                                                                                                                   >  2:28824/1‑251 (MQ=255)
                  ggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGaca                                                                                                                                                                                                                             <  1:112457/251‑1 (MQ=255)
                     aTGGGGGTGCTGCTGGCGTACATGTTCATTGGCCGTGGTAACGCTAAACAGCCTGCGGGCGCTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCACGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGAATGGCGTGTTCACGCTAACTATACTGGGCGTTGGCCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACAcc                                                                                                                                                                                                                           >  2:326352/1‑250 (MQ=255)
                       gggCGTGCTGCTGGCGGACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGTGGTGCGGCAATCATCCACTTACTGGGGGGTATCCACGAAATCTACTTCCCGGATGTGCTGATGAATCCGCGTCTGAGCCTCGCAGGCATCCTCGGCGGTATGACGGTCGTGTGCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTGGCTGTACTGGCGATGACCCCaaa                                                                                                                                                                                                                        <  2:467050/251‑1 (MQ=255)
                                        aCATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCt                                                                                                                                                                                                        <  2:465284/250‑1 (MQ=255)
                                                 ttGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCCCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGGACTGGCGATGACACCAAAAGGGGCTTACTTCGCTAACATCTCgg                                                                                                                                                                                              >  1:374708/1‑251 (MQ=255)
                                                   ggttggtgTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTGCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCTCAGTCATCCTCGGCGTTATGACTGGCGTGTTCACGCTGACTATTCTGGGCGGTGGTCTGGTTTCGCCGGCATCTTCGG‑TTGTATCCTTGCTGTACTGGCGATGACATCAAAAGGTGCTTACTTCGCTAACATTGCggt                                                                                                                                                                                             <  2:360224/243‑2 (MQ=255)
                                                      cGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGgtgt                                                                                                                                                                                          <  1:443489/250‑1 (MQ=255)
                                                           tAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGACATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCACAAGGTGCTTACTTCGCTAACATCGCGGTTGTGTTTgc                                                                                                                                                                                    >  1:296848/1‑251 (MQ=255)
                                                           tAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATACTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAACAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgc                                                                                                                                                                                    >  1:168512/1‑251 (MQ=255)
                                                            aGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATACGCGTCTGATACTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGCTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgcg                                                                                                                                                                                   >  1:336726/1‑251 (MQ=255)
                                                              cgcTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGTCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCCTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgcggc                                                                                                                                                                                 <  1:40526/251‑1 (MQ=255)
                                                               gcttaaCAGTCTGCGGGCGGGGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGTGTGTGTGTgcggcg                                                                                                                                                                                <  2:494685/247‑1 (MQ=255)
                                                                  aaaCAGTCTGCGGGCGGTGCGGCATTCATCCACTTCCTGTGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTAGCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTT‑TATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTCCTTCGCTAACATCGCGGGTGTGTGTGCGGCGAt                                                                                                                                                                              <  1:527506/250‑1 (MQ=255)
                                                                       ggtggggggcgggtcGGCAATCCTTCACTTCCTGGGGGGGTGCCACGAAATCCACCTTCCGGATGGGCTGATGAATCCGCGTCTGGTCCTCGGAGTTATCCTCGGCGGTATGACCGGCGGGTTCCCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTTCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGGGCGGCGATGGCTGt                                                                                                                                                                        <  2:16831/237‑1 (MQ=255)
                                                                        tCTGCGAGCGGTGCGGCGATCATCCACTTCCTGGGGGGTATACACGAAATCTACTTCCCGTATGTGCTGCTGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGCTGAAACCCAAAGGTGCTTACTTCGCTAACATTGCGGGTGTGTGTGCGGCGATGGCTGtt                                                                                                                                                                       >  2:38612/1‑250 (MQ=255)
                                                                         cTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGCATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTCTCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCACAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGtcc                                                                                                                                                                      >  2:227876/1‑250 (MQ=255)
                                                                             gggcggTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGGTGGCTGTCTCCtt                                                                                                                                                                  >  1:349898/1‑251 (MQ=255)
                                                                              ggcggTGCGGCAATCAGCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTc                                                                                                                                                                 <  1:387953/251‑1 (MQ=255)
                                                                                  gTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTg                                                                                                                                                             <  1:192544/251‑1 (MQ=255)
                                                                                       gCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCCGTCTCCTTCGTTGTctct                                                                                                                                                        >  1:106106/1‑251 (MQ=255)
                                                                                       gCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCGCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTGCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTctct                                                                                                                                                        <  2:523344/251‑1 (MQ=255)
                                                                                         gatcatcCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGAATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTAGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACAGCACAAGGGGCTTACTTCGTTTACATCGCGGGTGGGTGTGCGGCGATGGCGGTCTCCCTCGGTTTCTCTGc                                                                                                                                                      >  2:137881/2‑251 (MQ=255)
                                                                                          atcatcCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGCCTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTTCTTCGCTCACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCt                                                                                                                                                     >  1:30441/1‑251 (MQ=255)
                                                                                                    tCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTg                                                                                                                                            <  1:473766/250‑1 (MQ=255)
                                                                                                       tGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCAGCGTCTGATCCTCGCAGTCATCCTCGGCGGTCTGCCTGGCGTGTTCCCGCTGACTATCCTGGGAGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGAAACCAAAAGGTGCTTACTTCGCTCACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGaaa                                                                                                                                         >  1:255902/1‑250 (MQ=255)
                                                                                                         gggggTATCCACGACATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAAcc                                                                                                                                      >  1:423536/1‑251 (MQ=255)
                                                                                                                  cACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTg                                                                                                                             >  1:271058/1‑251 (MQ=255)
                                                                                                                    cGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACGGGCGTGTTCACGCTGACTATCCTGGGCTGTGGTCTTGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACTCCAAAAGGTGCTTACTGCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTTTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGaa                                                                                                                           <  1:489563/251‑1 (MQ=255)
                                                                                                                       aaTCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGaa                                                                                                                           >  1:532836/1‑248 (MQ=255)
                                                                                                                            aCTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACGGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAaagaagaa                                                                                                                    <  2:405215/250‑1 (MQ=255)
                                                                                                                               tCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAgat                                                                                                                 <  1:458895/250‑1 (MQ=255)
                                                                                                                                cccGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAgatga                                                                                                               <  1:482955/251‑1 (MQ=255)
                                                                                                                                    tATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATAtt                                                                                                           >  1:303981/1‑251 (MQ=255)
                                                                                                                                         gCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAag                                                                                                       <  1:479203/250‑1 (MQ=255)
                                                                                                                                          cTGATGAATCCGCGACTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAagc                                                                                                      >  1:90421/1‑250 (MQ=255)
                                                                                                                                               gaATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAAc                                                                                                <  1:149899/251‑1 (MQ=255)
                                                                                                                                                 aTCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCCAAAGGTGCTTACTTCGCTAACATCGCGGGTGGGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACtc                                                                                              >  1:512663/1‑251 (MQ=255)
                                                                                                                                                      cgTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTcgtcgt                                                                                         >  2:416874/1‑251 (MQ=255)
                                                                                                                                                            aTCCTCGCAGTCCTCCTAGGCGGTAAGACTGGCGTGTTCACGCTGACTATCCTGGTAGGGGGTCTGGTTTCTTCGGCATCACCGGCTT‑TATCCTTGCTGTACTGGCGATGACCCCAAAAGGTGCTTACTTCGCTAACCTCGCGGGCGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGc                                                                                     >  2:543423/1‑249 (MQ=255)
                                                                                                                                                            aTCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCCCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCac                                                                                   <  2:303981/251‑2 (MQ=255)
                                                                                                                                                            aTCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATACTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGGGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGc                                                                                     >  1:345450/1‑249 (MQ=255)
                                                                                                                                                               cTCGCAGTCATCCTCGGCGGTACGACTGGCGTGTTCAAGCGGACTTTCCTGGGCGGTGGTCTGGTTTCTCCGGGATCTCCGG‑TTTTATCCTTGCCGTACTGGCGATGACACCAACAGGTGATTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGGCGTATGCAGGAc                                                                                >  2:389110/1‑251 (MQ=255)
                                                                                                                                                                  gCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATg                                                                             <  1:416874/251‑1 (MQ=255)
                                                                                                                                                                     gTCATCCTCGGCGGTATGACTGGCGTGTTCTCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTATCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGaaa                                                                          <  1:336442/251‑1 (MQ=255)
                                                                                                                                                                        aTCCTCGGCGGTATGACTGGCGTGTTCACGCTGCCTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCt                                                                       >  1:100864/1‑251 (MQ=255)
                                                                                                                                                                             cggcggTATGACTGGCGTGTTAACGGTGGCTATCATGGGCGGTGGGCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGGTGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCCGCCAAGTGAAAGAAGAAGATTATATTGAAACAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTc                                                                  >  2:527900/1‑251 (MQ=255)
                                                                                                                                                                                   tATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAgg                                                            >  2:163353/1‑251 (MQ=255)
                                                                                                                                                                                     tGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCCAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGgcg                                                          >  2:250125/1‑251 (MQ=255)
                                                                                                                                                                                      gACTGGAGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATACTTGCTGTACTGGCGATGACACACAAAGGTGCTTACTTCGCTAACATCGCGGGTGAGTGTGCGGCGATGGCTGTCTCCTTCGTAGTCTCTGCTTTTTTGATGAAAACCAGCAACGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGacgc                                                         >  1:549135/1‑251 (MQ=255)
                                                                                                                                                                                         tGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCAtc                                                      >  1:397553/1‑251 (MQ=255)
                                                                                                                                                                                          ggCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCAtct                                                     <  2:168176/251‑1 (MQ=255)
                                                                                                                                                                                                     gCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTctgc                                          >  2:373302/1‑251 (MQ=255)
                                                                                                                                                                                                        gACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTctgct                                         <  2:439741/249‑1 (MQ=255)
                                                                                                                                                                                                        gACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTgg                                       <  2:397553/251‑1 (MQ=255)
                                                                                                                                                                                                            aTCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGAt                                   >  2:88547/1‑251 (MQ=255)
                                                                                                                                                                                                                    cGGTGGTCTGGTTTCTACGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTa                            >  2:252936/1‑250 (MQ=255)
                                                                                                                                                                                                                      gtggtCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACg                         <  2:297654/251‑1 (MQ=255)
                                                                                                                                                                                                                      gtggtCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATATCCGCTGTCTGCTGGCGATGTGACTAACg                         <  1:544567/251‑1 (MQ=255)
                                                                                                                                                                                                                                 tctcCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACAACAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCATTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAACGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCac              >  1:61369/1‑251 (MQ=255)
                                                                                                                                                                                                                                              ccggTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGTACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTAAAatc  <  1:527900/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
GAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTAAAATC  >  NC_000913/3772932‑3773420

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