Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,578,2870TC42.9% 56.5 / 16.3 42E147G (GAA→GGA) gntRd‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (23/1);  new base C (0/18);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 5.37e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTG  >  NC_000913/3578138‑3578429
                                                                                                                                                     |                                                                                                                                              
tgatgaTAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTc                                                                                                                                               >  2:87453/1‑151 (MQ=255)
          atgTTCGTCGAGACGTGCGCCGAGATAGGCAATGAGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACg                                                                                                                                        >  2:392103/3‑148 (MQ=255)
          aCGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACggg                                                                                                                                      <  1:785996/150‑1 (MQ=255)
                   gagaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGCTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGcc                                                                                                                            >  2:136763/1‑151 (MQ=255)
                      aCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACt                                                                                                                         >  1:990887/1‑151 (MQ=255)
                            gcCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatc                                                                                                                   >  2:351323/1‑151 (MQ=255)
                            gcCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCCCCTTCAAAGTTATCAAAACCGACGGAGATATCAAGGCATGGCGGCTTGCTGTCCATCAGTTCCACCACCGGAATACCCGACACTTCAatc                                                                                                                   >  2:627479/1‑151 (MQ=255)
                             ccGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatc                                                                                                                   >  1:1008638/1‑150 (MQ=255)
                                  aTAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCttt                                                                                                             >  2:413361/1‑151 (MQ=255)
                                     ggCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTa                                                                                                            >  1:979245/1‑149 (MQ=255)
                                        aggggggCGTTCCCCGCAAGAAAAATGGGCTGTGGTCTTCGGGCGTGCTGCTAAAAATTAAACAAACCCGACGGCAATACAAAGGCGGGGCGACTTGCTGTCCATCATTCCCCCCACGGGAATCCCCGCCACTTCAATCATCTTTa                                                                                                            <  2:121657/141‑1 (MQ=255)
                                           gtgGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTgcgc                                                                                                    >  2:33942/1‑151 (MQ=255)
                                           ggggagccgCCCGGGAAAAATAAGGGCAGTGCTTTTCGGGCGTGCTCCTTAAAAGTAATAAAACCCGAGGGGGATACCAAGGCAGGGCGTCTTGGCGCCAAGCGGTCCCCCCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTgc                                                                                                      <  2:850845/141‑1 (MQ=255)
                                           ggggCGATCCCCGGGACAAATAAGGGCAGGGGTCTTCTGGCGTGTTCCTCAAAATTTAAAAAACCCGACGGCAATACAAAGGCAGGGCTCCTTGCTTCCCATCATTCCCCCCGGGGAAATACCCGCCACTTCAATCATCTTTAAGGt                                                                                                        <  1:407232/145‑1 (MQ=255)
                                           ggggCGATCCCCGCGAGCAAAAAGGGCGGGGGTCATCGGGCGTGCTGCTTAAAATTTAACAAACCCGACGGCGATATCAAGGCAGGGCGCCTTGCTGTCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTgc                                                                                                      <  1:133324/147‑1 (MQ=255)
                                            ggccggccccggggagaAAGAATGGAAGTGTCCATGGGGGGGGCTGCTTCAAAGTAAACAAACCCGACGGCGATATCAAGGCGTGGCAACTCGCTCCCCATCAGTCCCCCCGGGGAAACCCCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                                   <  1:44939/138‑1 (MQ=255)
                                             ggCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCgcgg                                                                                                  >  2:887497/1‑151 (MQ=255)
                                                             aataatGGCAGTGTTCATCTGGCGTTCTTCTAAAATTTTAACAAACCCGACGGCAATACCAGGGCAGGCCTTCTTGCTCTCCATCAGTCCCCCCACGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTAc                                                                                       <  1:233128/146‑1 (MQ=255)
                                                               taatGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCgg                                                                                 >  2:332834/1‑150 (MQ=255)
                                                               taatGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGt                                                                                >  1:745622/1‑151 (MQ=255)
                                                                atgggagttGGCCATCTGCCGTGCTCTTAAAAATTTAAAAAACCCGCCGGCGATATAAGGGCAGGGCACTTTGCTGCCCCTCATTCCCCCCGGGAAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCgg                                                                                 <  2:775415/141‑1 (MQ=255)
                                                                 atGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTg                                                                               >  1:389400/1‑150 (MQ=255)
                                                                  gggCAGGGGTCATCTGGCGTGCTGCTTAAAAGTTAAAAAACCCGGGGGCAATATCAGGGCATGGCGATTTGCTGCCCATCAGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAg                                                                             <  1:78260/150‑1 (MQ=255)
                                                                        tGGTCTTCTGGGGTGCTGCTTAAAAGTTAAAAAACCCGCCGGCAATATAAGGGCATGGCGATTTGCTTTCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                         <  1:351374/149‑1 (MQ=255)
                                                                               cTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCg                                                                >  2:826431/1‑151 (MQ=255)
                                                                                            tCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                      >  1:52494/1‑148 (MQ=255)
                                                                                              aaattaataaaccccgccgGCAATATCAAGGCATGGCGTTTTCCTGCCCATCGGTCCCCCCGCGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGga                                                     <  2:656365/134‑1 (MQ=255)
                                                                                                       aaccccgcggggaatatCAAGCCATGGCTTTTTGCGGTCCAGGAGTCCCACCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTc                                         <  1:650239/137‑1 (MQ=255)
                                                                                                         accggccGGCAATATCAAGGCATGCCGACTTGCTGTCCATCAGTCCCCCCGGGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAg                                      <  1:576727/145‑1 (MQ=255)
                                                                                                           ccGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGc                                    >  2:736585/1‑151 (MQ=255)
                                                                                                           ccGACGGCGATACAAAGGCAGGGCTACTTGTTCCCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAg                                      <  2:218802/149‑1 (MQ=255)
                                                                                                             gACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGt                                  >  1:356065/1‑151 (MQ=255)
                                                                                                                      atCAGGGCAGGGCGATTTGCTGCCCATCAGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTcc                         <  1:855714/151‑1 (MQ=255)
                                                                                                                           ggCATGGCGACTTGCTGCCCATCAGTCCCCCCACGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAttt                     <  1:174657/150‑1 (MQ=255)
                                                                                                                           ggCATGGCGACTTGCTGCCCATCAGTCCCACCAGGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAt                       <  1:1008462/148‑1 (MQ=255)
                                                                                                                             cATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCg                  >  2:16158/1‑151 (MQ=255)
                                                                                                                                ggCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGtt               >  2:970598/1‑151 (MQ=255)
                                                                                                                                   gACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtat            >  1:189005/1‑151 (MQ=255)
                                                                                                                                      ttgttgCCCATCAGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTcc                   <  1:161709/141‑1 (MQ=255)
                                                                                                                                             ccATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTg  >  2:393332/1‑151 (MQ=255)
                                                                                                                                             ccATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGCGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTg  >  2:680455/1‑151 (MQ=255)
                                                                                                                                              cATCAGTCCCCCCAGGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGt   <  1:591391/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TGATGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTG  >  NC_000913/3578138‑3578429

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