Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 703,479 C→G 31.3% R117P (CGC→CCC)  kdpE ← DNA‑binding response regulator in two‑component regulatory system with KdpD

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606703,4790C→G31.3% 6.0 ‑3.7 0.027350kdpEDNA‑binding response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  new base (9/0):  ref base (10/8):  total (19/8)
Fisher's exact test strand distribution p-value = 2.58e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AGCGGATCGGGCGCGGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCA  >  REL606/703446‑703514
                                 |                                   
aGCGGATCGGGCGCGGTGGTGGCAGAGTGGGGGGGt                                   >  1:6717386/1‑36 (MQ=37)
aGCGGATCGGGCGCGGTGGTGGCAGAGTGGCGGGGt                                   >  1:5490724/1‑36 (MQ=255)
    gATCGGGCGCGGTGGTGGCAGAGTGGCGGGGTAATg                               >  1:5836746/1‑36 (MQ=255)
    gATCGGGCGCGGTGGTGGCAGAGTGGCGGGGCAATg                               >  1:5260867/1‑36 (MQ=37)
     aTCGGGCGCGGTGGTGGCAGAGTGGCGGGGCAATgg                              >  1:1132515/1‑35 (MQ=37)
     aTCGGGCGCGGTGGTGGCAGAGTGGCGGCGGAATgg                              >  1:6472907/1‑35 (MQ=255)
      tCGGGCGCGGTGGTGGCAGAGTGGCGGGGTAAttgg                             >  1:6930682/1‑33 (MQ=37)
       cGGGCGCGGTGGTGGCAGAGTGGCGGCGTAATGCGa                            <  1:1730379/36‑1 (MQ=255)
        gggCGCGGTGGTGGCAGAGTGGCGGGGTTATGGGAc                           >  1:5779941/1‑36 (MQ=25)
        gggCGCGGTGGTGGCAGAGTGGCGGGGTAATGCggg                           >  1:4864720/1‑34 (MQ=38)
        gggCGCGGTGGTGGCAGAGTGGCGGCGCGATGTGAc                           >  1:4869788/1‑36 (MQ=25)
          gcgcGGTGGTGGCAGAGTGGCGGCGTAATGCGAcgc                         <  1:2887460/36‑1 (MQ=255)
          gcgcGGGGGTGGCAGAGTGGAGGGGTATTGCGaaga                         >  1:4870581/1‑33 (MQ=16)
           cgcgGTGGTGGCAGAGTGGCGGCGTAATGTGAcggg                        >  1:4479218/1‑34 (MQ=38)
              ggtggtGGCAGAGTGGCGGCGTAATGCGACgcgccg                     >  1:3130422/1‑34 (MQ=255)
                  gtggCAGAGTGGCGGCGTAATGCGACGCGCAGACgg                 <  1:2602743/36‑1 (MQ=255)
                   tggCAGAGTGGCGGCGTAATGCGACGCGCAGACggg                <  1:2768831/36‑1 (MQ=255)
                    ggCAGAGTGGCGGCGTAATTCGACGCGCAGACGGGc               >  1:1780081/1‑36 (MQ=255)
                    ggCAGAGTGGCGGCGTAATGCGACGCGGAGACgggg               >  1:2341882/1‑35 (MQ=255)
                    ggCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGc               <  1:6178926/36‑1 (MQ=255)
                     gCAGAGTGGCGGCGTAATGCGACGCGCAGACggggc              >  1:3367838/1‑34 (MQ=255)
                      cAGAGTGGCGGCGTAATGCGACGCGCAGACGGGGCt             >  1:5955187/1‑36 (MQ=255)
                       agagTGGCGGCGTAATGCGACGCGCAGACGGGCCTg            <  1:5385213/36‑1 (MQ=255)
                         agTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCa          <  1:3726028/36‑1 (MQ=255)
                           tGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAAt        >  1:7847894/1‑36 (MQ=255)
                           tGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAAt        <  1:5834108/36‑1 (MQ=255)
                           tGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAAt        >  1:4286666/1‑36 (MQ=255)
                            ggcggcGTAATGCGACGCGCAGACGGGCCTGCAAtt       <  1:6137095/36‑1 (MQ=255)
                            agcggcgTAATGCGACGCGCAGACGGGCCTGCAAtt       <  1:4676309/35‑1 (MQ=255)
                             gcggcgTAATGCGACGCGCAGACGGGCCTGCAATTc      >  1:2098670/1‑36 (MQ=255)
                              cggcgTAATGCGACGCGCAGACGGGCCTGCAATTCg     >  1:6113788/1‑36 (MQ=255)
                               ggcgTAATGCGACGCGCAGACGGGCCTGCAATTCGc    >  1:3857900/1‑36 (MQ=255)
                               ggcgTAATGCGACGCGCAGACGGGCCTGCAATTCGc    <  1:3940595/36‑1 (MQ=255)
                               ggcgTAATGCGACGCGCAGACGGGCCTGCAATTCGc    >  1:1160132/1‑36 (MQ=255)
                               ggcgTAATGCGACGCGCAGACGGGCCTGCAATTCGc    >  1:7285448/1‑36 (MQ=255)
                                gcgTCATGCGACGCGCAGACGGGCCTGCAATTCGcc   <  1:263514/36‑1 (MQ=255)
                                gcgTAATGCGACGCGCAGACGGGCCTGCAATTCGcc   >  1:2407503/1‑36 (MQ=255)
                                gcgTAATGCGACGCGCAGACGGGCCTGCAATTCGcc   >  1:6582022/1‑36 (MQ=255)
                                 cgTAATGCGACGCGCAGACGGGCCTGCAATTCGCCa  >  1:6912106/1‑36 (MQ=255)
                                 |                                   
AGCGGATCGGGCGCGGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCA  >  REL606/703446‑703514

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG