Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,671,600 C→T 72.0% P192L (CCG→CTG)  malQ → 4‑alpha‑glucanotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,671,6000CT72.0% 29.4 / 9.8 32P192L (CCG→CTG) malQ4‑alpha‑glucanotransferase
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base T (20/3);  total (24/8)
Fisher's exact test for biased strand distribution p-value = 2.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCTGGTGGATGTGGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATCCGGCAAATCCGGAGAGCGCCAGCCCATACAGCCCGTCTTCTCGCCGTTGGCTGAA  >  minE/2671533‑2671654
                                                                   |                                                      
aTGCTGGTGGATGTGGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATcc                                                        >  1:1206531/1‑68 (MQ=255)
aTGCTGGTGGATGTGGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATcc                                                        >  1:2765836/1‑68 (MQ=255)
aTGCTGGTGGATGTGGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATcc                                                        >  1:2447002/1‑68 (MQ=255)
aTGCTGGTGGATGTGGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATcc                                                        >  1:2136154/1‑68 (MQ=255)
aTGCTGGTGGATGTGGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATcc                                                        >  1:2028016/1‑68 (MQ=255)
aTGCTGGTGGATGTGGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATcc                                                        >  1:446717/1‑68 (MQ=255)
aTGCTGGTGGATGTGGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATcc                                                        >  1:120435/1‑68 (MQ=255)
     ggtggATGTGGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATCCGGCaaa                                                  <  1:3343184/69‑1 (MQ=255)
     ggtggATGTGGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATCCGGCaaa                                                  <  1:927127/69‑1 (MQ=255)
     ggtggATGTGGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATCCGGCaaa                                                  <  1:2073551/69‑1 (MQ=255)
     ggtggATGTGGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATCCGGCaaa                                                  <  1:1462228/69‑1 (MQ=255)
            gtgGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATCCGGCAAATCCGgag                                           <  1:1558114/69‑1 (MQ=255)
                                gTTTATTGGTCTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGGAAAGTGCCAGCCCGTACAGccc                       <  1:2123908/69‑1 (MQ=255)
                                gTTTATTGGTCTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGGAAAGTGCCAGCCCGTACAGccc                       <  1:30087/69‑1 (MQ=255)
                                  ttaTTGGTCTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGGAAAGTGCCAGCCCGTACAGccc                       <  1:1282253/65‑1 (MQ=255)
                                   cATTGGCCTGAACCCGATTCATGCGCTCTATCCGGCAAATCCGGAGAGCGCCAGCCCATACAGCCCGtc                    >  1:27020/1‑69 (MQ=255)
                                   cATTGGCCTGAACCCGATTCATGCGCTCTATCCGGCAAATCCGGAGAGCGCCAGCCCATACAGCCCGt                     >  1:1349442/1‑68 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:491205/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:606383/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:649323/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:758150/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:3335283/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:3330942/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:987190/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:759215/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:3311130/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:3113831/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:2283043/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:2276538/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:218584/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:1979817/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:1968167/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:1966171/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:1953555/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:1940829/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:1813485/2‑39 (MQ=255)
                                         tcTGAACCCGATTCATGCGCTCTATCTGGCAAATCCGgaa                                           >  1:1349451/2‑39 (MQ=255)
                                                     tCATGCGCTCTATCCGGCAAATCCGGAGAGCGCCAGCCCa                               >  1:921919/1‑40 (MQ=255)
                                                     tCATGCGCTCTATCCGGCAAATCCGGAGAGCGCCAGCCCATACAGCCCGTCTTCTCGCCGTTGGCTGaa  >  1:1473120/1‑69 (MQ=255)
                                                                   |                                                      
ATGCTGGTGGATGTGGCAAAACGTGGCGGGTCGTTCATTGGCCTGAACCCGATTCATGCGCTCTATCCGGCAAATCCGGAGAGCGCCAGCCCATACAGCCCGTCTTCTCGCCGTTGGCTGAA  >  minE/2671533‑2671654

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: