Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,298,777 G→T 40.5% R171R (CGC→CGA metF ← 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,298,7770GT40.5% 22.6 / 28.8 37R171R (CGC→CGAmetF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (15/7);  new base T (5/10);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 4.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGACATCGAAGAAGAACTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGC  >  minE/2298724‑2298839
                                                     |                                                              
tCGACATCGAAGAAGAACTGAGTAATGGCGCGGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAg                                                  >  1:918890/1‑68 (MQ=255)
tCGACATCGAAGAAGAACTGAGTAATGGCGCGGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAg                                                  >  1:2795630/1‑68 (MQ=255)
tCGACATCGAAGAAGAACTGAGTAATGGCGCGGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAg                                                  >  1:1415470/1‑68 (MQ=255)
tCGACATCGAAGAAGAACTGAGTAATGGCGCGGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAg                                                  >  1:1424857/1‑68 (MQ=255)
tCGACATCGAAGAAGAACTGAGTAATGGCGCGGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAg                                                  >  1:1621915/1‑68 (MQ=255)
                 cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt                                >  1:1081488/1‑69 (MQ=255)
                 cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt                                >  1:868425/1‑69 (MQ=255)
                 cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt                                >  1:2909590/1‑69 (MQ=255)
                 cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt                                >  1:2516599/1‑69 (MQ=255)
                 cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt                                >  1:2348330/1‑69 (MQ=255)
                 cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt                                >  1:1900685/1‑69 (MQ=255)
                 cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt                                <  1:1806242/69‑1 (MQ=255)
                 cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt                                >  1:1047157/1‑69 (MQ=255)
                 cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt                                >  1:1569972/1‑69 (MQ=255)
                 cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGACTGAGCGCt                                >  1:1213463/1‑69 (MQ=255)
                   gAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCttt                              >  1:2767563/1‑69 (MQ=255)
                   gAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCttt                              >  1:3226356/1‑69 (MQ=255)
                   gAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCttt                              >  1:3148501/1‑69 (MQ=255)
                   gAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCttt                              >  1:2107294/1‑69 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                       <  1:1631364/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                       <  1:1889088/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                       <  1:3318847/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                       <  1:1639675/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                       <  1:1989293/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                       <  1:3104833/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                       <  1:2086624/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                       <  1:2783955/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                       <  1:177740/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                       <  1:2202705/68‑1 (MQ=255)
                                 ttGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGa                <  1:1289613/69‑1 (MQ=255)
                                  tGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGa                <  1:2375722/68‑1 (MQ=255)
                                          gCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGa       <  1:2720106/69‑1 (MQ=255)
                                          gCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGa       <  1:236746/69‑1 (MQ=255)
                                          gCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGa       <  1:439144/69‑1 (MQ=255)
                                            aTCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGa       <  1:525237/67‑1 (MQ=255)
                                               cACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATAcgc  >  1:2820335/1‑69 (MQ=255)
                                               cACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATAcgc  >  1:1201815/1‑69 (MQ=255)
                                                     |                                                              
TCGACATCGAAGAAGAACTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGC  >  minE/2298724‑2298839

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