Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,473,322 G→A 42.1% P199P (CCG→CCA glmS → L‑glutamine:D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,473,3220GA42.1% 17.4 / 40.6 57P199P (CCG→CCAglmSL‑glutamine:D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (8/25);  new base A (4/20);  total (12/45)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGATTGGCCTGGGGATGGGCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAATCACTCGCCGTTCGGTAAA  >  minE/2473261‑2473387
                                                             |                                                                 
tGGTGATTGGCCTGGGTATGGGCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCAGTGActc                                                            >  1:589304/1‑67 (MQ=255)
tGGTGATTGGCCTGGGTATGGGCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCAGTGActc                                                            >  1:3149557/1‑67 (MQ=255)
tGGTGATTGGCCTGGGTATGGGCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCAGTGActc                                                            >  1:2992960/1‑67 (MQ=255)
tGGTGATTGGCCTGGGTATGGGCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCAGTGActc                                                            >  1:1980993/1‑67 (MQ=255)
 ggTGATTGGCCTGGGGATGGGCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCg                                                           >  1:694257/1‑69 (MQ=255)
 ggTGATTGGCCTGGGGATGGGCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCg                                                           >  1:197130/1‑69 (MQ=255)
        ggCCTGGGGATGGGCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCttt                                                    <  1:827277/69‑1 (MQ=255)
        ggCCTGGGGATGGGCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCttt                                                    <  1:2711515/69‑1 (MQ=255)
        ggCCTGGGGATGGGCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCttt                                                    <  1:2513998/69‑1 (MQ=255)
        ggCCTGGGGATGGGCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCttt                                                    <  1:2017011/69‑1 (MQ=255)
        ggCCTGGGGATGGGCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCttt                                                    <  1:174985/69‑1 (MQ=255)
                    ggCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGaa                                        <  1:1060343/69‑1 (MQ=255)
                    ggCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGaa                                        <  1:895617/69‑1 (MQ=255)
                    ggCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGaa                                        <  1:2994743/69‑1 (MQ=255)
                    ggCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGaa                                        <  1:2894976/69‑1 (MQ=255)
                      cGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAag                                       <  1:954878/68‑1 (MQ=255)
                      cGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAag                                       <  1:913148/68‑1 (MQ=255)
                      cGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAag                                       <  1:3270564/68‑1 (MQ=255)
                      cGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAag                                       <  1:2180323/68‑1 (MQ=255)
                      cGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAag                                       <  1:1178405/68‑1 (MQ=255)
                      cGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAag                                       <  1:3093915/68‑1 (MQ=255)
                          aaCTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGc                                  >  1:473081/1‑69 (MQ=255)
                              ttATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCt                                                >  1:1013143/1‑51 (MQ=255)
                              ttATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGat                               >  1:2325363/1‑68 (MQ=255)
                                     ttCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCTGa                       <  1:1305710/69‑1 (MQ=255)
                                     ttCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCGGa                       <  1:1799706/69‑1 (MQ=255)
                                     ttCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCGGa                       <  1:850359/69‑1 (MQ=255)
                                     ttCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCGGa                       <  1:735563/69‑1 (MQ=255)
                                     ttCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCGGa                       <  1:3147563/69‑1 (MQ=255)
                                     ttCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCGGa                       <  1:3055566/69‑1 (MQ=255)
                                     ttCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCGGa                       <  1:2081458/69‑1 (MQ=255)
                                     ttCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCGGa                       <  1:2137782/69‑1 (MQ=255)
                                     ttCTGACCAGCTGGCGCTGTTGCCGGTGACCCGGCGCTTTATCTTCCTTGAAGAGGGCGATATTGCGGa                       <  1:1629126/69‑1 (MQ=255)
                                         gtccAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:351404/66‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCTATATTGCGGAAAt                    <  1:3342961/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:96835/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:2757982/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:921628/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:2005128/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:1065610/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:749247/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:1090834/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:460925/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:3424931/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:2143396/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:3224163/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:2240800/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:1543691/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:1444273/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:1350524/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:1390042/68‑1 (MQ=255)
                                         gACCAGCTGGCGCTGTTGCCAGTGAATCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:1172987/68‑1 (MQ=255)
                                                        ttGCCAGTGACTCGTCGTTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAAt                    <  1:1720573/53‑1 (MQ=255)
                                                          gCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCgg                        >  1:1955627/1‑47 (MQ=255)
                                                          gCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCGGa                       <  1:2432381/48‑1 (MQ=255)
                                                          gCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAATCACTCGCCGTTCGGTaaa  >  1:3397595/1‑69 (MQ=255)
                                                          gCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAATCACTCGCCGTTCGGTaa   >  1:2328873/1‑68 (MQ=255)
                                                             |                                                                 
TGGTGATTGGCCTGGGGATGGGCGAAAACTTTATCGCTTCTGACCAGCTGGCGCTGTTGCCGGTGACCCGTCGCTTTATCTTCCTTGAAGAGGGCGATATTGCGGAAATCACTCGCCGTTCGGTAAA  >  minE/2473261‑2473387

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: