Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,578,2870TC50.0% 25.8 / 12.7 38E147G (GAA→GGA) gntRd‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (14/5);  new base C (0/19);  total (14/24)
Fisher's exact test for biased strand distribution p-value = 2.41e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGATGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAG  >  NC_000913/3578137‑3578435
                                                                                                                                                        |                                                                                                                                                     
ttGATGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGtt                                                                                                                                                       <  2:40529/151‑1 (MQ=255)
       tAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTccacca                                                                                                                                                 >  2:43722/1‑150 (MQ=255)
               tcgtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATa                                                                                                                                        >  1:82064/1‑150 (MQ=255)
               tcgtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATa                                                                                                                                        >  1:478143/1‑150 (MQ=255)
                      gaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCa                                                                                                                                 >  2:166902/1‑150 (MQ=255)
                                        catgggggCGAGCCCCGGGAGAAATAAGGGCAGTGGTCTTCTGGCGTGCTTCTTAAAATTTAAAAAACCCGACGGCGATAACAAGGCAGGGCGCCTTGCTTTCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAgg                                                                                                               <  1:113451/144‑1 (MQ=255)
                                        cagtgtgGCGATGCCCGCGACAATTAGGGGCGTTGCTCATCTGGCGTGCTGCTTCAAATTTATCAAACCCGCCGGCGATATCAAGGCATGCCGCTTTCCTTTCCATCAGCCCCCCCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                                                               <  1:376287/147‑1 (MQ=255)
                                        cagtggggCGATGCCCGCGAGAAATAAGGGGAGTGGTCATCGGGCGTGCTGCTACAAAGTTATAAAACCCGACGGAAATATCAAGGCATGGCGACTTGCTGTCCATCAGTCCCACCCCGGGAATCCCCGCCACTTCAATCATCTTTAAgg                                                                                                               <  2:15547/144‑1 (MQ=255)
                                        cAATGTGGCGATCCCCGCGAGAAATAATGGCAGTGGTCATCTGGCGTGCTGCTTAAAAGTTATAAAACCCGACGGCAATATCAAGGCATGGCGACTTGCTGTCCATCATTCCCCCCCGGGAAATACCCGCCACTTCAATCATCTTTAAgg                                                                                                               <  2:248339/150‑1 (MQ=255)
                                        aaaTGGGGCAACGCCCGGAAGAAAAAATGGCGGTGGTCATGTGGCTTGCTGCTTAAAATTTAACAAACCCGACGGCGATATAAGGGCATGGCTACTTGCTGTCCACCAGTTCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAgg                                                                                                               <  1:473266/149‑1 (MQ=255)
                                               gcgaccccggggAGCAATAATGGCAGTGTCCATCGGGGGTGCTGCTACAAAGTTACAAAACCCACGGCCGATAACAAGGCTTGGCGTTTGGCTGTCAATCAGTCCCCCCAGGGAAACCCCCGCCACTTCAATCATCTTTAAGGTGCgcgg                                                                                                        <  1:260138/139‑1 (MQ=255)
                                                      ccGCGAAAAATAAGGGCAGTGGTCCTCTGGGGTGCTCCTTAAAATTAAAAAAACCCGCCGGCAATACCAGGGCATGGCGTCTTGCTGCCCATCATTCCCCCCCCGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTggg                                                                                                <  1:508086/151‑1 (MQ=255)
                                                         cgAGAAATAATGGCGGGGGTCTTCGGGCGTGCTGCTTAAAAGTTATAAAACCCGAGGGGGATATCAAGGCATGGCGCTTTGCTGTCCATGGGTCCCCCCCGGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTa                                                                                              <  1:403496/150‑1 (MQ=255)
                                                           agaaataatGGCAGTGTTCATCGGGGGTGCTGCTTCAAAGTTAAAAAACCCGACGGCGATATAAGGGCATGGCGTCTTGCTGCCCATCAGTCCCACCACGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTggg                                                                                                <  2:507329/143‑1 (MQ=255)
                                                           aaaaaaaagggCGGTGGTCATCTGGCGTGCTGCTAAAAATTAACAAAACCCGAGGGCGATACAAGGGCAGGGCGATTTGCTCTCCACCAGTCCCACCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACg                                                                                            <  1:372050/141‑1 (MQ=255)
                                                           aGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAACCCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTCCCCCCCCGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACg                                                                                            <  2:236806/150‑1 (MQ=255)
                                                           aGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTggg                                                                                                >  1:507319/1‑146 (MQ=255)
                                                            gCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGtt                                                                                          >  2:104243/1‑151 (MQ=255)
                                                                        gTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa                                                                              >  2:129994/1‑151 (MQ=255)
                                                                        gTGGTCATCGGGCGTGCTGCTTAAAAGTTAAAAAACCCGCCGGCAATATCAGGGCAGGGCGCCTTGCTGTCCATCAGTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                               <  1:82584/150‑1 (MQ=255)
                                                                         tggtactCTGGCGTCCGTCTTCAAAGTTATCAACCCCGCCGCCATTATCAGGGCGTGGCGACTTTTTCTCTAACATTTCCACCGCGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                               <  1:466416/143‑1 (MQ=255)
                                                                          ggTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAgg                                                                            >  2:269755/1‑151 (MQ=255)
                                                                                      tgctgcTTAAAAGTTAAAAAACCCGACGGCGATATAAAGGCAGGGCTATTTGCTGCCCATCATTCCCCCCCCGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTc                                                                <  2:416366/151‑1 (MQ=255)
                                                                                            ttaaaaTTTAACAAACCCGGCGGCGATATCAAGGCATGGCGTCTTGCTGTCCATCATTCCCACCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGgag                                                          <  1:529048/148‑1 (MQ=255)
                                                                                             ttcaaGTTAACAAAACCCGACGGCGAAATCAAGGCAGGGCGATTTCCTGTCAATCTGTTCCACAGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGgag                                                          <  1:39004/148‑1 (MQ=255)
                                                                                                     aTCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAt                                                 >  2:468832/1‑151 (MQ=255)
                                                                                                      aaaaacCCGAGGGGAATTTCAAGGATGGGCGATTTGCTGTCCATCATTTCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAtt                                                <  2:385136/149‑1 (MQ=255)
                                                                                                              cGACGGCGATATCAAGGCAGGGCGATTTGCTGTCCATCAGTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCg                                         <  2:350087/151‑1 (MQ=255)
                                                                                                                 cGGCGATAAAAGGGCATGCCTTCTTGCTGCCCATCATTCCCCCCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTc                                      <  1:359798/151‑1 (MQ=255)
                                                                                                                       atatCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTg                                   >  2:420132/1‑149 (MQ=255)
                                                                                                                       atatCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTg                                   >  1:178744/1‑149 (MQ=255)
                                                                                                                        tatCAGGGCAGGGCGATTTGCTGCCCATCAGTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTc                                <  2:50470/151‑1 (MQ=255)
                                                                                                                                cATGGCGATTTGCTGCCCATCAGTCCCCCCAGGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTcc                         <  2:71517/150‑1 (MQ=255)
                                                                                                                                         ttgttgTCCATCGTTCCCACCGCGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAc               <  1:297672/151‑1 (MQ=255)
                                                                                                                                         ttGCTGTCCATCAGTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAc               <  1:290022/151‑1 (MQ=255)
                                                                                                                                         ttGCGGTCCACCAGTCCCACCCGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCGGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAc               <  2:188810/151‑1 (MQ=255)
                                                                                                                                              gTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAg          >  1:401585/1‑151 (MQ=255)
                                                                                                                                                     aGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCa   >  1:156620/1‑151 (MQ=255)
                                                                                                                                                      gTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAg  >  2:280317/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
TTGATGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAG  >  NC_000913/3578137‑3578435

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: