Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6610TG41.5% 29.2 / 20.7 55pseudogene (943/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base T (5/26);  new base G (21/1);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 4.05e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

CATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGG  >  NC_000913/3862518‑3862805
                                                                                                                                                |                                                                                                                                                 
cATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCggg                                                                                                                                             >  1:902311/1‑150 (MQ=255)
cATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCggg                                                                                                                                             >  2:802758/1‑150 (MQ=255)
   gTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATc                                                                                                                                         <  1:991787/151‑1 (MQ=255)
   gTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCgggga                                                                                                                                           >  1:18151/1‑147 (MQ=255)
     tGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCa                                                                                                                                        <  1:604068/150‑1 (MQ=255)
       cccACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCa                                                                                                                                    <  1:633/151‑1 (MQ=255)
       cccACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGTTGCggg                                                                                                                                             >  2:511030/1‑143 (MQ=255)
        ccaccaCACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGTggggaaca                                                                                                                                        >  1:47513/1‑142 (MQ=255)
         caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAg                                                                                                                                  <  2:157271/151‑1 (MQ=255)
         caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGc                                                                                                                                      >  2:1105383/1‑148 (MQ=255)
         caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGGTGCGGGGATCAGCa                                                                                                                                    >  1:351332/1‑149 (MQ=255)
          accacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAg                                                                                                                                  <  2:974883/150‑1 (MQ=255)
            cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGc                                                                                                                                 <  2:487864/149‑1 (MQ=255)
                  aaaGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTTAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGAACAGCAAACCCGcccc                                                                                                                          >  1:1000271/1‑148 (MQ=255)
                    agagATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTTGGGGAAAAAGCAAACCCCCCCCttt                                                                                                                       >  1:935567/1‑151 (MQ=255)
                             ttACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATAATCAAACCCGCCACCTTTTCCCGAtt                                                                                                              >  2:328574/1‑151 (MQ=255)
                              tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGGATTAGCAAGACCGCCCCTTTTTCCCGCttt                                                                                                             >  1:1036633/1‑151 (MQ=255)
                                cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATTAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCTGCATGGCGGCTAAAGGTGGGGGAAAAAGGAAACCCCCCCCCTTTCCCCCAtgt                                                                                                             >  1:444784/1‑147 (MQ=255)
                                cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGAACAGCAAACCCGCCACCTTTACCCCAtttt                                                                                                            >  1:1110584/1‑150 (MQ=255)
                                        aTTGAGCCCGCCAGCACGGGGTGCACCGCAAACAGCAACGGAAAAATGAACAGAAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCATGGCGTTTAAGGTGGCGGAAATCAGCACGCCCGCCACTTTTACCCGATTTTCCGGTg                                                                                                      <  2:419164/148‑1 (MQ=255)
                                           gcggCCCCCAGCAGGGCGCCTCCCGAAAAGAGCACCGGTAAAAGAAACAGGAAGGTAAATCCCGGCGGTGCGATAATCCCCACCAGAATGGGGTTAAGGGTTGCGGCAAACACCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAg                                                                                                <  2:334836/149‑1 (MQ=255)
                                              gccgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATAAGGTTGCGGGAATCAGCCAGGCCGCCACTTTTTCCCCACTTTTCGGTGTgg                                                                                                   >  2:10241/1‑145 (MQ=255)
                                                  ccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGAATCAGCAAGCCCGCCCCTTTTTCCCCCTTTTCCGGTGCGGcgggg                                                                                              >  2:213463/1‑142 (MQ=255)
                                                   cAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGAACAGCAAGCCCGCCACTTTTACCCGCTTTTTCGGTGCGGGAggg                                                                                              >  2:420850/1‑143 (MQ=255)
                                                          gCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAACCCCCCACTTTTACCCGATTTTTCCGTGCGGCAGTGAAGTACATCGc                                                                                   >  1:1003469/1‑149 (MQ=255)
                                                             tgtACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATAAGGGTGCGGGAAACAGCAAGCCCGCCCCCTTTACCCGCTTTTCCGGTGCGGCAGTGAAGTACAACGCtca                                                                                >  2:456970/1‑147 (MQ=255)
                                                                 ccGCAAACAGCAAGGGTAAAAGAAACGGGAAGGTAAATTCCGGGGTTTCGTTAATCCCCCCCAGCATGGCGGTTAGGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaaa                                                                           <  1:929259/150‑1 (MQ=255)
                                                                  cGCAAACAGCAACGGTAAAATAAACAGGAAGGTAAATCCCGGCGGTTCGGTAATTCCCCCCAGCAGGGGGTTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaaa                                                                           <  1:401511/149‑1 (MQ=255)
                                                                   gCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTCAAGGGGGGGGAATCAGCAAGCCCGCCCATTTTTCCCCGTATTCCCGGGCGGCAGGGCGGAACCTCGGTTAAGAAATg                                                                         >  1:280253/1‑150 (MQ=255)
                                                                      accAGCACGGGAGAAAGGAACAGAAAGATAAATCCCAGCGTTTCGGTATTCCCCCCCAGCGTGGGGTTTAGGGTGGGGGAAAGCAACAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATg                                                                         <  2:581352/145‑1 (MQ=255)
                                                                        cAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTTCGGTGCGGCAGTGAAGGACATCGCTAAAGAAATGCCCa                                                                     >  1:551992/1‑149 (MQ=255)
                                                                              cGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAGAAAAGCCCACGGCAcc                                                              >  1:347763/1‑150 (MQ=255)
                                                                                gTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGATGCGGGAAACAGCAAGCCCCCCACTTTTATCCGATTTTTCGGTGCGGCAGCGGAGTACATCGCCAAAAAAATTCCCACGGcccc                                                              >  2:214607/1‑148 (MQ=255)
                                                                                   aaaaGAACAGGAAGGAAAATTCCGGCGTTTGGGTAATCCCCACCAGCAGGGCGGTTAAGGTGGCGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGa                                                         <  1:132012/150‑1 (MQ=255)
                                                                                    aaTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTTCGGGGCGGGAGGGAAGTACCTCGCTAAAGAAAATCCCACGGCACCAAAGAtt                                                       >  2:299723/1‑151 (MQ=255)
                                                                                     aTGAACAGGAAGGTAAATTCCGGGGTTTCGGTATTCCCCACCAGCATGGCGTTTAGGGTTGCGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttt                                                      <  1:1013111/151‑1 (MQ=255)
                                                                                     aTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTATGCTTCCGGGATCAGCCAACCCCCCCCTTTTACCCGATTTTCCTGCGCGGCAGTGGAGTACAACGCTAAAAAAAAGCCCACCGGGCCAAAGAttt                                                      >  1:544434/1‑151 (MQ=255)
                                                                                      taaaCAGGAAGGAAAATTCCAGCGGTTCGGTAATCCCCACCAGCTTGGCGGTTAGGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttt                                                      <  1:921932/148‑1 (MQ=255)
                                                                                      tGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCCCCACTTTTACCCCATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAGAAAAGCCCACGGCACCAAAGAtttt                                                     >  2:248559/1‑151 (MQ=255)
                                                                                        aacggaaaGGAAATTACCGGCGGTGCGTTATTCCCCACCAGCGGGGCGATTAAGGGGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttttt                                                    <  1:1062072/144‑1 (MQ=255)
                                                                                        aaaaGGAAGAAAATTTCCAGCGGTTCGGTAATCCCCACCAGCAGGGGGGTAAGGGTGGCGGAATTCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttttt                                                    <  2:489789/147‑1 (MQ=255)
                                                                                          cAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTCGGGAATCAGCAAGCCCCCCACTTTTACCCGAATTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAAACCCACCGCACCAAAGAttttt                                                    >  2:284695/1‑148 (MQ=255)
                                                                                            ggAAGGTAAATCCCGGCGGTTCGGTAATCCCCCCCAGCAGGGGGGTTAAGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGa                                                  <  2:523712/148‑1 (MQ=255)
                                                                                              aaGGAAAATCCCAGCGGTTCGGTAATCCCCACCAGCATGGGGTTTAGGGTGGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTAc                                             <  2:551866/151‑1 (MQ=255)
                                                                                                 gTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTTCGGGGCGGCAGTGAAGTACATCGCTAAAAAAATGCCCACGGCCCCAAAGATTTTTGAGTTACCg                                           >  1:629831/1‑150 (MQ=255)
                                                                                                   aaaTTCCGGCGGTTCGGTAATCCCCACCAGCAGGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGc                                        <  1:26347/151‑1 (MQ=255)
                                                                                                   aaaTCCCAGCGGTTCGTTAATCCCCACCAGCAGGGCGGTTAGGGTGGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGc                                        <  2:546857/151‑1 (MQ=255)
                                                                                                             ggTTCGGTAACCCCCACCAGCGTGGCGGTTAAGGTGGCGGAAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCaaaa                               <  2:1116550/150‑1 (MQ=255)
                                                                                                                gcgatAATTCCCCCAAGCGGGGCGGTAAAGGTTGCGGAACCAACAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGc                           <  2:265275/147‑1 (MQ=255)
                                                                                                                  ggTATTCCCCCCCAGCGGGGCGGTAAGGTTTGGGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAAcc                             <  1:93146/147‑1 (MQ=255)
                                                                                                                         ccccccAGCATGGGGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTcc                      <  2:458776/147‑1 (MQ=255)
                                                                                                                         ccccccAGCATGGCGTTTAAGGTGGGGGAAACCAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGa                  <  1:315070/151‑1 (MQ=255)
                                                                                                                                gCAGGGCGGTTAAGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaacaac             <  1:802964/149‑1 (MQ=255)
                                                                                                                                 cttGGCGGTTAAGGTCGCGGAAAACAGCACCGCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACg            <  1:1026315/147‑1 (MQ=255)
                                                                                                                                         ttAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCgg  >  2:341562/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGG  >  NC_000913/3862518‑3862805

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: