Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,298,789 A→G 29.9% L167L (CTT→CTC metF ← 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,298,7890AG29.9% 43.7 / 21.0 50L167L (CTT→CTCmetF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (24/11);  new base G (5/10);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 2.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TCGACATCGAAGAAGAACTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGGAGATATCGAA  >  minE/2298724‑2298857
                                                                 |                                                                    
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:3226356/1‑69 (MQ=255)
                   gAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCttt                                                >  1:3148501/1‑69 (MQ=255)
                   gAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCttt                                                >  1:2107294/1‑69 (MQ=255)
                   gAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCttt                                                >  1:2767563/1‑69 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                         <  1:177740/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                         <  1:3318847/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                         <  1:3104833/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                         <  1:2783955/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                         <  1:2086624/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                         <  1:1639675/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                         <  1:1631364/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                         <  1:1889088/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                         <  1:1989293/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:439144/69‑1 (MQ=255)
                                          gCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGa                         <  1:236746/69‑1 (MQ=255)
                                          gCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGa                         <  1:2720106/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)
                                                        tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTcc                                        >  1:2730277/1‑40 (MQ=255)
                                                        tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCa               >  1:2231394/1‑65 (MQ=255)
                                                        tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCa               >  1:391303/1‑65 (MQ=255)
                                                        tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGga           <  1:2405833/69‑1 (MQ=255)
                                                        tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGga           <  1:262072/69‑1 (MQ=255)
                                                        tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGga           <  1:3355776/69‑1 (MQ=255)
                                                        tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGga           <  1:2472799/69‑1 (MQ=255)
                                                        tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTCTTCCGGGTGAACTTCCGGATAcgccg                  >  1:1154703/1‑61 (MQ=255)
                                                          tCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGgaga         >  1:878230/1‑69 (MQ=255)
                                                          tCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGgaga         >  1:2061507/1‑69 (MQ=255)
                                                          tCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGgaga         >  1:2391391/1‑69 (MQ=255)
                                                          tCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGgaga         >  1:2448704/1‑69 (MQ=255)
                                                          tCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGgaga         >  1:2326777/1‑69 (MQ=255)
                                                                 aaGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGGAGATATCGaa  <  1:1380333/69‑1 (MQ=255)
                                                                 |                                                                    
TCGACATCGAAGAAGAACTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGGAGATATCGAA  >  minE/2298724‑2298857

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: