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

TTGATGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGG‑CGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAG  >  NC_000913/3578137‑3578435
                                                                                                                                                        |                                                                                                                                                     
ttGATGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGG‑CGATATCAAGGCATGGCGACTTGCTGTCCATCAGtt                                                                                                                                                       <  2:40529/151‑1 (MQ=255)
       tAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGG‑CGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTccacca                                                                                                                                                 >  2:43722/1‑150 (MQ=255)
               tcgtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGG‑CGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATa                                                                                                                                        >  1:82064/1‑150 (MQ=255)
               tcgtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGG‑CGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATa                                                                                                                                        >  1:478143/1‑150 (MQ=255)
                      gaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGG‑CGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCa                                                                                                                                 >  2:166902/1‑150 (MQ=255)
                                        catgggggCGAGCCCCGGGAGAAATAAGGGCAGTGGTCTTCTGGCGTGCTTCTTAAAATTTAAAAAAC‑CCGACGG‑CGATAACAAGGCAGGGCGCCTTGCTTTCCATCAGTCCCCCCAGGGA‑AATCCCCGCCACTTCAATCATCTTTAAgg                                                                                                               <  1:113451/144‑1 (MQ=255)
                                        cagtgtgGCGATGCCCGCGACAATTAGGGGCGTTGCTCATCTGGCGTGCTGCTTCAAATTTATCAAAC‑CCGCCGG‑CGATATCAAGGCATGCCGCTTTCCTTTCCATCAGCCCCCCCACGGG‑AATACCCGCCACTTCAATCATCTTTAAgg                                                                                                               <  1:376287/147‑1 (MQ=255)
                                        cagtggggCGATGCCCGCGAGAAATAAGGGGAGTGGTCATCGGGCGTGCTGCTACAAAGTTATAAAAC‑CCGACGG‑AAATATCAAGGCATGGCGACTTGCTGTCCATCAGTCCCACCCCGGG‑AATCCCCGCCACTTCAATCATCTTTAAgg                                                                                                               <  2:15547/144‑1 (MQ=255)
                                        cAATGTGGCGATCCCCGCGAGAAATAATGGCAGTGGTCATCTGGCGTGCTGCTTAAAAGTTATAAAAC‑CCGACGG‑CAATATCAAGGCATGGCGACTTGCTGTCCATCA‑TTCCCCCCCGGGAAATACCCGCCACTTCAATCATCTTTAAgg                                                                                                               <  2:248339/150‑1 (MQ=255)
                                        aaaTGGGGCAACGCCCGGAAGAAAAAATGGCGGTGGTCATGTGGCTTGCTGCTTAAAATTTAACAAAC‑CCGACGG‑CGATATAAGGGCATGGCTACTTGCTGTCCACCAGTTCCCCCAGGGA‑AATCCCCGCCACTTCAATCATCTTTAAgg                                                                                                               <  1:473266/149‑1 (MQ=255)
                                               gcgaccccggggAGCAATAATGGCAGTGTCCATCGGGGGTGCTGCTACAAAGTTA‑CAAAA‑CCCACGGCCGATAACAAGGCTTGGCGTTTGGCTGTCAATCAGTCCCCCCAGGGA‑AACCCCCGCCACTTCAATCATCTTTAAGGTGCgcgg                                                                                                        <  1:260138/139‑1 (MQ=255)
                                                      ccGCGAAAAATAAGGGCAGTGGTCCTCTGGGGTGCTCCTTAAAATTAAAAAAAC‑CCGCCGG‑CAATACCAGGGCATGGCGTCTTGCTGCCCATCATTCCCCCCCCGGA‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTggg                                                                                                <  1:508086/151‑1 (MQ=255)
                                                         cgAGAAATAATGGCGGGGGTCTTCGGGCGTGCTGCTTAAAAGTTATAAAAC‑CCGAGGG‑GGATATCAAGGCATGGCGCTTTGCTGTCCATGGGTCCCCCCCGGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTa                                                                                              <  1:403496/150‑1 (MQ=255)
                                                           agaaataatGGCAGTGTTCATCGGGGGTGCTGCTTCAAAGTTAAAAAAC‑CCGACGG‑CGATATAAGGGCATGGCGTCTTGCTGCCCATCAGTCCCACCACGGA‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTggg                                                                                                <  2:507329/143‑1 (MQ=255)
                                                           aaaaaaaagggCGGTGGTCATCTGGCGTGCTGC‑TAAAAATTAACAAAACCCGAGGG‑CGATACAAGGGCAGGGCGATTTGCTCTCCACCAGTCCCACCGGGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACg                                                                                            <  1:372050/141‑1 (MQ=255)
                                                           aGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAC‑CCGACGG‑CGATATCAAGGCATGGCGACTTGCTGTCCATCAGTCCCCCCCCGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACg                                                                                            <  2:236806/150‑1 (MQ=255)
                                                           aGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGG‑CGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTggg                                                                                                >  1:507319/1‑146 (MQ=255)
                                                            gCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGG‑CGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGtt                                                                                          >  2:104243/1‑151 (MQ=255)
                                                                        gTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGG‑CGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa                                                                              >  2:129994/1‑151 (MQ=255)
                                                                        gTGGTCATCGGGCGTGCTGCTTAAAAGTTAAAAAAC‑CCGCCGG‑CAATATCAGGGCAGGGCGCCTTGCTGTCCATCAGTCCCCCCGGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                               <  1:82584/150‑1 (MQ=255)
                                                                         tggtactCTGGCGTCCGTCTTCAAAGTTATCAACC‑CCGCCGC‑CATTATCAGGGCGTGGCGACTTTTTCTCTAACATTTCCACCGCGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                               <  1:466416/143‑1 (MQ=255)
                                                                          ggTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGG‑CGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAgg                                                                            >  2:269755/1‑151 (MQ=255)
                                                                                      tgctgcTTAAAAGTTAAAAAAC‑CCGACGG‑CGATATAAAGGCAGGGCTATTTGCTGCCCATCATTCCCCCCCCGGA‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTc                                                                <  2:416366/151‑1 (MQ=255)
                                                                                            ttaaaaTTTAACAAAC‑CCGGCGG‑CGATATCAAGGCATGGCGTCTTGCTGTCCATCATTCCCACCAGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGgag                                                          <  1:529048/148‑1 (MQ=255)
                                                                                             ttcaaGTTAACAAAACCCGACGG‑CGAAATCAAGGCAGGGCGATTTCCTGTCAATCTGTTCCACAGGGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGgag                                                          <  1:39004/148‑1 (MQ=255)
                                                                                                     aTCAAAA‑CCGACGG‑CGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAt                                                 >  2:468832/1‑151 (MQ=255)
                                                                                                      aaaaac‑CCGAGGG‑GAATTTCAAGGATGGGCGATTTGCTGTCCATCATTTCCCCCAGGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAtt                                                <  2:385136/149‑1 (MQ=255)
                                                                                                              cGACGG‑CGATATCAAGGCAGGGCGATTTGCTGTCCATCAGTCCCCCCACGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCg                                         <  2:350087/151‑1 (MQ=255)
                                                                                                                 cGG‑CGATAAAAGGGCATGCCTTCTTGCTGCCCATCATTCCCCCCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTc                                      <  1:359798/151‑1 (MQ=255)
                                                                                                                       atatCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTg                                   >  2:420132/1‑149 (MQ=255)
                                                                                                                       atatCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTg                                   >  1:178744/1‑149 (MQ=255)
                                                                                                                        tatCAGGGCAGGGCGATTTGCTGCCCATCAGTCCCCCCACGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTc                                <  2:50470/151‑1 (MQ=255)
                                                                                                                                cATGGCGATTTGCTGCCCATCAGTCCCCCCAGGGA‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTcc                         <  2:71517/150‑1 (MQ=255)
                                                                                                                                         ttgttgTCCATCGTTCCCACCGCGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAc               <  1:297672/151‑1 (MQ=255)
                                                                                                                                         ttGCTGTCCATCAGTCCCCCCGGGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAc               <  1:290022/151‑1 (MQ=255)
                                                                                                                                         ttGCGGTCCACCAGTCCCACCCGGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCGGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAc               <  2:188810/151‑1 (MQ=255)
                                                                                                                                              gTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAg          >  1:401585/1‑151 (MQ=255)
                                                                                                                                                     aGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCa   >  1:156620/1‑151 (MQ=255)
                                                                                                                                                      gTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAg  >  2:280317/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
TTGATGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGG‑CGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAG  >  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: ‑