Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,269,273 A→G 45.8% I11I (ATT→ATC 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,2730AG45.8% 9.3 / 69.4 83I11I (ATT→ATCbirAbifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase and DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding
Reads supporting (aligned to +/- strand):  ref base A (18/27);  new base G (25/13);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 2.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAC  >  minE/2269216‑2269328
                                                         |                                                       
tCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAg                                              <  1:653288/69‑1 (MQ=255)
 cccAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGt                                             >  1:728670/1‑69 (MQ=255)
 cccAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGt                                             >  1:1777662/1‑69 (MQ=255)
 cccAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGt                                             >  1:952544/1‑69 (MQ=255)
 cccAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGt                                             >  1:2772770/1‑69 (MQ=255)
 cccAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGt                                             >  1:884603/1‑69 (MQ=255)
    aGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTgg                                           <  1:3392865/68‑1 (MQ=255)
    aGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTgg                                           <  1:42059/68‑1 (MQ=255)
    aGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTgg                                           <  1:2050086/68‑1 (MQ=255)
    aGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTgg                                           <  1:2287948/68‑1 (MQ=255)
    aGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTgg                                           <  1:2818487/68‑1 (MQ=255)
          tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                                      >  1:759568/1‑51 (MQ=255)
          tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                                      >  1:184015/1‑51 (MQ=255)
          tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                                      >  1:167284/1‑51 (MQ=255)
          tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                                      >  1:846617/1‑51 (MQ=255)
          tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                                      >  1:1575645/1‑51 (MQ=255)
          tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                                      >  1:1409837/1‑51 (MQ=255)
          tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                                      >  1:892279/1‑51 (MQ=255)
          tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                                      >  1:2879322/1‑51 (MQ=255)
            aCCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTgg                                           <  1:1087752/60‑1 (MQ=255)
                       cGCCAGAGTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCAt                         <  1:2271826/67‑1 (MQ=255)
                       cGCCAGAGTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCAt                         <  1:968333/67‑1 (MQ=255)
                       cGCCAGAGTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCAt                         <  1:327714/67‑1 (MQ=255)
                            gagTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTg                  <  1:1416126/69‑1 (MQ=255)
                            gagTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTg                  <  1:3281668/69‑1 (MQ=255)
                            gagTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTg                  <  1:2587799/69‑1 (MQ=255)
                            gagTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTg                  <  1:3212610/69‑1 (MQ=255)
                              gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                <  1:1016201/69‑1 (MQ=255)
                              gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                <  1:2327272/69‑1 (MQ=255)
                              gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                <  1:2331205/69‑1 (MQ=255)
                              gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                <  1:3083581/69‑1 (MQ=255)
                              gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                <  1:2577900/69‑1 (MQ=255)
                              gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                <  1:500443/69‑1 (MQ=255)
                              gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                <  1:2734258/69‑1 (MQ=255)
                              gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                <  1:896209/69‑1 (MQ=255)
                              gTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtg                >  1:3031970/1‑68 (MQ=255)
                              gTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                >  1:189070/1‑69 (MQ=255)
                              gTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                >  1:2788469/1‑69 (MQ=255)
                              gTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                >  1:421142/1‑69 (MQ=255)
                              gTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                <  1:2019777/69‑1 (MQ=255)
                               tGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATACTTCATGAAATTGtc                <  1:2097950/68‑1 (MQ=255)
                               tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                >  1:1986830/1‑68 (MQ=255)
                               tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                >  1:2133218/1‑68 (MQ=255)
                               tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                >  1:1239524/1‑68 (MQ=255)
                               tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                >  1:2280857/1‑68 (MQ=255)
                               tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                >  1:2860592/1‑68 (MQ=255)
                               tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                >  1:1070823/1‑68 (MQ=255)
                               tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc                >  1:257727/1‑68 (MQ=255)
                                 aaaTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTa                                 >  1:796557/1‑49 (MQ=255)
                                 aaaTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTa                                 >  1:1029879/1‑49 (MQ=255)
                                   aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:1104584/69‑1 (MQ=255)
                                   aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:1691266/69‑1 (MQ=255)
                                   aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:168927/69‑1 (MQ=255)
                                   aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:1152416/69‑1 (MQ=255)
                                   aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:191345/69‑1 (MQ=255)
                                   aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:3112750/69‑1 (MQ=255)
                                   aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:2206439/69‑1 (MQ=255)
                                   aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:2163263/69‑1 (MQ=255)
                                   aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:3108909/69‑1 (MQ=255)
                                    ttCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACg                                       >  1:777402/1‑40 (MQ=37)
                                    ttCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACg                                       >  1:2302347/1‑40 (MQ=37)
                                    ttCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGaagga                   >  1:2295394/1‑57 (MQ=255)
                                    ttCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGaagga                   >  1:3045782/1‑57 (MQ=255)
                                    ttCGCCATTCGCTAACAGCGCGATCAATTTCAGAGGCACGGTGTTATCCt                             >  1:1611131/1‑50 (MQ=25)
                                    ttCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:164792/68‑1 (MQ=255)
                                    ttCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:464960/68‑1 (MQ=255)
                                    ttCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:234338/68‑1 (MQ=255)
                                    ttCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:1773706/68‑1 (MQ=255)
                                          attCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCAt                         <  1:1562253/47‑1 (MQ=255)
                                             cGCTAAGAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:1393386/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:658930/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:982324/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:2039744/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:1018829/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:1184441/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:1202278/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:1343022/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:1457324/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:1663744/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:1691655/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:1748224/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:3431547/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:2083402/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:2287918/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:2456980/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:2691464/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:2885760/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:3251865/1‑68 (MQ=255)
                                             cGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  >  1:3368191/1‑68 (MQ=255)
                                             cGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAc  <  1:1088016/68‑1 (MQ=255)
                                                      ggCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc           <  1:3433220/50‑1 (MQ=255)
                                                         |                                                       
TCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCACTGCGCTCAC  >  minE/2269216‑2269328

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: