Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,503,523 A→C 5.9% R321R (CGT→CGG)  fryA ← putative PTS enzyme: Hpr, enzyme I and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,503,5230AC 5.9% 338.7 / 8.9 124R321R (CGT→CGG) fryAputative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base A (55/61);  new base C (4/4);  total (59/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

Only 100 of 125 total aligned reads displayed.

GGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTA  >  NC_000913/2503376‑2503671
                                                                                                                                                   |                                                                                                                                                    
ggCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACg                                                                                                                                                    >  1:1083208/1‑150 (MQ=255)
 gCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACgg                                                                                                                                                   <  1:1698707/150‑1 (MQ=255)
   cAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGc                                                                                                                                                 <  2:1151132/150‑1 (MQ=255)
    aGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCa                                                                                                                                                <  1:808506/150‑1 (MQ=255)
       tGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGATTGCCCGTTATGCCGGGCaacc                                                                                                                                             <  1:262598/150‑3 (MQ=255)
        gcgcGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaacct                                                                                                                                            <  2:2101249/150‑4 (MQ=255)
         cgcgGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGc                                                                                                                                           <  2:1715198/150‑1 (MQ=255)
             gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATACCGGGCAATCGCATcc                                                                                                                                       >  1:1640690/1‑150 (MQ=255)
             gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATcc                                                                                                                                       >  1:456391/1‑150 (MQ=255)
             gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATcc                                                                                                                                       <  1:1012319/150‑1 (MQ=255)
             gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATcc                                                                                                                                       <  2:454330/150‑1 (MQ=255)
               aGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAg                                                                                                                                     >  1:1739546/1‑150 (MQ=255)
                 ccTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCg                                                                                                                                   >  1:1142390/1‑150 (MQ=255)
                 ccTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCAAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCg                                                                                                                                   >  2:2065655/1‑150 (MQ=255)
                  cTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCGGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCgg                                                                                                                                  <  1:1344284/150‑1 (MQ=255)
                   ggctccTGTTGTTTTTGGCGGTTATCCGCCAGCGTCTGCGCGACCTGATACTAACCGCGCACCGCGTCATTCGGGTTAATCGCCCGCACGCCGCACTGGGCGGCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGc                                                                                                                                 <  1:1048639/146‑1 (MQ=255)
                    gcgcCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCa                                                                                                                                >  2:1590011/1‑150 (MQ=255)
                    gcgcCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCa                                                                                                                                >  2:1742404/1‑150 (MQ=255)
                      gcCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAg                                                                                                                              >  1:1477409/1‑150 (MQ=255)
                         tgttgtTATTGCAGATTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCGCACCGCGTCAATCGGGTTAATCGCCAGCACCCCGCACTGGGCGTCAAGCCCTGCAGGTTGCCCGCCATAACGGGCAATCGCATCCAGCGGCAAGACa                                                                                                                           >  2:1292453/1‑150 (MQ=255)
                           ttgttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACt                                                                                                                         <  1:1065847/150‑1 (MQ=255)
                             gttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTcc                                                                                                                       >  2:967367/1‑149 (MQ=255)
                             gttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCa                                                                                                                       >  2:677706/1‑150 (MQ=255)
                               ttttGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGt                                                                                                                     >  1:53730/1‑150 (MQ=255)
                               accagccGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGt                                                                                                                     <  2:2136095/144‑1 (MQ=255)
                                tttGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTa                                                                                                                    >  1:568727/1‑150 (MQ=255)
                                   ggCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCg                                                                                                                 >  2:834088/1‑150 (MQ=255)
                                        ttATCCGCCTGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATc                                                                                                            <  2:1941063/150‑1 (MQ=255)
                                           tCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGcc                                                                                                         >  2:726588/1‑150 (MQ=255)
                                             cGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGa                                                                                                       >  1:987435/1‑150 (MQ=255)
                                             cGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGACAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATACAGCGGCAAGCCACTCAGTACCGGGATCGCCgc                                                                                                       >  2:1126037/1‑149 (MQ=255)
                                              gCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAg                                                                                                      <  2:539285/150‑1 (MQ=255)
                                                cAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGc                                                                                                    <  2:158485/150‑1 (MQ=255)
                                                 aGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGTGTCAAGCACTGCAGGTTGCCAGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGcc                                                                                                   >  2:1799228/1‑149 (MQ=255)
                                                 aGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCa                                                                                                   >  2:824373/1‑150 (MQ=255)
                                                   cGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACg                                                                                                 >  2:891648/1‑150 (MQ=255)
                                                   cGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACg                                                                                                 <  2:703622/150‑1 (MQ=255)
                                                    gTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACgg                                                                                                <  2:836221/150‑1 (MQ=255)
                                                       tGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGcc                                                                                             <  2:1871400/150‑1 (MQ=255)
                                                        gcgcGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACGGGGCGTCAAGCACTGCAGGTGGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCa                                                                                            <  1:1185477/150‑1 (MQ=255)
                                                          gcgACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGa                                                                                          <  2:1453133/150‑1 (MQ=255)
                                                            gACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAAt                                                                                        <  1:1835904/150‑1 (MQ=255)
                                                              ccTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCa                                                                                      >  1:1606439/1‑150 (MQ=255)
                                                               cTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAg                                                                                     >  2:254645/1‑150 (MQ=255)
                                                               cTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAg                                                                                     <  2:241702/150‑1 (MQ=255)
                                                               cTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAg                                                                                     <  1:615188/150‑1 (MQ=255)
                                                                tGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCCCGGGCCAGCATCTgt                                                                                    >  2:348564/1‑150 (MQ=255)
                                                                tGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgt                                                                                    >  1:1250480/1‑150 (MQ=255)
                                                                    aataaCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGGCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTAt                                                                                <  1:185927/150‑1 (MQ=255)
                                                                     ataaCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGCTCGCCGCGGCCCGGGCCAGAATCAGTGTATg                                                                               >  2:233515/1‑150 (MQ=255)
                                                                     ataaCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATg                                                                               <  1:1904306/150‑1 (MQ=255)
                                                                      taaCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCACTCGCCTCCAGCGGCACGCCACTCAGTACCGGGATCGCCGAGGCGCGGGCCAGAATCAGTGTATGc                                                                              >  2:1144419/1‑150 (MQ=255)
                                                                      taaCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGc                                                                              >  2:959359/1‑150 (MQ=255)
                                                                       aaCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCCGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGCATCAGTCTATGCg                                                                             >  1:1386810/1‑150 (MQ=255)
                                                                        aCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGa                                                                            >  1:1454624/1‑150 (MQ=255)
                                                                            cTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCACTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTg                                                                        >  1:2165152/1‑150 (MQ=255)
                                                                               aCCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGACAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTgcgc                                                                     >  2:1470227/1‑150 (MQ=255)
                                                                                ccGCGTCAATCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCCTCCAGCGGCAAGCCACTCCGTACCTGGATCGCAGCGGCAAGGGGCAGAATCAGAGTATCCGCGGTGCGcc                                                                    >  1:1621821/1‑150 (MQ=255)
                                                                                  gcgTCATTCGGGTTAATCGCCAGCACGCCGGACTGTGCGGCAAGCACTGCAGTATGCCCGTCATAACGGGCAATCGCATCAAGCGGCAAGGCACTCAGTAGCGGGATCGCCGAGGCACGTGCCAGAATCAGTGTAAGCGAGGTGCGCCCg                                                                  <  2:332669/150‑1 (MQ=255)
                                                                                   cgTCATTCGGCGCACTCGCCAGCCCGCCGCACTGGGCGGCAAGCACAGCAGGTTGCCCGGCATAACGGGCAAGCGCACCCCGCGGCAAGCCACTCCGTACCGGGCGCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCgg                                                                 <  1:2027350/150‑1 (MQ=255)
                                                                                    gTCATTCGGGTTAATCGCCATCACGACGCACTGGGCGTCAAGCACTGAATGTTGCACGGCATAAAGGGCAATCGCATCCAGCGTCAAGACACTCATTACATCGATCGCCGCGGCACGGGCCAGAATCAGTTTATGCGATGTTCTACCGGt                                                                >  2:2151592/1‑150 (MQ=255)
                                                                                    gTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGt                                                                <  2:1739546/150‑1 (MQ=255)
                                                                                      cATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGttt                                                              <  2:839395/150‑1 (MQ=255)
                                                                                      cATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGGCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGttt                                                              <  1:601332/150‑1 (MQ=255)
                                                                                        ttCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATGACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTtc                                                            <  2:538291/150‑1 (MQ=255)
                                                                                             gctAATCGCCAGCCCGCCGCACTGGCCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAg                                                       <  1:292359/148‑1 (MQ=255)
                                                                                                aaTCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCCCGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGGTTTCTCCAGAAt                                                    >  2:1667466/1‑150 (MQ=255)
                                                                                                aaTCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAAt                                                    >  2:652611/1‑150 (MQ=255)
                                                                                                aaTCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAAt                                                    <  2:1420133/150‑1 (MQ=255)
                                                                                                 aTCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATc                                                   >  2:1197799/1‑150 (MQ=255)
                                                                                                  tCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCa                                                  <  2:2154242/150‑1 (MQ=255)
                                                                                                    gCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATg                                                >  2:1868178/1‑150 (MQ=255)
                                                                                                      cAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGcc                                              <  1:456816/150‑1 (MQ=255)
                                                                                                        gCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCg                                            <  2:762555/150‑1 (MQ=255)
                                                                                                                cACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGAttt                                    >  2:1234691/1‑150 (MQ=255)
                                                                                                                   tGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTc                                 <  1:1559885/150‑1 (MQ=255)
                                                                                                                    gggCGTCAAGCACTGCAGGTTGCCCGGCATACCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCa                                <  2:655526/150‑1 (MQ=255)
                                                                                                                      gCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCaaa                              <  2:368083/150‑1 (MQ=255)
                                                                                                                      gCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCaaa                              <  1:891648/150‑1 (MQ=255)
                                                                                                                         tCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCACCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCa                           <  2:2003775/150‑1 (MQ=255)
                                                                                                                          cAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCaa                          >  1:176477/1‑150 (MQ=255)
                                                                                                                             gCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGc                       <  2:182817/150‑1 (MQ=255)
                                                                                                                             gCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGc                       <  1:2087143/150‑1 (MQ=255)
                                                                                                                                cTGCAGGTTGCCCGGCATAACGGGCAATCGCGTCCAGCGGCAAGCCACTCAGTACCGGTATCGCCGAGGCACGTGCCAGAATCAGTGTATGCGAGGTGCGCCCGTTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCa                    <  2:290259/150‑1 (MQ=255)
                                                                                                                                  gCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaa                  >  2:1759198/1‑150 (MQ=255)
                                                                                                                                   cAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaaa                 <  2:616733/150‑1 (MQ=255)
                                                                                                                                    cgattGCCCGTTATGCCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaaaa                >  2:652368/4‑150 (MQ=255)
                                                                                                                                    cgattGCCCGTTATGCCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGACCGCAAGATTTTTCAAATCAAGGCTCaaaaa                >  2:1512326/4‑150 (MQ=255)
                                                                                                                                       ttGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTg             <  1:1700802/150‑1 (MQ=255)
                                                                                                                                       ttGCCCGGCATAACGGGAAATCGAATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTg             <  2:2056401/150‑1 (MQ=255)
                                                                                                                                         gCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCTAATCAAGGCTCAAAAACTGGc           >  1:634504/1‑150 (MQ=255)
                                                                                                                                             ggCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGACTCAAAAACTGGCTTgg       >  2:545582/1‑150 (MQ=255)
                                                                                                                                              gCATAACGGGCAATCGCTTCCAGGGGCAATCCACTCAGTACTGGGATAGCCGAGGCACGGGCCATAATCAGTGTATGCGAGGTTCGCCCGGTTTTCTCCATAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAACACTGGCTTggg      <  2:126934/150‑1 (MQ=255)
                                                                                                                                               cATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTgggg     <  1:824373/150‑1 (MQ=255)
                                                                                                                                               cATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTgggg     <  1:1217734/150‑1 (MQ=255)
                                                                                                                                               cATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAAGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTgggg     <  2:117970/150‑1 (MQ=255)
                                                                                                                                                aTAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGt    <  2:1454624/150‑1 (MQ=255)
                                                                                                                                                 tAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGtt   >  2:987687/1‑150 (MQ=255)
                                                                                                                                                 tAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTGTCAAATCAAGGCTCAAAACCGGGCTTgggggt   >  2:75600/1‑148 (MQ=255)
                                                                                                                                                  aaCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTa  >  1:2174497/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTA  >  NC_000913/2503376‑2503671

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