Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,225,152 A→T 54.6% G338G (GGA→GGT yebU → predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,225,1520AT54.6% 5.6 / 45.8 44G338G (GGA→GGTyebUpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (11/9);  new base T (11/13);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGC  >  minE/1225095‑1225204
                                                         |                                                    
cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTTAAAGATCGCGAAGCCGGTCAAATTCGTCa                                           <  1:1433232/69‑1 (MQ=255)
cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                                           <  1:2915736/69‑1 (MQ=255)
cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                                           <  1:2892243/69‑1 (MQ=255)
cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                                           <  1:2877694/69‑1 (MQ=255)
cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                                           <  1:1210863/69‑1 (MQ=255)
cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                                           <  1:1326113/69‑1 (MQ=255)
cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                                           <  1:1405204/69‑1 (MQ=255)
cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                                           <  1:280742/69‑1 (MQ=255)
cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                                           <  1:2718581/69‑1 (MQ=255)
cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                                           <  1:2574778/69‑1 (MQ=255)
cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                                           <  1:367174/69‑1 (MQ=255)
cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                                           <  1:2290201/69‑1 (MQ=255)
cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                                           <  1:374922/69‑1 (MQ=255)
           aaGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACa                                 <  1:28168/68‑1 (MQ=255)
           aaGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACa                                 <  1:2521756/68‑1 (MQ=255)
             gTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGgt                              <  1:660482/69‑1 (MQ=255)
              tCGGTAATTTTCCGTTCAGCCCTGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGgtg                             <  1:1165110/69‑1 (MQ=255)
              tCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGgtg                             <  1:1915852/69‑1 (MQ=255)
              tCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGgtg                             <  1:2592483/69‑1 (MQ=255)
                   aaTTTTCCGTTCAGCCCGGTGAAAGAACGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTgg                         >  1:2437888/1‑68 (MQ=255)
                   aaTTTTCCGTTCAGCCCGGTGAAAGAACGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTgg                         >  1:791233/1‑68 (MQ=255)
                   aaTTTTCCGTTCAGCCCGGTGAAAGAACGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTgg                         >  1:3209771/1‑68 (MQ=255)
                   aaTTTTCCGTTCAGCCCGGTGAAAGAACGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTgg                         >  1:1071348/1‑68 (MQ=255)
                   aaTATTCCGTTCAGCCAGGTGAAAGAACGCGAAGCTGGACAAATTCGTCAGGCGGCTAAAGGTGTTgg                         >  1:1638171/1‑68 (MQ=255)
                     tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg                          >  1:985518/1‑65 (MQ=255)
                     tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg                          >  1:744287/1‑65 (MQ=255)
                     tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg                          >  1:3419487/1‑65 (MQ=255)
                     tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg                          >  1:840729/1‑65 (MQ=255)
                     tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg                          >  1:878697/1‑65 (MQ=255)
                     tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg                          >  1:913334/1‑65 (MQ=255)
                     tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg                          >  1:102977/1‑65 (MQ=255)
                     tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg                          >  1:264233/1‑65 (MQ=255)
                     tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg                          >  1:2546941/1‑65 (MQ=255)
                     tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg                          >  1:2309988/1‑65 (MQ=255)
                     tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg                          >  1:2067294/1‑65 (MQ=255)
                                   cGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGaaa        <  1:1032347/69‑1 (MQ=255)
                                       gAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAAcc     >  1:3328435/1‑68 (MQ=255)
                                       gAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCt    >  1:1650519/1‑69 (MQ=255)
                                       gAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCt    >  1:970882/1‑69 (MQ=255)
                                       gAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTAAAGGTGTTGGCTTAAACTGGGATGAAAACCt    >  1:203490/1‑69 (MQ=255)
                                        aaaGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAgg                               >  1:2043127/1‑41 (MQ=255)
                                        aaaGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAgg                               >  1:1700362/1‑41 (MQ=255)
                                          aGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTgc  <  1:1411220/68‑1 (MQ=255)
                                          aGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTgc  <  1:1105498/68‑1 (MQ=255)
                                                         |                                                    
CCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGC  >  minE/1225095‑1225204

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