Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,608,350 G→A 32.7% A84A (GCC→GCT nikR ← DNA‑binding transcriptional regulator, Ni‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,608,3500GA32.7% 120.1 / 72.5 137A84A (GCC→GCTnikRDNA‑binding transcriptional regulator, Ni‑binding
Reads supporting (aligned to +/- strand):  ref base G (18/74);  new base A (24/21);  total (42/95)
Fisher's exact test for biased strand distribution p-value = 1.31e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

CATGTCACCTTTCAACACGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGGAGACAATGCGGCTGGCTAAGTC  >  minE/2608287‑2608406
                                                               |                                                        
cATGTCACCTTTCAACACGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGa                                                      >  1:482842/1‑68 (MQ=255)
    tCACCTTTCAACACGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGgag                                                 >  1:1403351/1‑69 (MQ=255)
           tCAACACGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCgt                                          >  1:976270/1‑69 (MQ=255)
           tCAACACGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCgt                                          >  1:937469/1‑69 (MQ=255)
           tCAACACGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCg                                           >  1:672569/1‑68 (MQ=255)
           tCAACACGGCGATTTCCAGACAGTCGTCGTGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCgt                                          >  1:2371917/1‑69 (MQ=255)
           tCAACACGGCGATTTCCAGACAGTCGTCGTGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCgt                                          >  1:63967/1‑69 (MQ=255)
           tCAACACGGCGATTTCCAGACAGTCGTCGTGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCgt                                          >  1:2196958/1‑69 (MQ=255)
           tCAACACGGCGATTTCCAGACAGTCGTCGTGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCgt                                          >  1:1962500/1‑69 (MQ=255)
           tCAACACGGCGATTTCCAGACAGTCGTCGTGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCgt                                          >  1:1665106/1‑69 (MQ=255)
           tCAACACGGCGATTTCCAGACAGTCGTCGTGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCgt                                          >  1:114437/1‑69 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:1917526/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:2410845/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:2437349/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:2664431/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:2722017/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:2198654/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:2728845/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:3178656/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:2048962/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:3199509/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:1172553/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:3328935/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:1874796/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:338830/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:344498/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:531742/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:782643/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:657747/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:1083721/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:152634/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:44578/69‑1 (MQ=255)
                 cGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgat                                    <  1:371080/69‑1 (MQ=255)
                   gCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:195223/69‑1 (MQ=255)
                   gCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:101448/69‑1 (MQ=255)
                   gCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:2555585/69‑1 (MQ=255)
                   gCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:1415552/69‑1 (MQ=255)
                   gCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:333752/69‑1 (MQ=255)
                   gCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:2739746/69‑1 (MQ=255)
                   gCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:21343/69‑1 (MQ=255)
                   gCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:1235241/69‑1 (MQ=255)
                   gCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:2066473/69‑1 (MQ=255)
                   gCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:429486/69‑1 (MQ=255)
                   gCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:3282415/69‑1 (MQ=255)
                   gCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:1590093/69‑1 (MQ=255)
                    cGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:765458/68‑1 (MQ=255)
                      aTTTCCAGACAGTCGTCGTTGTTGATGTTCACATGCAGCGTGGCGACGGAGAGGTCGTGGttat                                    <  1:635422/64‑1 (MQ=255)
                              aCAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATTCTGGGTGga                       <  1:2103093/69‑1 (MQ=255)
                              aCAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGga                       <  1:2739458/69‑1 (MQ=255)
                              aCAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGga                       <  1:53113/69‑1 (MQ=255)
                              aCAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGga                       <  1:1034208/69‑1 (MQ=255)
                              aCAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGga                       <  1:3434987/69‑1 (MQ=255)
                              aCAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGga                       <  1:1541643/69‑1 (MQ=255)
                              aCAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGga                       <  1:2136495/69‑1 (MQ=255)
                              aCAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGga                       <  1:1834952/69‑1 (MQ=255)
                              aCAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGga                       <  1:1994753/69‑1 (MQ=255)
                              aCAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTCGGTGga                       <  1:440036/69‑1 (MQ=255)
                                 gtcgtcGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGGAGAc                    <  1:888066/69‑1 (MQ=255)
                                    gtcgtGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACAAt                 >  1:6961/1‑69 (MQ=255)
                                    gtcgtGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACAAt                 >  1:1490373/1‑69 (MQ=255)
                                    gtcgtGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACAAt                 >  1:149822/1‑69 (MQ=255)
                                    gtcgtGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACAAt                 >  1:825004/1‑69 (MQ=255)
                                    gtcgtGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACAAt                 >  1:85565/1‑69 (MQ=255)
                                    gtcgtGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACAAt                 >  1:1645435/1‑69 (MQ=255)
                                    gtcgtGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACAAt                 >  1:1759075/1‑69 (MQ=255)
                                    gtcgtGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACAAt                 >  1:912637/1‑69 (MQ=255)
                                    gtcgtGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACAAt                 >  1:2083835/1‑69 (MQ=255)
                                    gtcgtGATTAATGTGCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACAAt                 >  1:2818169/1‑69 (MQ=255)
                                          gTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGGAGACAATGCggct           >  1:2281802/1‑69 (MQ=255)
                                          gTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGGAGACAATGCggct           >  1:1918544/1‑69 (MQ=255)
                                          gTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGGAGACAATGCggct           >  1:116872/1‑69 (MQ=255)
                                          gTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGGAGACAATGCggct           >  1:47114/1‑69 (MQ=255)
                                            tGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGGAGACAATGCggctg          >  1:3143892/1‑68 (MQ=255)
                                            tGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGGAGACAATGCggctg          >  1:2138068/1‑68 (MQ=255)
                                            tGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGGAGACAATGCggctg          >  1:2969121/1‑68 (MQ=255)
                                             gATGTGCACATGCAGCGTGGCGACGGAGAGGTCgtggt                                       <  1:17314/38‑1 (MQ=255)
                                             gATGTGCACATGCAGCGTGGCGACGGAGAGGTCgtggt                                       <  1:2801272/38‑1 (MQ=255)
                                             gATGTGCACATGCAGCGTGGCGACGGAGAGGTCgtggt                                       <  1:236473/38‑1 (MQ=255)
                                             gATGTGCACATGCAGCGTGGCGACGGAGAGGTCgtggt                                       <  1:2802021/38‑1 (MQ=255)
                                             gATGTGCACATGCAGCGTGGCGACGGAGAGGTCgtggt                                       <  1:1546594/38‑1 (MQ=255)
                                             gATGTGCACATGCAGCGTGGCGACGGAGAGGTCgtggt                                       <  1:3410160/38‑1 (MQ=255)
                                             gATGTGCACATGCAGCGTGGCGACGGAGAGGTCgtggt                                       <  1:3018197/38‑1 (MQ=255)
                                             gATGTGCACATGCAGCGTGGCGACGGAGAGGTCgtggt                                       <  1:2567313/38‑1 (MQ=255)
                                             gATGTGCACATGCAGCGTGGCGACGGAGAGGTCgtggt                                       <  1:128624/38‑1 (MQ=255)
                                             gATGTGCACATGCAGCGTGGCGACGGAGAGGTCgtggt                                       <  1:1465645/38‑1 (MQ=255)
                                                 tgCACATGCAGAGTAGCAACTGAGAGATCTTGGTGATGATGCTTGGTGGAGACa                   <  1:1260820/54‑1 (MQ=255)
                                                 tgCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACa                   <  1:973421/54‑1 (MQ=37)
                                                 tgCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACa                   <  1:1349078/54‑1 (MQ=37)
                                                 tgCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACa                   <  1:1619454/54‑1 (MQ=37)
                                                 tgCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACa                   <  1:1406988/54‑1 (MQ=37)
                                                 tgCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACa                   <  1:603000/54‑1 (MQ=37)
                                                 tgCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACa                   <  1:1502119/54‑1 (MQ=37)
                                                 tgCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACa                   <  1:1831189/54‑1 (MQ=37)
                                                 tgCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACa                   <  1:2201790/54‑1 (MQ=37)
                                                 tgCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACa                   <  1:2323890/54‑1 (MQ=37)
                                                  gCACATGCAGAGTAGCAACTGAGAGATCGTGGTGATGATGCTGGGTGGAGACa                   <  1:1036915/53‑1 (MQ=25)
                                                   cacaTGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGGAGACAATGCGGCTGGCTAAGTc  >  1:1074587/1‑69 (MQ=255)
                                                   cacaTGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGGAGACAATGCGGCTGGCTAAGTc  >  1:1980633/1‑69 (MQ=255)
                                                   cacaTGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGGAGACAATGCGGCTGGCTAAGTc  >  1:945277/1‑69 (MQ=255)
                                                    acaTGCAGCGTGGCGACGGAGAGGTCGTGGtgatga                                  <  1:507178/36‑1 (MQ=255)
                                                               |                                                        
CATGTCACCTTTCAACACGGCGATTTCCAGACAGTCGTCGTGGTTGATGTGCACATGCAGCGTGGCGACGGAGAGGTCGTGGTGATGATGCTGGGTGGAGACAATGCGGCTGGCTAAGTC  >  minE/2608287‑2608406

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: