Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6680AG46.9% 20.7 / 13.3 32I312I (ATT→ATCglvCputative PTS enzyme II component GlvC
Reads supporting (aligned to +/- strand):  ref base A (2/15);  new base G (15/0);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 2.42e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAA  >  NC_000913/3862428‑3862913
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                             
aaTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCaaac                                                                                                                                                                                                                                                    >  1:691646/1‑249 (MQ=255)
   ggCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGATATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAAGTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGccc                                                                                                                                                                                                                                                  <  1:806825/250‑1 (MQ=255)
      gATTTGGGTTAGCCTCAAATCCGCGGGTTTGCTGAACATCGGGATCCTGTTTTGCGGTAAAACCTAGTCAATCAGGCCTCCCCCCATGGTGCCCCCCACACCAAAGAGGTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAAGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGcc                                                                                                                                                                                                                                              <  2:663452/250‑1 (MQ=255)
          tGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGTCCGCCACtt                                                                                                                                                                                                                                         <  1:141753/250‑1 (MQ=255)
                           cgcgtgtttgtaaaccATCGGTACCCATTTTTGCGTAAAAACCTAGACAACCAGCCCCCCCCCATTGTGGCCCCCCACACAAAAGAAATACATTAGGGGCGTCATTGAGGCCGCCACCGGGGCGTGTCCCGCAAACAGCAGGGTGAAAATGAACGGGAAGGTAAATTCCAGGGGTGCGGAAATTCCCCCTAGCATGGGGGTTAGGTTTGCGTAAATAAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTg                                                                                                                                                                                                                       <  1:513923/236‑1 (MQ=255)
                                  ttGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCCGCAAGCCCCCCACTTTTTCCCGCttttttgggggggcagg                                                                                                                                                                                                                >  1:607653/1‑238 (MQ=255)
                                     cTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGCGGGGATCAGCAAGCACGCCACTTTTTCCCGCttttttggggctggaggg                                                                                                                                                                                                               >  2:300500/1‑235 (MQ=255)
                                      tGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGATCAGCCAGCCCCCCACCTTTTTCCCGTTTTTcgtggggggcgggatg                                                                                                                                                                                                             >  1:73864/1‑236 (MQ=255)
                                       gAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGGTGCGGGGAGCAGCAAGGCCCCCACCTTTTCCCCCATTTTctgtgggggggggg                                                                                                                                                                                                               >  2:48957/1‑235 (MQ=255)
                                            ccggaaTCCAGTTTTGCGGTAAAACCTATCCAACCACCCCCCCCCCCATGTTCCCCCCCCCCCCAAAGAGATCCATTCCGGCCGACATGGAGGCCGCCAGCACGGCGTGTACCGAAAACAGCAAGGGTGAAAGGAACAGGAAGGTAAATCCCGGCGTTTGGGTAATTCCCCCCAGCATGGCGGTTAAGGTTGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAt                                                                                                                                                                                                       <  1:300500/245‑1 (MQ=255)
                                                 aTCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGGATTAGAAAACCCCCCCCTTTTTACCGATTTTCCCGTGCGGCCGGGAAGTACATAGCtt                                                                                                                                                                                                  >  1:429649/1‑249 (MQ=255)
                                                                aaaaCCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGTTTCGGGGAACAGCAAACCCCCCCCTTTTATCCGCTTTTCCGGGGGGGGCGGGGAGTACATCGCTAAAGaaaaacccccggc                                                                                                                                                                                  >  1:667204/1‑241 (MQ=255)
                                                                          caataagccccccccccATTTTGCCCCCCACACAAAAGAGATACATCACGCGCGAAATAGAGCCCCCCACCTGGGCGTAAACCCAACACAGCAAGGGTAAAAAGAACAGAAAGGAAAATCCCAGCGGTTGGGTAATCCCCACCAGCATGGGGGTTAAGGTGGGGGGAAACAACAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAtt                                                                                                                                                                     <  1:43018/237‑1 (MQ=255)
                                                                                                                               cgaaaTTGAGGCCGCCAGCACGGCGTGTACCGAAAACAGCAACGGTAAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATCCCCCCCAGCGTGGGGTTTAGGGTGGGGGGAATCAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCg                                                                                                               <  1:441009/247‑1 (MQ=255)
                                                                                                                               cGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGTTGCGGGGATCAGCAAAACCCCCCCCTTTACCCCATTTTTCGGGGCGGCAGGGAAGTACATCGCTAAAGAAAAACCCCCGGCACCCAAAAATTTTTAGTTACCGGGCCGGGGAAAACCGCCCTCCCGGaaaaaacatttcaacg                                                                                                               >  1:584311/1‑237 (MQ=255)
                                                                                                                                gacatggggcccgccAGCCCGGCTTGTACCGCAAACAGAAGCGTGAAAAGAAACGGGAGGGAAAAGCCCAGCGGTTGGATAATTCCACCAGGCATGGGTGTTAGGGTGGGGGAAACCAGAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCgg                                                                                                              <  1:718609/241‑1 (MQ=255)
                                                                                                                                   atgggggCCCCCAGCACGGCGTGTACCGAAAACAGCAACGGGAAAATAAACAGGAAGGAAACCTCCGGGGGTTCGGTAATCCCCACCAGCATGGGGGTAAGGGTGGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTc                                                                                                           <  1:792948/246‑1 (MQ=255)
                                                                                                                                          ccgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGCGGGGATCAGCCAAACCCCCCCCTTTAACCCACTTTTCCGGGCGGCAGGGAAGGAAATCGCTAAAAAAAAGCCCCCGGCACCCAAAAATTTTTAGTTAACCGGCAGGGGAAAAACACCCTCCGGGGACAACGATTTTTGCGGCTccgggccc                                                                                                    >  1:475638/1‑244 (MQ=255)
                                                                                                                                                             cGCAAACAGCAATGGAGAAATAAACGGGAAGGTAAATTCCAGCGGTACGTTAATCCCCCCCAGCAGGGCGTTAAGGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAt                                                                                 <  1:115803/251‑1 (MQ=255)
                                                                                                                                                                        aaCGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGCTAAAGGTGTGGGGAACAGGAAAGCCGGCACTTTTTACCGATTTTTCGGGGCGGCAGGGGAGGACATCGCTAAAGAAAATGCCACGGCAACAAAGAATTTTTGGTTATCGTGCAGGGAAAAAACCGCTTCCGGGAACCACGAATTCAGCGGACCGGGGCTCGAAATGAACTCTTGAAGAtaaggcgcccc                                                                       >  1:778615/1‑241 (MQ=255)
                                                                                                                                                                         aCGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGATCAGCAAACCCGCCACTTTTTCCCCATTTTTCGGGGGGGCCGGGAAATACCTCGCTAAAGAAAAGCCCACCGCACCAAAAATTTTTGAGTTTACGTGCAGGGCAAAAACGCCCTCCGGGAAAAAAGATTTCAACGGGCCGGCGCTCCAACTGAACTCTTGCAGAtagctcgcccac                                                                      >  1:325157/1‑240 (MQ=255)
                                                                                                                                                                          agggtaaaaggAACAGAAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCAGGGGGGTAAAGGTTGGGGGAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                                      <  1:67058/241‑1 (MQ=255)
                                                                                                                                                                             tGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATAAGGATGCGGGGAACAGCAAGCCCGCCACTTTTTCCCGATTTTCCGGGGCGGCAGGGAAGGACAACGCTCAAGAAATACCCACGGCCCCCAAGATTTTTGAGTTAACGTGGAGGGGAAAACCGCCCTCCGGGAACAACGATTTCAGCGGCGCGGCGCTCCAAATGAACTCTCTCACATGGTGCGcccccgacc                                                                  >  1:59860/1‑245 (MQ=255)
                                                                                                                                                                              gtaaaagaaaagggAAGGTAAACTCCAGGGGTTCGTTAATCCCCAACAGCAGGGGGGTTAGGGGGGGGGAAAGCAGAAGCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTAc                                                                   <  2:667204/239‑1 (MQ=255)
                                                                                                                                                                                             aGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGATCAGCAAGCCCGCCCCTTTTACCCGATTTTTCCGTGCGGCAGTGGAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTTAGTTACCGTGCAAGGCAAAACCCCCCTCCGGGGACAACGATTTCAACGGCTCGGCGCTCAAACTGAACTCTTTCTGATACTGCGCCCCGTACATCtataagacgcct                                                    >  1:434041/1‑240 (MQ=255)
                                                                                                                                                                                                 aaatcccggCGGTTCGTTAATACCCCCCAGCTGGGCGTTTAGGGTGGGGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCaa                                                <  1:362015/243‑1 (MQ=255)
                                                                                                                                                                                                  aaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGACCGCCACTTTTACCCGCTTTTTCGGTGCGGCAGGGAAGGACATCGCTAAAGAAATACCCACCGCACCAAAGATTTTTTAGTTTCCGTGCAGGGCAAAAACGCCTTCCGGGGACAACGATTTCAACGGCTCGGGGCGCAAACTCAACCCTTGCAGGTGCTGCGCCCCGTACATCTTAATTCCGCCtttcac                                               >  1:814026/1‑247 (MQ=255)
                                                                                                                                                                                                        ggcgttacgttAATCCCCACCAGCAGGGCGGTTAAGGTTGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGc                                         <  1:382568/241‑1 (MQ=255)
                                                                                                                                                                                                           ggtgcggaattcccccccaGCAGGGCGGTTAGGGTGGTGGGAACAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCgg                                      <  1:430565/235‑1 (MQ=255)
                                                                                                                                                                                                           ggTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCCCTTTTACCCGATTTTTCGGGGCGGCAGTGGAGTACATCGCTTAAGAAAATCCCCCGGCACCAAAGATTTTTTAGTTACCGGGCAGGGCAAAACCGCCTTCCGGGAACAACGAtt                                                                                                                     >  1:37196/1‑172 (MQ=255)
                                                                                                                                                                                                                                      ggTTAAGGATGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCCGGGAAGTACATCCCTAAAGACATACCCACGGCCACAAAGATTTTTGAGTTACCGTGAAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCCGATGCTGCGCACAGCACATCTGAATGCCCCCTTCAACACATGCCGGACCAAAAATAAACTGTCCGTAGATgg            >  1:112203/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   cGGGGATCAGCAAGCAAGCCACTTTTACACGATTTTACGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCAACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACGGCCTTCCGGGAACCACGATTTCAGCGGCTCGGCGCGCACACTGACCTCTTGCCCATGCTGCGCCCACTACCGCTGCATACCGCCTTCAACAGCGGCCCGCCCAAAGATCCACTGTCCGTAGATGCAgtgtggtaa  >  2:264940/1‑244 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                             
AATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAA  >  NC_000913/3862428‑3862913

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