Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,103,531 C→T 46.3% T74T (ACG→ACA agaR ← DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,103,5310CT46.3% 5.8 / 48.3 61T74T (ACG→ACAagaRDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (12/20);  new base T (9/19);  total (22/39)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

TCATCGCGGTGTTCAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAG  >  minE/2103476‑2103599
                                                       |                                                                    
tCATCGCGGTGTTCAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAg                                                          >  1:2640539/1‑68 (MQ=255)
   tCGCGGTGTTCAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATc                                                       <  1:978597/68‑1 (MQ=255)
   tCGCGGTGTTCAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATc                                                       <  1:2838524/68‑1 (MQ=255)
    cgcgGTGTTCAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATc                                                       <  1:1127829/67‑1 (MQ=255)
             cAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATGGCAGATCAACGCGccac                                             <  1:2960914/68‑1 (MQ=255)
             cAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGccac                                             <  1:2554341/68‑1 (MQ=255)
             cAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGccac                                             <  1:906091/68‑1 (MQ=255)
             cAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGccac                                             <  1:51036/68‑1 (MQ=255)
             cAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGccac                                             <  1:3431649/68‑1 (MQ=255)
             cAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGccac                                             <  1:3175162/68‑1 (MQ=255)
             cAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGccac                                             <  1:1747773/68‑1 (MQ=255)
             cAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGccac                                             <  1:1314987/68‑1 (MQ=255)
                 gcgcTTTTATCTTCCACTGATGGCTCGACTGACGGCGTTGTGTTATCGCAGATc                                                       >  1:329943/1‑54 (MQ=255)
                 gcgcTTTTATCTTCCACTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATc                                                       >  1:1474577/1‑54 (MQ=255)
                 gcgcTTTTATCTTCCACTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATc                                                       >  1:11400/1‑54 (MQ=255)
                 gcgcTTTTATCTTCCACTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATc                                                       >  1:947485/1‑54 (MQ=255)
                 gcgcTTTTATCTTCCACTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATc                                                       >  1:819688/1‑54 (MQ=255)
                 gcgcTTTTATCTTCCACTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATc                                                       >  1:3429462/1‑54 (MQ=255)
                 gcgcTTTTATCTTCCACTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATc                                                       >  1:2747276/1‑54 (MQ=255)
                 gcgcTTTTATCTTCCACTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATc                                                       >  1:1839729/1‑54 (MQ=255)
                 gcgcTTTTATCTTCCACTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATc                                                       >  1:1740327/1‑54 (MQ=255)
                      tttATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGgcacg                                   <  1:3284981/69‑1 (MQ=255)
                      tttATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGgcacg                                   <  1:490366/69‑1 (MQ=255)
                      tttATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGgcacg                                   <  1:2269678/69‑1 (MQ=255)
                      tttATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGgcacg                                   <  1:607664/69‑1 (MQ=255)
                      tttATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGgcacg                                   <  1:698568/69‑1 (MQ=255)
                        tATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGgcacg                                   >  1:1359964/1‑67 (MQ=255)
                        tATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGgcacg                                   >  1:2444082/1‑67 (MQ=255)
                        tATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGgcacg                                   >  1:824196/1‑67 (MQ=255)
                        tATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGgcacg                                   >  1:861892/1‑67 (MQ=255)
                         aTCTTCCACTGATGGCTCGACTGACGGCGTCGTTCTATc                                                              >  1:2599739/1‑39 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:611837/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:1147484/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:1143759/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:910935/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:401605/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:1417193/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:3193230/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:3138508/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:124988/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:2662335/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:2657014/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:2336000/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:2396972/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:2096116/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:2499782/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:1988615/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:1849616/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:1692852/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTTGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATc                         <  1:1278662/68‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATcc                        <  1:783908/69‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATcc                        <  1:1220414/69‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATcc                        <  1:1646644/69‑1 (MQ=255)
                                 cTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATcc                        <  1:2825327/69‑1 (MQ=255)
                                       gCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATcccc                      >  1:1576453/1‑65 (MQ=255)
                                       gCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCt                  >  1:2114390/1‑69 (MQ=255)
                                       gCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCt                  >  1:55411/1‑69 (MQ=255)
                                       gCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCt                  >  1:1485277/1‑69 (MQ=255)
                                       gCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCt                  >  1:1836988/1‑69 (MQ=255)
                                       gCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCt                  >  1:3104072/1‑69 (MQ=255)
                                       gCTCGACTGACGGCGTAGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCt                  >  1:2498658/1‑69 (MQ=255)
                                                       cgtGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAg  <  1:2136204/69‑1 (MQ=255)
                                                       cgtGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAg  <  1:1183633/69‑1 (MQ=255)
                                                       cgtGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAg  <  1:1529681/69‑1 (MQ=255)
                                                       cgtGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAg  <  1:719789/69‑1 (MQ=255)
                                                       cgtGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAg  <  1:739501/69‑1 (MQ=255)
                                                       cgtGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAg  <  1:2609856/69‑1 (MQ=255)
                                                       cgtGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAg  <  1:1503393/69‑1 (MQ=255)
                                                       cgtGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAg  <  1:191945/69‑1 (MQ=255)
                                                       cgtGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAg  <  1:309346/69‑1 (MQ=255)
                                                       cgtGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAg  <  1:1226455/69‑1 (MQ=255)
                                                       cgtGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAg  <  1:161861/69‑1 (MQ=255)
                                                       cgtGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAg  <  1:2675958/69‑1 (MQ=255)
                                                       |                                                                    
TCATCGCGGTGTTCAGTGCGCTTTTATCTTCCACTGATGGCTCGACTGACGGCGTCGTGCTATCGCAGATCAACGCGCCACCATAGGCACGCACAGCGATCCCCTGCTTTTCCAGAAACGCCAG  >  minE/2103476‑2103599

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: