Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,712,033 G→A 13.6% M328I (ATG→ATA gidA → glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,712,0330GA13.6% 119.0 / 9.7 44M328I (ATG→ATAgidAglucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base G (18/20);  new base A (4/2);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

ATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGATCCTC  >  W3110S.gb/3711963‑3712100
                                                                      |                                                                   
aTGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATg                                                                     <  1:2731331/71‑1 (MQ=255)
aTGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATg                                                                     >  1:2557779/1‑71 (MQ=255)
aTGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATg                                                                     <  1:2168636/71‑1 (MQ=255)
 tGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATg                                                                     >  1:2142641/1‑70 (MQ=255)
   aaaTTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATACAg                                                                  >  1:844342/1‑71 (MQ=255)
     aTTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAggg                                                                <  1:353741/71‑1 (MQ=255)
      tttATCCGAACGGTATCACCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATACAggg                                                                >  1:741205/1‑70 (MQ=255)
                cGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAg                                                                  >  1:2673726/1‑58 (MQ=255)
                       tcCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAAcgcg                                                   >  1:261164/1‑66 (MQ=255)
                       tcCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAAcgcg                                                   >  1:2588299/1‑66 (MQ=255)
                         caccaGCCTTCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATACAGGGGATGGAAAACGCGAAGATCg                                            <  1:562605/71‑1 (MQ=255)
                          accaTCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGAt                                              <  1:139612/68‑1 (MQ=255)
                            caGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGc                                         >  1:1118948/1‑71 (MQ=255)
                            caGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGc                                         >  1:293488/1‑71 (MQ=255)
                            caGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGc                                         >  1:2544943/1‑71 (MQ=255)
                            caGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGc                                         >  1:1785681/1‑71 (MQ=255)
                             aGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCg                                        >  1:1530619/1‑71 (MQ=255)
                             aGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCg                                        <  1:662479/71‑1 (MQ=255)
                                 tGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCg                                    >  1:1071157/1‑71 (MQ=255)
                                 tGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCATCCg                                    >  1:709130/1‑71 (MQ=255)
                                  gCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATg                                                                     <  1:1631593/37‑1 (MQ=255)
                                   ccGTTCGATGTGCAGATGCAAATCGTCCGCTCTATACAGGGGATGGAAAACGCGAAGATCGTGCGTCCggg                                  >  1:2630739/1‑71 (MQ=255)
                                         gATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGc                            <  1:1493362/71‑1 (MQ=255)
                                          aTGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGcc                           <  1:1245105/71‑1 (MQ=255)
                                             tgCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCAtt                        <  1:2719860/71‑1 (MQ=255)
                                                aGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAg                     <  1:990882/71‑1 (MQ=255)
                                                aGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAg                     <  1:234595/71‑1 (MQ=255)
                                                aGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAg                     <  1:926438/71‑1 (MQ=255)
                                                aGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAg                     <  1:2136671/71‑1 (MQ=255)
                                                 gATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAg                     >  1:1574332/1‑70 (MQ=255)
                                                   tGCAAATCGTCCGCTCTTTGCAGGGGATGGAAAACGCGAAGAt                                              <  1:1998100/43‑1 (MQ=255)
                                                   tGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGAt                                              <  1:2621849/43‑1 (MQ=255)
                                                    gCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTAt                  >  1:175606/1‑70 (MQ=255)
                                                      aaaTCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc               <  1:1632698/71‑1 (MQ=255)
                                                       aaTCGTCCGCTCTATACAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTAt                  <  1:332897/67‑1 (MQ=255)
                                                        aTCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCg                                        <  1:2656770/44‑1 (MQ=255)
                                                        aTCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCg                                        <  1:1581078/44‑1 (MQ=255)
                                                        aTCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCg                                        <  1:1988397/44‑1 (MQ=255)
                                                        aTCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGActt             >  1:2609402/1‑71 (MQ=255)
                                                        aTCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGActt             >  1:258362/1‑71 (MQ=255)
                                                        aTCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGActt             >  1:2140957/1‑71 (MQ=255)
                                                          cGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGActtct           >  1:1643034/1‑71 (MQ=255)
                                                           gTCCGCTCTATACAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTa                               >  1:1699367/1‑50 (MQ=255)
                                                              cGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCg        <  1:887471/70‑1 (MQ=255)
                                                              cGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa       <  1:317619/71‑1 (MQ=255)
                                                              cGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa       >  1:26479/1‑71 (MQ=255)
                                                                ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGt                                 <  1:184125/43‑1 (MQ=255)
                                                                   tATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc               <  1:589507/58‑1 (MQ=255)
                                                                   tATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGATCCTc  >  1:1811092/1‑71 (MQ=255)
                                                                      |                                                                   
ATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGATCCTC  >  W3110S.gb/3711963‑3712100

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