Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,070,862 G→T 100% R228L (CGC→CTC)  ygjG → putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,070,8620GT75.5% 69.9 / 34.9 49R228L (CGC→CTC) ygjGputrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (12/0);  new base T (0/37);  total (12/37)
Fisher's exact test for biased strand distribution p-value = 1.08e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAA  >  minE/2070796‑2070886
                                                                  |                        
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCGGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:1695816/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTTATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:1657444/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:1082522/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:663569/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:625355/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:51375/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:472340/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:1996317/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:1940500/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:1675723/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:1559886/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:1553879/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:1419634/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:1406760/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:1309433/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:1165951/1‑67 (MQ=255)
gCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATgcg                          >  1:1166179/1‑67 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCTTCATGTACCGTTTGTCAATATCGAAGCCATGCTCACGGcgc                  <  1:542017/67‑3 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACTGctc                  <  1:607027/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1998149/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:212515/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:204115/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:2115943/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:2153514/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:2193832/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:266169/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:315874/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1201359/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:495766/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:554033/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1038677/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:68117/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:87131/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:92339/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:951814/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:965900/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:19803/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1423737/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1399540/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1573570/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1576088/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1617391/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1390205/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1369961/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:16866/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1317355/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1753042/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1812812/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1268184/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1868025/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:182089/67‑1 (MQ=255)
        tttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1856558/67‑1 (MQ=255)
         ttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:1641125/66‑1 (MQ=255)
         ttATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGctc                  <  1:2156382/66‑1 (MQ=255)
                                         ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaa  >  1:237429/1‑50 (MQ=255)
                                         ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaa  >  1:1834486/1‑50 (MQ=255)
                                         ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaa  >  1:630661/1‑50 (MQ=255)
                                         ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaa  >  1:269248/1‑50 (MQ=255)
                                         ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaa  >  1:169685/1‑50 (MQ=255)
                                         ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaa  >  1:708122/1‑50 (MQ=255)
                                         ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaa  >  1:752690/1‑50 (MQ=255)
                                         ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaa  >  1:839526/1‑50 (MQ=255)
                                         ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaa  >  1:1411295/1‑50 (MQ=255)
                                         ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaa  >  1:1871972/1‑50 (MQ=255)
                                         ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaa  >  1:306304/1‑50 (MQ=255)
                                         ccGTATGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaa  >  1:294742/1‑50 (MQ=255)
                                                                  |                        
GCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAA  >  minE/2070796‑2070886

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: