Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,225,125 C→T 63.2% F329F (TTC→TTT yebU → predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,225,1250CT63.2% 18.3 / 20.8 38F329F (TTC→TTTyebUpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (7/7);  new base T (11/13);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGG  >  minE/1225081‑1225181
                                            |                                                        
cccGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGc                                   >  1:2643875/1‑68 (MQ=255)
cccGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGc                                   >  1:3302868/1‑68 (MQ=255)
         cccGCGCCCAAATACAAAGTCGGTAATTTTCCGTTCAGccc                                                     <  1:1566297/41‑1 (MQ=255)
              cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTTAAAGATCGCGAAGCCGGTCAAATTCGTCa                    <  1:1433232/69‑1 (MQ=255)
              cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                    <  1:367174/69‑1 (MQ=255)
              cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                    <  1:2892243/69‑1 (MQ=255)
              cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                    <  1:2877694/69‑1 (MQ=255)
              cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                    <  1:280742/69‑1 (MQ=255)
              cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                    <  1:2718581/69‑1 (MQ=255)
              cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                    <  1:2915736/69‑1 (MQ=255)
              cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                    <  1:2574778/69‑1 (MQ=255)
              cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                    <  1:2290201/69‑1 (MQ=255)
              cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                    <  1:374922/69‑1 (MQ=255)
              cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                    <  1:1405204/69‑1 (MQ=255)
              cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                    <  1:1326113/69‑1 (MQ=255)
              cccaAAATACAAAGTGGGCAATTTCCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCa                    <  1:1210863/69‑1 (MQ=255)
                         aaGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACa          <  1:2521756/68‑1 (MQ=255)
                         aaGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACa          <  1:28168/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:1071348/1‑68 (MQ=255)
                                 aaTTTTCCGTTCAGCCCGGTGAAAGAACGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTgg  >  1:3209771/1‑68 (MQ=255)
                                 aaTTTTCCGTTCAGCCCGGTGAAAGAACGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTgg  >  1:791233/1‑68 (MQ=255)
                                 aaTTTTCCGTTCAGCCCGGTGAAAGAACGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTgg  >  1:2437888/1‑68 (MQ=255)
                                 aaTATTCCGTTCAGCCAGGTGAAAGAACGCGAAGCTGGACAAATTCGTCAGGCGGCTAAAGGTGTTgg  >  1:1638171/1‑68 (MQ=255)
                                   tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg   >  1:264233/1‑65 (MQ=255)
                                   tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg   >  1:2546941/1‑65 (MQ=255)
                                   tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg   >  1:3419487/1‑65 (MQ=255)
                                   tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg   >  1:2309988/1‑65 (MQ=255)
                                   tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg   >  1:2067294/1‑65 (MQ=255)
                                   tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg   >  1:744287/1‑65 (MQ=255)
                                   tttcCCGTTTAGCCCGGTGAAAGATCGCGAAGCCGGTCAAATTCGTCAGGCGGCTGCAAGTGTTg   >  1:102977/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:985518/1‑65 (MQ=255)
                                            |                                                        
CCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGG  >  minE/1225081‑1225181

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: