Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,940,706 G→A 27.8% Q42Q (CAG→CAA ygfZ → predicted folate‑dependent regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,940,7060GA27.8% 93.6 / 34.9 86Q42Q (CAG→CAAygfZpredicted folate‑dependent regulatory protein
Reads supporting (aligned to +/- strand):  ref base G (36/26);  new base A (18/6);  total (54/32)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCA  >  minE/1940649‑1940765
                                                         |                                                           
gCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAc                                                  >  1:1234808/1‑69 (MQ=255)
gCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAc                                                  >  1:1521967/1‑69 (MQ=255)
gCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAc                                                  >  1:144865/1‑69 (MQ=255)
gCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAc                                                  >  1:1527790/1‑69 (MQ=255)
gCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAc                                                  >  1:432845/1‑69 (MQ=255)
gCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAc                                                  >  1:1317635/1‑69 (MQ=255)
gCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAc                                                  >  1:3351994/1‑69 (MQ=255)
gCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAc                                                  >  1:3074080/1‑69 (MQ=255)
gCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAc                                                  >  1:762774/1‑69 (MQ=255)
gCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGAAGGGTCa                                                        >  1:2578065/1‑63 (MQ=255)
gCTTGATGACTGGGCGCTTGCCACCATTAATGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAc                                                  >  1:2516048/1‑69 (MQ=255)
      tgaCTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGa                                            >  1:2257208/1‑69 (MQ=255)
           gggCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCa                                       >  1:2513820/1‑69 (MQ=255)
               gcTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   <  1:1878518/69‑1 (MQ=255)
               gcTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   <  1:3186942/69‑1 (MQ=255)
               gcTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   <  1:2543151/69‑1 (MQ=255)
               gcTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   <  1:1741650/69‑1 (MQ=255)
               gcTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   <  1:2690024/69‑1 (MQ=255)
               gcTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   <  1:2817538/69‑1 (MQ=255)
                  tgacaccaTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGAt                                >  1:3338139/4‑69 (MQ=255)
                  tcccaccaTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGAt                                >  1:3228132/3‑69 (MQ=255)
                  tGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGAt                                >  1:722400/1‑69 (MQ=255)
                  tGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGAt                                >  1:60604/1‑69 (MQ=255)
                  tGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGAt                                >  1:137717/1‑69 (MQ=255)
                  tGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGAt                                >  1:1339892/1‑69 (MQ=255)
                  tGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGAt                                >  1:2090013/1‑69 (MQ=255)
                  tGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGAt                                >  1:131452/1‑69 (MQ=255)
                       ccaTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   <  1:2152518/61‑1 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGcc                                    >  1:1082707/1‑58 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGcc                                    >  1:2301574/1‑58 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   >  1:1945061/1‑59 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   >  1:1478287/1‑59 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   >  1:885405/1‑59 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   >  1:3003436/1‑59 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   >  1:1179592/1‑59 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   >  1:2498157/1‑59 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   >  1:2399095/1‑59 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   >  1:1728813/1‑59 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   >  1:1902454/1‑59 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   >  1:346992/1‑59 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   >  1:1898600/1‑59 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   >  1:1877705/1‑59 (MQ=255)
                         aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   >  1:3344080/1‑59 (MQ=255)
                          ttACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCagaag                        <  1:1162579/69‑1 (MQ=255)
                               ggCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAg                   <  1:1396840/69‑1 (MQ=255)
                               ggCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAg                   <  1:2672796/69‑1 (MQ=255)
                               ggCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAg                   <  1:116570/69‑1 (MQ=255)
                               ggCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAg                   <  1:497588/69‑1 (MQ=255)
                               ggCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAg                   <  1:177899/69‑1 (MQ=255)
                               ggCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAg                   <  1:1258583/69‑1 (MQ=255)
                               ggCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAg                   <  1:557450/69‑1 (MQ=255)
                               ggCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAg                   <  1:3275531/69‑1 (MQ=255)
                               ggCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAg                   <  1:2191651/69‑1 (MQ=255)
                                gcgcGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAg                   <  1:1554521/68‑1 (MQ=255)
                                  gcgGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCa                 <  1:2319841/68‑1 (MQ=255)
                                  gcgGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCa                 <  1:3143996/68‑1 (MQ=255)
                                  gcgGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCa                 <  1:2853765/68‑1 (MQ=255)
                                  gcgGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCa                 <  1:2703698/68‑1 (MQ=255)
                                  gcgGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCa                 <  1:1964884/68‑1 (MQ=255)
                                  gcgGACAGCGAAAAATATATGCAGGGTCAGGTGACACCAGATGTCAGCCAGATGGCAGAAGATCAGCa                 <  1:2197716/68‑1 (MQ=255)
                                   cgGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCa                 <  1:2148002/67‑1 (MQ=255)
                                       cAGCGAAAAATATATGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCt           <  1:1101832/69‑1 (MQ=255)
                                       cAGCGAAAAATATATGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCt           <  1:872869/69‑1 (MQ=255)
                                       cAGCGAAAAATATATGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCt           <  1:2523523/69‑1 (MQ=255)
                                       cAGCGAAAAATATATGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCt           <  1:798537/69‑1 (MQ=255)
                                       cAGCGAAAAATATATGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCt           <  1:1146191/69‑1 (MQ=255)
                                       cAGCGAAAAATATATGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCt           <  1:1299593/69‑1 (MQ=255)
                                            aaaaaTATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                   <  1:396915/40‑1 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTcgccgc     >  1:3123724/1‑66 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGccc   >  1:401785/1‑68 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:12029/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:969133/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:2997910/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:2590218/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:7788/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:77759/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:2442312/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:1179301/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:2044883/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:1032034/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:1468270/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:45782/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:1682378/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:2788958/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:2129011/1‑69 (MQ=255)
                                                atatatGCAAGGCCAGGTGACAGCAGATGTCAGCCAGATGACTGAAGATCAGCACCTGCTCGCCGCCCa  >  1:1973466/1‑69 (MQ=255)
                                                         |                                                           
GCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCA  >  minE/1940649‑1940765

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: