Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,387 C→T 36.0% intergenic (‑193/‑344) mtlA ← / → yibI fused mannitol‑specific PTS enzyme IIABC components/predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,3870CT36.0% 48.3 / 40.3 72intergenic (‑193/‑344)mtlA/yibIfused mannitol‑specific PTS enzyme IIABC components/predicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (31/15);  new base T (14/12);  total (45/27)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGACGGCCTTTTGTGTGATATCAATCACATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAA  >  minE/2569324‑2569454
                                                               |                                                                   
tCGACGGCCTTTTGTGTGATATCAATCACATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTtc                                                                 >  1:295474/1‑68 (MQ=255)
tCGACGGCCTTTTGTGTGATATCAATCACATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTtc                                                                 >  1:1240241/1‑68 (MQ=255)
              tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                   <  1:1263072/68‑1 (MQ=255)
              tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                   <  1:2990831/68‑1 (MQ=255)
              tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                   <  1:3158782/68‑1 (MQ=255)
              tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                   <  1:3162857/68‑1 (MQ=255)
              tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                   <  1:1422686/68‑1 (MQ=255)
              tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                   <  1:2707026/68‑1 (MQ=255)
              tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                   <  1:627967/68‑1 (MQ=255)
              tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                   <  1:1225041/68‑1 (MQ=255)
              tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                   <  1:227148/68‑1 (MQ=255)
              tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                   <  1:108492/68‑1 (MQ=255)
              tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                   <  1:2134542/68‑1 (MQ=255)
                    atcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACaa                                                          >  1:47797/1‑55 (MQ=255)
                    atcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                             >  1:672710/1‑68 (MQ=255)
                    atcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                             >  1:14034/1‑68 (MQ=255)
                    atcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                             >  1:707349/1‑68 (MQ=255)
                    atcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                             >  1:1450704/1‑68 (MQ=255)
                    atcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                             >  1:177219/1‑68 (MQ=255)
                    atcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                             >  1:2189671/1‑68 (MQ=255)
                     tcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                             >  1:52824/1‑67 (MQ=255)
                     tcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                             >  1:719218/1‑67 (MQ=255)
                             aTAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGcaac                                    <  1:2052887/68‑1 (MQ=255)
                             aTAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGcaac                                    <  1:3414286/68‑1 (MQ=255)
                                aaTATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACg                                                      <  1:2487985/47‑1 (MQ=255)
                                aaTATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACg                                                      <  1:704003/47‑1 (MQ=255)
                                aaTATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACg                                                      <  1:89151/47‑1 (MQ=255)
                                aaTATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACg                                                      <  1:932580/47‑1 (MQ=255)
                                aaTATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACg                                                      <  1:944231/47‑1 (MQ=255)
                                aaTATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACg                                                      <  1:1016466/47‑1 (MQ=255)
                                aaTATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACg                                                      <  1:1750456/47‑1 (MQ=255)
                                aaTATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACg                                                      <  1:1135078/47‑1 (MQ=255)
                                aaTATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACg                                                      <  1:119985/47‑1 (MQ=255)
                                 atatAGGGGGT‑‑AGAGGTAATTGACGTGACGTTCATCACAAAACg                                                      <  1:2874033/44‑1 (MQ=38)
                                   ataGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGcaac                                    <  1:1562275/62‑1 (MQ=255)
                                    taGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                   <  1:830985/46‑1 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:1961062/1‑52 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:1108059/1‑52 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:818856/1‑52 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:1108178/1‑52 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:1177982/1‑52 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:582951/1‑52 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:1348522/1‑52 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:451795/1‑52 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:3378825/1‑52 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:3297281/1‑52 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:1922197/1‑52 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:2154126/1‑52 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:2783331/1‑52 (MQ=38)
                                      gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                           >  1:284309/1‑52 (MQ=38)
                                            agagagGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGAtttt                      >  1:4866/1‑67 (MQ=255)
                                            agagagGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTc                     >  1:336799/1‑68 (MQ=255)
                                                 ggTAATTTACGTTACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaa               >  1:546808/1‑69 (MQ=255)
                                                 ggTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaa               >  1:1585401/1‑69 (MQ=255)
                                                 ggTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaa               <  1:2773487/69‑1 (MQ=255)
                                                 ggTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaa                >  1:656360/1‑68 (MQ=255)
                                                  gTAATTGACGTGCCGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaaa              >  1:3186254/1‑69 (MQ=255)
                                                  gTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaaa              >  1:1526874/1‑69 (MQ=255)
                                                  gTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaaa              >  1:101484/1‑69 (MQ=255)
                                                  gTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaaa              >  1:1608997/1‑69 (MQ=255)
                                                  gTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaaa              >  1:3173288/1‑69 (MQ=255)
                                                  gTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaa               <  1:552662/68‑1 (MQ=255)
                                                    aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt            >  1:3209167/1‑69 (MQ=255)
                                                    aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt            >  1:2047654/1‑69 (MQ=255)
                                                    aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt            >  1:39578/1‑69 (MQ=255)
                                                    aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt            >  1:35644/1‑69 (MQ=255)
                                                    aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt            >  1:2960755/1‑69 (MQ=255)
                                                    aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt            >  1:1418759/1‑69 (MQ=255)
                                                    aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt            >  1:1731960/1‑69 (MQ=255)
                                                    aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt            >  1:3430206/1‑69 (MQ=255)
                                                    aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATATTCGCTAAAAAt            >  1:2948409/1‑69 (MQ=255)
                                                      ttGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATg           >  1:167361/1‑68 (MQ=255)
                                                       tgacgtgacgTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGc         <  1:2952990/69‑1 (MQ=255)
                                                       tgacgtgacgTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGc         >  1:31724/1‑69 (MQ=255)
                                                        gacgtgacgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGc         >  1:92958/1‑68 (MQ=255)
                                                        gacgtgacgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCa        >  1:1319858/1‑69 (MQ=255)
                                                        gacgtgacgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCa        >  1:2494515/1‑69 (MQ=255)
                                                              acgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTaaa  <  1:693114/69‑1 (MQ=255)
                                                              acgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTaaa  >  1:1779690/1‑69 (MQ=255)
                                                              acgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTaaa  <  1:57489/69‑1 (MQ=255)
                                                              acgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTaaa  >  1:1383932/1‑69 (MQ=255)
                                                              acgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTaaa  <  1:1071301/69‑1 (MQ=255)
                                                              acgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCAGCGATTTTCGCTAAAAATGGCATTTaaa  <  1:406137/69‑1 (MQ=255)
                                                               |                                                                   
TCGACGGCCTTTTGTGTGATATCAATCACATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAA  >  minE/2569324‑2569454

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: