Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 473,506 G→A 53.1% K356K (AAG→AAA mngA → fused 2‑O‑a‑mannosyl‑D‑glycerate specific PTS enzyme IIABC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE473,5060GA53.1% ‑0.0 / 69.1 49K356K (AAG→AAAmngAfused 2‑O‑a‑mannosyl‑D‑glycerate specific PTS enzyme IIABC components
Reads supporting (aligned to +/- strand):  ref base G (12/11);  new base A (17/9);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

TCGCAAATCTTCGGGCTACAAGATCTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCGGCCTATACCGCCTATTCTCTGG  >  minE/473441‑473573
                                                                 |                                                                   
tCGCAAATCTTCGGGCTACAAGATCTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAGCTg                                                                  <  1:3301730/69‑1 (MQ=255)
      aTCTTCGGGCTACAAGATCTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAGCTGGGcggcgg                                                          >  1:1254100/1‑71 (MQ=255)
      aTCTTCGGGCTACAAGATCTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAGCTGGGcggcgg                                                          >  1:3782419/1‑71 (MQ=255)
        cTTCGGGCTACAAGATCTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGgc                                                        >  1:3474064/1‑71 (MQ=255)
           cGGGCTACAAGATCTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCggg                                                         <  1:3238711/67‑1 (MQ=255)
           cGGGCTACAAGATCTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGgc                                                        <  1:2132885/68‑1 (MQ=255)
               cTACAAGATCTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCgg                                                 <  1:2384453/71‑1 (MQ=255)
               cTACAAGATCTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGGCTGCTCgg                                                 <  1:1083786/71‑1 (MQ=255)
                  cAAGATCTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAAt                                              >  1:1648425/1‑71 (MQ=255)
                    aGATCTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGggcggc                                                            >  1:2245640/1‑55 (MQ=255)
                     gATCTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAAtttt                                           <  1:764585/71‑1 (MQ=255)
                        cTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAATTTTGAt                                        >  1:1403180/1‑71 (MQ=255)
                           tttAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAGCt                                                                   >  1:3618411/1‑41 (MQ=255)
                           tttAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAAtttt                                           >  1:3460979/1‑65 (MQ=255)
                           tttAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAAtttt                                           >  1:480020/1‑65 (MQ=255)
                           tttAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAATTTTGAt                                        <  1:3978929/68‑1 (MQ=255)
                           tttAATGAAGAAAACTCCTGGC‑‑TGGATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAATTTTGAt                                        <  1:3372004/66‑1 (MQ=255)
                             tAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAATTTTGAt                                        >  1:2660182/1‑66 (MQ=255)
                             tAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGggcggc                                                            >  1:1631794/1‑46 (MQ=255)
                             aaaTGAAGAAAACTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGgctgc                                                     >  1:2715238/2‑53 (MQ=255)
                               aTGAAGAAAACTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAATTTTGATGGTAcc                                  <  1:2582037/70‑1 (MQ=255)
                                tGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGgctg                                                      >  1:1526277/1‑49 (MQ=255)
                                tGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGgctg                                                      >  1:1531125/1‑49 (MQ=255)
                                tGAAGAAAACTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAAtttt                                           <  1:3216055/60‑1 (MQ=255)
                                 gaagaaAACTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAATTTTGAt                                        <  1:1899664/62‑1 (MQ=255)
                                    gaaAACTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAATTTTGa                                         >  1:3438981/1‑58 (MQ=255)
                                    gaaAACTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAATTTTg                                          >  1:803846/1‑57 (MQ=255)
                                     aaaaCTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAAtt                                             <  1:918132/53‑1 (MQ=255)
                                     aaaaCTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAAtt                                             <  1:2657710/53‑1 (MQ=255)
                                       aaCTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGaa                                               >  1:3382522/1‑49 (MQ=255)
                                       aaCTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGaa                                               >  1:798374/1‑49 (MQ=255)
                                         cTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAATTTTGATg                                       >  1:1319474/1‑55 (MQ=255)
                                         cTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCgg                       <  1:3028821/71‑1 (MQ=255)
                                         cTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCgg                       <  1:2732039/71‑1 (MQ=255)
                                         cTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAATTTTg                                          >  1:739767/1‑52 (MQ=255)
                                         cTCCTGGCTGTGGATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCgg                       >  1:1849790/1‑71 (MQ=255)
                                            cTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCGGcc                     <  1:431439/70‑1 (MQ=255)
                                            cTGGCTGTGGATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCGGcc                     <  1:3381861/70‑1 (MQ=255)
                                            cTGGCTGTGGATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCGGcc                     <  1:2919703/70‑1 (MQ=255)
                                            cTGGCTGTGGATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCGGcc                     >  1:269703/1‑70 (MQ=255)
                                                   tgGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAAttt                                            <  1:688319/40‑1 (MQ=255)
                                                   tgGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAATTTTGATGGTAcc                                  <  1:2460106/50‑1 (MQ=255)
                                                     gATGTACCGCAAACTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCGGCCTATACCGcc            >  1:3151170/1‑70 (MQ=255)
                                                       tGTACCGCAAACTGGGCGGCGGGCTGCTCGGAATTTTGATgg                                      >  1:844027/1‑42 (MQ=255)
                                                         tACCGCAAACTGGGCGGCGGGCTGCTCGGAATTTTGa                                         >  1:2463525/1‑37 (MQ=255)
                                                          aCCGCAAACTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCGGCCTATACCGCCTATtct      >  1:3616388/1‑71 (MQ=255)
                                                           ccGCAAGCTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCGGCCTATACCGCCTATtctc     >  1:3812150/1‑71 (MQ=255)
                                                           ccGCAAACTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCGGCCTATACCGCCTATtctc     >  1:3262884/1‑71 (MQ=255)
                                                              cAAGCTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCGGCCTATACCGCCTATTCTCTgg  >  1:2787645/1‑71 (MQ=255)
                                                                 |                                                                   
TCGCAAATCTTCGGGCTACAAGATCTGTTTAATGAAGAAAACTCCTGGCTGTGGATGTACCGCAAGCTGGGCGGCGGGCTGCTCGGAATTTTGATGGTACCGGTGCTCGCGGCCTATACCGCCTATTCTCTGG  >  minE/473441‑473573

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: