Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,960,751 C→T 33.3% R83H (CGT→CAT)  kdgK ← ketodeoxygluconokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,960,7510CT33.3% 36.6 / 26.7 33R83H (CGT→CAT) kdgKketodeoxygluconokinase
Reads supporting (aligned to +/- strand):  ref base C (14/8);  new base T (5/6);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

CGTTCCGCCAGTAGTAGAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGCATCTGCTGG  >  W3110S.gb/3960683‑3960819
                                                                    |                                                                    
cGTTCCGCCAGTAGTAGAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACg                                                                     <  1:3329535/70‑1 (MQ=255)
       ccAGTAGTAGAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGG‑‑GACGGTTTTCCATc                                                           >  1:993765/1‑71 (MQ=255)
           tagtagAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCg                                                         >  1:1282926/1‑71 (MQ=255)
             gtagAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGtt                                                       <  1:581509/71‑1 (MQ=255)
             gtagAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGtt                                                       >  1:2970232/1‑71 (MQ=255)
              tagAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTg                                                      <  1:1329072/71‑1 (MQ=255)
                gAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATcc                                                          <  1:3708730/65‑1 (MQ=255)
                         cTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAg                                           <  1:3494319/71‑1 (MQ=255)
                             ccGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATc                                       <  1:734058/71‑1 (MQ=255)
                               ggTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGTATCaa                                     >  1:1096964/1‑71 (MQ=255)
                                  gCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAg                                           >  1:3329152/1‑62 (MQ=255)
                                   cTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGt                                  <  1:1018481/70‑1 (MQ=255)
                                   cTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGtt                                 <  1:4030783/71‑1 (MQ=255)
                                        ggTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTcc                                                              >  1:3785188/1‑37 (MQ=255)
                                        ggTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAgg                                               >  1:498565/1‑52 (MQ=255)
                                          tttCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAggg                                              >  1:3910925/1‑51 (MQ=255)
                                            tCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGt                         >  1:4152162/1‑70 (MQ=255)
                                            tCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTg                        >  1:498883/1‑71 (MQ=255)
                                            tCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTg                        >  1:4146846/1‑71 (MQ=255)
                                              aaTGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCa                                                 >  1:897007/1‑44 (MQ=255)
                                                  tagtaAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGaa                                            <  1:2848816/45‑1 (MQ=255)
                                                     taAATGCCCGGCAGACGGTTTTCCATCCGTTGGGTc                                                  <  1:2867318/36‑1 (MQ=255)
                                                     taAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAcgg                                              <  1:2381837/40‑1 (MQ=37)
                                                      aaaGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAggg                                              <  1:4153124/39‑1 (MQ=255)
                                                       aaGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTc               >  1:1944577/1‑69 (MQ=255)
                                                          gCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGCa          >  1:3776354/1‑71 (MQ=255)
                                                             cGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGt                                          <  1:243649/36‑1 (MQ=255)
                                                             cGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGt                                          <  1:249222/36‑1 (MQ=255)
                                                              ggCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCg                          >  1:1714166/1‑51 (MQ=255)
                                                               gCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGc                       >  1:3627810/1‑53 (MQ=255)
                                                                cAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGCATctgct    >  1:3488060/1‑71 (MQ=255)
                                                                 aGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGCATctgctg   >  1:3659152/1‑71 (MQ=255)
                                                                  gACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGCATCTGCTgg  >  1:798804/1‑71 (MQ=255)
                                                                    |                                                                    
CGTTCCGCCAGTAGTAGAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGCATCTGCTGG  >  W3110S.gb/3960683‑3960819

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