Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,073,250 T→C 18.3% intergenic (+88/‑96) ebgR → / → ebgA DNA‑binding transcriptional repressor/cryptic beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,073,2500TC18.3% 97.7 / 12.6 65intergenic (+88/‑96)ebgR/ebgADNA‑binding transcriptional repressor/cryptic beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (32/21);  new base C (5/7);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCGTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCACCATG  >  minE/2073187‑2073312
                                                               |                                                              
tCGTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTAc                                                           >  1:2315892/1‑69 (MQ=255)
tCGTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTAc                                                           >  1:1999557/1‑69 (MQ=255)
tCGTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTAc                                                           >  1:797543/1‑69 (MQ=255)
tCGTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTAc                                                           >  1:3102719/1‑69 (MQ=255)
 cGTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTACt                                                          >  1:1144448/1‑69 (MQ=255)
 cGTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTACt                                                          >  1:2859674/1‑69 (MQ=255)
 cGTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTACt                                                          >  1:2093650/1‑69 (MQ=255)
 cGTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTACt                                                          >  1:571934/1‑69 (MQ=255)
 cGTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTACt                                                          >  1:1680212/1‑69 (MQ=255)
 cGTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTACt                                                          >  1:1326664/1‑69 (MQ=255)
  gTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACTCGCTACtt                                                         >  1:1720355/1‑69 (MQ=255)
  gTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACTCGCTACtt                                                         >  1:1770594/1‑69 (MQ=255)
  gTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACTCGCTACtt                                                         >  1:395825/1‑69 (MQ=255)
  gTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACTCGCTACtt                                                         >  1:2935008/1‑69 (MQ=255)
  gTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACTCGCTACtt                                                         >  1:2879464/1‑69 (MQ=255)
              cAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAAt                                             >  1:2183685/1‑69 (MQ=255)
                 tgctgcCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAAtttt                                          >  1:884086/1‑69 (MQ=255)
                   ctgcCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCt                                                                >  1:574271/1‑45 (MQ=255)
                     gcCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTa                                      <  1:2777609/69‑1 (MQ=255)
                           tATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTc                                 >  1:426732/1‑68 (MQ=255)
                           tATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTc                                 >  1:1297965/1‑68 (MQ=255)
                                        cAGTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTacg                    >  1:2729570/1‑67 (MQ=255)
                                        cAGTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTaca                    >  1:1152121/1‑68 (MQ=255)
                                        cAGTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTaca                    >  1:602436/1‑68 (MQ=255)
                                        cAGTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTaca                    >  1:2275213/1‑68 (MQ=255)
                                        cAGTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTaca                    >  1:1628693/1‑68 (MQ=255)
                                        cAGTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTaca                    >  1:1771031/1‑68 (MQ=255)
                                          gTTAAAGTAAATGCATTTACTCGCTACTTTTTAGTAAAATTTTTACTAAACTCCCCAGCAATTAcaca                  <  1:1337072/68‑1 (MQ=38)
                                          gTTAAAGTAAATGCATTTACTCGCTACTTTTTAGTAAAATTTTTACTAAACTCCCCAGCAATTAcaca                  <  1:773769/68‑1 (MQ=38)
                                          gTTAAAGTAAATGCATTTACTCGCTACTTTTTAGTAAAATTTTTACTAAACTCCCCAGCAATTAcaca                  <  1:760921/68‑1 (MQ=38)
                                          gTTAAAGTAAATGCATTTACTCGCTACTTTTTAGTAAAATTTTTACTAAACTCCCCAGCAATTAcaca                  <  1:2372645/68‑1 (MQ=38)
                                          gTTAAAGTAAATGCATTTACTCGCTACTTTTTAGTAAAATTTTTACTAAACTCCCCAGCAATTAcaca                  <  1:630672/68‑1 (MQ=38)
                                          gTTAAAGTAAATGCATTTACTCGCTACTTTTTAGTAAAATTTTTACTAAACTCCCCAGCAATTAcaca                  <  1:2537033/68‑1 (MQ=38)
                                          gTTAAAGTAAATGCATTTACTCGCTACTTTTTAGTAAAATTTTTACTAAACTCCCCAGCAATTAcaca                  <  1:3365436/68‑1 (MQ=38)
                                          gTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACaa                 <  1:2188648/69‑1 (MQ=255)
                                          gTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACaa                 <  1:1964091/69‑1 (MQ=255)
                                          gTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACaa                 <  1:2704347/69‑1 (MQ=255)
                                                     tGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:231047/69‑1 (MQ=255)
                                                     tGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:663916/69‑1 (MQ=255)
                                                     tGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:3056787/69‑1 (MQ=255)
                                                     tGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:1681674/69‑1 (MQ=255)
                                                     tGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:3127683/69‑1 (MQ=255)
                                                       cATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAcc     >  1:2242128/1‑68 (MQ=255)
                                                       cATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAcc     >  1:1956150/1‑68 (MQ=255)
                                                         tttACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCACCATg  <  1:852266/69‑1 (MQ=255)
                                                         tttACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCACCATg  <  1:285467/69‑1 (MQ=255)
                                                         tttACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCACCATg  <  1:245293/69‑1 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGc                           >  1:2730033/1‑42 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGc                           >  1:3390983/1‑42 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCa       >  1:1790293/1‑62 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      >  1:2184991/1‑63 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      >  1:2737254/1‑63 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:3168874/63‑1 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:529805/63‑1 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      >  1:1703972/1‑63 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:1852829/63‑1 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:2504529/63‑1 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      >  1:3097334/1‑63 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:2920470/63‑1 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:2187490/63‑1 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:1545085/63‑1 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      >  1:2145850/1‑63 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      >  1:2737841/1‑63 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACAATCAc      >  1:44432/1‑63 (MQ=255)
                                                           tACCTGCTACTTTTTAGTAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:867079/62‑1 (MQ=255)
                                                            aCCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCAc      <  1:2094120/62‑1 (MQ=255)
                                                               |                                                              
TCGTTTTTGCCGGACAGTGCTGCCGTTTATTTTCGTGATCCAGTTAAAGTAAATGCATTTACCTGCTACTTTTTAGTAAAAATTTTACTAAACTCCCCAGCAATTACACAAACTACCATCACCATG  >  minE/2073187‑2073312

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: