Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,513,641 G→A 29.0% L80L (CTG→TTG)  ubiX ← 3‑octaprenyl‑4‑hydroxybenzoate carboxy‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,513,6410GA29.0% 113.9 / 49.7 110L80L (CTG→TTG) ubiX3‑octaprenyl‑4‑hydroxybenzoate carboxy‑lyase
Reads supporting (aligned to +/- strand):  ref base G (26/52);  new base A (15/17);  total (41/69)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 110 total aligned reads displayed.

CCATCAGTATAGCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATCGTGCGTGACATCGGCT  >  minE/1513574‑1513705
                                                                   |                                                                
ccATCAGTATAGCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGc                                                                 <  1:1421418/69‑1 (MQ=255)
ccATCAGTATAGCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGc                                                                 <  1:2991838/69‑1 (MQ=255)
  aTCAGTATAGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGt                                                               <  1:1627560/69‑1 (MQ=255)
  aTCAGTATAGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGt                                                               <  1:203301/69‑1 (MQ=255)
  aTCAGTATAGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGt                                                               <  1:391275/69‑1 (MQ=255)
  aTCAGTATAGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGt                                                               <  1:373813/69‑1 (MQ=255)
  aTCAGTATAGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGt                                                               <  1:359905/69‑1 (MQ=255)
  aTCAGTATAGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGt                                                               <  1:2059725/69‑1 (MQ=255)
  aTCAGTATAGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGt                                                               <  1:2272120/69‑1 (MQ=255)
  aTCAGTATAGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGt                                                               <  1:3133391/69‑1 (MQ=255)
   tCAGTATAGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGt                                                               <  1:2131682/68‑1 (MQ=255)
   tCAGTATAGCTATGGACAATGCCGGACAGGGTTATGATTGAACAGGGTAAAATCACCATCCCCAACGt                                                               <  1:2920938/68‑1 (MQ=255)
     aGTATAGCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCt                                                             >  1:2026013/1‑68 (MQ=255)
     aGTATAGCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCt                                                             >  1:3279795/1‑68 (MQ=255)
     aGTATAGCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCt                                                             >  1:2280/1‑68 (MQ=255)
     aGTATAGCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCt                                                             >  1:172669/1‑68 (MQ=255)
         taGCTGTGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGaa                                                         >  1:1649911/1‑68 (MQ=255)
         taGCTATGGACAATGCCGTACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGaa                                                         >  1:1090709/1‑68 (MQ=255)
         taGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGaa                                                         >  1:231500/1‑68 (MQ=255)
         taGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGaa                                                         >  1:2249111/1‑68 (MQ=255)
         taGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGaa                                                         >  1:1985818/1‑68 (MQ=255)
         taGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGaa                                                         >  1:2634726/1‑68 (MQ=255)
         taGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGaa                                                         >  1:1446463/1‑68 (MQ=255)
         taGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGaa                                                         >  1:2051487/1‑68 (MQ=255)
         taGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGaa                                                         >  1:3254714/1‑68 (MQ=255)
         taGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGaa                                                         >  1:221358/1‑68 (MQ=255)
         taGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGaa                                                         >  1:1059742/1‑68 (MQ=255)
         taGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGa                                                          >  1:1614697/1‑67 (MQ=255)
         taGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGa                                                          >  1:1502128/1‑67 (MQ=255)
         taGCTATGGACAATGCCGGACAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAACGTCTGGa                                                          >  1:1288362/1‑67 (MQ=255)
           gatATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGa                                                      >  1:2304641/3‑69 (MQ=255)
           gCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAg                                                       <  1:2191377/68‑1 (MQ=255)
           gCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAg                                                       <  1:1950340/68‑1 (MQ=255)
           gCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAg                                                       <  1:2633444/68‑1 (MQ=255)
           gCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAg                                                       <  1:1726828/68‑1 (MQ=255)
           gCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAg                                                       <  1:3076012/68‑1 (MQ=255)
           gCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAg                                                       <  1:1012754/68‑1 (MQ=255)
           gCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGa                                                      >  1:2378911/1‑69 (MQ=255)
           gCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGa                                                      >  1:1658271/1‑69 (MQ=255)
           gCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGa                                                      >  1:1441631/1‑69 (MQ=255)
           gCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGa                                                      >  1:3425001/1‑69 (MQ=255)
           gCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACAATCCCCAGCGTCTGGAAAg                                                       <  1:1040317/68‑1 (MQ=255)
                 gACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAg                                                       >  1:2088917/1‑62 (MQ=255)
                 gACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAAcc                                                   >  1:2164679/1‑66 (MQ=255)
                 gACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAAcc                                                   >  1:542232/1‑66 (MQ=255)
                 gACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAAcc                                                   >  1:1785766/1‑66 (MQ=255)
                 gACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAAcc                                                   >  1:2595631/1‑66 (MQ=255)
                 gACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAAcc                                                   >  1:1105780/1‑66 (MQ=255)
                 gACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAAcc                                                   >  1:1355533/1‑66 (MQ=255)
                 gACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAAcc                                                   >  1:3045857/1‑66 (MQ=255)
                    aaTGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAaga                                             <  1:3299660/69‑1 (MQ=255)
                    aaTGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAaga                                             <  1:2577630/69‑1 (MQ=255)
                    aaTGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAaga                                             <  1:1349378/69‑1 (MQ=255)
                    aaTGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAaga                                             <  1:2775353/69‑1 (MQ=255)
                    aaTGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAaga                                             <  1:404864/69‑1 (MQ=255)
                       gCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCTTCTGGAAAGAACCGGAAGAGAt                                          <  1:1129244/69‑1 (MQ=255)
                       gCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGAt                                          <  1:3038980/69‑1 (MQ=255)
                       gCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGAt                                          <  1:3137907/69‑1 (MQ=255)
                       gCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGAt                                          <  1:1337440/69‑1 (MQ=255)
                       gCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGAt                                          <  1:1993897/69‑1 (MQ=255)
                       gCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGAt                                          <  1:1506435/69‑1 (MQ=255)
                       gCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGAt                                          <  1:2488692/69‑1 (MQ=255)
                          ggAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCt                                       <  1:682057/69‑1 (MQ=255)
                          ggAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCt                                       <  1:1430321/69‑1 (MQ=255)
                          ggAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCt                                       <  1:423617/69‑1 (MQ=255)
                                ggTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTggcggc                                 <  1:1353928/69‑1 (MQ=255)
                                  ttttGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGAt                                          <  1:3334732/58‑1 (MQ=255)
                                    ttgattgaACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAat                              >  1:695429/1‑68 (MQ=255)
                                    ttgattgaACAGGGTAAAATCACCATCCCCAACGt                                                               <  1:1856922/35‑1 (MQ=255)
                                         tgaACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATcgcg                        >  1:2498917/1‑69 (MQ=255)
                                         tgaACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATcgcg                        >  1:819904/1‑69 (MQ=255)
                                             cAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTggcggc                                 <  1:342500/56‑1 (MQ=255)
                                             cAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTggcggc                                 <  1:283275/56‑1 (MQ=255)
                                             cAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTggcggc                                 <  1:893563/56‑1 (MQ=255)
                                                ggTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATcg                 <  1:746105/69‑1 (MQ=255)
                                                ggTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATcg                 <  1:907101/69‑1 (MQ=255)
                                                ggTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATcg                 <  1:186717/69‑1 (MQ=255)
                                                ggTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATcg                 <  1:2936477/69‑1 (MQ=255)
                                                ggTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATcg                 <  1:2671915/69‑1 (MQ=255)
                                                 gTAAAATCACCATCCCCAGCGTCTTGAAAGAACCGGAAGAGAt                                          <  1:546498/43‑1 (MQ=255)
                                                 gTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAaga                                             <  1:3278017/40‑1 (MQ=255)
                                                 gTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAaga                                             <  1:1991408/40‑1 (MQ=255)
                                                 gTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAaga                                             <  1:475397/40‑1 (MQ=255)
                                                 gTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAaga                                             <  1:2744915/40‑1 (MQ=255)
                                                 gTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAaga                                             <  1:2842626/40‑1 (MQ=255)
                                                 gTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAaga                                             <  1:364786/40‑1 (MQ=255)
                                                 gTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAaga                                             <  1:3248210/40‑1 (MQ=255)
                                                               cccAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATCGTGCGTGACATCGGCt  <  1:1096944/69‑1 (MQ=255)
                                                               cccAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATCGTGCGTGACATCGGCt  >  1:2879787/1‑69 (MQ=255)
                                                               cccAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATCGTGCGTGACATCGGCt  <  1:1176502/69‑1 (MQ=255)
                                                               cccAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATCGTGCGTGACATCGGCt  <  1:637240/69‑1 (MQ=255)
                                                               cccAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATCGTGCGTGACATCGGCt  >  1:694729/1‑69 (MQ=255)
                                                               cccAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATCGTGCGTGACATCGGCt  <  1:183540/69‑1 (MQ=255)
                                                               cccAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATCGTGCGTGACATCGGCt  <  1:1211195/69‑1 (MQ=255)
                                                               cccAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATCGTGCGTGACATCGGCt  <  1:2048627/69‑1 (MQ=255)
                                                               cccAACGTCTGGAAAGAACCGGAAGAAATGCTGGCAGCAATATCGCGCGCATCGTGTGTCACATCAGCt  <  1:288011/69‑1 (MQ=255)
                                                               cccAACGTCTGGAAAGAACCGGAAGAAATGCTGGCAGCAATATCGCGCGCATCGTGTGTCACATCAGCt  <  1:265750/69‑1 (MQ=255)
                                                               cccAACGTCTGGAAAGAACCGGAAGAAATGCTGGCAGCAATATCGCGCGCATCGTGTGTCACATCAGCt  <  1:2905959/69‑1 (MQ=255)
                                                               cccAACGTCTGGAAAGAACCGGAAGAAATGCTGGCAGCAATATCGCGCGCATCGTGTGTCACATCAGCt  <  1:2644384/69‑1 (MQ=255)
                                                               cccAACGTCTGGAAAGAACCGGAAGAAATGCTGGCAGCAATATCGCGCGCATCGTGTGTCACATCAGCt  <  1:2153445/69‑1 (MQ=255)
                                                                   |                                                                
CCATCAGTATAGCTATGGACAATGCCGGAAAGGGTTTTGATTGAACAGGGTAAAATCACCATCCCCAGCGTCTGGAAAGAACCGGAAGAGATGCTGGCGGCAATATCGCGCGCATCGTGCGTGACATCGGCT  >  minE/1513574‑1513705

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: