Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,826,6290TC26.9% 59.0 / 31.3 53Y80H (TAT→CAT) srlEglucitol/sorbitol‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base T (20/18);  new base C (5/9);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-02
Rejected as consensus: Frequency below/above cutoff threshold.

TACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAAAGAAGAAAACATCACTGTAGTAGGTGATGCG  >  NC_000913/2826490‑2826760
                                                                                                                                           |                                                                                                                                   
tACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCGGGCTATTGACGGAGTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCAtctctc                                                                                                                                  >  2:76206/1‑143 (MQ=255)
    cGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGACGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGaa                                                                                                                              >  1:839545/1‑143 (MQ=255)
     ggCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGaaa                                                                                                                             <  1:325738/143‑1 (MQ=255)
        cTGCGATTGTTGACAAACTGGCACAGCGTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCGATCCGAAacga                                                                                                                          <  1:346381/143‑1 (MQ=255)
              ttgttgACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTAtt                                                                                                                    >  1:520905/1‑143 (MQ=255)
              ttgttgACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTAtt                                                                                                                    >  1:520492/1‑143 (MQ=255)
              ttgttgACAAACTGGCACAGCGTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTAtt                                                                                                                    <  2:541127/143‑1 (MQ=255)
                    aCAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAAGTGGTGTCGCGGTAATCGACTGTGGAGACACATTCCGCTGCGGCATCTATCCGACACGACGTATGCccacc                                                                                                              >  2:623861/1‑143 (MQ=255)
                     cagaCTGGCGCAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCCATCCGAAACGACGTATTCccacca                                                                                                             <  2:887060/140‑1 (MQ=255)
                        aCTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCCATCCGAAACGACGTATTCCCACCATTa                                                                                                          >  1:736434/1‑143 (MQ=255)
                            gCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCAGTAatat                                                                                                      >  1:412615/1‑143 (MQ=255)
                             cccaGCTTCCTGGCGGGCAGGCAATAGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGGGGCGGCACATTACGCTGCGGCATCCATCCGAAACGACGTATTCCCACCATTAATATc                                                                                                     <  2:361276/141‑1 (MQ=255)
                              acaGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATcc                                                                                                    >  1:749941/1‑143 (MQ=255)
                              acaGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCGGCGCCATCTATCCGAAACGCCGTAGTCCCACACTTCATATcc                                                                                                    >  2:77722/1‑143 (MQ=255)
                                       tggctggcAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACgg                                                                                           >  2:26102/1‑143 (MQ=255)
                                         gctggcAGGCTGTTGTCGGAGTCAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGCGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGc                                                                                         <  2:485796/143‑1 (MQ=255)
                                           gcgcAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCCATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCaa                                                                                       <  1:90525/141‑1 (MQ=255)
                                                   tATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAAAGACGCAGTGTCACCATTAATATCCACTCGACGGGCAAgtccggtc                                                                               >  2:315715/1‑143 (MQ=255)
                                                     ttGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCCATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAgtccggcccc                                                                             >  2:517988/1‑142 (MQ=255)
                                                     ttGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGATGTCGCGGTAATCGACTGTGGCGGCACAATACGCCGCCGCATCTATCAGAAAACACGTATTCCACCCCTCAATCACCCCCCCACGCGCAAgtcccgcccg                                                                             >  2:746496/1‑143 (MQ=255)
                                                       gACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCt                                                                           >  1:854297/1‑143 (MQ=255)
                                                             tttAACGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCa                                                                     <  1:322846/143‑1 (MQ=255)
                                                                aaaGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTa                                                                  <  1:651821/143‑1 (MQ=255)
                                                                 aagaagGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCAGCTATCCGAAACGACGTATTCCCACCATTAATCTCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTAc                                                                 >  2:202097/1‑143 (MQ=255)
                                                                      ggTGGACCCGGGGAGGCGGAACGTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATtgt                                                            <  2:412615/143‑1 (MQ=255)
                                                                      ggTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCCATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATtgt                                                            >  1:798503/1‑143 (MQ=255)
                                                                       gttacgccGCGGGGGCGGAAATTGGTGTCGCGGTAATCGACTGTGCCGGCCCATGACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATtgtg                                                           <  2:546098/137‑1 (MQ=255)
                                                                         gAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGa                                                         >  2:419437/1‑143 (MQ=255)
                                                                          agcccGCGGAGGCGGAAATGGGTGTCGCGGTAAGCGACTGTGGCGGCACATTACGCTGCGGCAGCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGaa                                                        <  2:350763/141‑1 (MQ=255)
                                                                          aaCTCGCGGAGGCGGAAATTGGTGTCGCGGTAATGGACTGTGGCGGGACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGaa                                                        <  1:41034/143‑1 (MQ=255)
                                                                          aaCCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCCATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGaa                                                        <  1:816010/143‑1 (MQ=255)
                                                                            cccGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGa                                                      <  1:483421/143‑1 (MQ=255)
                                                                             ccGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGat                                                     >  1:14559/1‑143 (MQ=255)
                                                                                   aGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCCATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTAt                                               <  1:778840/143‑1 (MQ=255)
                                                                                       ggAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACAATTAGACACTACTCGACGGGCTAGTCCGGTCCGCTGGCGCCGTCCGCTGTGGAAGATATGTATGtcg                                           >  1:907115/1‑142 (MQ=255)
                                                                                       ggAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCAGCTATCCGAAACGACGTATTCCCACCATGAATATCCACTCGACGGGCAAGTCCGGTCCGCAGGCGCAGAAAACTGGGGAAGGTAATTATgagt                                           >  2:661080/1‑140 (MQ=255)
                                                                                         aaaTTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCCATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACGTTGTGGAAGATATTTATGTctct                                         >  1:97341/1‑143 (MQ=255)
                                                                                                tgtCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGCGCACGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTac                                  >  2:54049/1‑142 (MQ=255)
                                                                                                tgtCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACAGGCAAGTCCGGTCCGCTGGCGCAGTACATTGCGGAAGATATATATGTCGCTCGCGTaa                                  >  2:117065/1‑143 (MQ=255)
                                                                                                  tCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACGTTGTGGAAGATATTTATGTCTCTCGCGTaaaa                                >  2:244248/1‑143 (MQ=255)
                                                                                                   cgcgGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATGCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAaag                               <  2:576829/143‑1 (MQ=255)
                                                                                                    gcgGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCCATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAaaga                              <  2:255487/143‑1 (MQ=255)
                                                                                                     cgtcaaTCGGCTGTGGCGGCACACTACGCTGCGCCATCCATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAaagaa                             <  2:539821/139‑1 (MQ=255)
                                                                                                     cgGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAaagaa                             <  1:854401/143‑1 (MQ=255)
                                                                                                               cTGTCGCGGCACATTACGCTGCGGCATCCATCCGAAACGACCTATTCCCACCATTAATATCCACTCGACCGCCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGCTATTTATGTCTCTGGCGTAAAAGAAGAAAACATCa                   <  1:449127/143‑1 (MQ=255)
                                                                                                                  tgGCGCCACAGTACCCTCCGGCCTCTATCCGAAACGCCGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAAAGAAGAAAACATCACTg                <  2:49004/143‑1 (MQ=255)
                                                                                                                   ggcggcTCATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATGGCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAAAGAAGAAAACATCACTgt               <  2:105719/143‑1 (MQ=255)
                                                                                                                    gcggcACATTACGCTGCGGCATCTATCCGAAACGACGGATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGGCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAAAGAAGAAAACATCACTgta              <  2:108898/143‑1 (MQ=255)
                                                                                                                    gcggcACATTACGCTGCGGCATCCATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACTGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAAAGAAGAAAACATCACTgta              <  2:736434/143‑1 (MQ=255)
                                                                                                                     cggcACATTACGCTGCGGCATCTATCCGAAACGACGTCTTCCCACCATTCATATCCACTCGATGGACAAGTCCGGTCCGCTGGCGCAGTACATCGTGGAAGATATTTATGTCCCGGGCGTCACAGAAGAACCCATCACTgtac             >  2:162818/1‑142 (MQ=255)
                                                                                                                     cgacacaTTACGCTGCGGCATCTATCCGACACTACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGGCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAAAGAAGAAAACATCACTgtag             <  1:409516/140‑1 (MQ=255)
                                                                                                                               cGCTGGCGCATCTATCCGTAACGACGTATTCCCACTATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAAAGAAGAAAACATCACTGTAGTAGGTGATGc   <  2:749941/143‑1 (MQ=255)
                                                                                                                               cGCTGCGGCATCTATCCGAAACGACGGATGCCCACCATTAATATCCGCTCGACGGGCAAGGCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAAAGAAGAAAACATCACTGTAGTAGGTGATGc   <  1:517988/143‑1 (MQ=255)
                                                                                                                                gCTGCGGCATCCATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAAAGAAGAAAACATCACTGAAGTCGGTGGTGCg  >  2:877679/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                   
TACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATTTATGTCTCTGGCGTAAAAGAAGAAAACATCACTGTAGTAGGTGATGCG  >  NC_000913/2826490‑2826760

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