Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,743,6570AC45.8% 68.4 / 18.6 49T309P (ACC→CCC) yiaK2,3‑diketo‑L‑gulonate reductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (23/3);  new base C (0/22);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 8.48e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAA  >  NC_000913/3743507‑3743800
                                                                                                                                                           |                                                                                                                                                
aTTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTtacta                                                                                                                                                  <  2:701164/151‑1 (MQ=255)
 tttCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTtactac                                                                                                                                                 >  2:688216/1‑151 (MQ=255)
        aaTTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCctgctg                                                                                                                                         <  2:432153/151‑1 (MQ=255)
            tttATTGCCATTGACGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGCGCTGACGAAAATCAGGCCATTTGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCg                                                                                                                                     >  2:658299/1‑151 (MQ=255)
                  gCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAc                                                                                                                                <  1:905960/150‑1 (MQ=255)
                    cATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAccgc                                                                                                                             >  1:71571/1‑151 (MQ=255)
                     aTTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAccgcc                                                                                                                            <  2:41655/151‑1 (MQ=255)
                     aTTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAccgcc                                                                                                                            >  1:315445/1‑151 (MQ=255)
                                     ttACCGACGCCCCCTCCCGTGATCCCAAGCTGAAACACATCATGGTTTCCTTTACCAGCGCCGAGCGTCCTGACAAAAACCAGCCCATTCGCTTCCCCGGCCAAGAATTTCCTCCCCGGTTGGCCGAAAACCGCCGTAACGGCATCACt                                                                                                              <  1:657486/149‑1 (MQ=255)
                                        tCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTtg                                                                                                          >  1:314626/1‑150 (MQ=255)
                                            cGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGAc                                                                                                     >  1:636292/1‑151 (MQ=255)
                                               cccccccCGGGATCCCAGGCTGAAACGCATCATGGATTTCGTACCTGGCGCCGGGGGTGCGGCAAAAACCCAGCCCATCCGCTTCCCCGGCAATGATTTTCCTCCCCGGGGGCCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGc                                                                                                  <  2:144329/150‑1 (MQ=255)
                                                ccccccCGGGATCCCAAGCTGAACCGCATCAGGGATTACGTTACTAGGGCGGGGGGTGCTGACAAAAACCAGGCCTTCCGCTTCCCCGCCCATAATTTTACTCCCCGGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCg                                                                                                 <  1:297931/151‑1 (MQ=255)
                                                ccccccCGCGATGCAAGGCTCCAACGCATAAGGGATTATTTTACTAGTGCCGAGCGTGATGACAAAAACCAGGCCTTCCTCTTCCCGGGCTATAATTTACCTCCCCTGCGGGCAAAAACCCGCCGTAACGGCATCACTGTTGATGACAg                                                                                                   <  2:352997/149‑1 (MQ=255)
                                                    cccGGGATGCAACCTGCAACGCCTCATTGTATACGTTTCCTGGCCCCGGGGGTGCTACCAAAAACCCGCCAGTCCGCTTCCCCGTCTTTAATTTTCCTCCCCGGCTGGCAAAAACCCGCCGTAACGGCATCACTGTTGATGACAGCGtgtg                                                                                             <  2:276511/151‑1 (MQ=255)
                                                        cgATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGcc                                                                                         >  1:753689/1‑151 (MQ=255)
                                                         gaggCCAAGCTGCAACGTATCGTGGTTTCCTTTACTAGTGCCGAGGGTGCTAAAAAAAATCAGGCCATTCCTTTCCCCGCCCATGAATTTACTCCCCTGCTGCCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTggg                                                                                           <  1:508551/145‑1 (MQ=255)
                                                               aaGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATcc                                                                                  >  1:813523/1‑151 (MQ=255)
                                                               aaGCGGCACCGCATCAGGGTTTTTTTTACTGGGGCCGAGCGTGCTGACAAAAATCGGCCCTTTCGTTTCCCGGGCAATAATTTTACCCCCCGGCGGCCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCaa                                                                                       <  1:368887/146‑1 (MQ=255)
                                                                   tgcacccaattatGGATTACTTTCCTGGTCCCGGGGGTGCTGACAAAAACCAGCCCTTTCCCTTCCCCGGCCATGATTTTCCTCCCCTGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCa                                                                                 <  1:621116/137‑1 (MQ=255)
                                                                    gCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCg                                                                             >  2:847520/1‑151 (MQ=255)
                                                                      accgcgTCATGTATTACTTCACGAGTGGCGGGGGTGCTGAAGAAATGCAGGCTATCCCTTACCGCGGCCATGATTTGCCGCCCTGGCTGGCGAAAAGCCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGtt                                                                           <  2:45077/149‑1 (MQ=255)
                                                                       aCGCATCATGGATTCCTTTACTAGTGCCGGGCGTGCTGACAAAAACCAGGCCTTTCCCTTCCCCGGCCATGATTTTACTCCCCTGCGGCCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTa                                                                          <  1:34493/151‑1 (MQ=255)
                                                                       aCGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGtt                                                                           >  2:737409/1‑150 (MQ=255)
                                                                           atcatGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATg                                                                        >  1:575681/1‑149 (MQ=255)
                                                                            tcatGGATTCCTTTCCTAGTGCCGAGCGTGCTGACAAAAACCAGGCCATTCGCTCCCCCGGCCATGAATTTACTCCCCTGCTGGCCAAAACCCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATga                                                                       <  1:349081/149‑1 (MQ=255)
                                                                                ggATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATgag                                                                      >  1:687475/1‑146 (MQ=255)
                                                                                 gATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGAt                                                                 >  1:275774/1‑150 (MQ=255)
                                                                                 gATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGAt                                                                 >  2:306607/1‑150 (MQ=255)
                                                                                         ttctcgtcgCCGAGCGTCCCAACGAAATGCAGCCCTTCCTTTCACCCGCCCATAAATTTCCTCCCCGGAAGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTc                                                            <  1:808735/140‑1 (MQ=255)
                                                                                            cTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAtt                                                      >  2:977611/1‑151 (MQ=255)
                                                                                                  cggggCGTGCGGAAAAAAACCGGGCCTTCCGCTCCCCCGCCCATGATTTTACTCCCCGGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAtt                                                      <  1:754538/141‑1 (MQ=255)
                                                                                                    gAGCGTGCTGACGAAAACCAGGCCATTCGCTCCCCCGGCCATAATTTTACCCCCCGGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGa                                                <  1:176647/149‑1 (MQ=255)
                                                                                                    gAGCGGGCTGCCAAAACGCGGGCCATCCTCCTCCCCGCCAAAGATTTTGCTCCCCGGCTCGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATg                                                         <  2:412334/140‑1 (MQ=255)
                                                                                                      gCGTGCTAACAAAAACCAGCCCTTTCGCTTCCCCGCCAATAATTTTACCCCCCGGCGGCCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACata                                            <  2:286471/151‑1 (MQ=255)
                                                                                                        gtgcaaaaaaaaaCCGGGCCATCCCCTTCCCCGCCCAAAATTTTCCTCCCCTGCGGGACAAAACCCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAttt                                                     <  2:880966/131‑1 (MQ=255)
                                                                                                          gCTGACAAAAATCAGGCCTTCCCTTTCCCCGGCAATGAATTTACTCCCCTGTTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATcg                                         <  1:847146/150‑1 (MQ=255)
                                                                                                               cGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGc                                   >  1:782951/1‑151 (MQ=255)
                                                                                                                 aaaaGCAGGCAATTCGTTCCCCCGCCAATAATTTCACTCCCTTGCTGCCGAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACata                                            <  2:434808/140‑1 (MQ=255)
                                                                                                                 aaaaACAGCCCATCCCCTTCCCCGGAATTTATTTTCCTCCCCTCCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCa                                     <  2:882248/147‑1 (MQ=255)
                                                                                                                  aaaCCAGCCCATCCGCTTCCCCGGCAATAATTTTCCCCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCa                                <  2:416539/151‑1 (MQ=255)
                                                                                                                      cAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATcc                            >  2:718684/1‑151 (MQ=255)
                                                                                                                               tCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGtt                    >  2:1040587/1‑151 (MQ=255)
                                                                                                                                ctctcccccGCCCATGATTTTACTCCCCTGGGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGtt                    <  1:12786/144‑1 (MQ=255)
                                                                                                                                 gCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTg                  >  1:1019003/1‑151 (MQ=255)
                                                                                                                                  ctcccccGGCCTTAAATTTACCCCCGTGGTGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGcc                       <  2:739363/141‑1 (MQ=255)
                                                                                                                                  cTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGc                 >  2:1071039/1‑151 (MQ=255)
                                                                                                                                  cTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAAAGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGAAGCCCAATCTGTGCCGTTTg                  >  2:143275/1‑150 (MQ=255)
                                                                                                                                     aCCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGcc                >  2:562704/1‑149 (MQ=255)
                                                                                                                                         cggcccgggttttgactccCCTGCGGCCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCccgcc            <  2:743371/137‑1 (MQ=255)
                                                                                                                                                   tttACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGTCCGCCGCCATTGaaa  >  1:806183/1‑150 (MQ=255)
                                                                                                                                                    ttACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGaaa  >  1:364143/1‑149 (MQ=255)
                                                                                                                                                           |                                                                                                                                                
ATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAA  >  NC_000913/3743507‑3743800

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