Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,298,781 T→C 40.5% K170R (AAA→AGA)  metF ← 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,298,7810TC40.5% 10.3 / 23.3 37K170R (AAA→AGA) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (15/7);  new base C (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

TCGACATCGAAGAAGAACTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGGA  >  minE/2298724‑2298848
                                                         |                                                                   
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:1639675/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:177740/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                <  1:1889088/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                <  1:1631364/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                <  1:2202705/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                <  1:1989293/68‑1 (MQ=255)
                          ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc                                <  1:2086624/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:2720106/69‑1 (MQ=255)
                                          gCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGa                <  1:236746/69‑1 (MQ=255)
                                            aTCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGa                <  1:525237/67‑1 (MQ=255)
                                               cACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATAcgc           >  1:1201815/1‑69 (MQ=255)
                                               cACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATAcgc           >  1:2820335/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:3355776/69‑1 (MQ=255)
                                                        tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGga  <  1:262072/69‑1 (MQ=255)
                                                        tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGga  <  1:2472799/69‑1 (MQ=255)
                                                        tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGga  <  1:2405833/69‑1 (MQ=255)
                                                        tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTCTTCCGGGTGAACTTCCGGATAcgccg         >  1:1154703/1‑61 (MQ=255)
                                                         |                                                                   
TCGACATCGAAGAAGAACTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGGA  >  minE/2298724‑2298848

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: