Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,830,554 A→G 50.0% S126P (TCA→CCA)  ygbI ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,830,5540AG50.0% ‑1.0 / 95.0 104S126P (TCA→CCA) ygbIpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (30/22);  new base G (20/32);  total (50/54)
Fisher's exact test for biased strand distribution p-value = 7.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

TGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACAACTT  >  minE/1830489‑1830621
                                                                 |                                                                   
tgtgCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTc                                                                  <  1:1568662/69‑1 (MQ=255)
 gtgCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTcc                                                                 <  1:1470244/69‑1 (MQ=255)
 gtgCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTcc                                                                 <  1:3127902/69‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:2269570/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:2423443/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:1800320/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:2571043/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:2383126/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:2636586/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:1981067/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:3243689/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:407232/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:583813/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:800953/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:2131038/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:940037/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:1093662/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:962571/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCAtgt                                                             <  1:2070439/68‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCAtgtg                                                            <  1:2148164/69‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCAtgtg                                                            <  1:237981/69‑1 (MQ=255)
      aTTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCAtgtg                                                            <  1:922773/69‑1 (MQ=255)
       ttACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCAtgtg                                                            <  1:2839428/68‑1 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:2633376/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:1044390/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:2270210/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:2864448/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:331916/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:2208481/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:62768/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:858379/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:985416/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:1940684/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:1510162/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:1921332/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:1250176/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:2057701/1‑67 (MQ=255)
         aCTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGa                                                           >  1:129023/1‑67 (MQ=255)
                                  aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:390025/69‑1 (MQ=255)
                                  aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:378320/69‑1 (MQ=255)
                                  aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:1420424/69‑1 (MQ=255)
                                  aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:2094056/69‑1 (MQ=255)
                                  aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:1457322/69‑1 (MQ=255)
                                  aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:3090915/69‑1 (MQ=255)
                                  aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:1147116/69‑1 (MQ=255)
                                  aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:1508819/69‑1 (MQ=255)
                                  aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:2750711/69‑1 (MQ=255)
                                  aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:2695364/69‑1 (MQ=255)
                                  aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:828222/69‑1 (MQ=255)
                                  aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:977110/69‑1 (MQ=255)
                                  aTAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:2206125/69‑1 (MQ=255)
                                   tAACGAAGTCGTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTg                                <  1:2077556/68‑1 (MQ=255)
                                    aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTgg                               >  1:2488163/1‑68 (MQ=255)
                                    aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTgg                               >  1:1865076/1‑68 (MQ=255)
                                    aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                              >  1:1048544/1‑69 (MQ=255)
                                    aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                              >  1:883826/1‑69 (MQ=255)
                                    aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                              >  1:929666/1‑69 (MQ=255)
                                    aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                              >  1:49881/1‑69 (MQ=255)
                                    aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                              >  1:878650/1‑69 (MQ=255)
                                    aaCGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGt                              >  1:666752/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:2174044/1‑69 (MQ=255)
                                             gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc                     >  1:380608/1‑69 (MQ=255)
                                             gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc                     >  1:3430954/1‑69 (MQ=255)
                                             gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc                     >  1:3423471/1‑69 (MQ=255)
                                             gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc                     >  1:205491/1‑69 (MQ=255)
                                             gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc                     >  1:3188365/1‑69 (MQ=255)
                                             gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc                     >  1:2366124/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:2750659/1‑69 (MQ=255)
                                             gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc                     >  1:2748468/1‑69 (MQ=255)
                                             gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc                     >  1:1825498/1‑69 (MQ=255)
                                             gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc                     >  1:2426018/1‑69 (MQ=255)
                                             gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc                     >  1:1576200/1‑69 (MQ=255)
                                             gTTAGTAACCACTGTGAGTGGCTCCATGTGAGTCAGATGCTGTGCAATCGCTAACGTGGTCGTTCCCGc                     >  1:1736090/1‑69 (MQ=255)
                                                tgtgACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCa                    >  1:922828/1‑67 (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:821871/68‑1 (MQ=255)
                                                   gACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCa                <  1:2460809/68‑1 (MQ=255)
                                                             agtGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACaa     >  1:2069238/1‑69 (MQ=255)
                                                             agtGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTACCGCATCCAGATAGATACaa     >  1:663421/1‑69 (MQ=255)
                                                                gACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATGCAACtt  <  1:734880/69‑1 (MQ=255)
                                                                gACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACAACtt  <  1:145479/69‑1 (MQ=255)
                                                                gACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACAACtt  <  1:2046756/69‑1 (MQ=255)
                                                                gACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACAACtt  <  1:2392363/69‑1 (MQ=255)
                                                                gACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACAACtt  <  1:2858850/69‑1 (MQ=255)
                                                                gACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACAACtt  <  1:930781/69‑1 (MQ=255)
                                                                gACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACAACtt  <  1:2862697/69‑1 (MQ=255)
                                                                gACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACAACtt  <  1:927512/69‑1 (MQ=255)
                                                                gACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACAACtt  <  1:3286890/69‑1 (MQ=255)
                                                                gACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACAACtt  <  1:572598/69‑1 (MQ=255)
                                                                gACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACAACtt  <  1:1406915/69‑1 (MQ=255)
                                                                gACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACAACtt  <  1:787424/69‑1 (MQ=255)
                                                                 |                                                                   
TGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAGTGACTCCATGTGAATCAGATGCTGTGCTATCGCTAACGTGGTCGTTCCCGCATCCAGATAGATACAACTT  >  minE/1830489‑1830621

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: