Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,269,117 A→T 25.7% P63P (CCT→CCA 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,1170AT25.7% 157.3 / 57.3 105P63P (CCT→CCAbirAbifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase and DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding
Reads supporting (aligned to +/- strand):  ref base A (44/34);  new base T (21/6);  total (65/40)
Fisher's exact test for biased strand distribution p-value = 6.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 105 total aligned reads displayed.

ACACTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTG  >  minE/2269064‑2269181
                                                     |                                                                
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                   <  1:764333/69‑1 (MQ=255)
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                   <  1:106467/69‑1 (MQ=255)
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                   <  1:3200387/69‑1 (MQ=255)
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                   <  1:1266194/69‑1 (MQ=255)
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                   <  1:1241138/69‑1 (MQ=255)
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                   <  1:3056011/69‑1 (MQ=255)
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                   <  1:50261/69‑1 (MQ=255)
acacTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                   <  1:3097860/69‑1 (MQ=255)
   cTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTAt                                                >  1:740311/1‑69 (MQ=255)
   cTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTAt                                                >  1:3245272/1‑69 (MQ=255)
   cTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTAt                                                >  1:1911981/1‑69 (MQ=255)
     aCCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATcc                                              >  1:549349/1‑69 (MQ=255)
     aCCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATcc                                              >  1:2750272/1‑69 (MQ=255)
     aCCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATcc                                              >  1:1281618/1‑69 (MQ=255)
     aCCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATcc                                              >  1:2616405/1‑69 (MQ=255)
     aCCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATcc                                              >  1:2880623/1‑69 (MQ=255)
             ccAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTg                                                   <  1:1321519/56‑1 (MQ=255)
              cAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCgg                                     <  1:533760/69‑1 (MQ=255)
              cAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCgg                                     <  1:2076700/69‑1 (MQ=255)
              cAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCgg                                     <  1:2129430/69‑1 (MQ=255)
              cAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCgg                                     <  1:1223618/69‑1 (MQ=255)
              cAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCgg                                     <  1:1171572/69‑1 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTa                                          >  1:1463163/1‑54 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                             >  1:428389/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                             >  1:1036552/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                             >  1:528529/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:1845793/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                             >  1:1178682/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                             >  1:736102/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                             >  1:1597465/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                             >  1:1739309/1‑67 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                              >  1:2010183/1‑66 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                              >  1:2996610/1‑66 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                              >  1:846025/1‑66 (MQ=255)
                        aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                              >  1:1204265/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:711975/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                          >  1:1803552/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                          >  1:2300763/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:1126352/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                          >  1:1738614/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                          >  1:304790/1‑66 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACa                        >  1:499037/1‑68 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACa                        >  1:2760100/1‑68 (MQ=255)
                            tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACa                        >  1:3352118/1‑68 (MQ=255)
                             cTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                      >  1:3210307/1‑69 (MQ=255)
                             cTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                      >  1:188621/1‑69 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                          <  1:3259458/64‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                          <  1:668437/64‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                          <  1:2641626/64‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                          <  1:1356086/64‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                          <  1:1528698/64‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                          <  1:1536519/64‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                          <  1:2642416/64‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                      <  1:1911044/68‑1 (MQ=255)
                              tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCa                     >  1:2929981/1‑69 (MQ=255)
                                  tAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                          >  1:2008346/1‑60 (MQ=255)
                                  tAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGc                 >  1:2064127/1‑69 (MQ=255)
                                    gCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGcc                <  1:3328072/68‑1 (MQ=255)
                                    gCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGcc                <  1:1402068/68‑1 (MQ=255)
                                    gCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGcc                <  1:1971956/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:1202521/1‑69 (MQ=255)
                                         aaGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTc          >  1:2290585/1‑69 (MQ=255)
                                         aaGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTc          >  1:1504353/1‑69 (MQ=255)
                                         aaGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTc          >  1:1727879/1‑69 (MQ=255)
                                            tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc        <  1:905108/68‑1 (MQ=255)
                                            tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc        <  1:2702973/68‑1 (MQ=255)
                                            tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc        <  1:897581/68‑1 (MQ=255)
                                            tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc        <  1:67164/68‑1 (MQ=255)
                                            tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc        <  1:2955934/68‑1 (MQ=255)
                                            tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc        <  1:2450390/68‑1 (MQ=255)
                                             aaCTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGc      <  1:1361341/69‑1 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTTTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg  <  1:2620037/69‑1 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACAt                       >  1:1318542/1‑48 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                      >  1:464187/1‑49 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                      >  1:3043532/1‑49 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                      >  1:2566918/1‑49 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                      >  1:889826/1‑49 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                      >  1:335700/1‑49 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                      >  1:1998358/1‑49 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                      >  1:2153685/1‑49 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg  <  1:1363424/69‑1 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg  <  1:2024964/69‑1 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg  <  1:3013579/69‑1 (MQ=255)
                                                 ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg  <  1:2727903/69‑1 (MQ=255)
                                                  gATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg  <  1:899145/68‑1 (MQ=255)
                                                     |                                                                
ACACTACCGCCATCCAGCTGACCCAATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTG  >  minE/2269064‑2269181

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: