Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,781,284 C→A 38.7% T753T (ACC→ACA metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,781,2840CA38.7% 90.7 / 91.3 111T753T (ACC→ACAmetHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base C (45/23);  new base A (27/16);  total (72/39)
Fisher's exact test for biased strand distribution p-value = 8.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 111 total aligned reads displayed.

CCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACA  >  minE/2781218‑2781345
                                                                  |                                                             
ccTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCTCCACCg                                                              >  1:3226619/1‑68 (MQ=255)
ccTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCg                                                              >  1:1886705/1‑68 (MQ=255)
ccTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCg                                                              >  1:2731142/1‑68 (MQ=255)
ccTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCg                                                              >  1:2399435/1‑68 (MQ=255)
ccTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCg                                                              >  1:2948826/1‑68 (MQ=255)
ccTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCg                                                              >  1:2100121/1‑68 (MQ=255)
ccTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCg                                                              >  1:3153762/1‑68 (MQ=255)
ccTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCg                                                              >  1:2014725/1‑68 (MQ=255)
ccTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCg                                                              >  1:3309781/1‑68 (MQ=255)
ccTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCg                                                              >  1:1498104/1‑68 (MQ=255)
ccTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCg                                                              >  1:663323/1‑68 (MQ=255)
ccTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCg                                                              >  1:1342321/1‑68 (MQ=255)
ccTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCg                                                              >  1:1298314/1‑68 (MQ=255)
  tCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTg                                                            <  1:2300927/68‑1 (MQ=255)
  tCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTg                                                            <  1:2866063/68‑1 (MQ=255)
  tCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTg                                                            <  1:3163603/68‑1 (MQ=255)
  tCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTg                                                            <  1:1274379/68‑1 (MQ=255)
  tCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTg                                                            <  1:3392662/68‑1 (MQ=255)
  tCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTg                                                            <  1:892201/68‑1 (MQ=255)
  tCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTg                                                            <  1:1482334/68‑1 (MQ=255)
  tCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTg                                                            <  1:2900070/68‑1 (MQ=255)
                                    cAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGt                         >  1:2917619/1‑69 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCTCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:543048/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:2008655/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:3114169/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:3018241/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:3152992/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:2375192/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:3333860/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:2803622/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:1478327/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:710277/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:2629330/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:1989583/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:798481/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:3293811/69‑1 (MQ=255)
                                     aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                        <  1:3288766/69‑1 (MQ=255)
                                        aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                     >  1:868427/1‑69 (MQ=255)
                                        aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                     >  1:2109498/1‑69 (MQ=255)
                                        aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                     >  1:2583204/1‑69 (MQ=255)
                                        aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                     >  1:525118/1‑69 (MQ=255)
                                        aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                     >  1:2029106/1‑69 (MQ=255)
                                        aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                     >  1:2376366/1‑69 (MQ=255)
                                        aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                     >  1:3252308/1‑69 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGAc                                          >  1:584638/1‑40 (MQ=37)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCgg                                     >  1:3041493/1‑45 (MQ=38)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAatat                            >  1:1061750/1‑54 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTtggc                     >  1:1409059/1‑60 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:913321/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:922153/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:3395826/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:284385/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:613774/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:2798400/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:70020/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:2797802/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:1600425/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:1644986/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:1540479/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:1720866/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:180587/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:1873263/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:1881839/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:2133585/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:2291267/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:2316646/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:2570231/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:1111979/1‑68 (MQ=255)
                                                cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc              >  1:1367746/1‑68 (MQ=255)
                                                    aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa         >  1:793766/1‑69 (MQ=255)
                                                    aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa         >  1:2634497/1‑69 (MQ=255)
                                                    aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa         >  1:2619079/1‑69 (MQ=255)
                                                    aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa         >  1:1817206/1‑69 (MQ=255)
                                                    aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa         >  1:1253976/1‑69 (MQ=255)
                                                    aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa         >  1:2019161/1‑69 (MQ=255)
                                                    aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa         >  1:2972453/1‑69 (MQ=255)
                                                    aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa         >  1:3139874/1‑69 (MQ=255)
                                                    aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa         >  1:141455/1‑69 (MQ=255)
                                                    aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa         >  1:2992226/1‑69 (MQ=255)
                                                       gTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGTTAAAAATATCGTTGGTGTGGTGCTGCAATGt      <  1:613811/69‑1 (MQ=255)
                                                       gTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGTTAAAAATATCGTTGGTGTGGTGCTGCAATGt      <  1:144999/69‑1 (MQ=255)
                                                       gTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGTTAAAAATATCGTTGGTGTGGTGCTGCAATGt      <  1:672394/69‑1 (MQ=255)
                                                       gTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGTTAAAAATATCGTTGGTGTGGTGCTGCAATGt      <  1:2664519/69‑1 (MQ=255)
                                                       gTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGTTAAAAATATCGTTGGTGTGGTGCTGCAATGt      <  1:3022979/69‑1 (MQ=255)
                                                       gTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGTTGCAATGt      <  1:1233763/69‑1 (MQ=255)
                                                       gTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGt      <  1:593151/69‑1 (MQ=255)
                                                       gTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGt      <  1:1631860/69‑1 (MQ=255)
                                                       gTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGt      <  1:2472238/69‑1 (MQ=255)
                                                        tGATCGCCAGCGTGAAGGGCGACGTCCACGACTTCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGt      <  1:799566/68‑1 (MQ=255)
                                                        tGATCGCCACCGTGAAGGGCGACGTCCACGACATCGTTAAAAATATCGTTGGTGTGGTGCTGCAATGt      <  1:2097143/68‑1 (MQ=255)
                                                          aTCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaac   <  1:1705154/69‑1 (MQ=255)
                                                          aTCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaac   <  1:2061081/69‑1 (MQ=255)
                                                          aTCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaac   <  1:2249387/69‑1 (MQ=255)
                                                          aTCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGAGGTGCTGCAATGTaac   <  1:2536984/69‑1 (MQ=255)
                                                            cgcgACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaaca  >  1:3211761/1‑68 (MQ=255)
                                                            cGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                     >  1:1393344/1‑49 (MQ=255)
                                                            cGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaaca  >  1:869536/1‑68 (MQ=255)
                                                            cGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaaca  >  1:2332191/1‑68 (MQ=255)
                                                            cGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaaca  >  1:1534029/1‑68 (MQ=255)
                                                            cGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaaca  >  1:1006713/1‑68 (MQ=255)
                                                                  |                                                             
CCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACA  >  minE/2781218‑2781345

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: