Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 726,253 A→G 13.6% 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,2530AG13.6% 107.8 / 10.1 44R402R (CGT→CGCkdpBpotassium translocating ATPase, subunit B
Reads supporting (aligned to +/- strand):  ref base A (21/17);  new base G (3/3);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGC  >  W3110S.gb/726184‑726322
                                                                     |                                                                     
aCGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGcgac                                                                      <  1:2677725/71‑1 (MQ=255)
aCGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGcgac                                                                      <  1:2348342/71‑1 (MQ=255)
 cgcgCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGcgacg                                                                     <  1:3277577/71‑1 (MQ=255)
  gcgcAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGcgacg                                                                     >  1:1415457/1‑70 (MQ=255)
          tGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATc                                                            >  1:3737354/1‑71 (MQ=255)
          tGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATc                                                            >  1:3414943/1‑71 (MQ=255)
              cGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGGCGAATGGCATCGa                                                          >  1:1052319/1‑69 (MQ=255)
              cGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGGCGAATGGCATCGa                                                          >  1:2093296/1‑69 (MQ=255)
              cGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGaa                                                                   >  1:2336142/1‑60 (MQ=255)
               gACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGa                                                          >  1:1637829/1‑68 (MQ=255)
                 cTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGcga                                                                       >  1:2237713/1‑53 (MQ=255)
                      tGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCtt                                                 <  1:1138293/70‑1 (MQ=255)
                         tCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACg                                             <  1:999565/71‑1 (MQ=255)
                         tCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACg                                             <  1:3989250/71‑1 (MQ=255)
                           aaCATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATgg                                                                >  1:1996515/1‑50 (MQ=255)
                           aaCATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATgg                                                                >  1:3626864/1‑50 (MQ=255)
                           aaCATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATgg                                                                >  1:1307117/1‑50 (MQ=255)
                           aaCATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATgg                                                                >  1:3021232/1‑50 (MQ=255)
                           aaCATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATgg                                                                >  1:2898862/1‑50 (MQ=255)
                           aaCATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGa                                           <  1:346873/71‑1 (MQ=255)
                                 ggTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGa                                                          >  1:2911739/1‑50 (MQ=255)
                                  gTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCg                                    <  1:3287304/71‑1 (MQ=255)
                                   tAGGGAAGTGACCACCGTTAGCCTCAACATGGcgacg                                                                     <  1:1762673/37‑1 (MQ=255)
                                        aaGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTc                              >  1:3190526/1‑71 (MQ=255)
                                          gTGACCACCGTTAGCCTCAACATGTCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGa                            <  1:2741017/71‑1 (MQ=255)
                                          gTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATg                                      <  1:3576491/61‑1 (MQ=255)
                                          gTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATg                                      <  1:1323233/61‑1 (MQ=255)
                                          gTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGa                            <  1:613558/71‑1 (MQ=255)
                                              ccaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGtt                        <  1:2486213/71‑1 (MQ=255)
                                              ccaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGtt                        <  1:277353/71‑1 (MQ=255)
                                                accGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCtt                                                 <  1:1732334/44‑1 (MQ=255)
                                                 ccGTTAGCCTCAACATGGCGGCGAATGGCATCGTCAGAACCTTTACGGATCATGCGGTTGTc                              <  1:4081203/62‑1 (MQ=255)
                                                 ccGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTc                              <  1:565867/62‑1 (MQ=255)
                                                  cGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTa                                               >  1:2366258/1‑44 (MQ=255)
                                                     tAGCCTCAACATGGCGACGAATGGCATCGACAGAAcc                                                   >  1:209904/1‑37 (MQ=255)
                                                     tAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCg                 >  1:2638691/1‑71 (MQ=255)
                                                      aGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCa                                        >  1:2456932/1‑47 (MQ=255)
                                                         cTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCTGTTGTCGATGTTGATCCCGCTCa             <  1:1034346/71‑1 (MQ=255)
                                                         cTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCa             <  1:3602889/71‑1 (MQ=255)
                                                         cTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCa             <  1:3376828/71‑1 (MQ=255)
                                                            aaCATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGa                            >  1:2945352/1‑53 (MQ=255)
                                                            aaCATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATc                    >  1:3932642/1‑61 (MQ=255)
                                                            aaCATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATc                    >  1:318895/1‑61 (MQ=255)
                                                            aaCATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATc                    >  1:60240/1‑61 (MQ=255)
                                                            aaCATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATc                    >  1:1856830/1‑61 (MQ=255)
                                                            aaCATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATc           <  1:1819527/70‑1 (MQ=255)
                                                              cATGGCGACGAATGGCATCGACAGAACCTTTACGGATCAt                                       <  1:1593566/40‑1 (MQ=255)
                                                               aTGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGc       <  1:391128/71‑1 (MQ=255)
                                                                 ggCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCtt     >  1:3289703/1‑71 (MQ=255)
                                                                  gCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCttt    >  1:524207/1‑71 (MQ=255)
                                                                   cgacgaATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTg   <  1:3438763/71‑1 (MQ=255)
                                                                    gacgaATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTgc  >  1:4004379/1‑71 (MQ=255)
                                                                     |                                                                     
ACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGC  >  W3110S.gb/726184‑726322

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: