Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 227,251 A→G 23.1% K15K (AAA→AAG crl → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE227,2510AG23.1% 41.7 / 10.7 26K15K (AAA→AAGcrlDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (11/9);  new base G (3/3);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

TGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT  >  minE/227185‑227319
                                                                  |                                                                    
tGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTa                                                                  >  1:2098637/1‑71 (MQ=255)
   aTCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCg                                                               >  1:902363/1‑71 (MQ=255)
   aTCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCg                                                               <  1:806920/71‑1 (MQ=255)
     cacaACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCa                                                             >  1:3232931/1‑71 (MQ=255)
      acaacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCa                                                             <  1:1061945/70‑1 (MQ=255)
       caacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCAc                                                            >  1:534566/1‑70 (MQ=255)
               gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCACTAGGCCCg                                                    >  1:1869500/1‑70 (MQ=255)
               gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCa                                                             <  1:1749346/61‑1 (MQ=255)
                aTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGta                                                  <  1:2390814/71‑1 (MQ=255)
                       tGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGt                                           >  1:1508111/1‑71 (MQ=255)
                        gACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGt                                           >  1:1147237/1‑70 (MQ=255)
                         aCGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTg                                          <  1:485637/70‑1 (MQ=255)
                            ttACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCACTAGGCCCGTATATTCGTGaa                                        <  1:1462802/69‑1 (MQ=255)
                               ccGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTaa                                   <  1:1617735/71‑1 (MQ=255)
                                cGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGACCGTATATTCGTGAAGGTAAg                                  <  1:1492310/71‑1 (MQ=255)
                                      gACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGaaa                             <  1:3453113/70‑1 (MQ=255)
                                        cacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCa                              >  1:1613836/1‑67 (MQ=255)
                                             gAAGAGCAGATTGATCAAAAAGTTTACCGCACTAGGCCCGTATATTCGTg                                          >  1:2989335/1‑50 (MQ=255)
                                             gAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATc                     >  1:1558451/1‑71 (MQ=255)
                                              aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATc                     >  1:2136965/1‑70 (MQ=255)
                                               agagCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGa                   >  1:1975703/1‑71 (MQ=255)
                                                   cAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCt               >  1:863557/1‑71 (MQ=255)
                                                    aGATTGATCAAAAAGTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCtt              <  1:388054/71‑1 (MQ=255)
                                                     gATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCttt             <  1:3254169/71‑1 (MQ=255)
                                                     gATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCttt             <  1:1445585/71‑1 (MQ=255)
                                                     gATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCttt             >  1:3751720/1‑71 (MQ=255)
                                                             aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTg     <  1:1551411/71‑1 (MQ=255)
                                                              aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTg     >  1:1786757/1‑70 (MQ=255)
                                                               aaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCaa                             <  1:2821710/45‑1 (MQ=255)
                                                               aaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTc   >  1:3615053/1‑71 (MQ=255)
                                                               aaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTc   >  1:536249/1‑71 (MQ=255)
                                                               aaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTc   >  1:1801330/1‑71 (MQ=255)
                                                                aaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCt  >  1:3658840/1‑71 (MQ=255)
                                                                  aTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGAt                         <  1:1602857/46‑1 (MQ=255)
                                                                  |                                                                    
TGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT  >  minE/227185‑227319

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: