Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,503,522 A→G 5.8% Y322H (TAT→CAT)  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,5220AG 5.8% 311.3 / 7.4 103Y322H (TAT→CAT) fryAputative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base A (42/55);  new base G (2/4);  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 = 4.17e-01

Only 100 of 106 total aligned reads displayed.

AGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTT  >  NC_000913/2503375‑2503670
                                                                                                                                                   |                                                                                                                                                    
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)
        tgagcgTCAGCCTGCGCCGGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaac                                                                                                                                             <  2:247357/146‑2 (MQ=255)
        tGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTCTCCGCCAGCGTCTGCGCGACCTGATAA‑AACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATc                                                                                                                                            <  2:519013/149‑1 (MQ=255)
        tGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaac                                                                                                                                             <  1:1449303/149‑2 (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)
         gcgcGGCAGCCTGCGCAGGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACAGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGAGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATc                                                                                                                                            <  1:69743/149‑1 (MQ=255)
         gcgcAGCAGCCTGCGCCTGTTGTTTTTGGCGTGTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGGTGCCCGGCATAACGGGCAATc                                                                                                                                            <  1:877435/149‑1 (MQ=255)
           tcgGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACATGATAATAACCTCTAACCGCGTCATTCGGGTTAATCGCAAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGAATAACGTGCAATctc                                                                                                                                          >  2:576656/2‑147 (MQ=255)
           gcgGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGc                                                                                                                                          >  2:570190/1‑149 (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)
                    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: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)
                                                                         aCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCg                                                                            >  1:378951/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                                                                        <  1:553204/149‑1 (MQ=255)
                                                                             cTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGt                                                                        <  1:668527/149‑1 (MQ=255)
                                                                             cTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGt                                                                        <  2:776812/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)
                                                                                             ggTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTcc                                                        >  2:1108980/1‑149 (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)
                                                                                                        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)
                                                                                                                                       gTTGCCAGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGACAGAATCAGTGTATGCGAGGTGCGCCAGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAAc              >  2:837441/1‑149 (MQ=255)
                                                                                                                                       attGCCCGTTATGCCGGGCAATCGCATCCACCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCACGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaaaaa              >  2:588456/2‑148 (MQ=255)
                                                                                                                                          gCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTgg           >  1:1553908/1‑149 (MQ=255)
                                                                                                                                           cccGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAACGCTCAAAAACTGGc          <  2:378951/149‑1 (MQ=255)
                                                                                                                                            ccGGCATAACGGGCAATCGCATCCAGCGGCAAGCCCCTCAGTACCGGGCTCGCCGAGGCACGGGCCAGCATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCCAATCACGGCTCAAAAACTGGCt         >  1:989389/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTT  >  NC_000913/2503375‑2503670

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: ‑