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,1760CT97.0% 241.1 / ‑0.9 68coding (896/1914 nt)mtlAmannitol‑specific PTS enzyme IICBA component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/32);  minor base . (1/1);  total (34/34)
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.92e-01

GCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTA  >  NC_000913/3772935‑3773415
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
gCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCt                                                                                                                                                                                                                                       <  2:420266/251‑1 (MQ=255)
    aCCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTAc                                                                                                                                                                                                                                   <  1:266907/251‑1 (MQ=255)
    aCCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTAc                                                                                                                                                                                                                                   >  2:423989/1‑251 (MQ=255)
      ccaggtccaggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTg                                                                                                                                                                                                                                 >  1:176457/1‑251 (MQ=255)
        aggtccaggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGTTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTTTACTGGc                                                                                                                                                                                                                               <  1:503597/251‑1 (MQ=255)
        aggtccaggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGc                                                                                                                                                                                                                               >  2:360946/1‑251 (MQ=255)
        aggtccaggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGGGGTCTGGGTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGc                                                                                                                                                                                                                               >  2:170012/1‑251 (MQ=255)
         ggtccaggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTGTATCCTTGCTGTACTGGCg                                                                                                                                                                                                                              <  2:75883/251‑1 (MQ=255)
            ccaggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGG‑TATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATg                                                                                                                                                                                                                           <  1:308925/250‑1 (MQ=255)
                gtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGacac                                                                                                                                                                                                                       <  1:162109/251‑1 (MQ=255)
                    gggCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGGCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGAttccaccca                                                                                                                                                                                                                    >  2:436096/1‑242 (MQ=255)
                         tgctgcTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTg                                                                                                                                                                                                              >  1:38009/1‑251 (MQ=255)
                                 gCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATATGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTGCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTc                                                                                                                                                                                                      <  2:335096/251‑1 (MQ=255)
                                               tGGTCGTGGTAGCGCTAAACAGTCTGCGGTCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCgg                                                                                                                                                                                         <  2:550300/250‑1 (MQ=255)
                                                     tGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGgtgtgtg                                                                                                                                                                                  >  2:337201/1‑251 (MQ=255)
                                                      ggTAGCGCCAAACAGCCCGCGGGCGGTGCCGCAATCATCCCCTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGgtgtgtg                                                                                                                                                                                  <  2:574432/250‑1 (MQ=255)
                                                        tAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGGCTGGTTTCTCCGGCCTCTCCGGGTT‑TATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgc                                                                                                                                                                               >  2:132446/1‑251 (MQ=255)
                                                           cgcgaaccAGTCTGCGGGCGGTGCGGCAATCATCCACGTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCGCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgcggc                                                                                                                                                                            <  2:438210/244‑1 (MQ=255)
                                                           cgcgaaaCAGTCTGCGGGCGGTGCGGCAATCATCCGCTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTTGGCGGTGGTCTGGTTTCTCCGGCATCTGCGGGTT‑TATCCTTGCTGTACTGGCGGTGACACGAAAATGTGCTTACTTCGCTAACATCGCGGTTGTGTGTgcggc                                                                                                                                                                            <  1:92508/247‑1 (MQ=255)
                                                           cgcTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgcggc                                                                                                                                                                            >  1:224628/1‑251 (MQ=255)
                                                            gcTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGGGTGTGTgcggcg                                                                                                                                                                           >  1:158675/1‑251 (MQ=255)
                                                               aaaCAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGAt                                                                                                                                                                         <  2:200816/250‑1 (MQ=255)
                                                                        gcgggcgtTGCGGCAATCATCCACTGCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTcc                                                                                                                                                               <  2:223451/251‑1 (MQ=255)
                                                                           ggcggTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTc                                                                                                                                                            <  2:307560/251‑1 (MQ=255)
                                                                             cggTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGt                                                                                                                                                          <  1:433527/251‑1 (MQ=255)
                                                                                    gCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTctct                                                                                                                                                   <  2:143848/251‑1 (MQ=255)
                                                                                      aaTCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGc                                                                                                                                                 >  1:162398/1‑251 (MQ=255)
                                                                                         catcCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTAt                                                                                                                                              <  1:139940/251‑1 (MQ=255)
                                                                                         catcCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACGTCGCTAACATCGCGGGTGGGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTAt                                                                                                                                              <  1:58341/251‑1 (MQ=255)
                                                                                             cACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGCATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGGGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTAttttt                                                                                                                                          >  2:112360/1‑250 (MQ=255)
                                                                                             cACTTCCTGGGGGGTATCCACGAAATCTACTGCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTGCACGCTGACGATCCTGGGCGGGGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGAGGGCTGTCTCCTTCGTTGTCTCTGCTATTTTg                                                                                                                                          <  2:334161/251‑1 (MQ=255)
                                                                                                 tCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTg                                                                                                                                       <  1:239999/250‑1 (MQ=255)
                                                                                                         ggTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGc                                                                                                                              >  1:574356/1‑251 (MQ=255)
                                                                                                               cACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTg                                                                                                                        <  1:440848/251‑1 (MQ=255)
                                                                                                                   aaaTCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGaa                                                                                                                      >  1:321940/1‑249 (MQ=255)
                                                                                                                              ccGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAATCCAGCAAAGTGAAAGAAGAAgatgat                                                                                                         <  1:153846/251‑1 (MQ=255)
                                                                                                                                gTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGatat                                                                                                       <  1:170012/251‑1 (MQ=255)
                                                                                                                                  aTGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATAttt                                                                                                     >  2:171046/1‑250 (MQ=255)
                                                                                                                                  aTGTGCTGATGAATCAGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGCCTATCCTGGGCGGTGGTCTGGTTTCTCCGGCCTCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCACAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGCAAGAAGAAGATGATAttt                                                                                                     >  2:568454/1‑250 (MQ=255)
                                                                                                                                   tgtgCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTg                                                                                                     >  1:540511/1‑250 (MQ=255)
                                                                                                                                           tgaATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAgcagca                                                                                             >  2:459290/1‑250 (MQ=255)
                                                                                                                                            gaATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACGGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTGGTCTCTGCTATGTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAAc                                                                                           <  2:353788/251‑1 (MQ=255)
                                                                                                                                             aaTCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACt                                                                                          <  1:553466/251‑1 (MQ=255)
                                                                                                                                                ccGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACtcg                                                                                        <  1:423989/250‑1 (MQ=255)
                                                                                                                                                       tGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGAGGGCTGTCTCCTTCGTTGTCTCTGCCTTTTTTCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTAttc                                                                                >  2:391840/1‑249 (MQ=255)
                                                                                                                                                        gATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGc                                                                                >  1:144263/1‑250 (MQ=255)
                                                                                                                                                         aTCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTTCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGGGTGGGTGCGGCGATGGCTTTCTCCTTCGTTGTCTCCTCTATTTTGCTTAAAACCAGCAAAGTGAAAGAAGAAGATGATAATTAAGCAGCAAATCGTCGGATGc                                                                                >  2:471294/1‑249 (MQ=255)
                                                                                                                                                           ccTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAgg                                                                             >  1:148954/1‑250 (MQ=255)
                                                                                                                                                                 aGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGaa                                                                      <  2:176457/251‑1 (MQ=255)
                                                                                                                                                                 aGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGa                                                                       >  1:358014/1‑250 (MQ=255)
                                                                                                                                                                    cATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAg                                                                    <  1:30693/250‑1 (MQ=255)
                                                                                                                                                                     aTCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCt                                                                  >  1:237995/1‑251 (MQ=255)
                                                                                                                                                                     aTCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGGCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGGGAAAGAAGAAGATGATATTGAAGCAGCAACTGGTCGTATGCGGGACAGGAAAGCt                                                                  <  2:92086/251‑1 (MQ=255)
                                                                                                                                                                      tCCTCGGCGGTATGACTGGCGTGTTCAAGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGGGTGTGTGCGGCGATGGCTGTCTCCTTCGTTTTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGGTGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTg                                                                 >  2:156882/1‑251 (MQ=255)
                                                                                                                                                                        cTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTTTGGTTTCTCCGGCATCTTCGG‑TTTTTTCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGGGTGTGGGCGGCGATGGCTGTCTTCTTCGTTGTCTCTGCTATTTTGCTGAAAAACAGCAAAGTGAAAGCAGAAGATGATATTGAAGCCGGAACTCGTTGTATGCAGGACATGAAAGCTg                                                                 >  2:57595/1‑249 (MQ=255)
                                                                                                                                                                               gTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTaa                                                          >  1:239524/1‑249 (MQ=255)
                                                                                                                                                                                  tGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCCGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGgcg                                                     >  1:359212/1‑251 (MQ=255)
                                                                                                                                                                                        gCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAACGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCAtttc                                               >  2:507218/1‑248 (MQ=255)
                                                                                                                                                                                                aCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTct                                       <  1:217363/251‑1 (MQ=255)
                                                                                                                                                                                                         aTCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGAt                              >  2:175752/1‑251 (MQ=255)
                                                                                                                                                                                                           ccTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGAtg                             <  2:319198/250‑1 (MQ=255)
                                                                                                                                                                                                                gCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACt                        >  2:166993/1‑250 (MQ=255)
                                                                                                                                                                                                                    tggtCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACg                    >  1:459986/1‑250 (MQ=255)
                                                                                                                                                                                                                     ggtCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGa                   <  1:41046/250‑1 (MQ=255)
                                                                                                                                                                                                                          ggTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGGACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGa              <  2:415508/250‑1 (MQ=255)
                                                                                                                                                                                                                                 ccGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCacg        <  1:475232/249‑1 (MQ=255)
                                                                                                                                                                                                                                     cATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCacgtacgt   >  1:211417/1‑250 (MQ=255)
                                                                                                                                                                                                                                      aTCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCAcgtacgta  >  1:163684/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
GCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTA  >  NC_000913/3772935‑3773415

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