Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,269,258 G→A 21.9% N16N (AAC→AAT 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,2580GA21.9% 107.8 / 20.6 77N16N (AAC→AATbirAbifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase and DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base A (5/12);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTATTAATAGCCGCCCGGCTCATTCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAC  >  minE/2269190‑2269319
                                                                    |                                                             
tgtTTATTAATAGCCGCCCGGCTCATTCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCg                                                               >  1:826258/1‑69 (MQ=255)
            gccgccCGGCTCATTCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCt                                                         >  1:2343390/1‑63 (MQ=255)
            gccgccCGGCTCATTCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCt                                                         >  1:2060018/1‑63 (MQ=255)
            gccgccCGGCTCATTCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCt                                                         >  1:413533/1‑63 (MQ=255)
            gccgccCGGCTCATTCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCt                                                         >  1:721461/1‑63 (MQ=255)
            gccgccCGGCTCATTCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTCCGCt                                                         >  1:2962367/1‑63 (MQ=255)
                     cTCATTCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCt                                                         <  1:2185855/54‑1 (MQ=255)
                     cTCATTCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCt                                                         <  1:1084405/54‑1 (MQ=255)
                     cTCATTCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCt                                                         <  1:251263/54‑1 (MQ=255)
                     cTCATTCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCt                                                         <  1:200498/54‑1 (MQ=255)
                          tCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAg                                     <  1:653288/69‑1 (MQ=255)
                           cccAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGt                                    >  1:952544/1‑69 (MQ=255)
                           cccAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGt                                    >  1:2772770/1‑69 (MQ=255)
                           cccAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGt                                    >  1:1777662/1‑69 (MQ=255)
                           cccAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGt                                    >  1:728670/1‑69 (MQ=255)
                           cccAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGt                                    >  1:884603/1‑69 (MQ=255)
                              aGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTgg                                  <  1:42059/68‑1 (MQ=255)
                              aGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTgg                                  <  1:2287948/68‑1 (MQ=255)
                              aGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTgg                                  <  1:2818487/68‑1 (MQ=255)
                              aGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTgg                                  <  1:2050086/68‑1 (MQ=255)
                              aGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTgg                                  <  1:3392865/68‑1 (MQ=255)
                                    tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                             >  1:759568/1‑51 (MQ=255)
                                    tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                             >  1:1409837/1‑51 (MQ=255)
                                    tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                             >  1:2879322/1‑51 (MQ=255)
                                    tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                             >  1:846617/1‑51 (MQ=255)
                                    tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                             >  1:184015/1‑51 (MQ=255)
                                    tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                             >  1:1575645/1‑51 (MQ=255)
                                    tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                             >  1:892279/1‑51 (MQ=255)
                                    tCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGcaatc                                             >  1:167284/1‑51 (MQ=255)
                                      aCCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTgg                                  <  1:1087752/60‑1 (MQ=255)
                                                 cGCCAGAGTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCAt                <  1:327714/67‑1 (MQ=255)
                                                 cGCCAGAGTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCAt                <  1:968333/67‑1 (MQ=255)
                                                 cGCCAGAGTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCAt                <  1:2271826/67‑1 (MQ=255)
                                                      gagTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTg         <  1:1416126/69‑1 (MQ=255)
                                                      gagTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTg         <  1:3281668/69‑1 (MQ=255)
                                                      gagTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTg         <  1:3212610/69‑1 (MQ=255)
                                                      gagTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTg         <  1:2587799/69‑1 (MQ=255)
                                                        gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       <  1:1016201/69‑1 (MQ=255)
                                                        gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       <  1:896209/69‑1 (MQ=255)
                                                        gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       <  1:500443/69‑1 (MQ=255)
                                                        gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       <  1:3083581/69‑1 (MQ=255)
                                                        gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       <  1:2734258/69‑1 (MQ=255)
                                                        gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       <  1:2577900/69‑1 (MQ=255)
                                                        gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       <  1:2327272/69‑1 (MQ=255)
                                                        gTGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       <  1:2331205/69‑1 (MQ=255)
                                                        gTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtg       >  1:3031970/1‑68 (MQ=255)
                                                        gTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       >  1:2788469/1‑69 (MQ=255)
                                                        gTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       <  1:2019777/69‑1 (MQ=255)
                                                        gTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       >  1:421142/1‑69 (MQ=255)
                                                        gTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       >  1:189070/1‑69 (MQ=255)
                                                         tGAAATTCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACGGTGTTATACTTCATGAAATTGtc       <  1:2097950/68‑1 (MQ=255)
                                                         tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       >  1:257727/1‑68 (MQ=255)
                                                         tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       >  1:2280857/1‑68 (MQ=255)
                                                         tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       >  1:1070823/1‑68 (MQ=255)
                                                         tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       >  1:1239524/1‑68 (MQ=255)
                                                         tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       >  1:2133218/1‑68 (MQ=255)
                                                         tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       >  1:2860592/1‑68 (MQ=255)
                                                         tGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGtc       >  1:1986830/1‑68 (MQ=255)
                                                           aaaTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTa                        >  1:1029879/1‑49 (MQ=255)
                                                           aaaTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTa                        >  1:796557/1‑49 (MQ=255)
                                                             aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc  <  1:168927/69‑1 (MQ=255)
                                                             aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc  <  1:1104584/69‑1 (MQ=255)
                                                             aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc  <  1:1152416/69‑1 (MQ=255)
                                                             aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc  <  1:3112750/69‑1 (MQ=255)
                                                             aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc  <  1:1691266/69‑1 (MQ=255)
                                                             aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc  <  1:2206439/69‑1 (MQ=255)
                                                             aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc  <  1:3108909/69‑1 (MQ=255)
                                                             aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc  <  1:191345/69‑1 (MQ=255)
                                                             aTTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc  <  1:2163263/69‑1 (MQ=255)
                                                              ttCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACg                              >  1:2302347/1‑40 (MQ=37)
                                                              ttCGCCATTCGCTAACAGGGCGATCAATTTCAGTGGCACg                              >  1:777402/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:1773706/68‑1 (MQ=255)
                                                              ttCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAc  <  1:234338/68‑1 (MQ=255)
                                                                    |                                                             
TGTTTATTAATAGCCGCCCGGCTCATTCCCAGCGTTTCACCCAACTGCTCGCCAGAGTGAAATTCACCGTTCGCTAACAGGGCAATCAATTTCAGTGGCACGGTGTTATCCTTCATGAAATTGTCTCCAC  >  minE/2269190‑2269319

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: