Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,801,244 G→A 41.2% E234E (GAG→GAA yjbN → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,801,2440GA41.2% 40.5 / 68.4 97E234E (GAG→GAAyjbNtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (33/24);  new base A (12/28);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 7.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

CAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTAT  >  minE/2801179‑2801312
                                                                 |                                                                    
cAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGc                                                                    >  1:731375/1‑68 (MQ=255)
cAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGc                                                                    >  1:433315/1‑68 (MQ=255)
cAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGc                                                                    >  1:1128619/1‑68 (MQ=255)
cAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGc                                                                    >  1:263252/1‑68 (MQ=255)
cAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGc                                                                    >  1:2319883/1‑68 (MQ=255)
cAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGc                                                                    >  1:1896066/1‑68 (MQ=255)
cAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGc                                                                    >  1:1845820/1‑68 (MQ=255)
 aaCTGAAGCGTGACTTTCCGCATCTGACGATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGc                                                                    >  1:2432235/1‑67 (MQ=255)
 aaCTGAAGCGTGACTTTCCGCATCTGACGATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGc                                                                    >  1:2856267/1‑67 (MQ=255)
 aaCTGAAGCGTGACTTTCCGCATCTGACGATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGc                                                                    >  1:2561131/1‑67 (MQ=255)
 aaCTGAAGCGTGACTTTCCGCATCTGACGATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGc                                                                    >  1:2456336/1‑67 (MQ=255)
 aaCTGAAGCGTGACTTTCCGCATCTGACGATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGc                                                                    >  1:2357756/1‑67 (MQ=255)
 aaCTGAAGCGTGACTTTCCGCATCTGACGATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGc                                                                    >  1:2246121/1‑67 (MQ=255)
 aaCTGAAGCGTGACTTTCCGCATCTGACGATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGc                                                                    >  1:2185675/1‑67 (MQ=255)
 aaCTGAAGCGTGACTTTCCGCATCTGACGATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGc                                                                    >  1:1748677/1‑67 (MQ=255)
 aaCTGAAGCGTGACTTTCCGCATCTGACGATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGc                                                                    >  1:1447024/1‑67 (MQ=255)
 aaCTGAAGCGTGACTTTCCGCATCTGACGATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGc                                                                    >  1:126101/1‑67 (MQ=255)
 aaCTGAAGCGTGACTTTCCGCATCTGACGATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGc                                                                    >  1:1410415/1‑67 (MQ=255)
 aaCTGAAGCGTGACTTTCCGCATCTGACGATGTCGATTAACGGTGGTATCAAGTCGCAGGAAGAAGc                                                                    >  1:1292528/1‑67 (MQ=255)
             agtttCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                      <  1:2032732/67‑1 (MQ=255)
             aCTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGACGAGGCCAAACCACACCTGc                                                      <  1:3240817/69‑1 (MQ=255)
             aCTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCAGACCTGc                                                      <  1:3822/69‑1 (MQ=255)
             aCTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                      <  1:2800761/69‑1 (MQ=255)
             aCTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                      <  1:3031865/69‑1 (MQ=255)
             aCTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                      <  1:3093271/69‑1 (MQ=255)
             aCTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                      <  1:3182806/69‑1 (MQ=255)
             aCTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                      <  1:2897883/69‑1 (MQ=255)
             aCTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                      <  1:803688/69‑1 (MQ=255)
             aCTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                      <  1:346037/69‑1 (MQ=255)
             aCTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                      <  1:26848/69‑1 (MQ=255)
             aCTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                      <  1:1834428/69‑1 (MQ=255)
             aCTTTCCCCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                      <  1:3106574/69‑1 (MQ=255)
                 tCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAAc                                                   >  1:1161896/1‑68 (MQ=255)
                 tCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACa                                                  >  1:3398621/1‑69 (MQ=255)
                 tCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACa                                                  >  1:158404/1‑69 (MQ=255)
                 tCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACa                                                  >  1:2764706/1‑69 (MQ=255)
                 tCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACa                                                  >  1:2408513/1‑69 (MQ=255)
                 tCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACa                                                  >  1:1011511/1‑69 (MQ=255)
                 tCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACa                                                  >  1:822962/1‑69 (MQ=255)
                      aTCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatgg                                             <  1:2094421/69‑1 (MQ=255)
                       tCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatgga                                            >  1:2827747/1‑69 (MQ=255)
                       tCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatgga                                            >  1:93074/1‑69 (MQ=255)
                       tCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatgga                                            >  1:772427/1‑69 (MQ=255)
                       tCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatgga                                            >  1:2694028/1‑69 (MQ=255)
                        cTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatggat                                           <  1:434973/69‑1 (MQ=255)
                         tGACAATGTCGATTAACGGTGGTATCAAGTCGCTTGAAGAGGCCAAAGCACACCTGCAACATatggatg                                          >  1:240936/1‑69 (MQ=255)
                         tGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCaa                                                    <  1:119444/59‑1 (MQ=255)
                         tGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCaa                                                    <  1:2207426/59‑1 (MQ=255)
                         tGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACa                                                  >  1:1873651/1‑61 (MQ=255)
                         tGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatggatg                                          >  1:834554/1‑69 (MQ=255)
                         tGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatggatg                                          >  1:2105796/1‑69 (MQ=255)
                             gaTGTCGATTAACTGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:1077528/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:2927620/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:597985/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:783011/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:1072194/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:283690/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:691163/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:388549/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:1237930/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:897278/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:1347934/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:840794/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:1948630/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:2140405/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:2163416/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:2189121/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:2667214/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:2330269/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:2621933/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:2640373/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:116106/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:2509920/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:2504583/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:2492829/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:1188536/67‑1 (MQ=255)
                             gaTGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCAATTGCAACATATGGATGGCg                                       <  1:1886256/67‑1 (MQ=255)
                             gaTGCCGATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:3112910/67‑1 (MQ=255)
                                 tCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGca                                                             <  1:3067513/42‑1 (MQ=255)
                                 tCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGca                                                             <  1:849207/42‑1 (MQ=255)
                                   gATTAACGGTGGTATCAAGTCGCTGGAAGAAGCTAAAGCGCATTTGCAACATATGGATGGCg                                       <  1:38952/62‑1 (MQ=255)
                                       aaCGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTg                                     >  1:2083622/1‑60 (MQ=255)
                                       aaCGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCggg                            >  1:1148865/1‑69 (MQ=255)
                                       aaCGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCggg                            >  1:576831/1‑69 (MQ=255)
                                       aaCGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCggg                            >  1:671774/1‑69 (MQ=255)
                                       aaCGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCgg                             >  1:2535652/1‑68 (MQ=255)
                                       aaCGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCgg                             >  1:2847622/1‑68 (MQ=255)
                                                aTCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatg                                              <  1:788016/42‑1 (MQ=255)
                                                aTCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatg                                              <  1:2533365/42‑1 (MQ=255)
                                                aTCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatg                                              <  1:1672006/42‑1 (MQ=255)
                                                   aaGTCGCTGGAAGAGGCCAAAGCACACCTGGAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGt                  >  1:787526/1‑67 (MQ=255)
                                                   aaGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGt                  >  1:338220/1‑67 (MQ=255)
                                                   aaGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGt                  >  1:2177086/1‑67 (MQ=255)
                                                   aaGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGt                  >  1:1199609/1‑67 (MQ=255)
                                                   aaGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGt                  >  1:967010/1‑67 (MQ=255)
                                                   aaGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTa                 <  1:3195491/68‑1 (MQ=255)
                                                   aaGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTa                 <  1:1912438/68‑1 (MQ=255)
                                                                 ggCCAAAGCACACCTGCAACATATCGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTAt  >  1:1488714/1‑69 (MQ=255)
                                                                 |                                                                    
CAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTAT  >  minE/2801179‑2801312

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: