Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,343,955 A→C 21.5% S207S (TCA→TCC glnG → fused DNA‑binding response regulator in two‑component regulatory system with GlnL, nitrogen regulator I (NRI)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,343,9550AC21.5% 130.3 / 33.9 79S207S (TCA→TCCglnGfused DNA‑binding response regulator in two‑component regulatory system with GlnL, nitrogen regulator I (NRI)
Reads supporting (aligned to +/- strand):  ref base A (24/38);  new base C (15/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 2.97e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GTCCGCGCGCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTGAACA  >  minE/2343888‑2344023
                                                                   |                                                                    
gTCCGCGCGCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAg                                                                     >  1:2550343/1‑69 (MQ=255)
gTCCGCGCGCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAg                                                                     >  1:1871613/1‑69 (MQ=255)
gTCCGCGCGCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAg                                                                     >  1:3200797/1‑69 (MQ=255)
gTCCGCGCGCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAg                                                                     >  1:2191544/1‑69 (MQ=255)
gTCCGCGCGCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAg                                                                     >  1:2840404/1‑69 (MQ=255)
   cgcgcgCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGa                                                                    <  1:133782/67‑1 (MQ=255)
   cgcgcgCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGa                                                                    <  1:1390743/67‑1 (MQ=255)
   cgcgcgCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGa                                                                    <  1:472172/67‑1 (MQ=255)
   cgcgcgCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGa                                                                    <  1:1187467/67‑1 (MQ=255)
   cgcgcgCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGa                                                                    <  1:30027/67‑1 (MQ=255)
   cgcgcgCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGa                                                                    <  1:3045352/67‑1 (MQ=255)
   cgcgcgCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGa                                                                    <  1:896196/67‑1 (MQ=255)
   cgcgcgCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAAc                                                                  >  1:1344122/1‑69 (MQ=255)
   cgcgcgCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAAc                                                                  >  1:869382/1‑69 (MQ=255)
   cgcgcgCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAAc                                                                  >  1:3292539/1‑69 (MQ=255)
     cgcgcCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTg                                                                <  1:519497/69‑1 (MQ=255)
               cgcCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGcc                                                         >  1:411041/1‑66 (MQ=255)
                  cGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACgaga                                                   <  1:2655344/69‑1 (MQ=255)
                  cGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACgaga                                                   <  1:2839195/69‑1 (MQ=255)
                  cGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACgaga                                                   <  1:994434/69‑1 (MQ=255)
                  cGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACgaga                                                   <  1:3336976/69‑1 (MQ=255)
                      tATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAgg                                               >  1:416172/1‑69 (MQ=255)
                      tATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAgg                                               >  1:2466493/1‑69 (MQ=255)
                       aTCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGgc                                              >  1:2893396/1‑69 (MQ=255)
                       aTCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGgc                                              >  1:37460/1‑69 (MQ=255)
                            gcTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGtt                                         <  1:2125558/69‑1 (MQ=255)
                            gcTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGtt                                         <  1:2133681/69‑1 (MQ=255)
                            gcTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGtt                                         <  1:852066/69‑1 (MQ=255)
                                      gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                               <  1:3273378/69‑1 (MQ=255)
                                      gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                               <  1:3407160/69‑1 (MQ=255)
                                      gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                               <  1:29761/69‑1 (MQ=255)
                                      gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                               <  1:653244/69‑1 (MQ=255)
                                      gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                               <  1:2774305/69‑1 (MQ=255)
                                      gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                               <  1:1157286/69‑1 (MQ=255)
                                      gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                               <  1:2340768/69‑1 (MQ=255)
                                           tATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATAc                          <  1:2514510/69‑1 (MQ=255)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:1031823/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:1551141/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:1741576/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:1781850/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:1827942/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:3076302/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:1957089/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:2982820/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:2896671/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:2289714/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:2623740/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:2408691/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:2326230/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                         >  1:27543/1‑69 (MQ=25)
                                            aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAc                          >  1:2671945/1‑68 (MQ=25)
                                               ccAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCAtt                      >  1:118690/1‑69 (MQ=255)
                                               ccAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCAtt                      >  1:386046/1‑69 (MQ=255)
                                               ccAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCAtt                      >  1:3347594/1‑69 (MQ=255)
                                               ccAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCAtt                      >  1:647299/1‑69 (MQ=255)
                                               ccAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCAtt                      >  1:2798457/1‑69 (MQ=255)
                                               ccAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCAtt                      >  1:2963908/1‑69 (MQ=255)
                                                       tttGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAggg              >  1:679080/1‑69 (MQ=255)
                                                       tttGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAggg              >  1:955934/1‑69 (MQ=255)
                                                       tttGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAggg              >  1:208336/1‑69 (MQ=255)
                                                       attGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAggg              >  1:942335/2‑69 (MQ=255)
                                                        tcgATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATACCAt                       <  1:1782663/57‑1 (MQ=25)
                                                         tGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATACCAt                       <  1:2146601/58‑1 (MQ=25)
                                                                aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:1766870/68‑1 (MQ=255)
                                                                aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:1823854/68‑1 (MQ=255)
                                                                aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:631143/68‑1 (MQ=255)
                                                                aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:1322134/68‑1 (MQ=255)
                                                                aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:1199424/68‑1 (MQ=255)
                                                                aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:2214205/68‑1 (MQ=255)
                                                                aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:3314843/68‑1 (MQ=255)
                                                                aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:1843907/68‑1 (MQ=255)
                                                                aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:1869235/68‑1 (MQ=255)
                                                                aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:3211704/68‑1 (MQ=255)
                                                                aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:1940411/68‑1 (MQ=255)
                                                                aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:3115739/68‑1 (MQ=255)
                                                                aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:1967097/68‑1 (MQ=255)
                                                                aTCAGAACTGTATGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:464777/68‑1 (MQ=255)
                                                                 tCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg      <  1:2687079/67‑1 (MQ=255)
                                                                  cAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTGAAc   >  1:2478355/1‑69 (MQ=255)
                                                                   aGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTGAACa  >  1:1796866/1‑69 (MQ=255)
                                                                   |                                                                    
GTCCGCGCGCCAAAGCGCCGTTTATCGCGCTGAATATGGCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTGAACA  >  minE/2343888‑2344023

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: