Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,093,241 C→A 33.8% G207C (GGT→TGT)  gntT ← gluconate transporter, high‑affinity GNT I system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,093,2410CA33.8% 113.9 / 43.9 62G207C (GGT→TGT) gntTgluconate transporter, high‑affinity GNT I system
Reads supporting (aligned to +/- strand):  ref base C (27/14);  new base A (12/9);  total (39/23)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

AGGTCCAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTCGG  >  W3110S.gb/4093175‑4093310
                                                                  |                                                                     
agGTCCAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTcc                                                                   >  1:4619836/1‑71 (MQ=255)
agGTCCAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTcc                                                                   >  1:4840075/1‑71 (MQ=255)
  gTCCAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCgg                                                                 <  1:3914487/71‑1 (MQ=255)
  gTCCAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCgg                                                                 <  1:1703834/71‑1 (MQ=255)
     cAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGaa                                                               <  1:3425335/70‑1 (MQ=255)
     cAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGaa                                                               <  1:2101210/70‑1 (MQ=255)
     cAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAAt                                                              >  1:3097561/1‑71 (MQ=255)
     cAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAAt                                                              >  1:4627173/1‑71 (MQ=255)
      aGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAAtt                                                             <  1:2751576/71‑1 (MQ=255)
      aGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAAtt                                                             <  1:4680670/71‑1 (MQ=255)
       gACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTg                                                            <  1:3934119/71‑1 (MQ=255)
       gACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAATTg                                                            >  1:2931331/1‑71 (MQ=255)
           cTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCtt                                                        <  1:2999338/71‑1 (MQ=255)
           cTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCt                                                         <  1:670034/70‑1 (MQ=255)
             gACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTc                                                                    >  1:204591/1‑57 (MQ=255)
             gACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTc                                                                    >  1:4986590/1‑57 (MQ=255)
             gACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAATTGGCTTAt                                                      >  1:2604012/1‑71 (MQ=255)
              aCGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCTGGAATTGGCTTATc                                                     <  1:5121220/71‑1 (MQ=255)
                gCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTc                                                                    >  1:905586/1‑54 (MQ=255)
                   aaaGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAAt                                                              >  1:5134244/1‑57 (MQ=255)
                    aaGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTcc                                                                   >  1:5049005/1‑51 (MQ=255)
                     aGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAATTGGCTTATCGATACCt                                              >  1:4179592/1‑71 (MQ=255)
                      gCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAATTGGCTTATCGATACCtt                                             >  1:3338791/1‑71 (MQ=255)
                          ggCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCaa                                         <  1:22229/71‑1 (MQ=255)
                             aTCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAacac                                      >  1:4470631/1‑71 (MQ=255)
                              tctcTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCt                                                         <  1:4321719/51‑1 (MQ=255)
                              tctcTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACAcg                                     >  1:2280444/1‑71 (MQ=255)
                              tctcTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACAcg                                     >  1:296778/1‑71 (MQ=255)
                                tcttctTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTcc                                                                   >  1:3380668/1‑39 (MQ=255)
                                tcttctTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAAt                                                              >  1:3447776/1‑44 (MQ=255)
                                tcttctTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACAcgc                                    >  1:1391053/1‑70 (MQ=255)
                                tcttctACGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACAcgc                                    >  1:4904/1‑70 (MQ=255)
                                 cttcttCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGa                                                                <  1:3801310/41‑1 (MQ=255)
                                 cttcttCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAATTGGCTTATCGATACCTTTCAACAcgc                                    >  1:950647/1‑69 (MQ=255)
                                     ttcGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGa                                                   >  1:3502966/1‑50 (MQ=255)
                                     ttcGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGa                                                   >  1:1993464/1‑50 (MQ=255)
                                     ttcGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAATTg                                                            >  1:1752095/1‑41 (MQ=255)
                                     ttcGCTGAAGGTTTTCGCGCTGTAGAGAACTTCCGGAATTGGCTTATCGa                                                   >  1:3465587/1‑50 (MQ=255)
                                       cGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCg                               <  1:3851316/68‑1 (MQ=255)
                                        gCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTaa                            <  1:2004227/70‑1 (MQ=255)
                                            aaGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAaccgg                       <  1:3058780/71‑1 (MQ=255)
                                              ggTTTTCGCGCTGTAGAGACATTCCGGAATTGGCTTATCGAt                                                  >  1:4955531/1‑42 (MQ=255)
                                               gTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAaccggacc                    <  1:4207805/71‑1 (MQ=255)
                                               gTTTTCGCGCTGTAGAGACATTCCGGAATTGGCTTATCGa                                                   <  1:1214498/40‑1 (MQ=255)
                                               gTTTTCGCGCTGTAGAGACATTCCGGAATTGGCTTATCGa                                                   <  1:654214/40‑1 (MQ=255)
                                                  ttCGCGCTTTAGAGACCTTCCGGAATTGGCTTATc                                                     >  1:2233942/1‑35 (MQ=255)
                                                  ttCGCGCTGTAGAGACCTTCCGGAATTGGCTTATc                                                     >  1:828605/1‑35 (MQ=255)
                                                  ttCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGc                 >  1:1451243/1‑71 (MQ=255)
                                                  ttCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGc                 >  1:4214268/1‑71 (MQ=255)
                                                     gcgcTGTAGAGACATTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCgag              <  1:2707454/71‑1 (MQ=255)
                                                      cgcTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCgaga             <  1:2625757/71‑1 (MQ=255)
                                                         tGTAGAGACATTCCGGAATTGGCTTATCGATACCTTTCAACACGCGATCGTa                             <  1:4758766/52‑1 (MQ=255)
                                                           tAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCAc        >  1:2332543/1‑71 (MQ=255)
                                                            agagACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACg       >  1:4478037/1‑71 (MQ=255)
                                                             gagaCCTTCCGGAATTGGCTTATCGATACCTTTCaa                                         >  1:3777984/1‑36 (MQ=255)
                                                             gagaCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACgg      >  1:3598658/1‑71 (MQ=255)
                                                             gagaCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACgg      >  1:4555149/1‑71 (MQ=255)
                                                             gagaCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACg       <  1:2205894/70‑1 (MQ=255)
                                                             gagaCATTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAg                                 >  1:526731/1‑44 (MQ=255)
                                                               gaCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTc    >  1:4525278/1‑71 (MQ=255)
                                                               gaCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTc    >  1:1594330/1‑71 (MQ=255)
                                                               gaCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTc    >  1:1284793/1‑71 (MQ=255)
                                                                 ccTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTCgg  <  1:1888529/71‑1 (MQ=255)
                                                                 ccTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTCgg  <  1:1751004/71‑1 (MQ=255)
                                                                 ccTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTCgg  <  1:3870140/71‑1 (MQ=255)
                                                                  cTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCgaga             <  1:5092353/59‑1 (MQ=255)
                                                                  cTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTCgg  >  1:3044735/1‑70 (MQ=255)
                                                                  |                                                                     
AGGTCCAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTCGG  >  W3110S.gb/4093175‑4093310

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: