Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6640CG40.9% 12.7 / 24.9 44pseudogene (940/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 C (14/12);  new base G (4/14);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 6.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-01

GCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCG  >  NC_000913/3862515‑3862809
                                                                                                                                                         |                                                                                                                                                    
gCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCg                                                                                                                                                     <  1:505499/151‑1 (MQ=255)
  ccATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTgggg                                                                                                                                                    >  2:579138/1‑151 (MQ=255)
  ccATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCggg                                                                                                                                                   <  1:704818/151‑1 (MQ=255)
      gTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATc                                                                                                                                              <  1:643397/151‑1 (MQ=255)
        tGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCggggaac                                                                                                                                               >  1:907707/1‑145 (MQ=255)
               cataccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATTAACAGTAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGccc                                                                                                                                    <  2:917933/148‑1 (MQ=255)
                     aaaGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGTTGCGGGGATCAGCAAAcccccc                                                                                                                                 >  2:96250/1‑148 (MQ=255)
                             aCATTACGGTCGACATGGAGGCCGCCAGCCCGGCGTGTACCGAAAACAGCAACGGTAAAATGAACAGGAAGGAAAATTCCGGCGGTTCGGTAATCCCCACCAGCAGGGCGGTTAAGGTGGCGGAAATCAGCAAGCCCGCCACTTTTACCCg                                                                                                                     <  1:817419/151‑1 (MQ=255)
                                   cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTATGGTTGCGGGAACAACAAAGCCCCCCCCTTTTCCCCAtttt                                                                                                                >  2:580412/1‑150 (MQ=255)
                                             tGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGAATCCGCAAACCCCCCACCTTTTCCCCCATTTTCGGTGCgg                                                                                                       >  2:985850/1‑149 (MQ=255)
                                              gAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCCAGCCCCCCCCTTTTACCCCGTTTTTCGGGGCgg                                                                                                       >  2:341911/1‑148 (MQ=255)
                                                    gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGGATCAAGAAAACCGCCCCCTTTTTCCCCTTTTTCGGGGCgggcgtgga                                                                                                >  2:189694/1‑142 (MQ=255)
                                                         cACGGCGTGTCCCGAAAACAGCAACGGAGAAATGAAAAGGAAGTTAAATCCCGGCGGTTCGGTAATCCCCACCAGCAGGGGGGTTAAGTTGGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGt                                                                                              <  1:258138/146‑1 (MQ=255)
                                                         cACGGCGTGCACCGCAAACACCAACGGGAAAAAGAACAGAAAGGTAATTCCCAGCGGTTCGAAAATCCCCCCCACCGGGGCGTTTAGGGTGGGGGAAATAAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAt                                                                                          <  2:1011179/150‑1 (MQ=255)
                                                           cGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTTGGGGAAAAGCAAGACCACCCCTTTTACCCCCTTTTTCGGGTTGGCAGTGaa                                                                                                >  1:960286/1‑142 (MQ=255)
                                                            ggCGTGTCCCGCAAACAGCAACGGTAAAATGAACAGGAAGGTAAATTCCGGCGGTTGGGTAATCCCCACCAGCATGGCGTTAAAGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATc                                                                                         <  1:847722/148‑1 (MQ=255)
                                                              cGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGATTGGGGGCATCAGTCAACCCCCCACTTTTTCCCGATTTTCCGGTGAGGGAGTGAAgcatatagc                                                                                       >  2:807086/1‑140 (MQ=255)
                                                                     cGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGGAGTACATCGCTAAAGaaa                                                                               >  1:430362/1‑149 (MQ=255)
                                                                        aaaCAGCACGGGTAAAGGAAAAAGAAAGGAAATTCCCGGGGGTGGGGTATTCCCCACAAGCATGGGGGTTAGGGTGGGGGAAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAt                                                                              <  2:398031/147‑1 (MQ=255)
                                                                        aaaCACCAACGGTAAAATGAACAGAAAGATAAATCCCAGCGTTTCGGTAATTCCCCCCAGCTGGGCGTTTAGGTTGGGGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGa                                                                                 <  1:494494/144‑1 (MQ=255)
                                                                                   gtaaaagaaaCAGGAAGGTAAATACATGCGGTGCGGTAACCCCCCCCGGCTGGGTGTTAAGGGTGGGGGAAATCAACAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCaa                                                                <  1:1007800/142‑1 (MQ=255)
                                                                                      aaaTGAACAGGAAGGTAAATCCCGGCGGTTGGGTATTCCCCCCCAGCAGGGCGGTAAGGGTGGGGGAAATGAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCaaa                                                               <  1:532008/148‑1 (MQ=255)
                                                                                          gAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGGAACAGCCAGGCCGCCACTTTTAACCGCTTTTCCGGTGCGGGAGTGAAGTGAATAGCTAAAGAAAATCCCCCAGCGCCAAAGAttttt                                                        >  1:569745/1‑151 (MQ=255)
                                                                                            aaggggaGGAAAAATCCGGGGGTTTGGATAATCCCCACCAGCTTGCCGTTTAGGGTTGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttttt                                                        <  2:75236/144‑1 (MQ=255)
                                                                                            aCAGGAAGGAAAATCCCAGGGTTTCGGTAATCCCCCCCAGCAGGGCGGTTAGGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGa                                                      <  1:317374/151‑1 (MQ=255)
                                                                                               ggAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCCGTGAAGTACATCGCTAAAGAAATTCCCACGGCACCAAAGATTTTTTAGtt                                                   >  1:682175/1‑151 (MQ=255)
                                                                                                gAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAACCCCCCACTTTTACCCGATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGtt                                                   >  1:102364/1‑150 (MQ=255)
                                                                                                 aaGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATAAGGTTTCGGGGAACAGCAAGCCCGCCACTTTTACCCGCTTTTCCGGTGGGGGAGTGGAGTACATCGCCAAAAAAATGCCCACGGCGCCCAAGATTTTTGAGtt                                                   >  2:34911/1‑149 (MQ=255)
                                                                                                 aaGGTAAATCCCGGGGGTTGGGTAATTCCCCCCAGCATGGCGTTTAAGGTTGCGGAAATCACCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTa                                                  <  1:299302/150‑1 (MQ=255)
                                                                                                      aaaTTCCGGCGTTTCGGTAATCCCCCCCAGCATGGCGTTAAGGGTGGGGGAAATCAGAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTg                                             <  1:846909/150‑1 (MQ=255)
                                                                                                      aaaTTCCGGCGGTGCGGTAATCCCCCCCAGCGTGGGGGTTAAGGTGGGGGAATCAAGCAGCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGc                                            <  2:281800/151‑1 (MQ=255)
                                                                                                      aaaTCCCGGCGGTTTGGAAATCCCCACCAGCAGGGGGGTTAGGGTGGGGGAAATCAAAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTg                                             <  2:559978/150‑1 (MQ=255)
                                                                                                       aaTTCCAGCGGTTCGATAACCCCCACCAGAGGGGCGTTTAGGGTGGGGGAAATCAGAAGCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCa                                           <  2:967037/151‑1 (MQ=255)
                                                                                                             cAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGGGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAggg                                        >  1:339097/1‑149 (MQ=255)
                                                                                                                  ggttggtaaaccccccccAGCAGGGTGGTTAAGGTGGCGGAAACCAGAACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAggg                                        <  2:179668/133‑1 (MQ=255)
                                                                                                                    tCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAACAAACCCGCCCCTTTTACCCGATTTTTCGGGGCGGCAGTGAAGTACATCGCTAAAAAAAATCCCACGGCACCAAAGATTTTTTAGTTACCGTGCCGGGCaaaaaccc                               >  2:767222/1‑147 (MQ=255)
                                                                                                                         tAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAGCCCGCCCCTTTTTCCCGATTTTTCGGTGCGGCAGGGAAGTACATCGCTAAAGAAAATCCCACGGCACCAAAGATTTTTTAGTTACCGTGCAGGGCAAAACCGCCtt                            >  1:743825/1‑150 (MQ=255)
                                                                                                                              cccACCAGCGTGGGGGTTAAGGTGGGGGAAAGCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCt                             <  2:519566/144‑1 (MQ=255)
                                                                                                                                        tGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGAttt            >  1:692847/1‑151 (MQ=255)
                                                                                                                                         ggCGGTTAGGGTGGGGGGAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGAtt             <  1:432012/149‑1 (MQ=255)
                                                                                                                                          gCGGTTAAGGTGGCGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCa          <  1:826020/151‑1 (MQ=255)
                                                                                                                                                 tagggTGGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCgg      <  2:566566/147‑1 (MQ=255)
                                                                                                                                                  aaGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCgg      <  1:34954/149‑1 (MQ=255)
                                                                                                                                                    ggTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCg  <  1:677228/151‑1 (MQ=255)
                                                                                                                                                     gTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCg  >  1:475004/1‑150 (MQ=255)
                                                                                                                                                         |                                                                                                                                                    
GCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCG  >  NC_000913/3862515‑3862809

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