Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,830,549 C→G 57.8% L127L (CTG→CTC ygbI ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,830,5490CG57.8% 39.0 / 93.4 90L127L (CTG→CTCygbIpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (28/10);  new base G (20/32);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 1.28e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCA  >  minE/1830481‑1830607
                                                                    |                                                          
tgAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtc                                                            >  1:1328707/1‑69 (MQ=255)
tgAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtc                                                            >  1:1162703/1‑69 (MQ=255)
    taatTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtcagt                                                         <  1:2339182/68‑1 (MQ=255)
    taatTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtcagt                                                         <  1:3342900/68‑1 (MQ=255)
    taatTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtcagt                                                         <  1:2854415/68‑1 (MQ=255)
    taatTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtcagt                                                         <  1:1934094/68‑1 (MQ=255)
        tgtgCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTc                                                    <  1:1568662/69‑1 (MQ=255)
         gtgCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTcc                                                   <  1:1470244/69‑1 (MQ=255)
         gtgCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTcc                                                   <  1:3127902/69‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2070439/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2481601/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2571043/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2423443/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:1981067/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2383126/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2636586/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2131038/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:1800320/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:2269570/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:800953/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:3243689/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:940037/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:962571/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                               <  1:1093662/68‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCAtgtg                                              <  1:997508/69‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCAtgtg                                              <  1:237981/69‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCAtgtg                                              <  1:922773/69‑1 (MQ=255)
              aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCAtgtg                                              <  1:2148164/69‑1 (MQ=255)
               ttACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCAtgtg                                              <  1:2839428/68‑1 (MQ=255)
                 aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:858379/1‑67 (MQ=255)
                 aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:1044390/1‑67 (MQ=255)
                 aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:2270210/1‑67 (MQ=255)
                 aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:985416/1‑67 (MQ=255)
                 aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:2633376/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:129023/1‑67 (MQ=255)
                 aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:1510162/1‑67 (MQ=255)
                 aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:1921332/1‑67 (MQ=255)
                 aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:2208481/1‑67 (MQ=255)
                 aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:1940684/1‑67 (MQ=255)
                 aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:2057701/1‑67 (MQ=255)
                 aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                             >  1:1250176/1‑67 (MQ=255)
                                          aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:390025/69‑1 (MQ=255)
                                          aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:2206125/69‑1 (MQ=255)
                                          aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:828222/69‑1 (MQ=255)
                                          aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:3090915/69‑1 (MQ=255)
                                          aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:1508819/69‑1 (MQ=255)
                                          aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:1457322/69‑1 (MQ=255)
                                          aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:2750711/69‑1 (MQ=255)
                                          aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:977110/69‑1 (MQ=255)
                                          aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:2695364/69‑1 (MQ=255)
                                          aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:378320/69‑1 (MQ=255)
                                          aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:1147116/69‑1 (MQ=255)
                                          aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:2094056/69‑1 (MQ=255)
                                          aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                  <  1:1420424/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:666752/1‑69 (MQ=255)
                                            aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                >  1:49881/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)
                                            aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                >  1:929666/1‑69 (MQ=255)
                                            aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                >  1:1048544/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCga       >  1:961298/1‑68 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:787005/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:1364505/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:1324251/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:3423471/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:380608/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:2366124/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:237317/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:2174044/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:2400759/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:2426018/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:205491/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:2748468/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:2750659/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:1825498/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:1736090/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:3188365/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:1576200/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:3430954/1‑69 (MQ=255)
                                                     gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc       >  1:1461198/1‑69 (MQ=255)
                                                        tgtgACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCa      >  1:922828/1‑67 (MQ=255)
                                                        tgtgACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCa      >  1:664110/1‑67 (MQ=255)
                                                        tgtgACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCa      >  1:1419728/1‑67 (MQ=255)
                                                        tgtgACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCa      >  1:1027947/1‑67 (MQ=255)
                                                           gACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCa  <  1:821871/68‑1 (MQ=255)
                                                           gACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCa  <  1:2460809/68‑1 (MQ=255)
                                                                    |                                                          
TGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCA  >  minE/1830481‑1830607

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