Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 75,694 A→G 16.8% L287L (TTA→TTG 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,6940AG16.8% 101.5 / 11.0 65L287L (TTA→TTGmraWS‑adenosyl‑dependent methyltransferase activity on membrane‑located substrates
Reads supporting (aligned to +/- strand):  ref base A (14/40);  new base G (8/3);  total (22/43)
Fisher's exact test for biased strand distribution p-value = 4.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

GTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAG  >  minE/75628‑75762
                                                                  |                                                                    
gTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAAt                                                                    >  1:1147992/1‑69 (MQ=255)
                gagCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAagaaga                                                     <  1:2586418/68‑1 (MQ=255)
                gagCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAagaaga                                                     <  1:1348120/68‑1 (MQ=255)
                gagCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAagaaga                                                     <  1:656917/68‑1 (MQ=255)
                gagCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAagaaga                                                     <  1:3297905/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:2311719/44‑1 (MQ=255)
                                    tGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCgaag                                                         <  1:2191269/44‑1 (MQ=255)
                                    tGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCgaag                                                         <  1:2808839/44‑1 (MQ=255)
                                    tGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCgaag                                                         <  1:1596012/44‑1 (MQ=255)
                                    tGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCgaag                                                         <  1:134295/44‑1 (MQ=255)
                                                   aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt                 >  1:1622714/1‑69 (MQ=255)
                                                   aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt                 >  1:3345369/1‑69 (MQ=255)
                                                   aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt                 >  1:3181273/1‑69 (MQ=255)
                                                   aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt                 >  1:433616/1‑69 (MQ=255)
                                                   aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt                 >  1:1883939/1‑69 (MQ=255)
                                                   aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt                 >  1:2564957/1‑69 (MQ=255)
                                                   aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt                 >  1:2545293/1‑69 (MQ=255)
                                                   aGCACTAGGTAAGTTGATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt                 >  1:2141932/1‑69 (MQ=255)
                                                      aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg              <  1:2803122/69‑1 (MQ=255)
                                                      aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg              <  1:2675053/69‑1 (MQ=255)
                                                      aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg              <  1:677758/69‑1 (MQ=255)
                                                      aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg              <  1:772072/69‑1 (MQ=255)
                                                      aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg              <  1:1120939/69‑1 (MQ=255)
                                                      aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg              <  1:2214039/69‑1 (MQ=255)
                                                      aCTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCg              <  1:1124964/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:1782616/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:2914045/1‑69 (MQ=255)
                                                         aGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTAt           >  1:2396745/1‑69 (MQ=255)
                                                            cAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTg         <  1:1881734/68‑1 (MQ=255)
                                                            cAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTg         <  1:730085/68‑1 (MQ=255)
                                                            cAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTg         <  1:2836641/68‑1 (MQ=255)
                                                            cAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTg         <  1:1278982/68‑1 (MQ=255)
                                                            cAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTg         <  1:619105/68‑1 (MQ=255)
                                                            cAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTg         <  1:1907849/68‑1 (MQ=255)
                                                                ttAATGCCGGGCGAAGAAGAGGTGGGTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    <  1:2475714/69‑1 (MQ=255)
                                                                ttAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    <  1:2334676/69‑1 (MQ=255)
                                                                ttAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    <  1:2401896/69‑1 (MQ=255)
                                                                ttAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    <  1:1373503/69‑1 (MQ=255)
                                                                ttAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    <  1:2737640/69‑1 (MQ=255)
                                                                ttAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    <  1:3381782/69‑1 (MQ=255)
                                                                ttAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    <  1:1722534/69‑1 (MQ=255)
                                                                ttAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    <  1:3234666/69‑1 (MQ=255)
                                                                ttAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    <  1:2259576/69‑1 (MQ=255)
                                                                ttAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    <  1:301406/69‑1 (MQ=255)
                                                                ttAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    <  1:2733675/69‑1 (MQ=255)
                                                                ttAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    <  1:1863514/69‑1 (MQ=255)
                                                                 tAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    >  1:843113/1‑68 (MQ=255)
                                                                 tAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    >  1:2675726/1‑68 (MQ=255)
                                                                 tAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    >  1:1442243/1‑68 (MQ=255)
                                                                 tAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    >  1:249619/1‑68 (MQ=255)
                                                                 tAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    >  1:2807230/1‑68 (MQ=255)
                                                                 tAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCACGTGCCCGTAGTTCAGTTCTGCGTATTGCagag    >  1:342775/1‑68 (MQ=255)
                                                                  aaTGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagagag  >  1:871729/1‑69 (MQ=255)
                                                                  aaTGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagagag  >  1:197927/1‑69 (MQ=255)
                                                                  aaTGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagagag  >  1:788339/1‑69 (MQ=255)
                                                                  aaTGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagagag  >  1:1712235/1‑69 (MQ=255)
                                                                  aaTGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCagagag  >  1:2780851/1‑69 (MQ=255)
                                                                  |                                                                    
GTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAG  >  minE/75628‑75762

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: