Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,798,470 C→A 35.5% V136V (GTG→GTT zur ← DNA‑binding transcriptional activator, Zn(II)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,798,4700CA35.5% 56.7 / 40.4 59V136V (GTG→GTTzurDNA‑binding transcriptional activator, Zn(II)‑binding
Reads supporting (aligned to +/- strand):  ref base C (19/19);  new base A (8/13);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCACGCCT  >  minE/2798406‑2798533
                                                                |                                                               
aCTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATcaca                                                              >  1:2306580/1‑68 (MQ=255)
   gTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTAt                                                          >  1:37843/1‑69 (MQ=255)
      cAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGc                                                        >  1:234023/1‑68 (MQ=255)
      cAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGc                                                        >  1:570870/1‑68 (MQ=255)
      cAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGc                                                        >  1:2570909/1‑68 (MQ=255)
      cAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGc                                                        >  1:2648012/1‑68 (MQ=255)
       aGGATGACGACACGCTTCTACTTCTACACATGCTGAACATAATCCATGTGCTTCAATAACATTATGc                                                        >  1:1616699/1‑67 (MQ=255)
       aGGATGACGACACGCTTCTACTTCTACACATGCTGAACATAATCCATGTGCTTCAATAACATTATGc                                                        >  1:2189260/1‑67 (MQ=255)
       aGGATGACGACACGCTTCTACTTCTACACATGCTGAACATAATCCATGTGCTTCAATAACATTATGc                                                        >  1:1032378/1‑67 (MQ=255)
       aGGATGACGACACGCTTCTACTTCTACACATGCTGAACATAATCCATGTGCTTCAATAACATTATGc                                                        >  1:1949700/1‑67 (MQ=255)
       aGGATGACGACACGCTTCTACTTCTACACATGCTGAACATAATCCATGTGCTTCAATAACATTATGc                                                        >  1:2788095/1‑67 (MQ=255)
       aGGATGACGACACGCTTCTACTTCTACACATGCTGAACATAATCCATGTGCTTCAATAACATTATGc                                                        >  1:2899810/1‑67 (MQ=255)
       aGGATGACGACACGCTTCTACTTCTACACATGCTGAACATAATCCATGTGCTTCAATAACATTATGc                                                        >  1:1595505/1‑67 (MQ=255)
       aGGATGACGACACGCTTCTACTTCTACACATGCTGAACATAATCCATGTGCTTCAATAACATTATGc                                                        >  1:533137/1‑67 (MQ=255)
               gacACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCa                                               <  1:1492052/68‑1 (MQ=255)
               gacACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCa                                               <  1:1054285/68‑1 (MQ=255)
                  acGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACAtt                                                            >  1:2846889/1‑52 (MQ=255)
                  acGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAAc                                            >  1:1912505/1‑68 (MQ=255)
                  acGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAAc                                            >  1:3397452/1‑68 (MQ=255)
                  acGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAAc                                            >  1:1430505/1‑68 (MQ=255)
                                    aTGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCg                         <  1:2437154/69‑1 (MQ=255)
                                    aTGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCg                         <  1:3128224/69‑1 (MQ=255)
                                    aTGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCg                         <  1:2427063/69‑1 (MQ=255)
                                    aTGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCg                         <  1:888447/69‑1 (MQ=255)
                                    aTGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCg                         <  1:107313/69‑1 (MQ=255)
                                    aTGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCg                         <  1:1142558/69‑1 (MQ=255)
                                    aTGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCg                         <  1:1487353/69‑1 (MQ=255)
                                    aTGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCg                         <  1:1818080/69‑1 (MQ=255)
                                     tGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGt                        <  1:964161/69‑1 (MQ=255)
                                     tGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGt                        <  1:221884/69‑1 (MQ=255)
                                     tGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGt                        <  1:735733/69‑1 (MQ=255)
                                     tGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGt                        <  1:1102876/69‑1 (MQ=255)
                                     tGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGt                        <  1:2593952/69‑1 (MQ=255)
                                     tGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGt                        <  1:2670344/69‑1 (MQ=255)
                                     tGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGt                        <  1:2903289/69‑1 (MQ=255)
                                                       tgCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTa                       <  1:3236023/52‑1 (MQ=255)
                                                       tgCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTATGCAt                  >  1:2316859/1‑57 (MQ=255)
                                                       tgCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTATGCAt                  >  1:1611181/1‑57 (MQ=255)
                                                       tgCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCAc      >  1:2737187/1‑69 (MQ=255)
                                                       tgCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCAc      >  1:2756029/1‑69 (MQ=255)
                                                       tgCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCAc      >  1:238258/1‑69 (MQ=255)
                                                       tgCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCAc      >  1:1955149/1‑69 (MQ=255)
                                                       tgCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCAc      >  1:2889580/1‑69 (MQ=255)
                                                       tgCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCAc      >  1:3411536/1‑69 (MQ=255)
                                                       tgCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCAc      >  1:1932033/1‑69 (MQ=255)
                                                       tgCTTCAATCACATTATGCCGCAGGGCAACCCCAttttt                                   >  1:2198804/1‑39 (MQ=39)
                                                           tCAATAACATTATGCCGCAGGGCAAAACCCATTTTTGCCTCCAGCGTATGCATGATGTCTTCCACGCCt  <  1:2747645/69‑1 (MQ=255)
                                                           tCAATAACATTATGCCGCAGGGCAAAACCCATTTTTGCCGCCAGCGTATGCATGATGTCTTCCACGCCt  <  1:2951887/69‑1 (MQ=255)
                                                           tCAATAACATTATGCCGCAGGGCAAAACCCATTTTTGCCGCCAGCGTATGCATGATGTCTTCCACGCCt  <  1:151979/69‑1 (MQ=255)
                                                           tCAATAACATTATGCCGCAGGGCAAAACCCATTTTTGCCGCCAGCGTATGCATGATGTCTTCCACGCCt  <  1:2698131/69‑1 (MQ=255)
                                                           tCAATAACATTATGCCGCAGGGCAAAACCCATTTTTGCCGCCAGCGTATGCATGATGTCTTCCACGCCt  <  1:2472415/69‑1 (MQ=255)
                                                           tCAATAACATTATGCCGCAGGGCAAAACCCATTTTTGCCGCCAGCGTATGCATGATGTCTTCCACGCCt  <  1:2718682/69‑1 (MQ=255)
                                                           tCAATAACATTATGCCGCAGGGCAAAACCCATTTTTGCCGCCAGCGTATGCATGATGTCTTCCACGCCt  <  1:1704740/69‑1 (MQ=255)
                                                           tCAATAACATTATGCCGCAGGGCAAAACCCATTTTTGCCGCCAGCGTATGCATGATGTCTTCCACGCCt  <  1:1133402/69‑1 (MQ=255)
                                                           tCAATAACATTATGCCGCAGGGCAAAACCCATTTTTGCCGCCAGCGTATGCATGATGTCTTCCACGCCt  <  1:2805915/69‑1 (MQ=255)
                                                           tCAATAACATTATGCCGCAGGGCAAAACCCATTTTTGCCGCCAGCGTATGCATGATGTCTTCCACGCCt  <  1:178094/69‑1 (MQ=255)
                                                           tCAATAACATTATGCCGCAGGGCAAAACCCATTTTTGCCGCCAGCGTATGCATGATGTCTTCCACGCCt  <  1:2352435/69‑1 (MQ=255)
                                                           tCAATAACATTATGCCGCAGGGCAAAACCCATTTTTGCCGCCAGCGTATGCATGATGTCTTCCACGCCt  <  1:922796/69‑1 (MQ=255)
                                                           tCAATAACATTATGCCGCAGGGCAAAACCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCACGCCt  <  1:812153/69‑1 (MQ=255)
                                                              aTCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTa                       <  1:2242640/45‑1 (MQ=255)
                                                                |                                                               
ACTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCACGCCT  >  minE/2798406‑2798533

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: