Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,268,951 T→G 7.6% I129S (ATC→AGC)  atoC → fused response regulator of ato operon, in two‑component system with AtoS: response regulator/sigma54 interaction protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,268,9510T→G 7.6% 3.2 0.0 ‑0.048386atoCfused response regulator of ato operon, in two‑component system with AtoS: response regulator/sigma54 interaction protein
Reads supporting (aligned to +/- strand):  new base (1/3):  ref base (16/28):  total (17/31)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.66e-01

GCTTTACAACTCCAGTCAATGAAAAAAGAGATCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGC  >  REL606/2268920‑2268986
                               |                                   
gcTTTACAACTCCAGTCAATGAAAAAAGAGATCCGt                                 >  1:5134587/1‑36 (MQ=255)
gcTTTACAACTCCAGTCAATGAAAAAAGAGATCCGt                                 <  1:6858959/36‑1 (MQ=255)
gcTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGt                                 <  1:2855928/36‑1 (MQ=255)
 cTTTACAACTCCAGTCAATGAAAAAAGAGATCCGTc                                <  1:7584317/36‑1 (MQ=255)
 cTTTACAACTCCAGTCAATGAAAAAAGAGATCCGTc                                <  1:7518758/36‑1 (MQ=255)
 cTTTACAACTCCAGTCAATGAAAAAAGAGATCCGTc                                >  1:4997435/1‑36 (MQ=255)
 cTTTACAACTCCAGTCAATGAAAAAAGAGATCCGTc                                <  1:3648916/36‑1 (MQ=255)
   ttACAACTCCAGTCAATGAAAAAAGAGATCCGTCAt                              <  1:3953356/36‑1 (MQ=255)
    tACAACTACAGTCAATGAAAAAAGAGATCCGTCATc                             <  1:6563631/36‑1 (MQ=255)
     aaaaCTCCAGTCAATGAAAAAAGAGATCCGTCATCt                            <  1:3541971/34‑1 (MQ=255)
      cAACTCCAGTCAATGAAAAAAGAGATCCGTCATCTg                           <  1:6600074/36‑1 (MQ=255)
      aaaCTCCAGTCAATGAAAAAAGAGATCCGTCATCTg                           <  1:4428655/35‑1 (MQ=255)
         cTCCAGTCAATGAAAAAAGAGATCCGTCATCTGCAc                        <  1:5061944/36‑1 (MQ=255)
         cTCCAGTCAATGAAAAAAGAGATCCGTCATCTGCAc                        >  1:6672624/1‑36 (MQ=255)
         cTCCAGTCAATGAAAAAAGAGATCCGTCATCTGCAc                        <  1:1534873/36‑1 (MQ=255)
          tCCAGTCAATGAAAAAAGAGATCCGTCATCTGCAcc                       <  1:2869726/36‑1 (MQ=255)
          tCCAGTCAATGAAAAAAGAGATCCGTCATCTGCAcc                       <  1:466063/36‑1 (MQ=255)
          tCCAGTCAATGAAAAAAGAGATCCGTCATCTGCAcc                       <  1:4690422/36‑1 (MQ=255)
           caaGTAAATGAAAACAGAGATCCGTCATCTGCACCa                      <  1:5328371/34‑1 (MQ=37)
             aGTCAATGAAAAAAGAGATCCGTCATCTGCACCAgg                    >  1:4713313/1‑36 (MQ=255)
              ggcAATGAAAAAAGAGATCCGTCATCTGCTCCAGGc                   <  1:2643270/34‑1 (MQ=37)
               tCAATGAAAAAAGAGATCCGTCATCTGCACCAGGCa                  <  1:2589775/36‑1 (MQ=255)
                cAATGAAAAAAGAGATCAGTCTTCTGCACCAGTCAc                 >  1:3062371/1‑36 (MQ=25)
                 aaTGAAAAAAGAGATCCGTCATCTGCACCAGGCACt                <  1:2934154/36‑1 (MQ=255)
                 aaTGAAAAAAGAGAGCCGTCATCTGCACCAGGCACt                <  1:6855083/36‑1 (MQ=255)
                 aaTGAAAAAAGACAGCCGTCATCTGCATCATGCACt                >  1:3316442/1‑36 (MQ=16)
                  aTGAAAAAAGAGATCCGTCATCTGCACCAGGCACTg               >  1:5195711/1‑36 (MQ=255)
                   tGAAAAAAGAGATCCGTCATCTGCACCAGGCACTGa              >  1:7543260/1‑36 (MQ=255)
                   tGAAAAAAGAGATCCGTCATCTGCACCAGGCACTGa              <  1:5673504/36‑1 (MQ=255)
                    gAAAAAAGAGATCCGTCATCTGCACCAGGCCCTGAg             <  1:60216/36‑1 (MQ=255)
                    gAAAAAAGAGATCCGTCATCTGCACCAGGCACTGAg             >  1:5885336/1‑36 (MQ=255)
                    gAAAAAAGAGATCCGTCATCTGCACCAGCCACTGAg             >  1:4732647/1‑36 (MQ=255)
                      aaaaaGAGATCCGTCATCTGCACCAGGCACTGAGCa           >  1:1679464/1‑36 (MQ=255)
                       aaaaGAGAGCCGTCATCTGCACCAGGCACTGAGCAc          <  1:2018964/36‑1 (MQ=255)
                        aaaGAGATCCGTCATCTGCACCAGGCACTGAGCAcc         <  1:844466/36‑1 (MQ=255)
                        aaaGAGATCCGTCATCTGCACCAGGCACTGAGCAcc         >  1:4116323/1‑36 (MQ=255)
                        aaaGAGATCCGTCATCTGAACAAGGCACTGAGCaac         >  1:4200895/1‑34 (MQ=21)
                         aaGAGATCCGTCATCTGCACCAGGCACTGAGCACCa        >  1:7565804/1‑36 (MQ=255)
                          agagATCCGTCATCTGCACCAGGCACTGAGCACCAg       <  1:5429319/36‑1 (MQ=255)
                          agagATCCGTCATCTGCACCAGGCACTGAGCACCAg       <  1:7071844/36‑1 (MQ=255)
                           gagaTCCGTCATCTGCACCAGGCACTGAGCACCAGc      <  1:2001850/36‑1 (MQ=255)
                           gagaTCCGTCATCTGCACCAGGCACTGAGCACCAGc      <  1:2082383/36‑1 (MQ=255)
                            agaTCCGTCATCTGCACCAGGCACTGAGCACCAGCt     >  1:4035467/1‑36 (MQ=255)
                            agaTCCGTCATCTGCACCAGGCACTGAGCACCAGCt     >  1:5540622/1‑36 (MQ=255)
                              aTCCGTCATCTGCACCAGGCACTGAGCACCAGCTgg   <  1:7223481/36‑1 (MQ=255)
                              aTCCGTCATCTGCACCAGGCACTGAGCACCAGCTgg   >  1:2842712/1‑36 (MQ=255)
                              aTCCGTCATCTGCACCAGGCACTGAGCACCAGCTgg   <  1:5590034/36‑1 (MQ=255)
                              aTCCGTCATCTGCACCAGGCACTGAGCACCAGCTgg   <  1:784062/36‑1 (MQ=255)
                               tCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGc  >  1:2280541/1‑36 (MQ=255)
                               |                                   
GCTTTACAACTCCAGTCAATGAAAAAAGAGATCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGC  >  REL606/2268920‑2268986

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG