Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,743,6690GA51.1% 30.0 / 66.7 47E313K (GAA→AAA) yiaK2,3‑diketo‑L‑gulonate reductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base G (20/3);  new base A (0/24);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 1.81e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.67e-08

TTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTT  >  NC_000913/3743519‑3743814
                                                                                                                                                      |                                                                                                                                                 
tttATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCg                                                                                                                                                   >  1:492858/1‑151 (MQ=255)
 ttATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCg                                                                                                                                                   >  2:459638/1‑150 (MQ=255)
   aTTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGaa                                                                                                                                                 <  2:463507/150‑1 (MQ=255)
   aTTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGaa                                                                                                                                                 >  1:624168/1‑150 (MQ=255)
      gCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAcc                                                                                                                                             >  2:807892/1‑151 (MQ=255)
         aTTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAccgcc                                                                                                                                          <  2:426492/151‑1 (MQ=255)
                 ggACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTATGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGc                                                                                                                                  >  2:679165/1‑151 (MQ=255)
                  gACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGc                                                                                                                                  >  2:958167/1‑150 (MQ=255)
                         ttATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGt                                                                                                                          >  2:158378/1‑151 (MQ=255)
                         tatccGACGCTCCCCCCCGCGATGCCAAGCTCCAACCCATCATGGTTTACGTTACTAGTGCGGGGCGTGCTGACAAAAATCAGGCCATTCGCTTCCCCGCCCATGAATTTCCTCCCCTGCGGGCAAAAAACCGCCGTAACGGCATCACt                                                                                                                            <  1:705375/145‑1 (MQ=255)
                          tATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGtt                                                                                                                         >  2:870914/1‑151 (MQ=255)
                             cgccGGCCCCCCCCGCGATGCCAAGCTGCAACGCATCATGTATTCCTTTACTGGTGCCGGGCGTGCTGACAAAAATCGGCCCTTTCGCTTCCCCGGCCATGATTTTACTCCCCTGCGGGCGAAAAACCGCCGTAACGGCATCACTGTtga                                                                                                                       <  1:749891/147‑1 (MQ=255)
                             cGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTtgat                                                                                                                      >  2:268035/1‑151 (MQ=255)
                              gacgtcccccccCGCGAGGCCAAGCTGCAACGCATCATGGTTTACTTTACTGGTGCGGAGCGTGCTGACAAAATTCAGGCCTTTCGCTTCCCCGGCCATAAATTTACTCCCCTGCGGCCCAAAAACCGCCGTAACGGCATCACTGt                                                                                                                          <  1:933012/140‑1 (MQ=255)
                                 ggcccccccCGCGATGCCAAGCTGCACCGCATCATGGATTACTTTACTAGTGCCGGGCGTGCTGACGAAAATCAGGCCATCCCCTTCCCCGGCAATGAATTTACTCCCCGGCTGCCCAAAAACCGCCGTAACGGCATCACTGTtgatga                                                                                                                    <  1:169493/146‑1 (MQ=255)
                                    ccccccCGCGATGCCAAGCTGCACCGAATAATGGATTACGTTACTAGGGCCGGGCGTGCTGACAAAAATCAGGCCATTCGCTTCCCCGGCAATGAATTTACTCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACa                                                                                                                  <  1:220177/148‑1 (MQ=255)
                                     ccccccGCGATGCAAAGCTGAACCGCATCTTGGATTCCGTTACTAGTGCCGGGCGTGCTGACAAAAACCAGGCCATCCGCTTCCCCGGCCATGATTTTCCTCCCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCgt                                                                                                              <  1:86996/148‑1 (MQ=255)
                                             gATCCCAAGCTGCAACCCTCCATGGTTTACGTTACTTGTGCCGAGCTTGCTGACAAAAACCAGGCCATTCGCTTCCCCGGCCATGAATTTCCTCCCCTCCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGc                                                                                                        <  1:817309/149‑1 (MQ=255)
                                                      cTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAgg                                                                                             >  2:540981/1‑151 (MQ=255)
                                                       tGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGc                                                                                            >  2:951992/1‑151 (MQ=255)
                                                          accGCATCATGGATTACTTTACTAGTGCCGGGGGTGCTGACAAAAACCAGGCCTTCCGCTTCCCCGCCAATAATTTTACTCCCCTGCTGGCGAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGtt                                                                                         <  2:569744/149‑1 (MQ=255)
                                                          aacacctcatGGATTACTTTACTGGTGCCGGGGGTGCTGACAAAAATCGGCCCCTCCGTTTCCCCGGCCATGTTTTTCCTCCCCTGTTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAgg                                                                                             <  2:926113/141‑1 (MQ=255)
                                                          aaCGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGtt                                                                                         >  2:85059/1‑151 (MQ=255)
                                                              catcatGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGCTGATGACAGCGTGTGGGCCAACATCCAGGCGTTAt                                                                                       >  1:794134/1‑149 (MQ=255)
                                                                           gTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCAt                                                                        >  2:280162/1‑151 (MQ=255)
                                                                            ttccTAGTGCCGGGCGTGCTGACAAAAACCAGCCCATTCCCTTCCCCGGCAATAAATTTCCTCCCCTGCTGCCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCa                                                                         <  2:204969/146‑1 (MQ=255)
                                                                               cTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAtt                                                                    >  2:483563/1‑151 (MQ=255)
                                                                                  gtccggggCGTGCTGAAAAAAAGCAGCCCTTCCGCTTCCCCGGCATTGATTTTACTCCCCTGTGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAtt                                                                    <  2:798314/141‑1 (MQ=255)
                                                                                     cggggCGTGCGGACAAAAACCGGGCCTTTCGCTTCCCCGGCCATAATTTTACTCCCCTGCTGGCGAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTg                                                                <  1:536727/145‑1 (MQ=255)
                                                                                     cggAGCGTGCTAACGAAAACCAGCCCATTCGTTTCCCCGGCAATAATTTTCCGCCCCTCCTGCCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATG‑GGAGATAAGTCATGATTTTTgg                                                               <  1:82154/147‑1 (MQ=255)
                                                                                     ccgggggTGCTGAAAAAAATCAGCCCATCCGCTTCCCCGCCCAGAAATTTACTACCCGGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTgg                                                               <  1:952360/144‑1 (MQ=255)
                                                                                     ccGGGCGTGCTAACAAAAACCGGCCCATCCCCTTCCCCGCCCATAATTTTACTCCCCTGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAt                                                                     <  1:897148/144‑1 (MQ=255)
                                                                                         gCGTGCGGACAAAAATCAGGCCTTTCCCTTCCCCGCCCATGAATTTACTCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACata                                                          <  1:967462/151‑1 (MQ=255)
                                                                                               tGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCgcgc                                                    >  1:35132/1‑151 (MQ=255)
                                                                                                    aaaaCCAGGCCTTTCCCTTCCCCGGCCATGATTTTACTCCCCGGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCgcgc                                                    <  2:623307/146‑1 (MQ=255)
                                                                                                    aaaaCCAGCCCTTCCTCTTCCCCGGCAATAATTTTACTCCCCGGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACa                                                            <  2:842202/138‑1 (MQ=255)
                                                                                                     aaaccggcccATTCGTTTCCCCGCCCATAATTTTCCTCCCCTGCTGCCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCa                                              <  1:908355/143‑1 (MQ=255)
                                                                                                           ggCCTTTCCCTTCCCCGGCAATAATTTTACTCCCCGGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAAt                                            <  2:186626/147‑1 (MQ=255)
                                                                                                               aTTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACAGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCg                                    >  1:223494/1‑151 (MQ=255)
                                                                                                                  cgtttccccGCCCATAATTTTCCTCCCCTGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGttt                                 <  2:954849/145‑1 (MQ=255)
                                                                                                                  cGCTTCCCCGGCCATAAATTTACTCCCCTGCCGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGttt                                 <  2:937638/151‑1 (MQ=255)
                                                                                                                        cccGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCccgc                           >  1:157918/1‑151 (MQ=255)
                                                                                                                             ccATAATTTTACTCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCa                      <  1:483351/151‑1 (MQ=255)
                                                                                                                                  aaTTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGa                  >  1:318135/1‑150 (MQ=255)
                                                                                                                                         tacCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGgcgc          >  2:171166/1‑151 (MQ=255)
                                                                                                                                           cccTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTc        >  2:781430/1‑151 (MQ=255)
                                                                                                                                           cccTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGgc            <  1:952115/147‑1 (MQ=255)
                                                                                                                                                ctgCCCAAAACCCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGa      <  2:868356/148‑1 (MQ=255)
                                                                                                                                                 tgGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGAtttt  >  2:497776/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
TTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTT  >  NC_000913/3743519‑3743814

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