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

GCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGTTGCGGG‑AATC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCG  >  NC_000913/3862515‑3862809
                                                                                                                                                         |                                                                                                                                                    
gCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGTTGCg                                                                                                                                                     <  1:505499/151‑1 (MQ=255)
  ccATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAAGGTTgggg                                                                                                                                                    >  2:579138/1‑151 (MQ=255)
  ccATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGTTGCggg                                                                                                                                                   <  1:704818/151‑1 (MQ=255)
      gTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGTTGCGGG‑AATc                                                                                                                                              <  1:643397/151‑1 (MQ=255)
        tGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGGTGCggggaac                                                                                                                                               >  1:907707/1‑145 (MQ=255)
               cataccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATTAACAGTAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGTTGCGGG‑AATC‑AGCAAGccc                                                                                                                                    <  2:917933/148‑1 (MQ=255)
                     aaaGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑GT‑AAAGTTGCGGG‑GATC‑AGCAAAcccccc                                                                                                                                 >  2:96250/1‑148 (MQ=255)
                             aCATTACGGTCGACATGGAGGCCGCCAGCCCGGCGTGTACCGAAAACAGCAACGGTAAAATGAACAGGAAGGAAAATT‑CCGGCGGTTCGG‑TAATCCCCACCAGCAGGGCGG‑TT‑AAGGTGGCGGA‑AATC‑AGCAAGCCCGCCA‑CTTTTACCCg                                                                                                                     <  1:817419/151‑1 (MQ=255)
                                   cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGGTTT‑ATGGTTGCGGG‑AACA‑ACAAAGCCCCCCC‑CTTTT‑CCCCAtttt                                                                                                                >  2:580412/1‑150 (MQ=255)
                                             tGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGGTGGGGG‑AATC‑CGCAAACCCCCCACCTTTTCCCCCATTTT‑CGGTGCgg                                                                                                       >  2:985850/1‑149 (MQ=255)
                                              gAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGGTGCGGG‑GATC‑AGCCAGCCCCCCC‑CTTTTACCCCGTTTTTCGGGGCgg                                                                                                       >  2:341911/1‑148 (MQ=255)
                                                    gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑GT‑AAGGGTGCGGG‑GATC‑AAGAAAACCGCCC‑CCTTTTTCCCCTTTTTCGGGGCgggcgtgga                                                                                                >  2:189694/1‑142 (MQ=255)
                                                         cACGGCGTGTCCCGAAAACAGCAACGGAGAAATGAAAAGGAAGTTAAATC‑CCGGCGGTTCGG‑TAATCCCCACCAGCAGGGGGG‑TT‑AAGTTGGCGGA‑AATC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGt                                                                                              <  1:258138/146‑1 (MQ=255)
                                                         cACGGCGTGCACCGCAAACACCAACGGGAAAAAGAACAGAAAGGTAATTC‑CCAGCGGTTCGA‑AAATCCCCCCCACCGGGGCGT‑TT‑AGGGTGGGGGA‑AATA‑AGCAACCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAt                                                                                          <  2:1011179/150‑1 (MQ=255)
                                                           cGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAAGGTTTGGGG‑AA‑A‑AGCAAGACCACCC‑CTTTTACCCCCTTTTTCGGGTTGGCAGTGaa                                                                                                >  1:960286/1‑142 (MQ=255)
                                                            ggCGTGTCCCGCAAACAGCAACGGTAAAATGAACAGGAAGGTAAATT‑CCGGCGGTTGGG‑TAATCCCCACCAGCATGGCGT‑TA‑AAGGTGGGGGG‑AATC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATc                                                                                         <  1:847722/148‑1 (MQ=255)
                                                              cGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGATTGGGGG‑CATC‑AGTCAACCCCCCA‑CTTTTTCCCGATTTTCCGGTGAGGGAGTGAAgcatatagc                                                                                       >  2:807086/1‑140 (MQ=255)
                                                                     cGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAAGTTGCGGG‑AATC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGGAGTACATCGCTAAAGaaa                                                                               >  1:430362/1‑149 (MQ=255)
                                                                        aaaCAGC‑ACGGGTAAAGGAAAAAGAAAGGAAATTCCCGGGGGTGGGG‑TATTCCCCACAAGCATGGGGG‑TT‑AGGGTGGGGGA‑AATC‑AGCAACCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAt                                                                              <  2:398031/147‑1 (MQ=255)
                                                                        aaaCACCAACGGTAAAATGAACAGAAAGATAAATC‑CCAGCGTTTCGG‑TAATTCCCCCCAGC‑TGGGCG‑TT‑TAGGTTGGGGGAAACC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGa                                                                                 <  1:494494/144‑1 (MQ=255)
                                                                                   gtaaaagaaaCAGGAAGGTAAATA‑CATGCGGTGCGG‑TAACCCCCCCCGGC‑TGGGTG‑TT‑AAGGGTGGGGGAAATC‑AACAGGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCaa                                                                <  1:1007800/142‑1 (MQ=255)
                                                                                      aaaTGAACAGGAAGGTAAATC‑CCGGCGGTTGGG‑TATTCCCCCCCAGCAGGGCGG‑TA‑AGGGTGGGGGA‑AATG‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCaaa                                                               <  1:532008/148‑1 (MQ=255)
                                                                                          gAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGGTGGGGG‑GAAC‑AGCCAGGCCGCCA‑CTTTTAACCGCTTTTCCGGTGCGGGAGTGAAGTGAATAGCTAAAGAAAATCCCCCAGCGCCAAAGAttttt                                                        >  1:569745/1‑151 (MQ=255)
                                                                                            aaggggaGGAAAAAT‑CCGGGGGTTTGGATAATCCCCACCAGCTTGCCGT‑TT‑AGGGTTGCGGG‑AACC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttttt                                                        <  2:75236/144‑1 (MQ=255)
                                                                                            aCAGGAAGGAAAATC‑CCAGGGTTTCGG‑TAATCCCCCCCAGCAGGGCGG‑TT‑AGGGTGGGGGG‑AATC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGa                                                      <  1:317374/151‑1 (MQ=255)
                                                                                               ggAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGTTGCGGG‑AATC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCCGTGAAGTACATCGCTAAAGAAATTCCCACGGCACCAAAGATTTTTTAGtt                                                   >  1:682175/1‑151 (MQ=255)
                                                                                                gAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGGTGCGGG‑AATC‑AGCAAACCCCCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGtt                                                   >  1:102364/1‑150 (MQ=255)
                                                                                                 aaGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑AT‑AAGGTTTCGGG‑GAAC‑AGCAAGCCCGCCA‑CTTTTACCCGCTTTTCCGGTGGGGGAGTGGAGTACATCGCCAAAAAAATGCCCACGGCGCCCAAGATTTTTGAGtt                                                   >  2:34911/1‑149 (MQ=255)
                                                                                                 aaGGTAAATC‑CCGGGGGTTGGG‑TAATTCCCCCCAGCATGGCGT‑TT‑AAGGTTGCGGA‑AATC‑ACCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTa                                                  <  1:299302/150‑1 (MQ=255)
                                                                                                      aaaTT‑CCGGCGTTTCGG‑TAATCCCCCCCAGCATGGC‑G‑TT‑AAGGGTGGGGGAAATC‑AGAAACCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTg                                             <  1:846909/150‑1 (MQ=255)
                                                                                                      aaaTT‑CCGGCGGTGCGG‑TAATCCCCCCCAGCGTGGGGG‑TT‑AAGG‑TGGGGG‑AATCAAGCAGCCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGc                                            <  2:281800/151‑1 (MQ=255)
                                                                                                      aaaTC‑CCGGCGGTTTGG‑AAATCCCCACCAGCAGGGGGG‑TT‑AGGGTGGGGGA‑AATC‑AAAAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTg                                             <  2:559978/150‑1 (MQ=255)
                                                                                                       aaTT‑CCAGCGGTTCGA‑TAACCCCCACCAGAGGGGCGT‑TT‑AGGGTGGGGGA‑AATC‑AGAAGCCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCa                                           <  2:967037/151‑1 (MQ=255)
                                                                                                             cAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGGTGCGGG‑AATC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGGGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAggg                                        >  1:339097/1‑149 (MQ=255)
                                                                                                                  ggttggtaaaccccccccAGCAGGGTGG‑TT‑AAGGTGGCGGA‑AACC‑AGAACGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAggg                                        <  2:179668/133‑1 (MQ=255)
                                                                                                                    tCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGGTGCGGG‑AATC‑AACAAACCCGCCC‑CTTTTACCCGATTTTTCGGGGCGGCAGTGAAGTACATCGCTAAAAAAAATCCCACGGCACCAAAGATTTTTTAGTTACCGTGCCGGGCaaaaaccc                               >  2:767222/1‑147 (MQ=255)
                                                                                                                         tAATTCCCACCAGCATGGCGG‑TT‑AAGGTTGCGGG‑GATC‑AGCAAGCCCGCCC‑CTTTTTCCCGATTTTTCGGTGCGGCAGGGAAGTACATCGCTAAAGAAAATCCCACGGCACCAAAGATTTTTTAGTTACCGTGCAGGGCAAAACCGCCtt                            >  1:743825/1‑150 (MQ=255)
                                                                                                                              cccACCAGCGTGGGGG‑TT‑AAGGTGGGGGA‑AAGC‑AGAAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCt                             <  2:519566/144‑1 (MQ=255)
                                                                                                                                        tGGCGG‑TT‑AAGGTTGCGGG‑AATC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGAttt            >  1:692847/1‑151 (MQ=255)
                                                                                                                                         ggCGG‑TT‑AGGGTGGGGGG‑AATC‑AGAAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGAtt             <  1:432012/149‑1 (MQ=255)
                                                                                                                                          gCGG‑TT‑AAGGTGGCGGA‑AACC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCa          <  1:826020/151‑1 (MQ=255)
                                                                                                                                                 tagggTGGCGGG‑AACC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCgg      <  2:566566/147‑1 (MQ=255)
                                                                                                                                                  aaGGTGGCGGG‑AATC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCgg      <  1:34954/149‑1 (MQ=255)
                                                                                                                                                    ggTTGGGGA‑AATC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCg  <  1:677228/151‑1 (MQ=255)
                                                                                                                                                     gTTGCGGG‑AATC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCg  >  1:475004/1‑150 (MQ=255)
                                                                                                                                                         |                                                                                                                                                    
GCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TT‑AAGGTTGCGGG‑AATC‑AGCAAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCG  >  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: ‑