Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,781,293 C→T 43.1% G756G (GGC→GGT 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,2930CT43.1% 34.2 / 77.3 102G756G (GGC→GGTmetHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base C (33/25);  new base T (27/17);  total (60/42)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

CAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGAA  >  minE/2781254‑2781353
                                       |                                                            
cAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGt                                 >  1:2917619/1‑69 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCTCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:543048/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:2629330/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:3114169/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:3152992/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:3018241/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:3288766/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:3293811/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:1478327/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:3333860/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:2803622/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:1989583/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:710277/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:798481/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:2008655/69‑1 (MQ=255)
 aaaaCCAACGGCAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:2375192/69‑1 (MQ=255)
    aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                             >  1:2109498/1‑69 (MQ=255)
    aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                             >  1:2583204/1‑69 (MQ=255)
    aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                             >  1:2376366/1‑69 (MQ=255)
    aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                             >  1:2029106/1‑69 (MQ=255)
    aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                             >  1:3252308/1‑69 (MQ=255)
    aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                             >  1:525118/1‑69 (MQ=255)
    aCCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTtggt                             >  1:868427/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:922153/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:913321/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:70020/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:613774/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:2472014/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:2797802/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:2798400/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:284385/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:3395826/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:1540479/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:180587/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:1720866/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:1367746/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:1873263/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:1881839/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:196597/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:1644986/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:1600425/1‑68 (MQ=255)
            cAAGATGGTGATCGCCACAGTGAAGGGTGATGTCCACGACATCGGTAAAAATATCGTTGGCGTGGtgc                      >  1:1111979/1‑68 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:2992226/1‑69 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:2987444/1‑69 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:2972453/1‑69 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:471327/1‑69 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:14074/1‑69 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:141455/1‑69 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:3139874/1‑69 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:2122956/1‑69 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:1817206/1‑69 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:1253976/1‑69 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:2634497/1‑69 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:2619079/1‑69 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:2019161/1‑69 (MQ=255)
                aTGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCaa                 >  1:793766/1‑69 (MQ=255)
                   gTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGTTAAAAATATCGTTGGTGTGGTGCTGCAATGt              <  1:2664519/69‑1 (MQ=255)
                   gTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGTTAAAAATATCGTTGGTGTGGTGCTGCAATGt              <  1:613811/69‑1 (MQ=255)
                   gTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGTTAAAAATATCGTTGGTGTGGTGCTGCAATGt              <  1:672394/69‑1 (MQ=255)
                   gTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGTTAAAAATATCGTTGGTGTGGTGCTGCAATGt              <  1:144999/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:2249387/69‑1 (MQ=255)
                      aTCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaac           <  1:2061081/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:2332191/1‑68 (MQ=255)
                        cGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaaca          >  1:707674/1‑68 (MQ=255)
                        cGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaaca          >  1:2759918/1‑68 (MQ=255)
                        cGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaaca          >  1:869536/1‑68 (MQ=255)
                        cGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaaca          >  1:1006713/1‑68 (MQ=255)
                        cGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTaaca          >  1:1534029/1‑68 (MQ=255)
                        cGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGATGCAATGTaaga          >  1:320635/1‑66 (MQ=255)
                               gTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGaa  <  1:1069672/69‑1 (MQ=255)
                               gTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGaa  <  1:2058141/69‑1 (MQ=255)
                               gTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGaa  <  1:100559/69‑1 (MQ=255)
                               gTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGaa  <  1:1818401/69‑1 (MQ=255)
                               gTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGaa  <  1:655137/69‑1 (MQ=255)
                               gTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGaa  >  1:2893915/1‑69 (MQ=255)
                               gTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGaa  <  1:1693776/69‑1 (MQ=255)
                               gTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGaa  <  1:3116707/69‑1 (MQ=255)
                               gTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGaa  <  1:1496048/69‑1 (MQ=255)
                               gTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGaa  <  1:3232827/69‑1 (MQ=255)
                               gTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGa   >  1:1017287/1‑68 (MQ=255)
                                  aaGGGTGATGTCCACGACATCGGTAAAAATATCGtt                                <  1:1542078/36‑1 (MQ=37)
                                  aaGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGt              <  1:2704432/54‑1 (MQ=255)
                                       |                                                            
CAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGAA  >  minE/2781254‑2781353

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: