Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,830,540 G→A 49.0% V130V (GTC→GTT ygbI ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,830,5400GA49.0% 2.5 / 95.5 106V130V (GTC→GTTygbIpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (42/12);  new base A (20/32);  total (62/44)
Fisher's exact test for biased strand distribution p-value = 6.49e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 115 total aligned reads displayed.

CCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCA  >  minE/1830472‑1830607
                                                                    |                                                                   
ccGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTtgtg                                                                     <  1:2744143/69‑1 (MQ=255)
ccGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTtgtg                                                                     <  1:2619182/69‑1 (MQ=255)
ccGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTtgtg                                                                     >  1:1669213/1‑69 (MQ=255)
        gtgAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACcaca                                                               >  1:1552156/1‑67 (MQ=255)
        gtgAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACcaca                                                               >  1:151194/1‑67 (MQ=255)
        gtgAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACcaca                                                               >  1:407810/1‑67 (MQ=255)
        gtgAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACcaca                                                               >  1:543079/1‑67 (MQ=255)
        gtgAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACcaca                                                               >  1:1505693/1‑67 (MQ=255)
        gtgAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACcaca                                                               >  1:3297628/1‑67 (MQ=255)
         tgAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtc                                                            >  1:1328707/1‑69 (MQ=255)
         tgAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtc                                                            >  1:1162703/1‑69 (MQ=255)
             taatTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtcagt                                                         <  1:2854415/68‑1 (MQ=255)
             taatTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtcagt                                                         <  1:1934094/68‑1 (MQ=255)
             taatTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtcagt                                                         <  1:3342900/68‑1 (MQ=255)
             taatTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtcagt                                                         <  1:2339182/68‑1 (MQ=255)
                 tgtgCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTc                                                    <  1:1568662/69‑1 (MQ=255)
                  gtgCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTcc                                                   <  1:3127902/69‑1 (MQ=255)
                  gtgCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTcc                                                   <  1:1470244/69‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:3243689/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2383126/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:1800320/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2070439/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2636586/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2131038/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2571043/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2269570/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2423443/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:1093662/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:940037/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:407232/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:562092/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:583813/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:962571/68‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCAtgtg                                              <  1:237981/69‑1 (MQ=255)
                       aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCAtgtg                                              <  1:922773/69‑1 (MQ=255)
                        ttACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCAtgtg                                              <  1:2839428/68‑1 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGAc                                                                   >  1:1765845/1‑45 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGAc                                                                   >  1:1164999/1‑45 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGAc                                                                   >  1:2180945/1‑45 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGAc                                                                   >  1:256515/1‑45 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGAc                                                                   >  1:915473/1‑45 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGAc                                                                   >  1:1455929/1‑45 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGAc                                                                   >  1:1726869/1‑45 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:985416/1‑67 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:1044390/1‑67 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:2633376/1‑67 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:858379/1‑67 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:62768/1‑67 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:2864448/1‑67 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:331916/1‑67 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:1510162/1‑67 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:2270210/1‑67 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:2208481/1‑67 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:1921332/1‑67 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:1940684/1‑67 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:2057701/1‑67 (MQ=255)
                          aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:129023/1‑67 (MQ=255)
                                                   aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:1420424/69‑1 (MQ=255)
                                                   aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:1508819/69‑1 (MQ=255)
                                                   aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:390025/69‑1 (MQ=255)
                                                   aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:378320/69‑1 (MQ=255)
                                                   aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:828222/69‑1 (MQ=255)
                                                   aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:3090915/69‑1 (MQ=255)
                                                   aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:2750711/69‑1 (MQ=255)
                                                   aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:977110/69‑1 (MQ=255)
                                                   aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:2094056/69‑1 (MQ=255)
                                                   aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:1147116/69‑1 (MQ=255)
                                                   aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:2206125/69‑1 (MQ=255)
                                                    tAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:2077556/68‑1 (MQ=255)
                                                     aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTgg                 >  1:1865076/1‑68 (MQ=255)
                                                     aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTgg                 >  1:2488163/1‑68 (MQ=255)
                                                     aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                >  1:929666/1‑69 (MQ=255)
                                                     aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                >  1:1048544/1‑69 (MQ=255)
                                                     aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                >  1:666752/1‑69 (MQ=255)
                                                     aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                >  1:878650/1‑69 (MQ=255)
                                                     aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                >  1:883826/1‑69 (MQ=255)
                                                     aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                >  1:1924814/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCga       >  1:961298/1‑68 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:1364505/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:787005/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:1324251/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:237317/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:1461198/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:2400759/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:1576200/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:380608/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:2426018/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:3430954/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:1736090/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:1825498/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:3188365/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:205491/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:2750659/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:2748468/1‑69 (MQ=255)
                                                              gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:2174044/1‑69 (MQ=255)
                                                                 tgtgACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCa      >  1:1027947/1‑67 (MQ=255)
                                                                 tgtgACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCa      >  1:664110/1‑67 (MQ=255)
                                                                 tgtgACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCa      >  1:1419728/1‑67 (MQ=255)
                                                                    gACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCa  <  1:2460809/68‑1 (MQ=255)
                                                                    gACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCa  <  1:821871/68‑1 (MQ=255)
                                                                    |                                                                   
CCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCA  >  minE/1830472‑1830607

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