Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 214,101 T→C 23.5% F154F (TTT→TTC yafS → predicted S‑adenosyl‑L‑methionine‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE214,1010TC23.5% 67.1 / 15.7 55F154F (TTT→TTCyafSpredicted S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (25/17);  new base C (3/10);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 2.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

TTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTACACGCCAGCCGTTTCCA  >  minE/214036‑214169
                                                                 |                                                                    
tttATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTAct                                                                   >  1:1154519/1‑69 (MQ=255)
tttATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTAct                                                                   >  1:1191634/1‑69 (MQ=255)
tttATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTAct                                                                   >  1:101516/1‑69 (MQ=255)
        aTTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGc                                                            <  1:3335259/68‑1 (MQ=255)
        aTTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGc                                                            <  1:768947/68‑1 (MQ=255)
        aTTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGc                                                            <  1:1391254/68‑1 (MQ=255)
        aTTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGc                                                            <  1:1403562/68‑1 (MQ=255)
        aTTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGc                                                            <  1:2668743/68‑1 (MQ=255)
          tACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGc                                                         <  1:2505278/69‑1 (MQ=255)
          tACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGc                                                         <  1:2544681/69‑1 (MQ=255)
           aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                                                        <  1:2116587/69‑1 (MQ=255)
           aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                                                        <  1:923972/69‑1 (MQ=255)
           aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                                                        <  1:2511165/69‑1 (MQ=255)
           aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                                                        <  1:2718424/69‑1 (MQ=255)
           aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                                                        <  1:2813443/69‑1 (MQ=255)
           aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                                                        <  1:2817348/69‑1 (MQ=255)
           aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                                                        <  1:3302819/69‑1 (MQ=255)
           aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                                                        <  1:843974/69‑1 (MQ=255)
           aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                                                        <  1:1743026/69‑1 (MQ=255)
           aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                                                        <  1:1613839/69‑1 (MQ=255)
               aaaCTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGc                                                            <  1:2954137/61‑1 (MQ=255)
                aaCTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTgg                                                   <  1:1531834/69‑1 (MQ=255)
                            tATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCtt                                       <  1:1096206/69‑1 (MQ=255)
                            tATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCCCTCtt                                       <  1:867913/69‑1 (MQ=255)
                              ttGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTttg                                     >  1:330919/1‑69 (MQ=255)
                                gcgcAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTg                                                    <  1:3252395/52‑1 (MQ=255)
                                gcgcAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTg                                                    <  1:3307413/52‑1 (MQ=255)
                                gcgcAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTg                                                    <  1:3066190/52‑1 (MQ=255)
                                gcgcAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTg                                                    <  1:71862/52‑1 (MQ=255)
                                gcgcAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTg                                                    <  1:1219454/52‑1 (MQ=255)
                                gcgcAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTg                                                    <  1:2979406/52‑1 (MQ=255)
                                       aCCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAAttt                             >  1:1952292/1‑68 (MQ=255)
                                       aCCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAAttt                             >  1:302534/1‑68 (MQ=255)
                                       aCCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAAttt                             >  1:498166/1‑68 (MQ=255)
                                       aCCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAAttt                             >  1:715271/1‑68 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGTTGGACTGGCTCTCTTTGTTGAATTTTg                           >  1:2149330/1‑57 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGctct                                          >  1:2084740/1‑42 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAAtttt                            >  1:1145772/1‑56 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTg                           >  1:902079/1‑57 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTg                           >  1:811345/1‑57 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTg                           >  1:2759913/1‑57 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTg                           >  1:216441/1‑57 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTacac                >  1:3306883/1‑68 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTacac                >  1:1066324/1‑68 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTacac                >  1:1920484/1‑68 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTacac                >  1:760627/1‑68 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTacac                >  1:705834/1‑68 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTacac                >  1:3336984/1‑68 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTacac                >  1:3035433/1‑68 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTacac                >  1:3191289/1‑68 (MQ=255)
                                                    aCAGCCGGATGTTTACTCTGATGAGGCAGCTGGACTg                                               >  1:3346807/1‑37 (MQ=255)
                                                    aCAGCCGGATGTTTACACTGATGCGGCAGCTGGACTGGctc                                           >  1:2375340/1‑41 (MQ=255)
                                                          ggATGTTCACTCTGATGCGGCAGCTGGACTGGctct                                          >  1:1849519/1‑36 (MQ=255)
                                                          ggATGTTCACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTACACGCCAGcc         >  1:231123/1‑69 (MQ=255)
                                                          ggATGTTCACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTACACGCCAGcc         >  1:2225953/1‑69 (MQ=255)
                                                               tttACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTACACGCCAGCCGTTTc    >  1:1873772/1‑69 (MQ=255)
                                                               tttACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTACACGCCAGCCGTTTc    >  1:3213127/1‑69 (MQ=255)
                                                                 tACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTACACGCCAGCCGTTTCCa  >  1:3018192/1‑69 (MQ=255)
                                                                 tACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTACACGCCAGCCGTTTCCa  >  1:1225284/1‑69 (MQ=255)
                                                                 |                                                                    
TTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTTTGAAGTGCTACACGCCAGCCGTTTCCA  >  minE/214036‑214169

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: