Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,722,403 G→A 65.9% intergenic (‑152/‑150) crp ← / → yhfA DNA‑binding transcriptional dual regulator/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,722,4030GA65.9% 82.8 / 65.6 114intergenic (‑152/‑150)crp/yhfADNA‑binding transcriptional dual regulator/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/14);  new base A (29/46);  total (54/60)
Fisher's exact test for biased strand distribution p-value = 1.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

AAGGGGCTATCAACTGTACTGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTCCTGGTGTCTCGCTTCAG  >  minE/2722340‑2722470
                                                               |                                                                   
aaGGGGCTATCAACTGTACTGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTAtt                                                                 <  1:1789184/68‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:2987332/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:1663509/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:1948900/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:1982234/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:2031030/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:2056944/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:2331643/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:2484715/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:251033/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:2649550/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:2517657/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:3397777/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:358576/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:3713/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:399689/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:551815/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:679161/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:710024/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:935127/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:986841/69‑1 (MQ=255)
  ggggCTATCAACTGTACTGCACGGTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACt                                                              <  1:1546612/69‑1 (MQ=255)
           aaCTGTACTGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATc                                                      >  1:2563186/1‑68 (MQ=255)
           aaCTGTACTGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATc                                                      >  1:3159764/1‑68 (MQ=255)
           aaCTGTACTGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATc                                                      >  1:1085225/1‑68 (MQ=255)
             cTGTACTGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATCCt                                                    <  1:2773581/68‑1 (MQ=255)
             cTGTACTGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATCCt                                                    <  1:2831998/68‑1 (MQ=255)
               gTACTGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATCCTGAc                                                 >  1:3357231/1‑69 (MQ=255)
               gTACTGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATCCTGAc                                                 >  1:1711441/1‑69 (MQ=255)
               gTACTGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATCCTGAc                                                 >  1:2922184/1‑69 (MQ=255)
                  cTGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGt                                              <  1:1598348/69‑1 (MQ=255)
                  cTGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGt                                              <  1:202315/69‑1 (MQ=255)
                   tGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGt                                              <  1:1323737/68‑1 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTa                                          >  1:956484/1‑67 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTa                                          >  1:1740807/1‑67 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTa                                          >  1:2599275/1‑67 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:2883145/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:2898545/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:3157599/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:2752711/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:3179071/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:3105634/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:3092179/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:3393633/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:2507093/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:1416604/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:1321802/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:1725968/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:971212/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:748165/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:1560250/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:630320/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:151340/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:2110988/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:606949/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:2157312/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:1146762/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:2332169/1‑68 (MQ=255)
                        ggTAATGTGGCGTCCTTTGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAg                                         >  1:2459951/1‑68 (MQ=255)
                            aTGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGt                                              <  1:1093741/59‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:3292994/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:592761/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:3328509/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:1211838/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:655104/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:3368618/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:1225071/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:1056569/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:1726230/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:1563100/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:2699895/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:2009771/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:2575183/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:1883501/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:1015329/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:2288717/69‑1 (MQ=255)
                                        ttGCATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:1842161/69‑1 (MQ=255)
                                        ttGCATAAATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:2821145/69‑1 (MQ=255)
                                           cATACATGCAGTACATCAATATATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTtgt                        <  1:1862710/66‑1 (MQ=255)
                                              aCATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTcc                  >  1:1520481/1‑69 (MQ=255)
                                              aCATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTcc                  >  1:1146626/1‑69 (MQ=255)
                                              aCATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTcc                  >  1:886695/1‑69 (MQ=255)
                                              aCATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTcc                  >  1:934011/1‑69 (MQ=255)
                                              aCATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTcc                  >  1:1057631/1‑69 (MQ=255)
                                              aCATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTcc                  >  1:952160/1‑69 (MQ=255)
                                              aCATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTcc                  >  1:1615336/1‑69 (MQ=255)
                                              aCATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTcc                  >  1:2136394/1‑69 (MQ=255)
                                              aCATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTcc                  >  1:2240476/1‑69 (MQ=255)
                                              aCATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTcc                  >  1:544161/1‑69 (MQ=255)
                                              aCATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTcc                  >  1:1375842/1‑69 (MQ=255)
                                              aCATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTcc                  >  1:3136168/1‑69 (MQ=255)
                                              aCATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTcc                  >  1:2991858/1‑69 (MQ=255)
                                               cATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGtctc                   <  1:2390147/67‑1 (MQ=255)
                                                 tGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTCCTgg               >  1:227808/1‑69 (MQ=255)
                                                 tGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTCCTgg               >  1:2757681/1‑69 (MQ=255)
                                                   cAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTCCTGgtg             >  1:885937/1‑69 (MQ=255)
                                                   cAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTCCTGgtg             >  1:2533752/1‑69 (MQ=255)
                                                              tGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTCCTGGTGTCTCGCTTCAg  <  1:3078869/69‑1 (MQ=255)
                                                              tGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTCCTGGTGTCTCGCTTCAg  <  1:1364334/69‑1 (MQ=255)
                                                               gTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTCCTGGTGTCTCGCTTCAg  <  1:1481875/68‑1 (MQ=255)
                                                               |                                                                   
AAGGGGCTATCAACTGTACTGCACGGTAATGTGACGTCCTTTGCATACATGCAGTACATCAATGTATTACTGTAGCATCCTGACTGTTTTAGCATAGCTTTCGCTTTGTGTCTCCTGGTGTCTCGCTTCAG  >  minE/2722340‑2722470

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: