Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,269,150 G→A 27.7% T52T (ACC→ACT 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,1500GA27.7% 107.6 / 38.6 97T52T (ACC→ACTbirAbifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase and DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding
Reads supporting (aligned to +/- strand):  ref base G (39/31);  new base A (20/7);  total (59/38)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGTTTATTAATAGCCGCCCGGCTCAT  >  minE/2269088‑2269215
                                                              |                                                                 
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:1036552/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:528529/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:1845793/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:1739309/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:428389/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:2420829/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:1597465/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:2517044/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:2656005/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:736102/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:1452602/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:2959785/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:1178682/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTaa                                                               >  1:299926/1‑67 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                                                                >  1:2996610/1‑66 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                                                                >  1:3072797/1‑66 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                                                                >  1:1204265/1‑66 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                                                                >  1:807507/1‑66 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                                                                >  1:846025/1‑66 (MQ=255)
aaTATCTTTTCGGCATTAAGTAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTa                                                                >  1:2010183/1‑66 (MQ=255)
    tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGTCTGTATCCTTTACCCGGAACGGTAAAGa                                                            >  1:1156613/1‑66 (MQ=255)
    tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            >  1:304790/1‑66 (MQ=255)
    tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            >  1:3385530/1‑66 (MQ=255)
    tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            >  1:120788/1‑66 (MQ=255)
    tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            >  1:711975/1‑66 (MQ=255)
    tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            >  1:1738614/1‑66 (MQ=255)
    tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            >  1:1126352/1‑66 (MQ=255)
    tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            >  1:1803552/1‑66 (MQ=255)
    tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            >  1:3171324/1‑66 (MQ=255)
    tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            >  1:2300763/1‑66 (MQ=255)
    tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACa                                                          >  1:499037/1‑68 (MQ=255)
    tCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACa                                                          >  1:295026/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:1528698/64‑1 (MQ=255)
      tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            <  1:2642416/64‑1 (MQ=255)
      tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            <  1:2641626/64‑1 (MQ=255)
      tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            <  1:668437/64‑1 (MQ=255)
      tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            <  1:1356086/64‑1 (MQ=255)
      tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            <  1:3259458/64‑1 (MQ=255)
      tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            <  1:1536519/64‑1 (MQ=255)
      tGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGa                                                            <  1:416781/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:934203/1‑69 (MQ=255)
          tAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGc                                                   >  1:2064127/1‑69 (MQ=255)
            gCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGcc                                                  <  1:3328072/68‑1 (MQ=255)
            gCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGcc                                                  <  1:1971956/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:1202521/1‑69 (MQ=255)
                 aaGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTc                                            >  1:1504353/1‑69 (MQ=255)
                 aaGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTc                                            >  1:1727879/1‑69 (MQ=255)
                 aaGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTc                                            >  1:2290585/1‑69 (MQ=255)
                    tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc                                          <  1:2702973/68‑1 (MQ=255)
                    tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc                                          <  1:2955934/68‑1 (MQ=255)
                    tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc                                          <  1:897581/68‑1 (MQ=255)
                    tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc                                          <  1:2450390/68‑1 (MQ=255)
                    tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc                                          <  1:905108/68‑1 (MQ=255)
                    tAACTGGATTGGCTCGGGCAAGCTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc                                          <  1:67164/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:335700/1‑49 (MQ=255)
                         ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                                        >  1:889826/1‑49 (MQ=255)
                         ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                                        >  1:464187/1‑49 (MQ=255)
                         ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                                        >  1:2566918/1‑49 (MQ=255)
                         ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                                        >  1:3043532/1‑49 (MQ=255)
                         ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                                        >  1:1998358/1‑49 (MQ=255)
                         ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATc                                                        >  1:2153685/1‑49 (MQ=255)
                         ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg                                    <  1:2727903/69‑1 (MQ=255)
                         ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg                                    <  1:1363424/69‑1 (MQ=255)
                         ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg                                    <  1:2024964/69‑1 (MQ=255)
                         ggATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg                                    <  1:3013579/69‑1 (MQ=255)
                          gATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAgtg                                    <  1:899145/68‑1 (MQ=255)
                                cTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGaa                              <  1:3400513/68‑1 (MQ=255)
                                cTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGaa                              <  1:3252024/68‑1 (MQ=255)
                                cTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGaa                              <  1:1125677/68‑1 (MQ=255)
                                cTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGaa                              <  1:2669189/68‑1 (MQ=255)
                                     gcaggcTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGtt                        <  1:2334599/69‑1 (MQ=255)
                                     gcaggcTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGtt                        <  1:1424113/69‑1 (MQ=255)
                                      caggcTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAAtgtg                          <  1:3404009/66‑1 (MQ=255)
                                          cTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGTttatt                    >  1:2801988/1‑68 (MQ=255)
                                          cTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGTttatt                    >  1:1293038/1‑68 (MQ=255)
                                          cTGTATCCTTTACCCGGAACAGTAAATACATCAACGCCCCAGTCAc                                          <  1:1584464/46‑1 (MQ=39)
                                               tCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGTTTATTAATAgc              >  1:245220/1‑69 (MQ=255)
                                               tCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGTTTATTAATAgc              >  1:1862554/1‑69 (MQ=255)
                                               tCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGTTTATTAATAgc              >  1:1374323/1‑69 (MQ=255)
                                                           aaCGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGTTTATTAATAGCCGCCCGGCTCAt  <  1:977683/69‑1 (MQ=255)
                                                           aaCGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGTTTATTAATAGCCGCCCGGCTCAt  <  1:1254040/69‑1 (MQ=255)
                                                           aaCGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGTTTATTAATAGCCGCCCGGCTCAt  <  1:1457450/69‑1 (MQ=255)
                                                           aaCGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGTTTATTAATAGCCGCCCGGCTCAt  <  1:1737957/69‑1 (MQ=255)
                                                              |                                                                 
AATATCTGTTTAGCATTAAGTAACTGGATAGGCTCAGGCAGGCTGTATCCTTTACCCGGAACGGTAAAGACATCAACGCCCCAGTCACGCAGTGTCTGAATGTGTTTATTAATAGCCGCCCGGCTCAT  >  minE/2269088‑2269215

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: