Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,293,320 C→G 66.7% R302R (CGG→CGC cbl ← DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,293,3200CG66.7% 46.5 / 37.1 51R302R (CGG→CGCcblDNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (10/7);  new base G (23/11);  total (33/18)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

TCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTA  >  minE/1293280‑1293383
                                        |                                                               
tCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCt                                      >  1:1474115/1‑68 (MQ=255)
tCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTg                                     >  1:1442042/1‑69 (MQ=255)
tCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTg                                     >  1:199931/1‑69 (MQ=255)
tCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTg                                     >  1:763880/1‑69 (MQ=255)
  tGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGc                                    >  1:3075710/1‑68 (MQ=255)
     tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                                        <  1:1717882/61‑1 (MQ=255)
     tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                                        <  1:2882802/61‑1 (MQ=255)
     tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                                        <  1:3253787/61‑1 (MQ=255)
     tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                                        <  1:2138765/61‑1 (MQ=255)
     tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                                        <  1:2385126/61‑1 (MQ=255)
     tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                                        <  1:1212435/61‑1 (MQ=255)
     tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                                        <  1:734395/61‑1 (MQ=255)
     tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                                        <  1:1551131/61‑1 (MQ=255)
     tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                                        <  1:838536/61‑1 (MQ=255)
     tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                                        <  1:1242511/61‑1 (MQ=255)
     tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                                        <  1:2662745/61‑1 (MQ=255)
              tactctTCACTGCTTTCCATCACCTGCCGCTTGATATCCTc                                                   >  1:2726557/3‑41 (MQ=255)
              tCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGt                          >  1:1468813/1‑66 (MQ=255)
              tCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGt                          >  1:2827880/1‑66 (MQ=255)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTAc                               >  1:89404/1‑52 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTAc                               >  1:2304063/1‑52 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:3194037/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:597326/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:667560/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:715427/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:721442/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:805333/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:2689226/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:2408775/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:2218375/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:1240379/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:1398729/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:1685529/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:1797702/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:2074863/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:2093034/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:2253072/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:2277630/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:233913/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:2423296/1‑53 (MQ=39)
                       cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATGACa                              >  1:3071955/1‑53 (MQ=38)
                       cTGTTTTCCATCACCTGGCGCTGGATATCCTCAACTGACAGTCCTGCATTACa                              >  1:3122992/1‑53 (MQ=38)
                               cATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTTCATTACAAAGTTCCAGAAAGCGCCAGACAt       <  1:254215/68‑1 (MQ=255)
                               cATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACAt       <  1:314052/68‑1 (MQ=255)
                               cATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACAt       <  1:3165548/68‑1 (MQ=255)
                                 tCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGt    <  1:2051130/69‑1 (MQ=255)
                                 tCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGt    <  1:1303972/69‑1 (MQ=255)
                                 tCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGt    <  1:2456183/69‑1 (MQ=255)
                                   aCCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTa  >  1:1820538/1‑69 (MQ=255)
                                   aCCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTACAGAAAGCGCCAGACATAGTTa  >  1:2197699/1‑69 (MQ=255)
                                    ccTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACAt       <  1:2977227/63‑1 (MQ=255)
                                        |                                                               
TCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTA  >  minE/1293280‑1293383

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: