Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,401 C→T 27.2% intergenic (‑207/‑330) 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,4010CT27.2% 111.1 / 37.1 95intergenic (‑207/‑330)mtlA/yibIfused mannitol‑specific PTS enzyme IIABC components/predicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (44/25);  new base T (14/12);  total (58/37)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

TGTGATATCAATCACATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAATGTGAGTAG  >  minE/2569338‑2569463
                                                               |                                                              
tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                            <  1:2134542/68‑1 (MQ=255)
tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                            <  1:108492/68‑1 (MQ=255)
tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                            <  1:227148/68‑1 (MQ=255)
tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                            <  1:2707026/68‑1 (MQ=255)
tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                            <  1:2990831/68‑1 (MQ=255)
tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                            <  1:627967/68‑1 (MQ=255)
tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                            <  1:1225041/68‑1 (MQ=255)
tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                            <  1:1263072/68‑1 (MQ=255)
tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                            <  1:1422686/68‑1 (MQ=255)
tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                            <  1:3158782/68‑1 (MQ=255)
tgtgATATCAGTCACATAAATGTAGGGGGTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCtt                                                            <  1:3162857/68‑1 (MQ=255)
      atcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                                      >  1:177219/1‑68 (MQ=255)
      atcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                                      >  1:2189671/1‑68 (MQ=255)
      atcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                                      >  1:1450704/1‑68 (MQ=255)
      atcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                                      >  1:14034/1‑68 (MQ=255)
      atcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                                      >  1:672710/1‑68 (MQ=255)
      atcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                                      >  1:707349/1‑68 (MQ=255)
       tcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                                      >  1:52824/1‑67 (MQ=255)
       tcaatcaCATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTaa                                                      >  1:719218/1‑67 (MQ=255)
               aTAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGcaac                                             <  1:3414286/68‑1 (MQ=255)
               aTAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGcaac                                             <  1:2052887/68‑1 (MQ=255)
                  aaTATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACg                                                               <  1:704003/47‑1 (MQ=255)
                  aaTATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACg                                                               <  1:2487985/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:119985/47‑1 (MQ=255)
                  aaTATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACg                                                               <  1:1135078/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:1348522/1‑52 (MQ=38)
                        gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                                    >  1:451795/1‑52 (MQ=38)
                        gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                                    >  1:284309/1‑52 (MQ=38)
                        gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                                    >  1:2783331/1‑52 (MQ=38)
                        gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                                    >  1:582951/1‑52 (MQ=38)
                        gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                                    >  1:1177982/1‑52 (MQ=38)
                        gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                                    >  1:1108178/1‑52 (MQ=38)
                        gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                                    >  1:2154126/1‑52 (MQ=38)
                        gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                                    >  1:1922197/1‑52 (MQ=38)
                        gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                                    >  1:1961062/1‑52 (MQ=38)
                        gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                                    >  1:3297281/1‑52 (MQ=38)
                        gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                                    >  1:818856/1‑52 (MQ=38)
                        gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                                    >  1:3378825/1‑52 (MQ=38)
                        gggggTAGAGAGGTAATTGACGTGATGCTCATCACAAAATGCTTGTGTaaaa                                                    >  1:1108059/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:2773487/69‑1 (MQ=255)
                                   ggTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaa                        >  1:1585401/1‑69 (MQ=255)
                                   ggTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaa                         >  1:656360/1‑68 (MQ=255)
                                    gTAATTGACGTGCCGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaaa                       >  1:3186254/1‑69 (MQ=255)
                                    gTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaaa                       >  1:3173288/1‑69 (MQ=255)
                                    gTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaaa                       >  1:101484/1‑69 (MQ=255)
                                    gTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaaa                       >  1:1526874/1‑69 (MQ=255)
                                    gTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaaa                       >  1:1608997/1‑69 (MQ=255)
                                    gTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTaaa                        <  1:552662/68‑1 (MQ=255)
                                      aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt                     >  1:2960755/1‑69 (MQ=255)
                                      aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt                     >  1:1418759/1‑69 (MQ=255)
                                      aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt                     >  1:39578/1‑69 (MQ=255)
                                      aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt                     >  1:3209167/1‑69 (MQ=255)
                                      aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt                     >  1:35644/1‑69 (MQ=255)
                                      aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt                     >  1:1731960/1‑69 (MQ=255)
                                      aaTTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAAt                     >  1:2047654/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:31724/1‑69 (MQ=255)
                                         tgacgtgacgTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGc                  <  1:2952990/69‑1 (MQ=255)
                                          gacgtgacgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGc                  >  1:92958/1‑68 (MQ=255)
                                          gacgtgacgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCa                 >  1:2494515/1‑69 (MQ=255)
                                          gacgtgacgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCa                 >  1:1319858/1‑69 (MQ=255)
                                                acgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTaaa           >  1:1779690/1‑69 (MQ=255)
                                                acgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTaaa           <  1:693114/69‑1 (MQ=255)
                                                acgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTaaa           <  1:1071301/69‑1 (MQ=255)
                                                acgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTaaa           <  1:57489/69‑1 (MQ=255)
                                                acgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTaaa           >  1:1383932/1‑69 (MQ=255)
                                                acgtTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCAGCGATTTTCGCTAAAAATGGCATTTaaa           <  1:406137/69‑1 (MQ=255)
                                                     catcaCAAAACGCTTGTGTAAAAAAGCAACAAggcggc                                     >  1:1443718/1‑38 (MQ=255)
                                                     catcaCAAAACGCTTGTGTAAAAAAGCAACAAggcggc                                     >  1:1960076/1‑38 (MQ=255)
                                                     catcaCAAAACGCTTGTGTAAAAAAGCAACAAggcggc                                     >  1:2224698/1‑38 (MQ=255)
                                                     catcaCAAAACGCTTGTGTAAAAAAGCAACAAggcg                                       >  1:2186992/1‑36 (MQ=255)
                                                        cacaAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAATGTGagta   >  1:1822291/1‑69 (MQ=255)
                                                        cacaAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAATGTGagta   >  1:2002460/1‑69 (MQ=255)
                                                        cacaAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAATGTGagta   >  1:240408/1‑69 (MQ=255)
                                                        cacaAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAATGTGagta   >  1:706844/1‑69 (MQ=255)
                                                        cacaAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAATGTGagta   >  1:2252111/1‑69 (MQ=255)
                                                        cacaAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAATGTGagt    <  1:1965206/68‑1 (MQ=255)
                                                        cacaAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAATGTGagt    <  1:3239816/68‑1 (MQ=255)
                                                        cacaAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAATGTGagt    <  1:2840199/68‑1 (MQ=255)
                                                        cacaAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAATGTGagt    <  1:3153701/68‑1 (MQ=255)
                                                        cacaAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAATGTGagt    <  1:375962/68‑1 (MQ=255)
                                                         acaAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAATGTGagtag  >  1:2447116/1‑69 (MQ=255)
                                                               |                                                              
TGTGATATCAATCACATAAATATAGGGGGTAGAGAGGTAATTGACGTGACGTTCATCACAAAACGCTTGTGTAAAAAAGCAACAAGGCGGCGATTTTCGCTAAAAATGGCATTTAAATGTGAGTAG  >  minE/2569338‑2569463

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: