Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 75,688 C→T 23.8% G285G (GGC→GGT mraW → S‑adenosyl‑dependent methyltransferase activity on membrane‑located substrates

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE75,6880CT23.8% 57.0 / 11.9 46G285G (GGC→GGTmraWS‑adenosyl‑dependent methyltransferase activity on membrane‑located substrates
Reads supporting (aligned to +/- strand):  ref base C (8/27);  new base T (8/3);  total (16/30)
Fisher's exact test for biased strand distribution p-value = 4.31e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

AGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTG  >  minE/75625‑75755
                                                               |                                                                   
aGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGtt                                                                <  1:1016251/69‑1 (MQ=255)
aGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGtt                                                                <  1:9262/69‑1 (MQ=255)
aGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGtt                                                                <  1:478176/69‑1 (MQ=255)
aGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGtt                                                                <  1:1998776/69‑1 (MQ=255)
aGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGGGGCCGTCAGCTGCGAGCACTAGGCAAGtt                                                                <  1:2599885/69‑1 (MQ=255)
   gTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAAt                                                             >  1:1147992/1‑69 (MQ=255)
                   gagCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAagaaga                                              <  1:656917/68‑1 (MQ=255)
                   gagCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAagaaga                                              <  1:1348120/68‑1 (MQ=255)
                   gagCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAagaaga                                              <  1:3297905/68‑1 (MQ=255)
                   gagCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAagaaga                                              <  1:2586418/68‑1 (MQ=255)
                   gagCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAagaaga                                              <  1:1668319/68‑1 (MQ=255)
                    agcagcTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAagaaga                                              <  1:2047449/67‑1 (MQ=255)
                                  cTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCt                              <  1:1021973/69‑1 (MQ=255)
                                  cTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCt                              <  1:3376307/69‑1 (MQ=255)
                                       tGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCgaag                                                  <  1:134295/44‑1 (MQ=255)
                                       tGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCgaag                                                  <  1:2191269/44‑1 (MQ=255)
                                       tGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCgaag                                                  <  1:2311719/44‑1 (MQ=255)
                                       tGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCgaag                                                  <  1:1596012/44‑1 (MQ=255)
                                       tGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCgaag                                                  <  1:2808839/44‑1 (MQ=255)
                                                      aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt          >  1:433616/1‑69 (MQ=255)
                                                      aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt          >  1:3345369/1‑69 (MQ=255)
                                                      aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt          >  1:3181273/1‑69 (MQ=255)
                                                      aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt          >  1:2141932/1‑69 (MQ=255)
                                                      aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt          >  1:1883939/1‑69 (MQ=255)
                                                      aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt          >  1:2545293/1‑69 (MQ=255)
                                                      aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt          >  1:2564957/1‑69 (MQ=255)
                                                      aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt          >  1:1622714/1‑69 (MQ=255)
                                                         aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg       <  1:1120939/69‑1 (MQ=255)
                                                         aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg       <  1:1124964/69‑1 (MQ=255)
                                                         aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg       <  1:677758/69‑1 (MQ=255)
                                                         aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg       <  1:2803122/69‑1 (MQ=255)
                                                         aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg       <  1:2675053/69‑1 (MQ=255)
                                                         aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg       <  1:2214039/69‑1 (MQ=255)
                                                         aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg       <  1:772072/69‑1 (MQ=255)
                                                         aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg       <  1:1983054/69‑1 (MQ=255)
                                                         aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg       <  1:1931477/69‑1 (MQ=255)
                                                           tAGGTAAGTTGTTGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTTCGTa     <  1:3413695/69‑1 (MQ=255)
                                                           tAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTa     <  1:741753/69‑1 (MQ=255)
                                                           tAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTa     <  1:3234868/69‑1 (MQ=255)
                                                            aGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTa     >  1:748634/1‑68 (MQ=255)
                                                            aGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTAt    >  1:2914045/1‑69 (MQ=255)
                                                            aGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTAt    >  1:611732/1‑69 (MQ=255)
                                                            aGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTAt    >  1:508863/1‑69 (MQ=255)
                                                            aGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTAt    >  1:1496130/1‑69 (MQ=255)
                                                            aGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTAt    >  1:2396745/1‑69 (MQ=255)
                                                            aGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTAt    >  1:1782616/1‑69 (MQ=255)
                                                               cAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTg  <  1:1278982/68‑1 (MQ=255)
                                                               cAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTg  <  1:619105/68‑1 (MQ=255)
                                                               cAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTg  <  1:2836641/68‑1 (MQ=255)
                                                               cAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTg  <  1:730085/68‑1 (MQ=255)
                                                               cAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTg  <  1:1881734/68‑1 (MQ=255)
                                                               cAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTg  <  1:1907849/68‑1 (MQ=255)
                                                               |                                                                   
AGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTG  >  minE/75625‑75755

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: