Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 726,253 A→G 10.1% R402R (CGT→CGC kdpB ← potassium translocating ATPase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb726,2530AG10.1% 182.3 / 11.0 69R402R (CGT→CGCkdpBpotassium translocating ATPase, subunit B
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base G (4/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCG  >  W3110S.gb/726185‑726323
                                                                    |                                                                      
cgcgCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAGCATGgcggcg                                                                      <  1:166174/71‑1 (MQ=255)
 gcgcAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGcgacga                                                                     <  1:1529027/71‑1 (MQ=255)
 gcgcAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGcgacg                                                                      >  1:994528/1‑70 (MQ=255)
  cgcAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGcgacg                                                                      <  1:3182338/69‑1 (MQ=255)
       ccTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGaa                                                                    >  1:3668138/1‑66 (MQ=255)
       ccTGATAGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCa                                                               >  1:908012/1‑71 (MQ=255)
          gATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATc                                                             >  1:2389451/1‑70 (MQ=255)
           aTCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATc                                                             <  1:1424627/69‑1 (MQ=255)
           aTCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATc                                                             <  1:782825/69‑1 (MQ=255)
                cTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGaa                                                      >  1:2688336/1‑71 (MQ=255)
                  ttttGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATc                                                             <  1:1003958/62‑1 (MQ=255)
                   tttGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCt                                                   >  1:29176/1‑71 (MQ=255)
                   tttGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCa                                                               >  1:1521776/1‑59 (MQ=255)
                   tttGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCt                                                   >  1:4028015/1‑71 (MQ=255)
                     tGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCtt                                                  <  1:763258/70‑1 (MQ=255)
                       aTCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGa                                                       <  1:4022131/63‑1 (MQ=255)
                          aaCATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGa                                                       <  1:1547868/60‑1 (MQ=255)
                          aaCATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGa                                                       <  1:3669580/60‑1 (MQ=255)
                          aaCATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGa                                                       <  1:4585924/60‑1 (MQ=255)
                           aCATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCtt                                                  <  1:87828/64‑1 (MQ=255)
                            cATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAAc                                                     >  1:3116116/1‑60 (MQ=255)
                            cATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAAc                                                     >  1:3767769/1‑60 (MQ=255)
                              tCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAAt                                                                   >  1:2671635/1‑44 (MQ=255)
                                 gTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCg                                     <  1:2931149/71‑1 (MQ=255)
                                  tAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCgg                                    <  1:1270417/71‑1 (MQ=255)
                                    gggAAGTGACCACCGTTAGCCTCAACATGGcgacg                                                                      <  1:2004516/35‑1 (MQ=255)
                                    gggAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGtt                                  >  1:3634674/1‑71 (MQ=255)
                                    gggAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGt                                   >  1:2467896/1‑70 (MQ=255)
                                      gAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGAc                                                          >  1:413722/1‑45 (MQ=255)
                                      gAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACa                                                         >  1:2123634/1‑46 (MQ=255)
                                      gAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGt                                >  1:2456960/1‑71 (MQ=255)
                                       aaGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTc                               <  1:2319998/71‑1 (MQ=255)
                                       aaGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTc                               >  1:1336518/1‑71 (MQ=255)
                                       aaGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTc                               <  1:1906890/71‑1 (MQ=255)
                                         gTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACgg                                             <  1:160739/55‑1 (MQ=255)
                                         gTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGa                             <  1:3544102/71‑1 (MQ=255)
                                          tGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGa                                                           >  1:4475492/1‑40 (MQ=255)
                                          tGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCg                                     <  1:1950382/62‑1 (MQ=255)
                                           gACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATg                           >  1:3879578/1‑71 (MQ=255)
                                            accaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATg                           <  1:3352103/70‑1 (MQ=255)
                                            accaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGt                          >  1:2336300/1‑71 (MQ=255)
                                             ccaccGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGa                             >  1:1976330/1‑67 (MQ=255)
                                             ccaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGtt                         <  1:4639163/71‑1 (MQ=255)
                                             ccaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGtt                         <  1:3574525/71‑1 (MQ=255)
                                              caccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCttt                                                 <  1:3104329/46‑1 (MQ=255)
                                              caccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACg                                              <  1:1339423/49‑1 (MQ=255)
                                              caccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTg                        <  1:2919420/71‑1 (MQ=255)
                                               accGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGa                       <  1:1148689/71‑1 (MQ=255)
                                               accGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGa                       <  1:627434/71‑1 (MQ=255)
                                               accGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGa                       <  1:3168383/71‑1 (MQ=255)
                                               accGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGa                       <  1:3530727/71‑1 (MQ=255)
                                               accGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGa                       <  1:3476653/71‑1 (MQ=255)
                                                 cGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCtt                                                  <  1:3187239/42‑1 (MQ=255)
                                                   ttAGCCTCAACATGGCGGCGAATGGCATCGACAGAAc                                                     >  1:404082/1‑37 (MQ=255)
                                                   ttAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATg                                       <  1:3310252/51‑1 (MQ=255)
                                                   ttAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGtt                         >  1:845094/1‑65 (MQ=255)
                                                     aGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGa                             >  1:2441072/1‑59 (MQ=255)
                                                        cTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATc                     >  1:239139/1‑64 (MQ=255)
                                                        cTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCACGCTCa              >  1:3013104/1‑71 (MQ=255)
                                                         tCAACATGGCGACGAATGGCATCGACAGAACCTTTa                                                >  1:293129/1‑36 (MQ=255)
                                                         tCAACATGGCGACGAATGGCATCGACAGAACCTTTa                                                >  1:4292461/1‑36 (MQ=255)
                                                         tCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCg                              <  1:4462608/54‑1 (MQ=255)
                                                         tCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGa                       <  1:3215322/61‑1 (MQ=255)
                                                         tCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGa                       <  1:2508572/61‑1 (MQ=255)
                                                          cAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATc                                          >  1:2798908/1‑41 (MQ=255)
                                                          cAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATc            >  1:841839/1‑71 (MQ=255)
                                                           aaCATGGCGACGAATGGCATCGACAGAACCTTTACGGa                                            >  1:1298096/1‑38 (MQ=255)
                                                           aaCATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGAt                            >  1:948725/1‑54 (MQ=255)
                                                           aaCATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGAt                            >  1:3436475/1‑54 (MQ=255)
                                                           aaCATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCACGCTCATcc           >  1:2849940/1‑71 (MQ=255)
                                                             cATGGCGACGAATGGCATCGACAGAACCTTTACgg                                             <  1:3015406/35‑1 (MQ=255)
                                                              aTGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCgg                                    <  1:1679382/43‑1 (MQ=255)
                                                               tGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCt       <  1:2644566/71‑1 (MQ=255)
                                                               tGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCt       <  1:4042206/71‑1 (MQ=255)
                                                                  cgacgaATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTg    <  1:2886631/71‑1 (MQ=255)
                                                                    acgaATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTgcg  <  1:716862/71‑1 (MQ=255)
                                                                    |                                                                      
CGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCG  >  W3110S.gb/726185‑726323

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: