Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,072,881 C→A 16.3% P235T (CCG→ACG)  ebgR → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,072,8810CA16.3% 169.7 / 14.3 61P235T (CCG→ACG) ebgRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (28/23);  new base A (6/4);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

GCGCGGCGGTTTTTCCAGTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAATTCATGAACGAGGCC  >  minE/2072814‑2072950
                                                                   |                                                                     
gcgcGGCGGTTTTTCCAGTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCg                                                                     <  1:1897882/70‑1 (MQ=255)
 cgcgGCGGTTTTTCCAGTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATcc                                                                      >  1:1254708/1‑68 (MQ=255)
 cgcgGCGGTTTTTCCAGTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGaa                                                                   <  1:1405886/71‑1 (MQ=255)
   cggcggTTTTTCCAGTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAgg                                                                 <  1:37087/71‑1 (MQ=255)
         tttttCCAGTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACtgt                                                           >  1:2871061/1‑71 (MQ=255)
          ttttCCAGTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACtgtt                                                          >  1:1300068/1‑71 (MQ=255)
            ttCCAGTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACtgtttg                                                        >  1:1442714/1‑71 (MQ=255)
             tCCAGTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATACGAAGGCACtgtttg                                                        <  1:3746635/70‑1 (MQ=255)
              ccAGTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACtgtttgtt                                                      >  1:1074498/1‑71 (MQ=255)
              ccAGTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACtgtttgt                                                       <  1:1719295/70‑1 (MQ=255)
              ccAGTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACtgtttg                                                        <  1:3159979/69‑1 (MQ=255)
                 gTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCt                                                   >  1:2865978/1‑71 (MQ=255)
                  ttCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATc                                                                       <  1:3605870/50‑1 (MQ=255)
                  ttCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCtt                                                  >  1:1758208/1‑71 (MQ=255)
                      tcgGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGa                                              >  1:673276/1‑71 (MQ=255)
                         ggTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTc                                           <  1:2872583/71‑1 (MQ=255)
                           ttATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGa                                                                    >  1:1974280/1‑44 (MQ=255)
                           ttATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGa                                                                    >  1:3357628/1‑44 (MQ=255)
                           ttATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCa                                         <  1:2427224/71‑1 (MQ=255)
                            tATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTc                                                 >  1:876130/1‑62 (MQ=255)
                              tGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTg                                      >  1:3130993/1‑71 (MQ=255)
                               gAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCt                                                   <  1:703837/57‑1 (MQ=255)
                                 aCTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACtgtttgt                                                       <  1:3252624/51‑1 (MQ=255)
                                 aCTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTa                                   <  1:1700328/71‑1 (MQ=255)
                                 aCTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTa                                   <  1:3648844/71‑1 (MQ=255)
                                   tGGCAAAACAAATGCTGGCGCGGGAAGACTATCCg                                                                     <  1:99526/35‑1 (MQ=255)
                                   tGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGaa                                                                   >  1:1453035/1‑37 (MQ=255)
                                   tGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGaa                                                                   >  1:2738552/1‑37 (MQ=255)
                                   tGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCa                                         >  1:3391260/1‑63 (MQ=255)
                                   tGGAAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATc                                 >  1:3465894/1‑71 (MQ=255)
                                     gCAAAACAAATGCTGGCGCGGGAAGACTATACGAAGGCACtgtttgtt                                                      >  1:1085380/1‑48 (MQ=255)
                                       aaaaCAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGc                              >  1:1298437/1‑70 (MQ=255)
                                         aaCAAATGCTGGCGCGGGAAGACTATCCGAAGGCACtgt                                                           >  1:1200843/1‑39 (MQ=255)
                                          aCAAATGCTGGCGCGGGAAGACTATCCGAAGGCACtgtttgt                                                       >  1:1629416/1‑42 (MQ=255)
                                          aCAAATGCTGGCGCGGGAAGACTATACGAAGGCACtgtttgtt                                                      >  1:3471521/1‑43 (MQ=255)
                                          aCAAATGCTGGCGCGGGAAGACTATACGAAGGCACtgtttgtt                                                      >  1:3216080/1‑43 (MQ=255)
                                          aCAAATGCTGGCGCGGGAAGACTATACGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTAc                          <  1:1622626/71‑1 (MQ=255)
                                           cAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGc                              >  1:2059113/1‑66 (MQ=255)
                                            aaaTGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCa                                         <  1:834378/54‑1 (MQ=255)
                                             aaTGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGc                       <  1:2419294/71‑1 (MQ=255)
                                              aTGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCg                      >  1:3952681/1‑71 (MQ=255)
                                              aTGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCg                      >  1:3477610/1‑71 (MQ=255)
                                              aTGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCg                      >  1:206600/1‑71 (MQ=255)
                                                 cTGGCGCGGGAAGACTATACGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCggg                    >  1:2611669/1‑70 (MQ=255)
                                                    gcgcGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCa                  <  1:3232696/69‑1 (MQ=255)
                                                    gcgcGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAAt                <  1:185058/71‑1 (MQ=255)
                                                    gcgcGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAAt                <  1:3453220/71‑1 (MQ=255)
                                                    gcgcGGGAAGACTATACGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCaa                 <  1:3284072/70‑1 (MQ=255)
                                                     cgcgGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAAtt               >  1:973466/1‑71 (MQ=255)
                                                      gcgGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGTCGTACTGCGGGCAAtt               <  1:88217/70‑1 (MQ=255)
                                                      gcgGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAAtt               >  1:1201872/1‑70 (MQ=255)
                                                      gcgGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAATTc              <  1:2240311/71‑1 (MQ=255)
                                                       cgGGAAGACTATACGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAATTCa             <  1:565971/71‑1 (MQ=255)
                                                        gggAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAATTCAt            >  1:1210588/1‑71 (MQ=255)
                                                         ggAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAATTCATg           <  1:1011626/71‑1 (MQ=255)
                                                            aGACTATACGAAGGCACTGTTTGTTGCTTCCGATTCCATTGc                                     >  1:1559072/1‑42 (MQ=255)
                                                            aGACTATACGAAGGCACTGTTTGTTGCTTCCGATTCCATTGc                                     >  1:1764436/1‑42 (MQ=255)
                                                              aCTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCg                             <  1:1108627/48‑1 (MQ=255)
                                                                tATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAATTCATGAACGa      >  1:2818459/1‑69 (MQ=255)
                                                                  tCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCTGCGTAc                          >  1:2166367/1‑47 (MQ=255)
                                                                  tCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAATTCATGAACGAgg    <  1:3265149/69‑1 (MQ=255)
                                                                  tCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAATTCATGAACGAgg    <  1:2245572/69‑1 (MQ=255)
                                                                  tCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAATTCATGAACGAGGcc  >  1:3563843/1‑71 (MQ=255)
                                                                   |                                                                     
GCGCGGCGGTTTTTCCAGTTCGTCGGGTTATGAACTGGCAAAACAAATGCTGGCGCGGGAAGACTATCCGAAGGCACTGTTTGTTGCTTCCGATTCCATTGCTATCGGCGTACTGCGGGCAATTCATGAACGAGGCC  >  minE/2072814‑2072950

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: