Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,503,523 A→C 6.8% 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 6.8% 249.2 / 5.5 103R321R (CGT→CGG) fryAputative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base A (41/54);  new base C (3/5);  total (44/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 107 total aligned reads displayed.

AGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTA  >  NC_000913/2503375‑2503671
                                                                                                                                                    |                                                                                                                                                    
aGGCCAGCTGCTCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCGTTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATaa                                                                                                                                                      >  2:1347999/1‑149 (MQ=255)
aGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATaa                                                                                                                                                      >  2:360433/1‑149 (MQ=255)
  gCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATACCg                                                                                                                                                    >  1:221557/1‑149 (MQ=255)
  gCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACg                                                                                                                                                    <  2:149544/149‑1 (MQ=255)
   ccAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCCCCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACgg                                                                                                                                                   >  1:796027/1‑149 (MQ=255)
        tGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTCTCCGCCAGCGTCTGCGCGACCTGATAA‑AACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATc                                                                                                                                             <  2:519013/149‑1 (MQ=255)
        tGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaac                                                                                                                                              <  1:509328/149‑2 (MQ=255)
        tGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCaag                                                                                                                                              >  1:718821/1‑148 (MQ=255)
        cgcgcGGCAGCCTGCGCCTGTTGTTGTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGCTAATCGCCAGCACGCCGCACTGGGCGTCTAGCACTGCAGGTTGCCCGTTATGCCGGGCaac                                                                                                                                              <  1:1066686/148‑2 (MQ=255)
         gcgcGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACAGGCAATc                                                                                                                                             <  1:785460/149‑1 (MQ=255)
         gcgcGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGCTTGCCCGGCATAACGGGCAATc                                                                                                                                             >  1:1149816/1‑149 (MQ=255)
         gcgcAGCAGCCTGCGCCTGTTGTTTTTGGCGTGTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGGTGCCCGGCATAACGGGCAATc                                                                                                                                             <  1:877435/149‑1 (MQ=255)
           tcgGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACATGATAATAACCTCTAACCGCGTCATTCGGGTTAATCGCAAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGAATAACGTGCAATctc                                                                                                                                           >  2:576656/2‑147 (MQ=255)
             ggCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCAt                                                                                                                                         >  1:1250175/1‑149 (MQ=255)
             ggCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCAt                                                                                                                                         <  1:544393/149‑1 (MQ=255)
              gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCCTTCGGGTTAATCGCCAGCACGCCGCACTGAGCGTCAAGCACTGCAGGTTGCCCGGCAGAACGGGCAATCGCATc                                                                                                                                        >  1:1391493/1‑149 (MQ=255)
              gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATc                                                                                                                                        <  2:808538/149‑1 (MQ=255)
                  ccTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGc                                                                                                                                    <  1:61256/149‑1 (MQ=255)
                   cTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCg                                                                                                                                   <  2:552123/149‑1 (MQ=255)
                   cTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCg                                                                                                                                   <  2:561106/149‑1 (MQ=255)
                    tGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCgg                                                                                                                                  <  2:81184/149‑1 (MQ=255)
                    tGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCgg                                                                                                                                  <  1:758295/149‑1 (MQ=255)
                       gcCTGTTGTTTTTGGCGTTTATCCGCCAGCGGCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCaa                                                                                                                               <  1:1565683/149‑1 (MQ=255)
                           gttgttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCCTTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCAGAACGGGCATTCGCATCCAGCGGCAAGCCa                                                                                                                           >  2:1014504/1‑149 (MQ=255)
                             tgttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACt                                                                                                                         <  1:11685/149‑1 (MQ=255)
                              gttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTc                                                                                                                        >  2:920887/1‑149 (MQ=255)
                              gttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTc                                                                                                                        <  2:499636/149‑1 (MQ=255)
                               tttttGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGTGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCATTCGCATCCAGCGGCAAGACCCTcc                                                                                                                       >  1:863501/1‑148 (MQ=255)
                               tttttGGCGTTTATCCGCCAGCGTCTGCGCGACATGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCCCACTGGGCGTCAAGAACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACtaa                                                                                                                       >  2:1186847/1‑147 (MQ=255)
                                 tttGGCGTTTAGCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAAGCGCATCCAGCGGCAAGCCACTCAGt                                                                                                                     <  2:530280/149‑1 (MQ=255)
                                  ttGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTa                                                                                                                    <  1:310843/149‑1 (MQ=255)
                                  ttGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTa                                                                                                                    <  1:1146283/149‑1 (MQ=255)
                                    ggCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTAcc                                                                                                                  >  2:1475496/1‑149 (MQ=255)
                                    ggCGTTAATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTAcc                                                                                                                  <  1:292065/149‑1 (MQ=255)
                                     gCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCg                                                                                                                 <  2:354465/149‑1 (MQ=255)
                                     gCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCg                                                                                                                 <  1:1068919/149‑1 (MQ=255)
                                       gTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCggg                                                                                                               >  2:579761/1‑149 (MQ=255)
                                         ttATCCGCCAGCGTCTGCGCGACCTGATAATACCCGCTCCCCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGAt                                                                                                             <  2:856647/149‑1 (MQ=255)
                                              ctccACCGTCTGCTCGACCTGATAATAACCGCTCACCGAGTCATTCGGGTTAATAGCAAGCACGCTGCACTGGTCGACAAGCACTTCACGTTGCACGCCATAACGTGAAATAGCATCCAGCCTCAAGCCACACAGTACCGGGATCGCCg                                                                                                        >  2:687065/3‑149 (MQ=255)
                                              cGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGACCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCg                                                                                                        >  2:1350803/1‑149 (MQ=255)
                                               gCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGa                                                                                                       <  1:516050/149‑1 (MQ=255)
                                               gCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGa                                                                                                       <  1:125559/149‑1 (MQ=255)
                                               gCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGa                                                                                                       >  2:31343/1‑149 (MQ=255)
                                                ccaggttCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGTGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCagag                                                                                                      <  2:126466/143‑1 (MQ=255)
                                                       cTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACggg                                                                                               <  1:1475496/149‑1 (MQ=255)
                                                         gcgcGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGcc                                                                                             <  1:1332147/149‑1 (MQ=255)
                                                          cgcgACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCa                                                                                            >  1:524360/1‑149 (MQ=255)
                                                           gcgACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCAAGCGGCAAGCCACTCAGAACCGGGACCGCCGAGGCACGGGCCAg                                                                                           >  1:1499528/1‑149 (MQ=255)
                                                              cccTGATAAAAACCGCTCACCGCCGCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGCTTGCCCGGCATACCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTCCCGGGATCGCCGAGGCACGGGCCAGAAt                                                                                        <  2:963132/148‑1 (MQ=255)
                                                                  gATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgt                                                                                    >  2:1503499/1‑149 (MQ=255)
                                                                  gATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTCGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgt                                                                                    <  2:390690/149‑1 (MQ=255)
                                                                   ataataACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgtg                                                                                   >  1:425498/1‑149 (MQ=255)
                                                                      ataaCCGCTCACCGAGTCATTCGAGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGATATCCCGGCATAACGGGCAATCGCATCCAGCGGAATGCCACTAAGTACCGGGATAGCCGAGGCACGCGCCAGAATCAGAGTAt                                                                                >  2:109769/1‑149 (MQ=255)
                                                                           cGCTCACCGCGGCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAg                                                                           <  1:301403/149‑1 (MQ=255)
                                                                            gCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAgg                                                                          <  1:920887/149‑1 (MQ=255)
                                                                            gCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGAGCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAgg                                                                          <  2:164485/149‑1 (MQ=255)
                                                                             cTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGt                                                                         <  2:776812/149‑1 (MQ=255)
                                                                             cTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGt                                                                         <  1:553204/149‑1 (MQ=255)
                                                                             cTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGt                                                                         <  1:668527/149‑1 (MQ=255)
                                                                                       cATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAAGGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGtt                                                               <  1:1347999/149‑1 (MQ=255)
                                                                                        aTTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCAGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCAAGAATCAGTGGATGCGCGGTGCGCCCGGttg                                                              >  2:465851/1‑148 (MQ=255)
                                                                                         ctCGGGTCACTCGCCAGCACGCCGCACTGGGCGGCAAGCACTGCAGGTTGCCCGGCATAACCGGCAATCGCATCCAGCGGCAAGCCACTCAGTCCCGGGAGCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGtttt                                                             <  1:207446/148‑1 (MQ=255)
                                                                                          ccGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTtc                                                            <  2:728999/148‑1 (MQ=255)
                                                                                             ggTTCATCGCCAGCACGCCGCACTGGGCGTCAAGCACGGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGAACCGGGATCGCCGAGGCACGGGCCAGAATCAGGGTATGCGAGGTGCGCCCGGGTTTCTcc                                                         <  1:227043/149‑1 (MQ=255)
                                                                                               ttACTCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAg                                                       <  1:1014504/149‑1 (MQ=255)
                                                                                               ttAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGGTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAg                                                       <  1:1493759/149‑1 (MQ=255)
                                                                                                tAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGa                                                      <  1:1119267/149‑1 (MQ=255)
                                                                                                tAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTATTCTCCAGa                                                      >  2:469567/1‑149 (MQ=255)
                                                                                                   tCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATc                                                   <  2:341491/149‑1 (MQ=255)
                                                                                                   tCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATc                                                   <  1:23311/149‑1 (MQ=255)
                                                                                                    cGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCa                                                  >  2:1106537/1‑149 (MQ=255)
                                                                                                     gCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCAt                                                 >  1:1224460/1‑149 (MQ=255)
                                                                                                      ccAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCCCTCAGTACCGGGATCGCCGAGGCACGGGCCAGATTCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATg                                                >  1:1352567/1‑149 (MQ=255)
                                                                                                       cAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGc                                               <  1:516686/149‑1 (MQ=255)
                                                                                                       cAGCACGCCGCACTGCGCGTCAAGCACTGCAGGTGGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGc                                               <  2:1391493/149‑1 (MQ=255)
                                                                                                        aGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGcc                                              <  2:1250175/149‑1 (MQ=255)
                                                                                                         gCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGccc                                             >  1:444243/1‑149 (MQ=255)
                                                                                                                  aCTTGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGAttt                                    <  1:1431813/149‑1 (MQ=255)
                                                                                                                  aCTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGAttt                                    >  2:12302/1‑149 (MQ=255)
                                                                                                                  aCTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGAttt                                    >  1:770096/1‑149 (MQ=255)
                                                                                                                     gggCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTc                                 <  1:826272/149‑1 (MQ=255)
                                                                                                                       gCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCaa                               <  2:56701/149‑1 (MQ=255)
                                                                                                                         gTCAAGCACTGCAGGTTGCCCGGCTTAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTGTTCTCCAGAATCATGCCCGCAAGATTTTTCAAAt                             >  1:692490/1‑149 (MQ=255)
                                                                                                                            aaGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCaa                          >  2:1139434/1‑149 (MQ=255)
                                                                                                                             aGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCTGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCacg                         >  1:248121/1‑147 (MQ=255)
                                                                                                                                aCTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCGGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCt                      >  2:919785/1‑149 (MQ=255)
                                                                                                                                   gCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaa                   >  2:964738/1‑149 (MQ=255)
                                                                                                                                    cAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACGGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGGGTTCTCCAGAATCCTGCCCGCCAGCTTTTTCATATCAAGGCTCcaa                  >  1:394342/1‑149 (MQ=255)
                                                                                                                                     aGGTTGCCCGGCAAAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTACAGAATCATGCCCGCAAGATTTTTCATATCAAGTCTCtaaa                 <  2:1214964/149‑1 (MQ=255)
                                                                                                                                      ggTTGGCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaaaa                <  1:570190/149‑1 (MQ=255)
                                                                                                                                      ggTTGCCCGGCATAACGGGCAATCGCATCCATCGGCAAGCCACTCAGTAGCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGACCGCAAGATTTTTCAAATCAAGGCTCaaaaa                <  1:576656/149‑1 (MQ=255)
                                                                                                                                       gTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCCCCCCGTACCGGGATCGCCGAGGCACGTGCCAGAATCAGTGTATGCGAGGTGCGCACGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATAAAGGCTCAAAAAc               <  2:1064733/149‑1 (MQ=255)
                                                                                                                                       gTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAAc               >  1:1491996/1‑149 (MQ=255)
                                                                                                                                       gTTGCCAGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGACAGAATCAGTGTATGCGAGGTGCGCCAGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAAc               >  2:837441/1‑149 (MQ=255)
                                                                                                                                       attGCCCGTTATGCCGGGCAATCGCATCCACCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCACGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaaaaa               >  2:588456/2‑148 (MQ=255)
                                                                                                                                        ttGCCCGTTATGCCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCCAGGCTCAAAAACt              >  1:836142/1‑149 (MQ=255)
                                                                                                                                          gCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTgg            >  1:1553908/1‑149 (MQ=255)
                                                                                                                                           cccGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAACGCTCAAAAACTGGc           <  2:378951/149‑1 (MQ=255)
                                                                                                                                            ccGGCATAACGGGCAATCGCATCCAGCGGCAAGCCCCTCAGTACCGGGCTCGCCGAGGCACGGGCCAGCATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCCAATCACGGCTCAAAAACTGGCt          >  1:989389/1‑149 (MQ=255)
                                                                                                                                                    aCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGACTCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTa  >  2:1110749/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTA  >  NC_000913/2503375‑2503671

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