Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 69,716 A→G 33.9% K203K (AAA→AAG leuO → DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE69,7160AG33.9% 35.5 / 26.9 53K203K (AAA→AAGleuODNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (20/15);  new base G (13/5);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAAAGATGAAATGGTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTT  >  minE/69652‑69783
                                                                |                                                                   
ttAAAGATGAAATGGTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACat                                                                  <  1:2439819/68‑1 (MQ=255)
ttAAAGATGAAATGGTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACat                                                                  <  1:2364906/68‑1 (MQ=255)
ttAAAGATGAAATGGTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACat                                                                  <  1:2136834/68‑1 (MQ=255)
ttAAAGATGAAATGGTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACat                                                                  <  1:639689/68‑1 (MQ=255)
 tAAAGATGAAATGGTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACat                                                                  <  1:2498100/67‑1 (MQ=255)
 tAAAGATGAAATGGTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACat                                                                  <  1:2761354/67‑1 (MQ=255)
   aaGATGAAATGGTGTTGGTAGCCAGCAAAAATCATCCAACAATTAAAGGCCCGTTACTGAAGCatgatg                                                              >  1:3009620/1‑69 (MQ=255)
   aaGATGAAATGGTGTTGGTAGCCAGCAAAAATCATCCAACAATTAAAGGCCCGTTACTGAAGCatgatg                                                              >  1:249559/1‑69 (MQ=255)
   aaGATGAAATGGTGTTGGTAGCCAGCAAAAATCATCCAACAATTAAAGGCCCGTTACTGAAGCatgatg                                                              >  1:3161825/1‑69 (MQ=255)
   aaGATGAAATGGTGTTGGTAGCCAGCAAAAATCATCCAACAATTAAAGGCCCGTTACTGAAGCatgatg                                                              >  1:2685336/1‑69 (MQ=255)
   aaGATGAAATGGTGTTGGTAGCCAGCAAAAATCATCCAACAATTAAAGGCCCGTTACTGAAGCatgatg                                                              >  1:2198002/1‑69 (MQ=255)
   aaGATGAAATGGTGTTGGTAGCCAGCAAAAATCATCCAACAATTAAAGGCCCGTTACTGAAGCatgatg                                                              >  1:2167865/1‑69 (MQ=255)
   aaGATGAAATGGTGTTGGTAGCCAGCAAAAATCATCCAACAATTAAAGGCCCGTTACTGAAGCatgatg                                                              >  1:294905/1‑69 (MQ=255)
   aaGATGAAATGGTGTTGGTAGCCAGCAAAAATCATCCAACAATTAAAGGCCCGTTACTGAAGCatgatg                                                              >  1:598044/1‑69 (MQ=255)
   aaGATGAAATGGTGTTGGTAGCCAGCAAAAATCATCCAACAATTAAAGGCCCGTTACTGAAGCatgatg                                                              >  1:1573352/1‑69 (MQ=255)
   aaGATGAAATGGTGTTGGTAGCCAGCAAAAATCATCCAACAATTAAAGGCCCGTTACTGAAGCatgatg                                                              >  1:1562572/1‑69 (MQ=255)
   aaGATGAAATGGTGTTGGTAGCCAGCAAAAATCATCCAACAATTAAAGGCCCGTTACTGAAGCatgatg                                                              >  1:644186/1‑69 (MQ=255)
   aaGATGAAATGGTGTTGGTAGCCAGCAAAAATCATCCAACAATTAAAGGCCCGTTACTGAAGCatgatg                                                              >  1:2821879/1‑69 (MQ=255)
   aaGATGAAATGGTGTTGGTAGCCAGCAAAAATCATCCAACAATTAAAGCCCGTTACTGAAGCatgatg                                                              >  1:167655/1‑68 (MQ=255)
    aGATGAAATGGTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACatgatg                                                              <  1:2195905/68‑1 (MQ=255)
     gATGAAATGGTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACatgatg                                                              <  1:3317844/67‑1 (MQ=255)
           aTGGTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATaa                                                       >  1:3010403/1‑68 (MQ=255)
              gTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGa                                                    >  1:1817104/1‑68 (MQ=255)
              gTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGa                                                    >  1:3176358/1‑68 (MQ=255)
              gTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGa                                                    >  1:1864339/1‑68 (MQ=255)
              gTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGa                                                    >  1:1065920/1‑68 (MQ=255)
                             aaaaaTCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGtt                                    >  1:2360186/1‑69 (MQ=255)
                             aaaaaTCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGtt                                    >  1:2581744/1‑69 (MQ=255)
                             aaaaaTCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGtt                                    >  1:1803788/1‑69 (MQ=255)
                                aaTCATCCAACAATTAAGGGCCCGTTACTGAAACatgatg                                                              <  1:2567748/40‑1 (MQ=255)
                                             ttAAAGGCCCGTTACTGAAGCATGATGTTTATAACGaac                                                  <  1:2645422/39‑1 (MQ=37)
                                             ttAAAGGCCCGTTACTGAAGCATGATGTTTATAACGaac                                                  <  1:2482094/39‑1 (MQ=37)
                                             ttAAAGGCCCGTTACTGAAGCATGATGTTTATAACGaac                                                  <  1:1172996/39‑1 (MQ=37)
                                             ttAAAGGCCCGTTACTGAAGCATGATGTTTATAACGaac                                                  <  1:75247/39‑1 (MQ=37)
                                             ttAAAGGCCCGTTACTGAAGCATGATGTTTATAACGaac                                                  <  1:2800523/39‑1 (MQ=37)
                                              tAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTcgc                   >  1:1143710/1‑69 (MQ=255)
                                                         tACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTc        <  1:2889279/69‑1 (MQ=255)
                                                         tACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTc        <  1:2462796/69‑1 (MQ=255)
                                                         tACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTc        <  1:378818/69‑1 (MQ=255)
                                                         tACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTc        <  1:733071/69‑1 (MQ=255)
                                                         tACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTc        <  1:1119801/69‑1 (MQ=255)
                                                          aCTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTc        <  1:2772782/68‑1 (MQ=255)
                                                             gAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTc                >  1:2327710/1‑57 (MQ=255)
                                                             gAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAAcc    >  1:2991041/1‑69 (MQ=255)
                                                             gAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAAcc    >  1:584202/1‑69 (MQ=255)
                                                             gAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAAcc    >  1:779945/1‑69 (MQ=255)
                                                             gAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAAcc    >  1:2187229/1‑69 (MQ=255)
                                                             gAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAAcc    >  1:2074097/1‑69 (MQ=255)
                                                             gAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAAcc    >  1:196680/1‑69 (MQ=255)
                                                             gAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAAcc    >  1:1357612/1‑69 (MQ=255)
                                                             gAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAAcc    >  1:1356436/1‑69 (MQ=255)
                                                             gAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAAcc    >  1:1133975/1‑69 (MQ=255)
                                                             gAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAAc     >  1:1789297/1‑68 (MQ=255)
                                                              aaaCATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAAcc    >  1:2196864/1‑68 (MQ=255)
                                                              aaaCATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCt   >  1:2546890/1‑69 (MQ=255)
                                                                aCATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCtt  >  1:274707/1‑68 (MQ=255)
                                                                aCATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCtt  >  1:2647031/1‑68 (MQ=255)
                                                                aCATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCtt  >  1:2596635/1‑68 (MQ=255)
                                                                aCATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCtt  >  1:2488477/1‑68 (MQ=255)
                                                                aCATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCtt  >  1:2135883/1‑68 (MQ=255)
                                                                aCATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCtt  >  1:1547314/1‑68 (MQ=255)
                                                                |                                                                   
TTAAAGATGAAATGGTGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTT  >  minE/69652‑69783

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: