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

TTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCT  >  NC_000913/3743508‑3743807
                                                                                                                                                     |                                                                                                                                                      
tttCACAAATTTTTATTGCCATTGAATTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTtactac                                                                                                                                                       >  1:403299/1‑151 (MQ=255)
        aTTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTgg                                                                                                                                               <  1:29006/151‑1 (MQ=255)
        aTTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTgg                                                                                                                                               >  1:750569/1‑151 (MQ=255)
          ttttATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGcc                                                                                                                                             >  1:326789/1‑151 (MQ=255)
            ttATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCg                                                                                                                                            <  2:227941/150‑1 (MQ=255)
             tATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGa                                                                                                                                           <  1:879152/150‑1 (MQ=255)
              aTTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGaaa                                                                                                                                         <  2:222613/151‑1 (MQ=255)
                  ccATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAccg                                                                                                                                     >  2:342770/1‑151 (MQ=255)
                     ttGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAccgccg                                                                                                                                  >  2:458450/1‑151 (MQ=255)
                       gAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCTGAAAACCGCCGt                                                                                                                                 >  1:86669/1‑150 (MQ=255)
                                        cGACGGCCCCCCCCGGGATGCAAAGCTACAACCCATCATGGTTTATGTTCTTGGTGCGGGGGGTGTTGACAAAAATCAGCCTTTTCGCTTCCCCGGCCATAAATTTACTCCCCTGCGGGCAAAAAACCGCCGTAACGGCATCACTGTtg                                                                                                                 <  2:661971/149‑1 (MQ=255)
                                            ggcccccccGGCGTTGCCAGCCTCCACCCCATCAGGAATTACTTTCTTGGTGCCGAGCGTCCTGACAAAAACCAGCCCTTCCGTTTCCCCGACCTAGATTTTACTCCCCTGCGGCCCAAAAACCGCCGTAACGGCATCACTGTTGATGACa                                                                                                           <  1:452725/148‑1 (MQ=255)
                                               ccccccCGCGATGCCAAGCTGCAACGCATCATGGTTTCCGTTACTAGTGCCGGGCGTGCTGACAAAAATCAGGCCATTCGCTTCCCCGGCCATGAATTTACTCCCCTGCGGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCg                                                                                                        <  2:302252/151‑1 (MQ=255)
                                                ccACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCgt                                                                                                       >  2:801275/1‑151 (MQ=255)
                                                 cACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCgtg                                                                                                      >  2:464711/1‑151 (MQ=255)
                                                  aCCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCgtg                                                                                                      >  1:877351/1‑150 (MQ=255)
                                                             cAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATc                                                                                          >  1:626251/1‑151 (MQ=255)
                                                              aaGCTACAACGCATCAGGGTTTCCTTTACTGGTGCGGGGCGTGCTGACAAAAACCGGGCCTTCCGCTTCCCCGGCCTTGATTTTACTCCCGTGCGGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATc                                                                                          <  1:423739/150‑1 (MQ=255)
                                                               aGCTGCAACGCATCAGGGATTACTTTCCTGGGGCCGGGCGTGCTGACAAAAACCAGGCCATCCGCTTCCCCGGCCATGATTTTACTCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCaaaa                                                                                            <  2:599129/147‑1 (MQ=255)
                                                                        gCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATg                                                                               >  2:170742/1‑151 (MQ=255)
                                                                          atcatGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATg                                                                               >  1:552294/1‑149 (MQ=255)
                                                                          atcatGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATg                                                                               >  1:755034/1‑149 (MQ=255)
                                                                              tgtaTTACTTTACTGGCCCCGAGCGTGCTAACAAAAATCACCCCTTCCGCTTCCCCGCCCATAATTTTACTCCCCTGCTGGCCAAACACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATg                                                                               <  2:61875/142‑1 (MQ=255)
                                                                               ggATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATgagg                                                                            >  1:529481/1‑147 (MQ=255)
                                                                                  ttcctttcctggTGCCGGGCGTGCTAAAAAAAACCGGGCCATCCGCTCCCCCGGCAATAATTTTCCTCCCCGGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATa                                                                       <  1:135926/138‑1 (MQ=255)
                                                                                       ttACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATg                                                                >  2:534683/1‑151 (MQ=255)
                                                                                         acggttgcggggCGTGCTGACAAAAAACGGCCCCTTCGCTCCCCCGGCCATAAATTTCCTCCCGTGCGGGCAGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAt                                                              <  2:185476/140‑1 (MQ=255)
                                                                                                cggggCGTCCTGAAAAAAACCAGCCCTTTCGCTTCCCCGGCAATAATTTTCCTCCCCTGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAtttt                                                           <  2:414608/143‑1 (MQ=255)
                                                                                                   aGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACat                                                    >  1:668168/1‑151 (MQ=255)
                                                                                                            aaaaaaaTCGGCCCTTTCCCTTCCCGGGCAATAAATTTCCTCCCCGGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAg                                           <  1:751167/148‑1 (MQ=255)
                                                                                                              aaaaaTCAGGCCTTCCGCTTCCCGGGCTATAATTTTCCTCCCCTGCTGCCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAt                                                              <  2:316004/129‑1 (MQ=255)
                                                                                                               aaaaccggccctttCGTTTCCCCGGCAATAAATTTACTCCCCGGCTGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATcg                                                <  1:746811/131‑1 (MQ=255)
                                                                                                                aactctggtatttCGCTTCCGCGGTCATTATTTTATTCCCTTCCTGGCTAAAACCCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCa                                            <  1:753943/135‑1 (MQ=255)
                                                                                                                          aTTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCg                             >  2:489578/1‑151 (MQ=255)
                                                                                                                           ttCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGt                            >  1:474913/1‑151 (MQ=255)
                                                                                                                                tcccccGCCCATGTTTTTCCTCCCCTGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGttt                          <  1:348066/145‑1 (MQ=255)
                                                                                                                                  ccccGGCCTTAATTTTACTCCCTAGCCGCCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGt                                 <  1:612361/138‑1 (MQ=255)
                                                                                                                                  ccccGGCCATGAATTTACTCCCCTGGGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCccg                     <  2:750662/150‑1 (MQ=255)
                                                                                                                                  aCCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGccc                      >  2:748026/1‑150 (MQ=255)
                                                                                                                                   cccGGCAATAATTTTACTCCCCTGCGGGCAGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCccgc                    <  2:877469/151‑1 (MQ=255)
                                                                                                                                          attagtttgcgacCCTGGTGGGCGAAAAGCGCCGAAACCGGATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCAtt             >  1:276156/12‑151 (MQ=255)
                                                                                                                                                ttACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGaaa         >  1:244529/1‑149 (MQ=255)
                                                                                                                                                     acCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGAGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCt  >  2:93128/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCT  >  NC_000913/3743508‑3743807

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