Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,343,988 T→C 31.0% T218T (ACT→ACC 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,9880TC31.0% 50.0 / 26.5 58T218T (ACT→ACCglnGfused DNA‑binding response regulator in two‑component regulatory system with GlnL, nitrogen regulator I (NRI)
Reads supporting (aligned to +/- strand):  ref base T (12/28);  new base C (15/3);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 2.04e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTGAACAGGCCGATGGCGGTACAT  >  minE/2343926‑2344040
                                                              |                                                    
gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                                                <  1:1157286/69‑1 (MQ=255)
gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                                                <  1:2340768/69‑1 (MQ=255)
gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                                                <  1:2774305/69‑1 (MQ=255)
gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                                                <  1:29761/69‑1 (MQ=255)
gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                                                <  1:3273378/69‑1 (MQ=255)
gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                                                <  1:653244/69‑1 (MQ=255)
gCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGcgcg                                                <  1:3407160/69‑1 (MQ=255)
     tATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATAc                                           <  1:2514510/69‑1 (MQ=255)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:1741576/1‑69 (MQ=25)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:27543/1‑69 (MQ=25)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:1031823/1‑69 (MQ=25)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:2289714/1‑69 (MQ=25)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:2896671/1‑69 (MQ=25)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:1957089/1‑69 (MQ=25)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:2326230/1‑69 (MQ=25)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:2982820/1‑69 (MQ=25)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:3076302/1‑69 (MQ=25)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:1827942/1‑69 (MQ=25)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:1551141/1‑69 (MQ=25)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:2408691/1‑69 (MQ=25)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:2623740/1‑69 (MQ=25)
      aTCCCGAAGGATTTGATCGAATCCGAACTGTTTGGTCACGAGAAAGGCGCATTTACCGGCGCGAATAcc                                          >  1:1781850/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:2798457/1‑69 (MQ=255)
         ccAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCAtt                                       >  1:2963908/1‑69 (MQ=255)
         ccAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCAtt                                       >  1:3347594/1‑69 (MQ=255)
         ccAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCAtt                                       >  1:647299/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:3115739/68‑1 (MQ=255)
                          aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg                       <  1:1199424/68‑1 (MQ=255)
                          aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg                       <  1:1322134/68‑1 (MQ=255)
                          aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg                       <  1:631143/68‑1 (MQ=255)
                          aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg                       <  1:1766870/68‑1 (MQ=255)
                          aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg                       <  1:3314843/68‑1 (MQ=255)
                          aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg                       <  1:1823854/68‑1 (MQ=255)
                          aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg                       <  1:3211704/68‑1 (MQ=255)
                          aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg                       <  1:1843907/68‑1 (MQ=255)
                          aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg                       <  1:1869235/68‑1 (MQ=255)
                          aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg                       <  1:1940411/68‑1 (MQ=255)
                          aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg                       <  1:1967097/68‑1 (MQ=255)
                          aTCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg                       <  1:2214205/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)
                                         tcACGAGAAAGGCGCATTTACCGGCGCGAATACCAt                                        <  1:1664700/35‑1 (MQ=37)
                                             gagaAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTg                       <  1:2688780/49‑1 (MQ=255)
                                              agaAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTGAACAGGCCGATGGCGGTACat  <  1:186061/69‑1 (MQ=255)
                                              agaAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTGAACAGGCCGATGGCGGTACat  <  1:1709739/69‑1 (MQ=255)
                                               gaAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTGAACAGGCCGATGGCGGTACat  <  1:221346/68‑1 (MQ=255)
                                                             cTGGCGCGAATACCATTCGTCAGGGGCGTTTTGAACAGGCCGATGGCGGTACat  <  1:3097647/54‑1 (MQ=255)
                                                              |                                                    
GCAGCTATCCCAAAAGATTTGATCGAATCAGAACTGTTTGGCCACGAGAAAGGCGCGTTTACTGGCGCGAATACCATTCGTCAGGGGCGTTTTGAACAGGCCGATGGCGGTACAT  >  minE/2343926‑2344040

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: