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

TTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCT  >  NC_000913/3743509‑3743807
                                                                                                                                                      |                                                                                                                                                      
ttCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTtactac                                                                                                                                                       <  2:688013/150‑1 (MQ=255)
ttCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTtactac                                                                                                                                                       <  1:730469/150‑1 (MQ=255)
   acaAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCctg                                                                                                                                                   >  2:15067/1‑151 (MQ=255)
      aaTTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCctgctg                                                                                                                                                <  1:59383/151‑1 (MQ=255)
         ttttATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGc                                                                                                                                              <  1:722801/150‑1 (MQ=255)
           ttATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCg                                                                                                                                            <  2:710649/150‑1 (MQ=255)
            tATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGa                                                                                                                                           <  1:686283/150‑1 (MQ=255)
              ttGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGaaa                                                                                                                                         >  2:462812/1‑150 (MQ=255)
               tGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAc                                                                                                                                       >  1:431363/1‑151 (MQ=255)
                gCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAcc                                                                                                                                      <  2:101213/151‑1 (MQ=255)
                  cATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAccgc                                                                                                                                    >  2:79356/1‑151 (MQ=255)
                   aTTGAAGTGGTCAAGATTATCGATGGTCCCACGCGCGATGCGAAGCTGCAAAGCATCATGGATGACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCTGTTCGCATAACAGGCGATGAAGATACTACCCTGATGGCGGAAAAccgc                                                                                                                                    >  2:214432/1‑150 (MQ=255)
                   aTTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAccgcc                                                                                                                                   >  2:380847/1‑151 (MQ=255)
                      gAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGt                                                                                                                                 >  2:342300/1‑150 (MQ=255)
                                   ttACCGACGTCCCCCCCCGCGACGCAAGGCTAAAACGCATCATGGTTTACTTTCCAAGTCCGGGGCGTGCGGACAAAAACCAGGCCATCCCCTCCCCCGGCAATAAATTTCCTCCCCTGCTGCCCAAAAACCGCCGTAACGGCATCACt                                                                                                                     <  1:335650/149‑1 (MQ=255)
                                      tCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTtg                                                                                                                 >  2:188302/1‑150 (MQ=255)
                                           ggcccccccCGCGATGCAAGGCGGAACCGCATCATGAATTCCTTTCCTAGGGCGGGGGGTGCTGACAAAAACCGGGCCTTCCGCTCCCCCGGCAATAATTTTCCTCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTtgatga                                                                                                             <  1:309819/146‑1 (MQ=255)
                                              ccccccCGCGATGCCAAGCTCCAACGCATCGTGTTTTCCTTTATTAGGGCCGGGCGTCTTGAAAAAAACCAGCCCATCCCCTTCCCCGCCCATAATTTTCCTCCCCTGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAg                                                                                                          <  1:83327/149‑1 (MQ=255)
                                              ccccccCGCGATGCAAACCTGCAACGCATCAGGGATTACTTTACTAGGGCGGGGCGTGCTGAAAAAAACCAGGCCATTCGCTTCCCCGGCCATAAATTTACTCCCCTGCGGGCGAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCg                                                                                                        <  1:233252/151‑1 (MQ=255)
                                              ccccccCGCGAGGCAAAGCGGAAACGCATCATGGTTTACTTTACTAGGGCGGGGCGTGCTAACAAAAACCGGCCTTTCCGCTCCCCCGCCAATAATTTTACTCCCCGCCGGGCGGAAAACCGCCGTAACGGCATCACTGTTGATGACAGc                                                                                                         <  2:225842/150‑1 (MQ=255)
                                                cccccGCGATGCCAAGCTGAAACGCATCTGGGTTTACGTTATGAGTGCCGGGGGTGCTGAAAAAAACCGGCCCATCCGCTTACCCGGCCATAATTTTACTACCCTGGTGGCAGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCgtg                                                                                                      <  2:579191/149‑1 (MQ=255)
                                                 ccccGCGATGCCAAGCTAAAACGCATCATGGATTTCTTTACTAGGGCGGGGGGTGCTGAAAAAAACCGGCCCATTCGCTTCCCCGCCCATAAATTTACTCCCCTGCGGGCGAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCgtgt                                                                                                     <  1:67371/150‑1 (MQ=255)
                                                     gcgATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGc                                                                                                 >  1:130379/1‑151 (MQ=255)
                                                         tGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCaa                                                                                              >  1:235470/1‑150 (MQ=255)
                                                           ccAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAAt                                                                                           >  1:194137/1‑151 (MQ=255)
                                                            cAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATc                                                                                          >  1:387520/1‑151 (MQ=255)
                                                                 tGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGc                                                                                     >  1:437182/1‑151 (MQ=255)
                                                                    accccgccatGTATTCTGTTACTGGGCCGGAGCGTCGTAACAAAAACCGGGCCTTTCCCTTCCCCGGCCATGATTTCATTCCCCCGGTGGCGAAACCCCGCCGTAACGGCATCACTGTTGATGACAGCGTGTggg                                                                                                  <  2:627618/128‑1 (MQ=255)
                                                                    accccataatGTATTACGTTATTGGTGCGGGGGGTCCTGAAAAAAACCGGCCCTTCCGCTTCCCCGGCCTTAATTTTCCCCCCCTGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCa                                                                                        <  1:342302/141‑1 (MQ=255)
                                                                     aCGCATCATGGATTGCTTTACGGGGCCCGAGCGTCCGGACAAAAATCAGGCCTTCCGCTTACCCGCCAATAATTTTACTACCCTGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTa                                                                                 <  2:481028/151‑1 (MQ=255)
                                                                             tgtaTTACTTTACTAGTGCGGGGCGTGCTGACAAAAATCAGGCCATTCGCTTCCCCGGCAATAAATTTCCTCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATgagg                                                                            <  1:25567/145‑1 (MQ=255)
                                                                              ggATTACTTTTCTAGTGCGGGGGGTGCTGACAAAATTCGGGCCTTCCCCTTCCCCGGCAATGATTTTCCCCCCCTGCTGGCGGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGgaga                                                                         <  2:686554/150‑1 (MQ=255)
                                                                                aTTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATa                                                                       >  1:487443/1‑150 (MQ=255)
                                                                                 ttactttactAGGGCCGAGCGTGCTGACAAAACACAGCCCTTCCCCTACCCCGCCCAGAATTTTCCTCCCCGGCTGGCAAAAACCCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATaa                                                                      <  1:177616/150‑1 (MQ=255)
                                                                                    cGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTc                                                                   >  1:126402/1‑150 (MQ=255)
                                                                                       tccTAGTGCCGGGCGTGCTGACAAAAACCAGCCCATTCCCTTCCCCGGCCATAATTTTACTCCCCTGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCAt                                                                 <  1:386012/147‑1 (MQ=255)
                                                                                                 gAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGAc                                                      >  1:627869/1‑150 (MQ=255)
                                                                                                 gAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGAc                                                      >  1:368422/1‑150 (MQ=255)
                                                                                                  aGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACat                                                    >  1:185348/1‑151 (MQ=255)
                                                                                                           aaaaaaaTCAGCCCTTTCCCTTCCCCGGCCATAATTTTACTCCCCTGCTGGCGAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCa                                            <  1:637410/147‑1 (MQ=255)
                                                                                                             gAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGcc                                         >  1:656929/1‑151 (MQ=255)
                                                                                                              aaaaCCAGCCCTTTCGTTTCCCCGGCCATGATTTTCCTCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGcc                                         <  2:575130/150‑1 (MQ=255)
                                                                                                                  tCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATc                                    >  2:58988/1‑151 (MQ=255)
                                                                                                                     gcccATTCGCTCCCCCGGCAATAAATTTACTCCCCGGCGGCCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAAt                                     <  1:607753/145‑1 (MQ=255)
                                                                                                                        ccttccGCTTCCCGGGCAATAATTTTACTCCCCGGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGc                               <  2:433341/146‑1 (MQ=255)
                                                                                                                               cTTCCCGGGCAATAAATTTCCTCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTg                         <  1:567466/150‑1 (MQ=255)
                                                                                                                               cTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGc                        >  2:235323/1‑151 (MQ=255)
                                                                                                                                   cccGGCCTTGTATTTACCCCCCTGCTGCCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGc                        <  1:690174/147‑1 (MQ=255)
                                                                                                                                   cccGGCAATGAATTTCCTCCCCTGCTGGCGAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGccc                      <  1:528915/149‑1 (MQ=255)
                                                                                                                                   cccGGCAATAATTTTCCCCCCCTGCTGGCAAAAAGCCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGttt                          <  2:711293/145‑1 (MQ=255)
                                                                                                                                   cccGCCAATAAATTTCCTCCCCGGGTGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGttt                          <  2:720910/145‑1 (MQ=255)
                                                                                                                                        ccATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCgccgcc                >  1:722558/1‑150 (MQ=255)
                                                                                                                                           agaattttCCTCCCCGGGTGGCCGAAACCCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCgccgcc                <  1:354791/143‑1 (MQ=255)
                                                                                                                                                  tactacCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAgg      >  1:603495/1‑151 (MQ=255)
                                                                                                                                                    ctacCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGgcg    >  2:165727/1‑151 (MQ=255)
                                                                                                                                                    ctacCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGgcg    >  1:490712/1‑151 (MQ=255)
                                                                                                                                                      acCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCt  >  2:555149/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
TTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCT  >  NC_000913/3743509‑3743807

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