Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,183,420 C→T 50.6% C1232C (TGC→TGT gltB → glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,183,4200CT50.6% ‑5.0 / 71.9 77C1232C (TGC→TGTgltBglutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base C (29/9);  new base T (10/29);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.12e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

GCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGT  >  minE/2183360‑2183487
                                                            |                                                                   
gCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGaaa                                                              >  1:1132087/1‑68 (MQ=255)
gCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGaaa                                                              >  1:2748602/1‑68 (MQ=255)
gCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGaaa                                                              >  1:1550754/1‑68 (MQ=255)
gCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGaaa                                                              >  1:1395957/1‑68 (MQ=255)
gCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGaa                                                               >  1:3257779/1‑67 (MQ=255)
   ggCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAaacaa                                                          <  1:1375244/69‑1 (MQ=255)
            gAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGt                                                  >  1:111614/1‑68 (MQ=255)
            gAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGt                                                  >  1:1894726/1‑68 (MQ=255)
            gAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGt                                                  >  1:1735917/1‑68 (MQ=255)
            gAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGt                                                  >  1:2191679/1‑68 (MQ=255)
            gAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGt                                                  >  1:1500274/1‑68 (MQ=255)
            gAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGt                                                  >  1:393442/1‑68 (MQ=255)
            gAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGt                                                  >  1:703498/1‑68 (MQ=255)
            gAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGt                                                  >  1:918324/1‑68 (MQ=255)
            gAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGt                                                  >  1:920822/1‑68 (MQ=255)
            gAAGCTGCTGGAAACTGCCGAGCCACATGCAGGTAAGTCACTCTACTGTACCGAAAACAACccgccg                                                   <  1:2639378/67‑1 (MQ=255)
            gAAGCTGCTGGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACccgccg                                                   <  1:330133/67‑1 (MQ=255)
            gAAGCTGCTGGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACccgccg                                                   <  1:357386/67‑1 (MQ=255)
            gAAGCTGCTGGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACccgccg                                                   <  1:2941575/67‑1 (MQ=255)
            gAAGCTGCTGGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACccgccg                                                   <  1:2082570/67‑1 (MQ=255)
            gAAGCTGCTGGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACccgccg                                                   <  1:2393067/67‑1 (MQ=255)
            gAAGCTGCTGGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACccgccg                                                   <  1:1834479/67‑1 (MQ=255)
                ctgctgGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGAt                                             <  1:982952/69‑1 (MQ=255)
                ctgctgGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGAt                                             <  1:3082223/69‑1 (MQ=255)
                ctgctgGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGAt                                             <  1:580048/69‑1 (MQ=255)
                ctgctgGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGAt                                             <  1:2610036/69‑1 (MQ=255)
                ctgctgGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGAt                                             <  1:1058424/69‑1 (MQ=255)
                ctgctgGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGAt                                             <  1:1719534/69‑1 (MQ=255)
                   ctgGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAaaca                                                           >  1:1240965/1‑52 (MQ=25)
                   ctgGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAAc                                          >  1:225756/1‑69 (MQ=255)
                   ctgGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAAc                                          >  1:261281/1‑69 (MQ=255)
                   ctgGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAAc                                          >  1:1505580/1‑69 (MQ=255)
                   ctgGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAAc                                          >  1:1395242/1‑69 (MQ=255)
                   ctgGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAAc                                          >  1:232890/1‑69 (MQ=255)
                   ctgGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAAc                                          >  1:1345464/1‑69 (MQ=255)
                   ctgGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAAc                                          >  1:514371/1‑69 (MQ=255)
                   ctgGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAAc                                          >  1:1784893/1‑69 (MQ=255)
                   ctgGAAACTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAAc                                          >  1:1856091/1‑69 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATACCGGc                                       <  1:3412125/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:351995/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:2780988/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:407125/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:460401/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:49821/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:3232499/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:1206617/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:3146145/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:589507/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:3009857/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:1924626/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:1993824/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:2043713/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:1679624/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:1720395/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:2464838/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:2360141/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:2843336/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:2160520/65‑1 (MQ=255)
                       attcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:1831284/65‑1 (MQ=255)
                        ttcTGCCGAGCCACATGCAGGTAAGGCACTCTACTGTACCGAAAACAACCCGCCGTTTGATAACGGc                                       <  1:790923/65‑1 (MQ=255)
                                  ccGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAAcgcg                           <  1:1563192/69‑1 (MQ=255)
                                  ccGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAAcgcg                           <  1:1675732/69‑1 (MQ=255)
                                              aaGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTtgctgc                 >  1:2227963/1‑67 (MQ=255)
                                                      ctACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTACTGAACGCGCAGTTGCTGCAACAGGCGaa       >  1:2690741/1‑69 (MQ=255)
                                                        aCTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAAc     >  1:2252992/1‑69 (MQ=255)
                                                        aCTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAAc     >  1:2480672/1‑69 (MQ=255)
                                                        aCTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAAc     >  1:2799157/1‑69 (MQ=255)
                                                        aCTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAAc     >  1:1295367/1‑69 (MQ=255)
                                                           gCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCg   >  1:364182/1‑68 (MQ=255)
                                                           gCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCg   >  1:1890099/1‑68 (MQ=255)
                                                           gCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGt  >  1:998330/1‑69 (MQ=255)
                                                           gCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGt  >  1:72274/1‑69 (MQ=255)
                                                           gCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGt  >  1:445664/1‑69 (MQ=255)
                                                           gCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGt  >  1:1316772/1‑69 (MQ=255)
                                                           gCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGt  >  1:3202338/1‑69 (MQ=255)
                                                           gCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGt  >  1:2870363/1‑69 (MQ=255)
                                                           gCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGt  >  1:2203469/1‑69 (MQ=255)
                                                            |                                                                   
GCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGT  >  minE/2183360‑2183487

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: