Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,566,476 T→C 40.7% A192A (GCA→GCG 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,4760TC40.7% 19.6 / 41.0 59A192A (GCA→GCGmtlDmannitol‑1‑phosphate dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (24/11);  new base C (9/15);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 3.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TTACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGG  >  minE/2566417‑2566541
                                                           |                                                                 
ttACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTg                                                          <  1:1145248/69‑1 (MQ=255)
ttACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTg                                                          <  1:890001/69‑1 (MQ=255)
 tACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGa                                                         >  1:2911975/1‑69 (MQ=255)
 tACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGa                                                         <  1:1982950/69‑1 (MQ=255)
 tACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGa                                                         >  1:2290288/1‑69 (MQ=255)
 tACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGa                                                         <  1:2433804/69‑1 (MQ=255)
     ttCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTg                                                     <  1:2675103/69‑1 (MQ=255)
               gCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATc                                            >  1:2665490/1‑68 (MQ=255)
               gCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATc                                            >  1:152780/1‑68 (MQ=255)
               gCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATc                                            >  1:1579838/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:2486076/50‑1 (MQ=255)
                                tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                             <  1:3125401/50‑1 (MQ=255)
                                tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                             <  1:3065205/50‑1 (MQ=255)
                                tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                             <  1:2706564/50‑1 (MQ=255)
                                tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                             <  1:743692/50‑1 (MQ=255)
                                tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                             <  1:100187/50‑1 (MQ=255)
                                tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                             <  1:754062/50‑1 (MQ=255)
                                tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                             <  1:1216927/50‑1 (MQ=255)
                                tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                             <  1:91391/50‑1 (MQ=255)
                                tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                             <  1:1393854/50‑1 (MQ=255)
                                tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                             <  1:1601123/50‑1 (MQ=255)
                                tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                             <  1:1430636/50‑1 (MQ=255)
                                tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt                                             <  1:1463242/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:1580575/68‑1 (MQ=255)
                                            tGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCg               <  1:2286244/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:316714/1‑69 (MQ=255)
                                                aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt          >  1:82009/1‑69 (MQ=255)
                                                aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt          >  1:34754/1‑69 (MQ=255)
                                                aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt          >  1:2955218/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:1259490/1‑69 (MQ=255)
                                                aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt          >  1:1610260/1‑69 (MQ=255)
                                                aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt          >  1:1778927/1‑69 (MQ=255)
                                                aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt          >  1:1914400/1‑69 (MQ=255)
                                                aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt          >  1:1919067/1‑69 (MQ=255)
                                                aTGTTAGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAAt                                      >  1:2635553/1‑41 (MQ=255)
                                                    tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTccg      >  1:1189838/1‑68 (MQ=255)
                                                    tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTccg      >  1:2399553/1‑68 (MQ=255)
                                                    tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa      >  1:3256386/1‑69 (MQ=255)
                                                    tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa      >  1:409354/1‑69 (MQ=255)
                                                    tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa      >  1:2529406/1‑69 (MQ=255)
                                                    tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa      >  1:181888/1‑69 (MQ=255)
                                                    tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa      >  1:3202237/1‑69 (MQ=255)
                                                    tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa      >  1:2589782/1‑69 (MQ=255)
                                                    tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa      >  1:3060920/1‑69 (MQ=255)
                                                        cAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCgg  >  1:1662977/1‑69 (MQ=255)
                                                        cAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCgg  >  1:2532689/1‑69 (MQ=255)
                                                           |                                                                 
TTACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGG  >  minE/2566417‑2566541

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: