Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 68,518 C→A 65.6% intergenic (‑70/‑590) leuL ← / → leuO leu operon leader peptide/DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE68,5180CA65.6% 42.8 / 47.9 61intergenic (‑70/‑590)leuL/leuOleu operon leader peptide/DNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (17/4);  new base A (22/18);  total (39/22)
Fisher's exact test for biased strand distribution p-value = 5.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTAAAATTAGCTAATTGTGCT  >  minE/68462‑68576
                                                        |                                                          
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:1839251/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:58748/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:2820995/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:2385522/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:2102740/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:295422/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:2519351/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:2785647/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:1809682/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:2757324/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:176523/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:912879/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:915301/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:1485148/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:1428900/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:3196556/69‑1 (MQ=255)
gATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:3340726/69‑1 (MQ=255)
 aTGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACg                                                <  1:2144026/68‑1 (MQ=255)
      tGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCACCCCTGACGCAAAAACAACGTcc                                             >  1:910699/1‑66 (MQ=255)
      tGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTccc                                            >  1:1917701/1‑67 (MQ=255)
          atatGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACg                                                >  1:2796199/1‑59 (MQ=255)
            atGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAg                                       >  1:74277/1‑66 (MQ=255)
            atGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGc                                      >  1:2621426/1‑67 (MQ=255)
            atGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGc                                      >  1:2986334/1‑67 (MQ=255)
            atGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGc                                      >  1:3376985/1‑67 (MQ=255)
            atGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGc                                      >  1:900105/1‑67 (MQ=255)
             tGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATaaa                                                    <  1:1079055/52‑1 (MQ=255)
             tGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATaaa                                                    >  1:542971/1‑52 (MQ=255)
             tGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATaaa                                                    >  1:717052/1‑52 (MQ=255)
             tGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATaaa                                                    >  1:1154076/1‑52 (MQ=255)
             tGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTg                                   <  1:486815/69‑1 (MQ=255)
                       gTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCAt                         >  1:801652/1‑69 (MQ=255)
                       gTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCAt                         >  1:303947/1‑69 (MQ=255)
                       gTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCAt                         >  1:3342424/1‑69 (MQ=255)
                        ttACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATc                        >  1:2655027/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCat                                       >  1:3105531/1‑47 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:497358/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:63770/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:36818/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:800151/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:868393/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:982330/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:3132536/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:1598495/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:1721695/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:1852527/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:1884093/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:2449580/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:2461862/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:2476925/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:2557701/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:2642834/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:2842997/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:3132557/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaa                   >  1:1080790/1‑68 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAATAAAAACGTCACGCCATTGCGAGTTCGGCATCCGTaaa                  >  1:1744236/1‑69 (MQ=255)
                              gATACAAAAACGGATGTCAACCCTGAAGCAAGAAAAACGTCACGCCATTGTGAGTTCGGCATCCGTaaa                  >  1:1777870/1‑69 (MQ=255)
                                     aaaCGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTAAAATTAGCt           <  1:2607044/69‑1 (MQ=255)
                                     aaaCGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTAAAATTAGCt           <  1:187118/69‑1 (MQ=255)
                                              tCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTAAAATTAGCTAATTGtgct  >  1:71136/1‑69 (MQ=255)
                                              tCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTAAAATTAGCTAATTGtgct  >  1:1875003/1‑69 (MQ=255)
                                                        |                                                          
GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTAAAATTAGCTAATTGTGCT  >  minE/68462‑68576

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: