Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 725,279 T→C 10.5% K586K (AAA→AAG kdpB ← K(+) transporting P‑type ATPase subunit KdpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913725,2790TC10.5% 103.3 / 4.0 47K586K (AAA→AAGkdpBK(+) transporting P‑type ATPase subunit KdpB
Reads supporting (aligned to +/- strand):  ref base T (14/25);  new base C (3/2);  total (17/30)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

GACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGC  >  NC_000913/725150‑725416
                                                                                                                                 |                                                                                                                                         
gACGATAGTCGAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGt                                                                                                                              <  1:784813/143‑1 (MQ=255)
   gagagtCAGGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCAGTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTgg                                                                                                                           <  2:180231/138‑1 (MQ=255)
   gATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGTAGTATTTCGCCACATCggcgg                                                                                                                           >  1:107765/1‑139 (MQ=255)
     tAATCGTGGCGATGAAGAGCACCGCACTGAGGATTGCGGAGTCGGGCAAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCa                                                                                                                         <  1:118787/143‑1 (MQ=255)
       ttCAGGGCGTGGAAAATCACCGCGCTCAGGATTGTGGACAGCAGGTTGTGCGTGGACATGATGTTCAGCGCATTCAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAAt                                                                                                                       <  2:786915/142‑1 (MQ=255)
       aTCAGGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATACAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAAt                                                                                                                       <  2:615795/143‑1 (MQ=255)
               atagAAAATCACAGCACTGAGGAATGCAGAGACGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGCTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAgg                                                                                                               <  1:16452/140‑1 (MQ=255)
                  gAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCCTTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTCTTTCGCCCCATCGTCGCCAATGCGGAAggtgg                                                                                                            >  2:101026/1‑143 (MQ=255)
                     gaTCACCGCACTGAGGAGTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCa                                                                                                         <  2:171423/142‑1 (MQ=255)
                       tCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATGTCGCCACATCGGTGGCAATGGTGACGGTGGTCAGc                                                                                                       >  1:428647/1‑143 (MQ=255)
                         cgcgtaCTGAGGACTGCGGAGTCGGGCGAATGCAGCCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGGCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGa                                                                                                     <  2:71293/138‑1 (MQ=255)
                                    gATTGCGGAGGCGGGCAAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCGGAAGGTGGTCAGCGAGCCACGGGtca                                                                                          <  1:878227/143‑1 (MQ=255)
                                      ttGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGtcatc                                                                                        <  1:189717/143‑1 (MQ=255)
                                      gggCGGAGTCGGGCGAATGCGGGCACATGATGTTCAGCGCATTTAACTGCGGATCCCAAAACTTTGATGTCGCTGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGtcatc                                                                                        <  1:909022/141‑1 (MQ=255)
                                         cgaaggcGGGCGAATGCAGGCACATGACGGTCAGCGCATGTAACTGCGGGCACGGTGCCGCGAATGCGGCCGGAAGAATGGCGCAGTATGGCGCCACAGCGCTGGCAAAGCTGACTGCGGTCAGCAAGCCAAGGGCCCTCAGc                                                                                     >  2:758319/7‑143 (MQ=255)
                                             gacGGGCGAAGGCGGGCACAGGAGGGTCAGCGCAGGTAAATGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCt                                                                                 <  1:255680/141‑1 (MQ=255)
                                                 ggCGAATGCAGGCACATGATGTTCGGCGCATTTAACTGCGGATATGTTGCCGCGAATGTCGCCGGAATACTGGCGACGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCCAGCCACGGGTCATCAGCATCTGttt                                                                             <  2:48258/143‑1 (MQ=255)
                                                 ggCGAATGCAGGCACATGATGTGCAGCGCATTTAACTGCGGATACGTTGCCTCGAATTCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGttt                                                                             <  2:359733/143‑1 (MQ=255)
                                                   ggAATGCATGCACATGATGTCCAACGCATTAAACTGCGAATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGc                                                                           <  2:915269/142‑1 (MQ=255)
                                                      aTGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCac                                                                        >  2:858092/1‑142 (MQ=255)
                                                         cAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAAtgt                                                                     <  1:311375/143‑1 (MQ=255)
                                                           gggacaTGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGGATTGCGCCACGTGGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGc                                                                   <  2:93951/140‑1 (MQ=255)
                                                                cgcatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTCGTCAGCGAGCCACGGGTCATCAGCATCTGCTTGCCAATGTGcaccac                                                              <  1:427516/140‑1 (MQ=255)
                                                                 ggatgTTAAGCGCATTGAACTGCGGATACGTGGCCGCCAATGCCGCCGGAACAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCaccacc                                                             <  1:310142/142‑1 (MQ=255)
                                                                  gatgTTCGGCGCAGTTAAGTGCGGATACGGTGCGGCTAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTGAGCGAGCCAGGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCt                                                            <  1:130768/143‑1 (MQ=255)
                                                                  gatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATACTGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCCTCTGTTTGCCAATGGGCACCAGCt                                                            >  2:915262/1‑143 (MQ=255)
                                                                  gatgTTCAGCGCATTTAACTGCGGATACGTTGCAGCGCATGCCACCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCCCGGGTCAGCAGCAGCTGGTGGCCAATGTGCACCACCt                                                            >  2:745303/1‑143 (MQ=255)
                                                                    tgTTCAGCGTATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGGGAGGGAGCAGGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCg                                                          <  1:909214/143‑1 (MQ=255)
                                                                      ttCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAATACTTCGCcatt                                                                                                                                  >  1:160479/1‑67 (MQ=255)
                                                                      ttCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAATACTTCGCcatt                                                                                                                                  <  2:160479/69‑3 (MQ=38)
                                                                                    cTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGt                                          >  1:89341/1‑143 (MQ=255)
                                                                                               ttGCCGCTAATGTCTCCGGAGTCGTGGCGATGTTCTTCTCGGCACCGTTGGCAATGCTGAAGGTGGCCAGCGATGCTCCGGTCACTAGCTTCTGTGTGCCAATGTGCACCACCTCGATCCACTTGGTCGGGTTAGAGTCgaga                               <  2:796179/143‑1 (MQ=255)
                                                                                                tGCCGCGACTGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGAt                              <  1:27958/143‑1 (MQ=255)
                                                                                                tGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGAt                              <  2:834012/143‑1 (MQ=255)
                                                                                                  ccGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCg                            >  1:130810/1‑143 (MQ=255)
                                                                                                   cgcgAAGGCCTCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGa                           <  2:426572/143‑1 (MQ=255)
                                                                                                    gcgAATGCCGCCGGAATCATGGCGAAGTAGTTCGCCACATCGTTCGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGAc                          <  2:628437/143‑1 (MQ=255)
                                                                                                            cgccgGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGc                  >  1:89465/1‑143 (MQ=255)
                                                                                                             cccttattaatGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTAAGCGAGCCACGGGGCATCAGCATCGGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGcc                 <  1:24907/136‑1 (MQ=255)
                                                                                                               ctgAATAATGGCGAATGAGCTCGCCACATCGTTGGCCTTGCAATATGTGCTCGTCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCg               <  2:399859/141‑1 (MQ=255)
                                                                                                               cgGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACATGGTCGGGTTAGAGTCGAGGACGACGATAATGccca               >  2:2835/1‑142 (MQ=255)
                                                                                                                    atgtGGCGAAATACTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCGGTTTGCCAATGTGCACCGCCTCGAGCGACTTGGTCGGGTTAGAGTCGAGCTCGGCGATATTGCCCGCctct          >  2:692046/4‑143 (MQ=255)
                                                                                                                      atGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTACCGTCGAGATCGACCATATTGCCCGCCTCttt        >  2:638392/1‑143 (MQ=255)
                                                                                                                       tGGCGAAATACTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGCCCCTATTGCCCGCCTCTTTc       >  1:772433/1‑143 (MQ=255)
                                                                                                                          cGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTcgcc    <  1:567218/143‑1 (MQ=255)
                                                                                                                            aaGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTcgccgc  >  1:222647/1‑143 (MQ=255)
                                                                                                                            aaGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTcgccgc  >  1:691772/1‑143 (MQ=255)
                                                                                                                                 |                                                                                                                                         
GACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGC  >  NC_000913/725150‑725416

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: