Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,269,128 G→A 24.6% L60L (CTG→TTG)  birA ← bifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase and DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,269,1280GA24.6% 139.9 / 37.3 109L60L (CTG→TTG) birAbifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase and DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding
Reads supporting (aligned to +/- strand):  ref base G (44/38);  new base A (21/6);  total (65/44)
Fisher's exact test for biased strand distribution p-value = 4.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

Only 100 of 112 total aligned reads displayed.

ACACTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGTT  >  minE/2269064‑2269193
                                                                |                                                                 
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                               <  1:106467/69‑1 (MQ=255)
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                               <  1:3200387/69‑1 (MQ=255)
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                               <  1:764333/69‑1 (MQ=255)
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                               <  1:50261/69‑1 (MQ=255)
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                               <  1:3097860/69‑1 (MQ=255)
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                               <  1:1241138/69‑1 (MQ=255)
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                               <  1:1266194/69‑1 (MQ=255)
   cTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTAt                                                            >  1:207912/1‑69 (MQ=255)
   cTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTAt                                                            >  1:1911981/1‑69 (MQ=255)
   cTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTAt                                                            >  1:3245272/1‑69 (MQ=255)
     aCCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATcc                                                          >  1:2880623/1‑69 (MQ=255)
     aCCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATcc                                                          >  1:549349/1‑69 (MQ=255)
     aCCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATcc                                                          >  1:1281618/1‑69 (MQ=255)
     aCCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATcc                                                          >  1:2616405/1‑69 (MQ=255)
     aCCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATcc                                                          >  1:2750272/1‑69 (MQ=255)
             ccAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                               <  1:1321519/56‑1 (MQ=255)
              cAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCgg                                                 <  1:1223618/69‑1 (MQ=255)
              cAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCgg                                                 <  1:1171572/69‑1 (MQ=255)
              cAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCgg                                                 <  1:2076700/69‑1 (MQ=255)
              cAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCgg                                                 <  1:2129430/69‑1 (MQ=255)
              cAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCgg                                                 <  1:533760/69‑1 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTa                                                      >  1:1463163/1‑54 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                         >  1:1845793/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                         >  1:428389/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                         >  1:2420829/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                         >  1:2517044/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                         >  1:2656005/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                         >  1:299926/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                         >  1:2959785/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                         >  1:736102/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                         >  1:1036552/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                         >  1:1178682/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                         >  1:1452602/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                         >  1:528529/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                         >  1:1739309/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                                          >  1:2010183/1‑66 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                                          >  1:846025/1‑66 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                                          >  1:2996610/1‑66 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                                          >  1:3072797/1‑66 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                                          >  1:807507/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGTCTGTATCCTTTACCCGGAACGGTAAAGa                                      >  1:1156613/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTTTATCCTTTACCAGGAAc                                              >  1:845339/1‑58 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      >  1:2300763/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      >  1:711975/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      >  1:3385530/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      >  1:3171324/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      >  1:120788/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      >  1:1738614/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      >  1:304790/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      >  1:1803552/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      >  1:1126352/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACa                                    >  1:2760100/1‑68 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACa                                    >  1:295026/1‑68 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACa                                    >  1:3352118/1‑68 (MQ=255)
                             cTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                  >  1:188621/1‑69 (MQ=255)
                             cTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                  >  1:3210307/1‑69 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      <  1:416781/64‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      <  1:2642416/64‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      <  1:668437/64‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      <  1:1356086/64‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      <  1:1536519/64‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      <  1:1528698/64‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                  <  1:1911044/68‑1 (MQ=255)
                                  tAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                      >  1:2008346/1‑60 (MQ=255)
                                  tAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGc                             >  1:2064127/1‑69 (MQ=255)
                                    gCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGcc                            <  1:1971956/68‑1 (MQ=255)
                                    gCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGcc                            <  1:3328072/68‑1 (MQ=255)
                                    gCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGcc                            <  1:1402068/68‑1 (MQ=255)
                                    gCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATACTTTACCCGGAACGGTAAAGACATCAACGcc                            <  1:3203999/68‑1 (MQ=255)
                                       ttAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTa                                          >  1:2733955/1‑51 (MQ=255)
                                       ttAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCa                         >  1:1029377/1‑68 (MQ=255)
                                         aaGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTc                      >  1:1504353/1‑69 (MQ=255)
                                         aaGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTc                      >  1:1727879/1‑69 (MQ=255)
                                         aaGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTc                      >  1:1202521/1‑69 (MQ=255)
                                            tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc                    <  1:67164/68‑1 (MQ=255)
                                            tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc                    <  1:897581/68‑1 (MQ=255)
                                            tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc                    <  1:905108/68‑1 (MQ=255)
                                            tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc                    <  1:2955934/68‑1 (MQ=255)
                                            tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc                    <  1:2702973/68‑1 (MQ=255)
                                            tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc                    <  1:2450390/68‑1 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTTTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg              <  1:2620037/69‑1 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACAt                                   >  1:1318542/1‑48 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                  >  1:335700/1‑49 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                  >  1:889826/1‑49 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                  >  1:464187/1‑49 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                  >  1:2153685/1‑49 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                  >  1:2566918/1‑49 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                  >  1:3043532/1‑49 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg              <  1:2727903/69‑1 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg              <  1:3013579/69‑1 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg              <  1:1363424/69‑1 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg              <  1:2024964/69‑1 (MQ=255)
                                                  gATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg              <  1:899145/68‑1 (MQ=255)
                                                        cTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGaa        <  1:3252024/68‑1 (MQ=255)
                                                        cTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGaa        <  1:2669189/68‑1 (MQ=255)
                                                        cTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGaa        <  1:1125677/68‑1 (MQ=255)
                                                        cTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGaa        <  1:3400513/68‑1 (MQ=255)
                                                             gcaggcTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGtt  <  1:1424113/69‑1 (MQ=255)
                                                             gcaggcTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGtt  <  1:2334599/69‑1 (MQ=255)
                                                              caggcTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAAtgtg    <  1:3404009/66‑1 (MQ=255)
                                                                |                                                                 
ACACTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGTT  >  minE/2269064‑2269193

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