Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6610TG39.5% 24.0 / 16.0 44pseudogene (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 (10/16);  new base G (15/2);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 1.55e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCG  >  NC_000913/3862513‑3862804
                                                                                                                                                    |                                                                                                                                               
ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGtt                                                                                                                                                <  1:175226/150‑1 (MQ=255)
ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGt                                                                                                                                                 <  1:704252/149‑1 (MQ=255)
ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGt                                                                                                                                                 <  2:309005/149‑1 (MQ=255)
 cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATTAGGtt                                                                                                                                                >  1:254239/1‑149 (MQ=255)
        gTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATc                                                                                                                                       <  1:813908/151‑1 (MQ=255)
             ccaccaCACCAAAGAGATACATTACGGTCGACATTAAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGTTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCggggaaaa                                                                                                                                      >  1:142197/1‑142 (MQ=255)
                 cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGccc                                                                                                                              <  2:721764/151‑1 (MQ=255)
                 cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAAccc                                                                                                                              >  1:77294/1‑151 (MQ=255)
                    accaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGTTGCGGGGATCAGCaaaaccccc                                                                                                                           >  2:682711/1‑144 (MQ=255)
                     ccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGGAACAGCAAACCCGcc                                                                                                                           >  1:37525/1‑150 (MQ=255)
                      caAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGGAACAGCAAAcccccccc                                                                                                                         >  1:429229/1‑149 (MQ=255)
                      caAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGCCCGCCAc                                                                                                                         >  1:94278/1‑151 (MQ=255)
                           agaTAAATTACGATCGACATGGGGGCCCCCAGCAGGGCGTGTACCGAAAACAGAAGCGGGGAAAGGAACAGGAAGGTAAATTCCGGCGGTTCGGTAATCCCCACCAGCATGGCGGTTAAGGTTGCGGAAACCAGCAAGCCCGCCACtttt                                                                                                                     <  2:9302/150‑1 (MQ=255)
                                   tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTTCGGGGATTAGCAAGGCCGCCACTTTTACCCGCttt                                                                                                            >  2:616208/1‑151 (MQ=255)
                                    aCGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGAACAGCAAGCCCGCCACCTTTACCCGAattt                                                                                                           >  1:800555/1‑151 (MQ=255)
                                          gACATTCAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGCGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGTTTCCGGGCCCAGCCAACCCCCCcccttttcccccttttcc                                                                                                         >  2:598701/1‑147 (MQ=255)
                                              ttGCGGCCGCCAGCACGGCGCGTACCGCAAACAGCAACGGTGTAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGCTGCGGGAATTTGCAAGTGCGCCCATTTTTCCCGACTTTc                                                                                                          >  2:57063/1‑142 (MQ=255)
                                               tggggCCGCCAGCACGGCGTTTCCCGCAAACAGCAAGGTTGAAAAGAAAAGGAAGGTAAATTCCGGCGGTTCGGTAATTCCCCCCAGCATGGCGGTTAGGGTGGCGGAAATCCGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGc                                                                                                 <  2:507419/147‑1 (MQ=255)
                                                    ccgcaACCACGGCGTGTACCGCAAACAGCAAGGGTAAAAGAAACGGGAAGTAAAATCCCAGGGGTTAGGTAACCCCCACCAGCTGGGCGGTTAAGGGGGCGGAAGTAAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGt                                                                                              <  2:642013/148‑1 (MQ=255)
                                                     cgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGTGGGGATCAGCAAACCCCCCCCTTTTAACCCATTTTTCGGTGCGGCAggg                                                                                             >  1:120420/1‑146 (MQ=255)
                                                             cGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAAACAGCAAGCCCCCCACTTTTACCCGATTTTCCGGCGCGGCAGGGGAGGACAACGc                                                                                  >  2:389297/1‑151 (MQ=255)
                                                             cGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGAATCAGCAAACCCGCCACTTTTTCCCCCTTTTCCGGTGCGGGAGGGGAGTa                                                                                        >  2:426750/1‑145 (MQ=255)
                                                              ggggtgtCCCGAAAACAGAAACGGTGAAATAAACAGGAAGTTAAATCCCGGCGTTTCGGTAATCCCCCCCAGCATGGCGGTTAGGGTTGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGt                                                                                         <  1:952968/140‑1 (MQ=255)
                                                                   gtACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAACCCGCCCCTTTTAACCGATTTTTCGGAGCGGCAGTGAAGTACa                                                                                      >  1:210803/1‑141 (MQ=255)
                                                                       cGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCCCTTTTACCCGATTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAGaaa                                                                          >  1:734888/1‑149 (MQ=255)
                                                                        gCAAACAGCAACGGTGAAATGAAAAGGAAGATAAATTCCAGGGGTTCGGTAATCCCCACCAGCAGGGCGTTTAAGGTTGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAt                                                                         <  1:487197/149‑1 (MQ=255)
                                                                          aaaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGGACATCGCTAAAGAAATAccc                                                                     >  1:173193/1‑151 (MQ=255)
                                                                          aaaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCCGGAAACAGGAAGCCCGCCACTTTTACCCGCTTTTCCGGTGCGGCAGGGAAGGACATCGCTAAAGAAAAGccc                                                                     >  2:865768/1‑151 (MQ=255)
                                                                          aaaCAGAACCGGTAAAATGAACGGAAAGAAAAATCCCAGCGGTTCGGTAATCCCCACCAGCAGGGCGATTATGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATg                                                                        <  1:198589/148‑1 (MQ=255)
                                                                           aaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTTCGGGAATCAGCAAGCCCGCCACTTTTTCCCGATTTTCCGGTGCGGCAGTGAAGGACATCGCTAAAGAAATGccc                                                                     >  1:937104/1‑150 (MQ=255)
                                                                               gaaccGGCGAATTGAACAGAGAGGTAAATCCCGCCTGCCCGTGATTCCCCCCCTGCTTGGCAGTTAAGGGGGCGGAAAGCAACAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACg                                                                  <  2:496530/145‑1 (MQ=255)
                                                                                cAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAAACAGCAAGCCCCCCACTTTTACCCGATTTTCCGGGGCGGCAGTGAAGTACATCGCTAAAAAAAAACCCACGGc                                                                >  2:144075/1‑150 (MQ=255)
                                                                                    ggTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGTTGCGGGAAGCAGCAAACCCGCCCCTTTTTCCCCACTTTCCGGTGAGGCAGTGAAGTAAAACGCTAAAAAAATGACCCCGGCAcc                                                             >  2:951788/1‑149 (MQ=255)
                                                                                               cAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGAACAGCAAACCCGCCACTTTTTCCCGATTTTTCGGTTGGGCAGCGAAGAACATCGCTAAAAGAAACCCCCCGGCACCAAAGAttttt                                                   >  1:599032/1‑148 (MQ=255)
                                                                                                aGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCCGGGGAGTAAATCCCCTAAGAAATGCCCACCGCACCCAAGATTTTTTAGt                                               >  1:626467/1‑151 (MQ=255)
                                                                                                   aaGGTAAATCCCAGGGGTTCGGTAATCCCCCCCAGCATGGGGGTTAAGGTTGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTg                                                  <  2:3470/145‑1 (MQ=255)
                                                                                                     ggTAAATCCCGGCGGTTCGGTAATCCCCCCCAGCTTGGCGTTTAGGGTGGCGGGAATCAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTa                                             <  1:266592/148‑1 (MQ=255)
                                                                                                        aaaTTCCGGCGGTTCGTTAATCCCCACCAGCATGGCGTTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGc                                       <  1:539026/151‑1 (MQ=255)
                                                                                                           ttCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAACCCGCCACTTTTACCCGATTTTTCGGGGCGGCAGTGAAGGACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAgg                                    >  1:786742/1‑151 (MQ=255)
                                                                                                             cgggggtttCGTTAATCCCCCCCAGCAGGGGGGTAAGGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAgg                                    <  2:473063/142‑1 (MQ=255)
                                                                                                             ccAGCGGTCGGGTAATCCCCACAAGCATGGGGGTTATGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGc                                  <  2:43448/151‑1 (MQ=255)
                                                                                                              cAGCGGTTCGGTAATCCCCACCAGCAGGGGGTTTAGGGTTGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAggg                                   <  1:138400/149‑1 (MQ=255)
                                                                                                                  ggTTCGGTAATTCCCACCAGCATGGCGGATAAGGGTGTGGGAATCAGCAAGACCCCCACTTTTAACCGATTTTCCGGGGCGGGAGTGTAGTACATCGCTAAAAAAATTGCCACCGCACCAAAGATTTTTGAATTACCGTGCAGGGCaaa                               >  1:363402/1‑149 (MQ=255)
                                                                                                                    ttCGGTAATCCCCCCCAGCATGGCGGTTAGGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGCGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAAcc                            <  1:682987/150‑1 (MQ=255)
                                                                                                                     tCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCCCGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGc                          >  2:991542/1‑151 (MQ=255)
                                                                                                                                caccaGGATGGGGGTTAAGGTTGGGGAAACCACAACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCg                    <  2:975423/146‑1 (MQ=255)
                                                                                                                                  ccaGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaaca              >  1:708114/1‑150 (MQ=255)
                                                                                                                                              ttAAGTTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCg  <  2:37517/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCG  >  NC_000913/3862513‑3862804

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