Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,780,693 G→A 26.8% A556A (GCG→GCA metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,780,6930GA26.8% 88.3 / 28.9 78A556A (GCG→GCAmetHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base G (46/11);  new base A (17/4);  total (63/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATC  >  minE/2780630‑2780745
                                                               |                                                    
caacaaCTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGAt                                                 >  1:578325/1‑69 (MQ=255)
caacaaCTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGAt                                                 >  1:293516/1‑69 (MQ=255)
   caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGAttcc                                              >  1:1532143/1‑67 (MQ=255)
   caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc                                              >  1:2275530/1‑69 (MQ=255)
   caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc                                              >  1:1630963/1‑69 (MQ=255)
   caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc                                              >  1:1578533/1‑69 (MQ=255)
   caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc                                              >  1:3146101/1‑69 (MQ=255)
   caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc                                              >  1:1987524/1‑69 (MQ=255)
   caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc                                              >  1:526380/1‑69 (MQ=255)
   caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc                                              >  1:570884/1‑69 (MQ=255)
   caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc                                              >  1:2004006/1‑69 (MQ=255)
   caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc                                              >  1:1090750/1‑69 (MQ=255)
   caatTACGCACAGGACTTTATCGGCGCGTATGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc                                              >  1:1145402/1‑69 (MQ=255)
     acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg                                            >  1:1955838/1‑69 (MQ=255)
     acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg                                            >  1:2106008/1‑69 (MQ=255)
     acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg                                            >  1:2333671/1‑69 (MQ=255)
     acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg                                            >  1:1008746/1‑69 (MQ=255)
     acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg                                            >  1:1879428/1‑69 (MQ=255)
     acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg                                            >  1:2988190/1‑69 (MQ=255)
     acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg                                            >  1:1482367/1‑69 (MQ=255)
        aCGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTgcta                                                >  1:2318153/1‑59 (MQ=255)
        aCGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGAtt                                                >  1:167108/1‑62 (MQ=255)
        aCGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGAtt                                                >  1:2217370/1‑62 (MQ=255)
        aCGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGAtt                                                >  1:2625434/1‑62 (MQ=255)
        aCGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGAtt                                                >  1:3424644/1‑62 (MQ=255)
           cgcATGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcggcgg                                        <  1:1706355/67‑1 (MQ=255)
           cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcggcgg                                        <  1:2877195/67‑1 (MQ=255)
           cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcggcgg                                        <  1:1147613/67‑1 (MQ=255)
           cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcggcgg                                        <  1:3168474/67‑1 (MQ=255)
           cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcggcgg                                        <  1:646857/67‑1 (MQ=255)
           cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGgcggcg                                      <  1:1993706/69‑1 (MQ=255)
           cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGgcggcg                                      <  1:1215262/69‑1 (MQ=255)
           cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGgcggcg                                      <  1:306019/69‑1 (MQ=255)
            gcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGgcggcg                                      <  1:2142107/68‑1 (MQ=255)
                 aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                >  1:2864585/1‑69 (MQ=255)
                 aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                >  1:2511895/1‑69 (MQ=255)
                 aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                >  1:2342857/1‑69 (MQ=255)
                 aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                >  1:1044215/1‑69 (MQ=255)
                 aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                >  1:1076592/1‑69 (MQ=255)
                 aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                >  1:1318247/1‑69 (MQ=255)
                  cTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa                                >  1:633209/1‑68 (MQ=255)
                            gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCcggcgg                                        >  1:1524191/1‑50 (MQ=255)
                            gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc                     >  1:3140688/1‑69 (MQ=255)
                            gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc                     >  1:3175312/1‑69 (MQ=255)
                            gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc                     >  1:2073586/1‑69 (MQ=255)
                            gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc                     >  1:1694831/1‑69 (MQ=255)
                            gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc                     >  1:370677/1‑69 (MQ=255)
                              gtgtGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGgcggcg                                      <  1:2772742/50‑1 (MQ=255)
                                   aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGGGTATCTAACGTTTCTTTCTCGTTcc              >  1:224672/1‑69 (MQ=255)
                                   aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc              >  1:3426869/1‑69 (MQ=255)
                                   aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc              >  1:1930512/1‑69 (MQ=255)
                                   aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc              >  1:2925311/1‑69 (MQ=255)
                                   aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc              >  1:19044/1‑69 (MQ=255)
                                   aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc              >  1:2654285/1‑69 (MQ=255)
                                   aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc              >  1:259632/1‑69 (MQ=255)
                                   aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc              >  1:1262110/1‑69 (MQ=255)
                                   aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc              >  1:1649691/1‑69 (MQ=255)
                                         tCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCa        <  1:554944/69‑1 (MQ=255)
                                           aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCt                                 >  1:2839176/1‑42 (MQ=255)
                                           aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAAc      >  1:2115145/1‑69 (MQ=255)
                                           aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAAc      >  1:1553536/1‑69 (MQ=255)
                                           aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAAc      >  1:2388940/1‑69 (MQ=255)
                                           aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAAc      >  1:3007848/1‑69 (MQ=255)
                                               gtgAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGTGTCGTTCTCGTTTCGTGGCAACGATc  <  1:646523/67‑1 (MQ=255)
                                               gtgAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGTGTCGTTCTCGTTTCGTGGCAACGATc  <  1:1229258/67‑1 (MQ=255)
                                               gtgAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGTGTCGTTCTCGTTTCGTGGCAACGATc  <  1:517945/67‑1 (MQ=255)
                                               gtgAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGTGTCGTTCTCGTTTCGTGGCAACGATc  <  1:2103134/67‑1 (MQ=255)
                                                      gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg     >  1:699103/1‑59 (MQ=255)
                                                      gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg     >  1:2383541/1‑59 (MQ=255)
                                                      gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg     >  1:1123011/1‑59 (MQ=255)
                                                      gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg     >  1:2461658/1‑59 (MQ=255)
                                                      gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg     >  1:376749/1‑59 (MQ=255)
                                                      gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg     >  1:268059/1‑59 (MQ=255)
                                                      gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg     >  1:198392/1‑59 (MQ=255)
                                                      gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg     >  1:3299863/1‑59 (MQ=255)
                                                      gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg     >  1:3255811/1‑59 (MQ=255)
                                                      gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg     >  1:3253766/1‑59 (MQ=255)
                                                      gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg     >  1:3028480/1‑59 (MQ=255)
                                                               |                                                    
CAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATC  >  minE/2780630‑2780745

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: