Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 214,077 C→A 27.0% T146T (ACC→ACA 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,0770CA27.0% 70.5 / 21.9 48T146T (ACC→ACAyafSpredicted S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (13/22);  new base A (3/10);  total (16/32)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TAGTGGCTTCAATCCCATCAGTTTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTT  >  minE/214014‑214142
                                                               |                                                                 
tAGTGGCTTCAATCCCATCAGTTTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGc                                                               <  1:3106568/68‑1 (MQ=255)
 aGTGGCTTCAATCCCATCAGTTTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGccc                                                             >  1:246945/1‑69 (MQ=255)
 aGTGGCTTCAATCCCATCAGTTTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGccc                                                             >  1:3325202/1‑69 (MQ=255)
     gCTTCAATCCCATCAGTTTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAATCCTCGCCCtat                                                          >  1:2568948/1‑68 (MQ=255)
     gCTTCAATCCCATCAGTTTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGcc                                                              <  1:397752/64‑1 (MQ=255)
     gCTTCAATCCCATCAGTTTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCtat                                                          >  1:2514039/1‑68 (MQ=255)
     gCTTCAATCCCATCAGTTTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCtat                                                          >  1:2055209/1‑68 (MQ=255)
         cAATCCCATCAGTTTAATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAg                                                     >  1:2107394/1‑69 (MQ=255)
         cAATCCCATCAGTTTAATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAg                                                     >  1:171163/1‑69 (MQ=255)
         cAATCCCATCAGTTTAATGGGATTACGCAAACTTGTGCCGGTATTGAGCAAAACATCGCGCTATAACa                                                      >  1:2583379/1‑68 (MQ=255)
            tCCCATCAGTTTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCg                                                  <  1:3172292/69‑1 (MQ=255)
            tCCCATCAGTTTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCg                                                  <  1:2418470/69‑1 (MQ=255)
            tCCCATCAGTTTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCg                                                  <  1:3228889/69‑1 (MQ=255)
                      tttATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTAct                                        >  1:101516/1‑69 (MQ=255)
                      tttATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTAct                                        >  1:1154519/1‑69 (MQ=255)
                      tttATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTAct                                        >  1:1191634/1‑69 (MQ=255)
                              aTTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGc                                 <  1:3335259/68‑1 (MQ=255)
                              aTTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGc                                 <  1:2668743/68‑1 (MQ=255)
                              aTTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGc                                 <  1:1391254/68‑1 (MQ=255)
                              aTTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGc                                 <  1:1403562/68‑1 (MQ=255)
                              aTTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGc                                 <  1:768947/68‑1 (MQ=255)
                                tACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGc                              <  1:2505278/69‑1 (MQ=255)
                                tACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGc                              <  1:2544681/69‑1 (MQ=255)
                                 aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                             <  1:843974/69‑1 (MQ=255)
                                 aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                             <  1:3302819/69‑1 (MQ=255)
                                 aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                             <  1:2718424/69‑1 (MQ=255)
                                 aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                             <  1:923972/69‑1 (MQ=255)
                                 aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                             <  1:2817348/69‑1 (MQ=255)
                                 aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                             <  1:2813443/69‑1 (MQ=255)
                                 aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                             <  1:2511165/69‑1 (MQ=255)
                                 aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                             <  1:1613839/69‑1 (MQ=255)
                                 aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                             <  1:2116587/69‑1 (MQ=255)
                                 aCGCAAACTTGTGCCGGTATTGCGCAAAACATCGCCCTATAACAGCCGGATGTTCACTCTGATGCGGCa                             <  1:1743026/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:3307413/52‑1 (MQ=255)
                                                      gcgcAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTg                         <  1:3252395/52‑1 (MQ=255)
                                                      gcgcAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTg                         <  1:71862/52‑1 (MQ=255)
                                                      gcgcAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTg                         <  1:3066190/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:498166/1‑68 (MQ=255)
                                                             aCCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAAttt  >  1:715271/1‑68 (MQ=255)
                                                             aCCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAAttt  >  1:302534/1‑68 (MQ=255)
                                                               |                                                                 
TAGTGGCTTCAATCCCATCAGTTTTATGGGATTACGCAAACTTGTGCCGGTATTGCGCAAAACCTCGCCCTATAACAGCCGGATGTTTACTCTGATGCGGCAGCTGGACTGGCTCTCTTTGTTGAATTT  >  minE/214014‑214142

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: