Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,780,711 A→T 20.0% V562V (GTA→GTT metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,780,7110AT20.0% 90.9 / 16.0 50V562V (GTA→GTTmetHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base A (32/8);  new base T (5/5);  total (37/13)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTC  >  minE/2780647‑2780777
                                                                |                                                                  
aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                                                >  1:1076592/1‑69 (MQ=255)
aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                                                >  1:2864585/1‑69 (MQ=255)
aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                                                >  1:1044215/1‑69 (MQ=255)
aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                                                >  1:2511895/1‑69 (MQ=255)
aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                                                >  1:2342857/1‑69 (MQ=255)
aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                                                >  1:1318247/1‑69 (MQ=255)
 cTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                                                >  1:633209/1‑68 (MQ=255)
           gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc                                                     >  1:3140688/1‑69 (MQ=255)
           gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc                                                     >  1:370677/1‑69 (MQ=255)
           gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc                                                     >  1:2073586/1‑69 (MQ=255)
           gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc                                                     >  1:3175312/1‑69 (MQ=255)
           gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc                                                     >  1:1694831/1‑69 (MQ=255)
                  aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGGGTATCTAACGTTTCTTTCTCGTTcc                                              >  1:224672/1‑69 (MQ=255)
                  aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc                                              >  1:1930512/1‑69 (MQ=255)
                  aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc                                              >  1:19044/1‑69 (MQ=255)
                  aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc                                              >  1:3426869/1‑69 (MQ=255)
                  aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc                                              >  1:1262110/1‑69 (MQ=255)
                  aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc                                              >  1:2925311/1‑69 (MQ=255)
                  aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc                                              >  1:1649691/1‑69 (MQ=255)
                  aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc                                              >  1:259632/1‑69 (MQ=255)
                  aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc                                              >  1:2654285/1‑69 (MQ=255)
                        tCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCa                                        <  1:554944/69‑1 (MQ=255)
                          aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCt                                                                 >  1:2839176/1‑42 (MQ=255)
                          aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAAc                                      >  1:1553536/1‑69 (MQ=255)
                          aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAAc                                      >  1:3007848/1‑69 (MQ=255)
                          aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAAc                                      >  1:2388940/1‑69 (MQ=255)
                          aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAAc                                      >  1:2115145/1‑69 (MQ=255)
                              gtgAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGTGTCGTTCTCGTTTCGTGGCAACGATc                                  <  1:1229258/67‑1 (MQ=255)
                              gtgAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGTGTCGTTCTCGTTTCGTGGCAACGATc                                  <  1:517945/67‑1 (MQ=255)
                              gtgAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGTGTCGTTCTCGTTTCGTGGCAACGATc                                  <  1:2103134/67‑1 (MQ=255)
                              gtgAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGTGTCGTTCTCGTTTCGTGGCAACGATc                                  <  1:646523/67‑1 (MQ=255)
                                     gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg                                     >  1:699103/1‑59 (MQ=255)
                                     gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg                                     >  1:376749/1‑59 (MQ=255)
                                     gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg                                     >  1:268059/1‑59 (MQ=255)
                                     gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg                                     >  1:3299863/1‑59 (MQ=255)
                                     gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg                                     >  1:3255811/1‑59 (MQ=255)
                                     gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg                                     >  1:3253766/1‑59 (MQ=255)
                                     gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg                                     >  1:3028480/1‑59 (MQ=255)
                                     gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg                                     >  1:2461658/1‑59 (MQ=255)
                                     gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg                                     >  1:2383541/1‑59 (MQ=255)
                                     gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg                                     >  1:198392/1‑59 (MQ=255)
                                     gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg                                     >  1:1123011/1‑59 (MQ=255)
                                                           ggcgTTTCTAACGTGTCGTTCTCGTTTCGTGGCAACGATc                                  <  1:2005430/40‑1 (MQ=255)
                                                              gTATCTAACGTTTCTTTCTCGTTCCTTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc  <  1:3135273/69‑1 (MQ=255)
                                                              gTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc  <  1:3258733/69‑1 (MQ=255)
                                                              gTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc  <  1:2771015/69‑1 (MQ=255)
                                                              gTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc  <  1:1996381/69‑1 (MQ=255)
                                                              gTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc  <  1:68335/69‑1 (MQ=255)
                                                              gTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc  <  1:939471/69‑1 (MQ=255)
                                                              gTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc  <  1:992675/69‑1 (MQ=255)
                                                                |                                                                  
ACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTC  >  minE/2780647‑2780777

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: