Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,566,491 C→G 28.9% T187T (ACG→ACC mtlD ← mannitol‑1‑phosphate dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,566,4910CG28.9% 109.7 / 51.2 83T187T (ACG→ACCmtlDmannitol‑1‑phosphate dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (53/6);  new base G (9/15);  total (62/21)
Fisher's exact test for biased strand distribution p-value = 2.56e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

GCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTGCCGAAG  >  minE/2566432‑2566558
                                                           |                                                                   
gCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATc                                                             >  1:2665490/1‑68 (MQ=255)
gCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATc                                                             >  1:1579838/1‑68 (MQ=255)
gCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATc                                                             >  1:152780/1‑68 (MQ=255)
   aTCAGGTTGTCGGTTAACTCCACGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACa                                                         >  1:321522/1‑69 (MQ=255)
        gTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATcc                                                     >  1:1778188/1‑68 (MQ=255)
        gTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCa                                                    >  1:1792763/1‑69 (MQ=255)
                 tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:1463242/50‑1 (MQ=255)
                 tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:743692/50‑1 (MQ=255)
                 tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:754062/50‑1 (MQ=255)
                 tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:1601123/50‑1 (MQ=255)
                 tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:2486076/50‑1 (MQ=255)
                 tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:1430636/50‑1 (MQ=255)
                 tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:1393854/50‑1 (MQ=255)
                 tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:91391/50‑1 (MQ=255)
                 tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:2706564/50‑1 (MQ=255)
                 tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:1216927/50‑1 (MQ=255)
                 tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:3065205/50‑1 (MQ=255)
                 tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:3125401/50‑1 (MQ=255)
                 tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:100187/50‑1 (MQ=255)
                           ccTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:1438561/40‑1 (MQ=38)
                             tGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTAGCg                                <  1:2352130/68‑1 (MQ=255)
                             tGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCg                                <  1:516122/68‑1 (MQ=255)
                             tGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCg                                <  1:2566838/68‑1 (MQ=255)
                             tGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCg                                <  1:2286244/68‑1 (MQ=255)
                             tGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCg                                <  1:1580575/68‑1 (MQ=255)
                               ggATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCg                                <  1:291843/66‑1 (MQ=255)
                                gATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                                              <  1:952849/35‑1 (MQ=37)
                                 aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTa                                   >  1:3283925/1‑61 (MQ=255)
                                 aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt                           >  1:2955218/1‑69 (MQ=255)
                                 aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt                           >  1:1259490/1‑69 (MQ=255)
                                 aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt                           >  1:1076020/1‑69 (MQ=255)
                                 aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt                           >  1:2792618/1‑69 (MQ=255)
                                 aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt                           >  1:2893544/1‑69 (MQ=255)
                                 aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt                           >  1:82009/1‑69 (MQ=255)
                                 aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt                           >  1:1610260/1‑69 (MQ=255)
                                 aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt                           >  1:34754/1‑69 (MQ=255)
                                 aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt                           >  1:1919067/1‑69 (MQ=255)
                                 aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt                           >  1:1914400/1‑69 (MQ=255)
                                 aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt                           >  1:316714/1‑69 (MQ=255)
                                 aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt                           >  1:1778927/1‑69 (MQ=255)
                                 aTGTTAGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAAt                                                       >  1:2635553/1‑41 (MQ=255)
                                     tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTccg                       >  1:2399553/1‑68 (MQ=255)
                                     tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTccg                       >  1:1189838/1‑68 (MQ=255)
                                     tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa                       >  1:2589782/1‑69 (MQ=255)
                                     tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa                       >  1:3060920/1‑69 (MQ=255)
                                     tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa                       >  1:409354/1‑69 (MQ=255)
                                     tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa                       >  1:3202237/1‑69 (MQ=255)
                                     tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa                       >  1:2529406/1‑69 (MQ=255)
                                     tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa                       >  1:3256386/1‑69 (MQ=255)
                                     tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa                       >  1:181888/1‑69 (MQ=255)
                                         cAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCgg                   >  1:2532689/1‑69 (MQ=255)
                                         cAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCgg                   >  1:1662977/1‑69 (MQ=255)
                                             gcgcCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGa                                            >  1:539373/1‑40 (MQ=255)
                                              cgcCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATc                >  1:464729/1‑67 (MQ=255)
                                              cgcCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATc                >  1:658245/1‑67 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc                              >  1:3232922/1‑42 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc                              >  1:3272637/1‑42 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc                              >  1:331176/1‑42 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc                              >  1:1775047/1‑42 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc                              >  1:3211309/1‑42 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc                              >  1:3205733/1‑42 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc                              >  1:131942/1‑42 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc                              >  1:1954469/1‑42 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGAt                 >  1:780417/1‑55 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGt          >  1:2236224/1‑62 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa   >  1:500118/1‑69 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa   >  1:956632/1‑69 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa   >  1:325013/1‑69 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa   >  1:3249748/1‑69 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa   >  1:2669118/1‑69 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa   >  1:2663814/1‑69 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa   >  1:2194021/1‑69 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa   >  1:1746584/1‑69 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa   >  1:1713584/1‑69 (MQ=255)
                                                         tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa   >  1:1546786/1‑69 (MQ=255)
                                                          gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaag  >  1:3081363/1‑69 (MQ=255)
                                                          gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaag  >  1:2157509/1‑69 (MQ=255)
                                                          gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaag  >  1:1844802/1‑69 (MQ=255)
                                                          gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaag  >  1:1337802/1‑69 (MQ=255)
                                                          gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaag  >  1:1187491/1‑69 (MQ=255)
                                                          gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa   >  1:2930879/1‑68 (MQ=255)
                                                          gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa   >  1:1039215/1‑68 (MQ=255)
                                                          gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTACCAGCGGATCGTTAGTTgccgaag  >  1:1649968/1‑69 (MQ=255)
                                                           |                                                                   
GCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTGCCGAAG  >  minE/2566432‑2566558

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: