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

CATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGG  >  NC_000913/3862518‑3862805
                                                                                                                                                |                                                                                                                                                 
cATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCggg                                                                                                                                             >  1:902311/1‑150 (MQ=255)
cATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑AGGTGCggg                                                                                                                                             >  2:802758/1‑150 (MQ=255)
   gTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATc                                                                                                                                         <  1:991787/151‑1 (MQ=255)
   gTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCgggga                                                                                                                                           >  1:18151/1‑147 (MQ=255)
     tGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCa                                                                                                                                        <  1:604068/150‑1 (MQ=255)
       cccACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑a                                                                                                                                    <  1:633/151‑1 (MQ=255)
       cccACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTA‑AGTTGCggg                                                                                                                                             >  2:511030/1‑143 (MQ=255)
        ccaccaCACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGTggggaaca                                                                                                                                        >  1:47513/1‑142 (MQ=255)
         caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAg                                                                                                                                  <  2:157271/151‑1 (MQ=255)
         caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGGATCAGc                                                                                                                                      >  2:1105383/1‑148 (MQ=255)
         caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTA‑GGGTGCGGGGATCAGC‑a                                                                                                                                    >  1:351332/1‑149 (MQ=255)
          accacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAg                                                                                                                                  <  2:974883/150‑1 (MQ=255)
            cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAGc                                                                                                                                 <  2:487864/149‑1 (MQ=255)
                  aaaGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTTAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGGAACAGC‑AAACCCGcccc                                                                                                                          >  1:1000271/1‑148 (MQ=255)
                    agagATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTTGGGGAAAAAGC‑AAACCCCCCCCttt                                                                                                                       >  1:935567/1‑151 (MQ=255)
                             ttACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGGATAATC‑AAACCCGCCACCTTTTCCCGAtt                                                                                                              >  2:328574/1‑151 (MQ=255)
                              tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAA‑GGGTGCGGGGATTAGC‑AAGACCGCCCCTTTTTCCCGCttt                                                                                                             >  1:1036633/1‑151 (MQ=255)
                                cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATTAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCTGCATGGCGGCTAA‑AGGTGGGGGAAAAAGG‑AAACCCCCCCCCTTTCCCCCAtgt                                                                                                             >  1:444784/1‑147 (MQ=255)
                                cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGGAACAGC‑AAACCCGCCACCTTTACCCCAtttt                                                                                                            >  1:1110584/1‑150 (MQ=255)
                                        aTTGAGCCCGCCAGCACGGGGTGCACCGCAAACAGCAACGGAAAAATGAACAGAAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCATGGCGTTTAA‑GGTGGCGGAAATCAGC‑ACGCCCGCCACTTTTACCCGATTTTCCGGTg                                                                                                      <  2:419164/148‑1 (MQ=255)
                                           gcggCCCCCAGCAGGGCGCCTCCCGAAAAGAGCACCGGTAAAAGAAACAGGAAGGTAAATCCCGGCGGTGCGATAATCCCCACCAGAATGG‑GGTTAAGGGTTGCGGCAAACACC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAg                                                                                                <  2:334836/149‑1 (MQ=255)
                                              gccgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATAA‑GGTTGCGGGAATCAGC‑CAGGCCGCCACTTTTTCCCCACTTTTCGGTGTgg                                                                                                   >  2:10241/1‑145 (MQ=255)
                                                  ccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑AGGTGCGGGAATCAGC‑AAGCCCGCCCCTTTTTCCCCCTTTTCCGGTGCGGcgggg                                                                                              >  2:213463/1‑142 (MQ=255)
                                                   cAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGGAACAGC‑AAGCCCGCCACTTTTACCCGCTTTTTCGGTGCGGGAggg                                                                                              >  2:420850/1‑143 (MQ=255)
                                                          gCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGGATCAGC‑AAACCCCCCACTTTTACCCGATTTTTCCGTGCGGCAGTGAAGTACATCGc                                                                                   >  1:1003469/1‑149 (MQ=255)
                                                             tgtACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATAA‑GGGTGCGGGAAACAGC‑AAGCCCGCCCCCTTTACCCGCTTTTCCGGTGCGGCAGTGAAGTACAACGCtca                                                                                >  2:456970/1‑147 (MQ=255)
                                                                 ccGCAAACAGCAAGGGTAAAAGAAACGGGAAGGTAAATTCCGGGGTTTCGTTAATCCCCCCCAGCATGGCGGTTAG‑GGTTGCGGAAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaaa                                                                           <  1:929259/150‑1 (MQ=255)
                                                                  cGCAAACAGCAACGGTAAAATAAACAGGAAGGTAAATCCCGGCGGTTCGGTAATTCCCCCCAGCAGGGGGTTTAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaaa                                                                           <  1:401511/149‑1 (MQ=255)
                                                                   gCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTCA‑AGGGGGGGGAATCAGC‑AAGCCCGCCCATTTTTCCCCGTATTCCCGGGCGGCAGGGCGGAACCTCGGTTAAGAAATg                                                                         >  1:280253/1‑150 (MQ=255)
                                                                      accAGCACGGGAGAAAGGAACAGAAAGATAAATCCCAGCGTTTCGGTATTCCCCCCCAGCGTGGGGTTTAG‑GGTGGGGGAAAGCAAC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATg                                                                         <  2:581352/145‑1 (MQ=255)
                                                                        cAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTTCGGTGCGGCAGTGAAGGACATCGCTAAAGAAATGCCCa                                                                     >  1:551992/1‑149 (MQ=255)
                                                                              cGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGGATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAGAAAAGCCCACGGCAcc                                                              >  1:347763/1‑150 (MQ=255)
                                                                                gTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGATGCGGGAAACAGC‑AAGCCCCCCACTTTTATCCGATTTTTCGGTGCGGCAGCGGAGTACATCGCCAAAAAAATTCCCACGGcccc                                                              >  2:214607/1‑148 (MQ=255)
                                                                                   aaaaGAACAGGAAGGAAAATTCCGGCGTTTGGGTAATCCCCACCAGCAGGGCGGTTAA‑GGTGGCGGAAACCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGa                                                         <  1:132012/150‑1 (MQ=255)
                                                                                    aaTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTTCGGGGCGGGAGGGAAGTACCTCGCTAAAGAAAATCCCACGGCACCAAAGAtt                                                       >  2:299723/1‑151 (MQ=255)
                                                                                     aTGAACAGGAAGGTAAATTCCGGGGTTTCGGTATTCCCCACCAGCATGGCGTTTAG‑GGTTGCGGAAATCAGA‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttt                                                      <  1:1013111/151‑1 (MQ=255)
                                                                                     aTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTA‑TGCTTCCGGGATCAGC‑CAACCCCCCCCTTTTACCCGATTTTCCTGCGCGGCAGTGGAGTACAACGCTAAAAAAAAGCCCACCGGGCCAAAGAttt                                                      >  1:544434/1‑151 (MQ=255)
                                                                                      taaaCAGGAAGGAAAATTCCAGCGGTTCGGTAATCCCCACCAGCTTGGCGGTTAG‑GGTTGGGGAAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttt                                                      <  1:921932/148‑1 (MQ=255)
                                                                                      tGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGAATCAGC‑AAGCCCCCCACTTTTACCCCATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAGAAAAGCCCACGGCACCAAAGAtttt                                                     >  2:248559/1‑151 (MQ=255)
                                                                                        aacggaaaGGAAATTACCGGCGGTGCGTTATTCCCCACCAGCGGGGCGATTAA‑GGGGGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttttt                                                    <  1:1062072/144‑1 (MQ=255)
                                                                                        aaaaGGAAGAAAATTTCCAGCGGTTCGGTAATCCCCACCAGCAGGGGGGTAAG‑GGTGGCGGAATTCAGA‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttttt                                                    <  2:489789/147‑1 (MQ=255)
                                                                                          cAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑AGGTTCGGGAATCAGC‑AAGCCCCCCACTTTTACCCGAATTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAAACCCACCGCACCAAAGAttttt                                                    >  2:284695/1‑148 (MQ=255)
                                                                                            ggAAGGTAAATCCCGGCGGTTCGGTAATCCCCCCCAGCAGGGGGGTTAA‑GGTGGGGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGa                                                  <  2:523712/148‑1 (MQ=255)
                                                                                              aaGGAAAATCCCAGCGGTTCGGTAATCCCCACCAGCATGGGGTTTAG‑GGTGGCGGAAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTAc                                             <  2:551866/151‑1 (MQ=255)
                                                                                                 gTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTTCGGGGCGGCAGTGAAGTACATCGCTAAAAAAATGCCCACGGCCCCAAAGATTTTTGAGTTACCg                                           >  1:629831/1‑150 (MQ=255)
                                                                                                   aaaTTCCGGCGGTTCGGTAATCCCCACCAGCAGGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGc                                        <  1:26347/151‑1 (MQ=255)
                                                                                                   aaaTCCCAGCGGTTCGTTAATCCCCACCAGCAGGGCGGTTAG‑GGTGGCGGAAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGc                                        <  2:546857/151‑1 (MQ=255)
                                                                                                             ggTTCGGTAACCCCCACCAGCGTGGCGGTTAA‑GGTGGCGGAAATCAGC‑AACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCaaaa                               <  2:1116550/150‑1 (MQ=255)
                                                                                                                gcgatAATTCCCCCAAGCGGGGCGGTAAA‑GGTTGC‑GGAACCAACAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGc                           <  2:265275/147‑1 (MQ=255)
                                                                                                                  ggTATTCCCCCCCAGCGGGGCGGTAAG‑GTTTGGGGAAACCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAAcc                             <  1:93146/147‑1 (MQ=255)
                                                                                                                         ccccccAGCATGGGGGTTAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTcc                      <  2:458776/147‑1 (MQ=255)
                                                                                                                         ccccccAGCATGGCGTTTAA‑GGTGGGGGAAACCAGC‑AGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGa                  <  1:315070/151‑1 (MQ=255)
                                                                                                                                gCAGGGCGGTTAA‑GGTTGCGGAAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaacaac             <  1:802964/149‑1 (MQ=255)
                                                                                                                                 cttGGCGGTTAA‑GGTCGCGGAAAACAGC‑ACCGCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACg            <  1:1026315/147‑1 (MQ=255)
                                                                                                                                         ttAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCgg  >  2:341562/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGG  >  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: ‑