Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,105,113 A→T 35.3% V371V (GTA→GTT kbaZ → tagatose 6‑phosphate aldolase 1, kbaZ subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,105,1130AT35.3% 71.3 / 68.6 85V371V (GTA→GTTkbaZtagatose 6‑phosphate aldolase 1, kbaZ subunit
Reads supporting (aligned to +/- strand):  ref base A (23/32);  new base T (14/16);  total (37/48)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

CGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGTATCTGCCGCTGCAGTACGT  >  minE/2105052‑2105178
                                                             |                                                                 
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:96810/68‑1 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:2510964/68‑1 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:2685770/68‑1 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:2021700/68‑1 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:2836969/68‑1 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:751523/68‑1 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:1977068/68‑1 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:3124354/68‑1 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:1478582/68‑1 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:1356686/68‑1 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATc                                                            >  1:1588990/1‑69 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATc                                                            >  1:2380613/1‑69 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATc                                                            >  1:697935/1‑69 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATc                                                            >  1:1072511/1‑69 (MQ=255)
cGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATc                                                            >  1:1325379/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:962014/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:3016641/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:3039739/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:3081919/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:3083343/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:1795605/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:2917874/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:199240/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:1087931/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:1515907/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:376225/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:430670/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:1476830/1‑69 (MQ=255)
     cgcgTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCg                                                       >  1:1024919/1‑69 (MQ=255)
        gTGCGCTATTTCTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:622082/60‑1 (MQ=255)
        gTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGa                                                     <  1:1783497/68‑1 (MQ=255)
        gTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGa                                                     <  1:842721/68‑1 (MQ=255)
        gTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGa                                                     <  1:2751036/68‑1 (MQ=255)
        gTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGa                                                     <  1:1363170/68‑1 (MQ=255)
        gTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGa                                                     <  1:2503841/68‑1 (MQ=255)
             cTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCAcc                                               <  1:1223570/69‑1 (MQ=255)
             cTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCAcc                                               <  1:367021/69‑1 (MQ=255)
                     ggCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTccgc                                       <  1:372500/69‑1 (MQ=255)
                     ggCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTccgc                                       <  1:2664496/69‑1 (MQ=255)
                     ggCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTccgc                                       >  1:1503699/1‑69 (MQ=255)
                        cGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:1900242/44‑1 (MQ=255)
                          gACAGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAATCTGGCGGATTCACCAATTccgctgccg                                  <  1:1024517/69‑1 (MQ=255)
                             aGCGAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:619725/39‑1 (MQ=38)
                             aGCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:1458037/39‑1 (MQ=255)
                             aGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTccgcttccgctg                               <  1:441645/69‑1 (MQ=255)
                             aGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTccgctgccgctg                               <  1:2482234/69‑1 (MQ=255)
                             aGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTccgctgccgctg                               <  1:3300272/69‑1 (MQ=255)
                             aGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTccgctgccgctg                               <  1:1206245/69‑1 (MQ=255)
                             aGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTccgctgccgctg                               <  1:2009218/69‑1 (MQ=255)
                             aGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTccgctgccgctg                               <  1:2102521/69‑1 (MQ=255)
                              gCCAGATTGATGACGCTTTCGCTCATCTGGTTCGTAAt                                                             <  1:849284/38‑1 (MQ=255)
                              gCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGa                                                     <  1:1861754/46‑1 (MQ=255)
                                 aGATTGATTACGCTTTCGCTCTTCTGGTACGTAATCTGGCGGa                                                     <  1:338477/43‑1 (MQ=38)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGa                                                     >  1:2798253/1‑40 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             >  1:454355/1‑48 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             >  1:1477459/1‑48 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             >  1:563968/1‑48 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             <  1:57623/48‑1 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             >  1:40998/1‑48 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             >  1:3209841/1‑48 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             >  1:1650194/1‑48 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             <  1:312385/48‑1 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             <  1:3063292/48‑1 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             <  1:3059294/48‑1 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             >  1:1892222/1‑48 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             <  1:1337698/48‑1 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             >  1:2024552/1‑48 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             <  1:2127852/48‑1 (MQ=255)
                                    ttGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             >  1:1142170/1‑48 (MQ=255)
                                     tgatgaCGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCaa                                             <  1:994852/47‑1 (MQ=255)
                                     tgatgaCGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGcc                        >  1:3249313/1‑68 (MQ=255)
                                     tgatgaCGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGcc                        >  1:777078/1‑68 (MQ=255)
                                        tgaCGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGt                     <  1:2430833/68‑1 (MQ=255)
                                        tgaCGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGt                     <  1:1592975/68‑1 (MQ=255)
                                        tgaCGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGt                     <  1:553543/68‑1 (MQ=255)
                                        tgaCGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGt                     <  1:1108994/68‑1 (MQ=255)
                                        tgaCGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGt                     <  1:1901272/68‑1 (MQ=255)
                                         gaCGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGTa                    <  1:3274411/68‑1 (MQ=255)
                                         gaCGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGTa                    <  1:1816371/68‑1 (MQ=255)
                                         gaCGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGTa                    <  1:1494767/68‑1 (MQ=255)
                                                 cGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGTATCTGCCGc            >  1:1946727/1‑68 (MQ=255)
                                                 cGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGTATCTGCCGc            >  1:952031/1‑68 (MQ=255)
                                                     cATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGTATCTGCCGCTGCAg       >  1:222205/1‑69 (MQ=255)
                                                     cATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGTATCTGCCGCTGCAg       >  1:3046022/1‑69 (MQ=255)
                                                     cATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGTATCTGCCGCTGCAg       >  1:101541/1‑69 (MQ=255)
                                                             aCGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGTATCTGCCGCTGCAGTACGt  <  1:876237/66‑1 (MQ=255)
                                                             aCGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGTATCTGCCGCTGCAGTACGt  <  1:263167/66‑1 (MQ=255)
                                                             aCGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGTATCTGCCGCTGCAGTACGt  <  1:1583354/66‑1 (MQ=255)
                                                             |                                                                 
CGGATCGCGTGCGCTATTACTGGCCGGACAGCCAGATTGATGACGCTTTCGCTCATCTGGTACGTAATCTGGCGGATTCACCAATTCCGCTGCCGCTGATCAGCCAGTATCTGCCGCTGCAGTACGT  >  minE/2105052‑2105178

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: