Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,826,6270TC27.8% 58.5 / 33.7 55Y80H (TAT→CAT) srlEglucitol/sorbitol‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base T (21/18);  new base C (5/10);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-02
Rejected as consensus: Frequency below/above cutoff threshold.

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

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