Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,293,306 C→T 68.1% S307N (AGC→AAC)  cbl ← DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,293,3060CT68.1% 37.8 / 23.3 50S307N (AGC→AAC) cblDNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (10/6);  new base T (23/11);  total (33/17)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGT  >  minE/1293239‑1293359
                                                                   |                                                     
aTTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCt                                                      >  1:1011229/1‑69 (MQ=255)
aTTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCt                                                      >  1:2047359/1‑69 (MQ=255)
 ttCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCtt                                                     <  1:1065319/69‑1 (MQ=255)
 ttCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCtt                                                     <  1:993813/69‑1 (MQ=255)
          aGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCtt                                                     <  1:920450/60‑1 (MQ=255)
          aGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCtt                                                     <  1:170342/60‑1 (MQ=255)
          aGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCtt                                                     <  1:562001/60‑1 (MQ=255)
          aGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCtt                                                     <  1:3368016/60‑1 (MQ=255)
                                         tCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCt              >  1:1474115/1‑68 (MQ=255)
                                         tCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTg             >  1:1442042/1‑69 (MQ=255)
                                         tCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTg             >  1:763880/1‑69 (MQ=255)
                                         tCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTg             >  1:199931/1‑69 (MQ=255)
                                           tGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGc            >  1:3075710/1‑68 (MQ=255)
                                              tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                <  1:734395/61‑1 (MQ=255)
                                              tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                <  1:2385126/61‑1 (MQ=255)
                                              tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                <  1:1242511/61‑1 (MQ=255)
                                              tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                <  1:2882802/61‑1 (MQ=255)
                                              tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                <  1:3253787/61‑1 (MQ=255)
                                              tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                <  1:2662745/61‑1 (MQ=255)
                                              tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                <  1:838536/61‑1 (MQ=255)
                                              tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                <  1:1212435/61‑1 (MQ=255)
                                              tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                <  1:1717882/61‑1 (MQ=255)
                                              tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                <  1:2138765/61‑1 (MQ=255)
                                              tAATCAATTTCCTCTTCACTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTc                <  1:1551131/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:2304063/1‑52 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTAc       >  1:89404/1‑52 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:1685529/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:805333/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:2277630/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:1398729/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:721442/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:715427/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:667560/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:597326/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:1797702/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:1240379/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:2074863/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:3194037/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:2093034/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:2689226/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:2218375/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:2423296/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:2408775/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:2253072/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATTACa      >  1:233913/1‑53 (MQ=39)
                                                                cTGTTTTCCATCACCTGGCGCTTGATATCCTCAACTGACAGTCCTGCATGACa      >  1:3071955/1‑53 (MQ=38)
                                                                cTGTTTTCCATCACCTGGCGCTGGATATCCTCAACTGACAGTCCTGCATTACa      >  1:3122992/1‑53 (MQ=38)
                                                                   |                                                     
ATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGT  >  minE/1293239‑1293359

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: