Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,578,3010AC43.6% 85.7 / 14.9 55G142G (GGT→GGGgntRd‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (25/6);  new base C (0/24);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 2.42e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGATA  >  NC_000913/3578153‑3578447
                                                                                                                                                           |                                                                                                                                                  
cgtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATAcc                                                                                                                                                  <  2:1032336/151‑1 (MQ=255)
 gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATAccc                                                                                                                                                 >  1:4528/1‑151 (MQ=255)
 gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATAccc                                                                                                                                                 <  2:1114903/151‑1 (MQ=255)
    gagaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGcc                                                                                                                                              >  2:149913/1‑151 (MQ=255)
      gaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCa                                                                                                                                             <  1:633439/150‑1 (MQ=255)
                      ggaattggggggATGCCCGCGAGAAATAAGGGCGGTGTTCATCTGGCGTCGTGTTTAAAATTTATAAAACCCGACGGCGATATCAGGGCAGGGCGCCTTGCTGTCCATCAGTCCCACCAGGGAAATCCCCGCCACTTCAATCATCttt                                                                                                                               <  1:559619/137‑1 (MQ=255)
                      ggCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTa                                                                                                                              >  1:485358/1‑149 (MQ=255)
                       gCAATGTGGCGATGCCCGCGAGCAATATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGt                                                                                                                          >  1:384749/1‑151 (MQ=255)
                        aaagggggggATCCCCGGGAAAAATAAGGGCGGTGGTCTTGGGGGGGGCGGCTACAAATTTATAAAACCCGACGGCGATATCAAGGCATGGCTATTTGTTGTCCATCGGTTCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGt                                                                                                                          <  2:1057720/142‑1 (MQ=255)
                         aaTGGGGCGATGCCCGCAAAAAAAAAGGGCGGTGGTCACCGGGCGTGCTGCTCAAAATTTAAAAAACCCGCCGGCAATAAAAAGGCAGGGCTCTTTGCTGCCCATCAGTCCCCCCCGGGAAATCCCCGCCACTTCAATCATCTTTAAgg                                                                                                                           <  2:469651/149‑1 (MQ=255)
                           gggggCGATGCCGGCGAGAAATAGTGGAGGTGGTCATCGGGGTGCCTGCTAAAAAGTTACAAAACCCGACGGCGATATCAGGCCGTGGCGTCTTGCTGTCCATCATTCCCCCCAGGGAAACCCCCGCCACTTCAATCATCTTTAAGGTgc                                                                                                                        <  2:404166/147‑1 (MQ=255)
                              ggggATGCCCGGGAGAAATAAGGGCGGTGGTCATCTGGCGTTTTGCTCAAAAGTTAAAAAACCCGACGGCGATATCAAGGCAGGGCGCCTTCCTGTCCACCAGTCCCCCCACGGAAATACCCGCCACTTCAATCATCTTTAAGGTgc                                                                                                                        <  1:756390/144‑1 (MQ=255)
                                           aGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGt                                                                                                       >  1:253293/1‑151 (MQ=255)
                                               ataatGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCgg                                                                                                   >  1:451009/1‑151 (MQ=255)
                                                 aatGGCGGTGGTCTTCGGGCGTCCTCTTTAAAATTAACAAAACCCGCCGGCAATACAAAGGCAGGGCTCTTTGCTGCCCATCAGTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCgg                                                                                                   <  2:708115/149‑1 (MQ=255)
                                                  atGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTg                                                                                                 >  1:148000/1‑150 (MQ=255)
                                                  agggCGGTGGTCATCTGGTTTGCTTATTTAAGGTAACAAACGCCCGCGGCGATAACAAGGCAGGGCGTCTTGCTGTCCACCAGTTCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGt                                                                                                  <  1:726612/147‑1 (MQ=255)
                                                     ggaGTGGTTCATCTGGGGTGCTGTTCAAAATTTACAAAACCCGCCGGCGATATAAAGGCGTGGCGCTTTGCTGTCCATCAGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAg                                                                                               <  2:1007465/148‑1 (MQ=255)
                                                        gTGGTCTTCGGGCGTGCTCCTCAAAATTTAACAAACCCGACGAGAATACCAGGGCATGGCGTCTTGTTCCCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa                                                                                          <  2:607016/151‑1 (MQ=255)
                                                           gttcttCTGGCGTGCTGCTTAAAATTTATAAAACCCGCCGGCAATATAAAGGCATGGCGCCTTGTTGCCCACCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAg                                                                                         <  1:653334/145‑1 (MQ=255)
                                                           ggtcttggggCGTGCTGTTTAAAAGTTATAAAACCCGAGGGCAATACCAAGGCAGGGCGACTTGCTGCCCACCAGTCCCCCCGCGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                           <  1:429843/140‑1 (MQ=255)
                                                                tcgggggTGCTGTTTAAAATTTAAAAAACCCGGGGGCGATAACAAGGCAGGGCGATTTGCTTCCCATCAGTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGt                                                                                    <  1:464592/144‑1 (MQ=255)
                                                                   ggggTGCTGCTCAAATTTTAAAAAACCCGGCGGCGAAACAAGGGCAGGGCGTTTTCCTGTCCACCAGTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCg                                                                                     <  1:13839/143‑1 (MQ=255)
                                                                   ggCGTGCTTCTAAAAATTTACAAACCCCGAGGGCAATATCAGGGCATGGCGTCTTGCTGTCCATCATTCCCCCCCGGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                           <  1:742458/140‑1 (MQ=255)
                                                                   ggCGTGCTGCTTAAAATTAAAAAAACCCGCCGGCGATAACAAGGCATGGCGCTTTTCTGTCCATCGGTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTc                                                                                   <  1:251141/148‑1 (MQ=255)
                                                                      gTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATAtt                                                                             >  1:424612/1‑151 (MQ=255)
                                                                        gctgctTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTcc                                                                           >  1:652089/1‑151 (MQ=255)
                                                                           tgctaaaaaTTTACCAAACCCGACGGCGATATCAAGGCAGGGCGTCTTGCTGTCCATCAGTCCCCCCACGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTcc                                                                           <  1:277430/143‑1 (MQ=255)
                                                                            gcataaaattaaacAAACCCGAAGGCGAAATCAGGGCAGGGCGGTTTGCTCACCTTAGTTCCCCCCACGGAAAGCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCa                                                                          <  2:1109013/136‑1 (MQ=255)
                                                                              ttCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                                        >  1:574857/1‑149 (MQ=255)
                                                                                cAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGgag                                                                      >  2:1080674/1‑149 (MQ=255)
                                                                                 aaatttaaaaaacCCGAGGGCAATATCAAGGCTGGGCGCCTGGCTGTCCACCATTCCCCCCCGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGga                                                                       <  1:398417/138‑1 (MQ=255)
                                                                                 aaatttaaaaaacCCGACGGCGATATCAGGGCAGGGCGTTTTGCTGTCCATCAGTCCCACCAGGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTc                                                                            <  1:50606/133‑1 (MQ=255)
                                                                                      tATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATgg                                                               >  1:443155/1‑150 (MQ=255)
                                                                                       aacAAAACCGCCGGCAATATCAAGGCAGGGCGATTTGCGTCCTATCAGTCCCGCGAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATgg                                                               <  2:630797/147‑1 (MQ=255)
                                                                                                 gACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGt                                                    >  2:809026/1‑151 (MQ=255)
                                                                                                    ggCGATATCAAGGCATGGCGACTTGCTGTCTATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCt                                                 >  2:570392/1‑151 (MQ=255)
                                                                                                    ggCGATATCAAGGCATGGCGACTTGCTGCCCATCAGTCCCCCCGCGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCt                                                 <  1:993670/151‑1 (MQ=255)
                                                                                                     gCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCtt                                                >  2:159523/1‑151 (MQ=255)
                                                                                                       gATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGc                                              >  2:100875/1‑151 (MQ=255)
                                                                                                           tCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTcc                                           >  2:122503/1‑150 (MQ=255)
                                                                                                            cAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAt                                         >  1:17959/1‑151 (MQ=255)
                                                                                                             aaGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAtt                                        >  1:323218/1‑151 (MQ=255)
                                                                                                               ggcggggCGTCTTGCTGCCCATCAGTTCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAttt                                       <  1:548127/145‑1 (MQ=255)
                                                                                                                    ggCGTCTTCCTCTCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGt                                  <  2:452712/150‑1 (MQ=255)
                                                                                                                    ggCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGtt                                 >  2:353421/1‑151 (MQ=255)
                                                                                                                     gcccttctctccccATCATTCCCCCATCGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGt                                  <  1:771552/137‑1 (MQ=255)
                                                                                                                      ccttttGCTGCCCATCAGTCCCCCCCCGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGttt                                <  1:920442/146‑1 (MQ=255)
                                                                                                                               gTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGt                        >  2:1131524/1‑149 (MQ=255)
                                                                                                                                            ccaccaCGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCAt           >  1:64793/1‑151 (MQ=255)
                                                                                                                                              ccccgggggAACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCa               <  2:4528/140‑1 (MQ=255)
                                                                                                                                                   cGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTg     >  1:1133185/1‑151 (MQ=255)
                                                                                                                                                   cGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTg     >  1:1089166/1‑151 (MQ=255)
                                                                                                                                                       aaaTCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGAt   <  1:700068/149‑1 (MQ=255)
                                                                                                                                                        aaTACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGATa  >  1:344591/1‑150 (MQ=255)
                                                                                                                                                           |                                                                                                                                                  
CGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGATA  >  NC_000913/3578153‑3578447

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: