Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 721,604 G→A R44C (CGT→TGT)  kdpE ← DNA‑binding transcriptional activator KdpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913721,6040GA96.7% 90.0 / NA 30R44C (CGT→TGT) kdpEDNA‑binding transcriptional activator KdpE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/14);  minor base C (1/0);  total (16/14)
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.91e-01

GCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCC  >  NC_000913/721380‑721837
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
gcggcgTAATGCGACGCTCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAAGCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGtt                                                                                                                                                                                                                 <  2:21382/251‑1 (MQ=255)
                      cGGGCCTGCAATTCGCCACTGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCCCCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGAc                                                                                                                                                                                           >  2:182464/1‑251 (MQ=255)
                                      cAATGCCAAACGGCTTACTCCGATAATCATCCGCTACGGCATCCAGCGCGGAGATTTTGTCGCTCTCTTCGCTGAGTGCGGAAAGCACAATCACCGGCACCGCGCTCAACTGGCGCAGGACGCGGATAAAATCAATCCCATCACCATCGGGCAGGCCGAGATCGAGACTAATCAAATCAGGCTTACAGGTTGCCGCTTCCAGCAAGCAGCGTTGCAGCGTTTCGGCATCAAAGACGCGCATCCCGTccacc                                                                                                                                                                           >  1:106195/1‑247 (MQ=255)
                                          gCCAAACCGCGTACTCAGATACTCAACCGCTCCGCCATCCAGCCCGGCGAATTTCCCGCCATCTTCGCTCCGTGCGGCACGCACAATCACCGGCACCACGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGCGACCGAGAATAATCAAATCTGGCTTACAGGATGCCGCTTCCAGCAAGCCGCGTTGCCGCGTTTCGGCCTCAAAGCCGCGCATCCCGTCGCCCTccc                                                                                                                                                                       >  2:41488/1‑250 (MQ=255)
                                                       cTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTgcgc                                                                                                                                                            >  1:48013/1‑249 (MQ=255)
                                                          aGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGaa                                                                                                                                                        >  1:101483/1‑250 (MQ=255)
                                                               atcatcCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCGTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGcgac                                                                                                                                                  <  2:96577/251‑1 (MQ=255)
                                                                tcatcCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGaaa                                                                                                                                                       <  1:59267/245‑1 (MQ=255)
                                                                tcatcCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGaaa                                                                                                                                                       >  2:59267/1‑245 (MQ=255)
                                                                 cagccGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGcgacga                                                                                                                                                <  2:191422/248‑1 (MQ=255)
                                                                    ccGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAAt                                                                                                                                              >  1:210905/1‑250 (MQ=255)
                                                                             cATCCAGCGCTGCGTTTTTGTCGCTCTCTTCGCTTCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGTCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGGCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCa                                                                                                                                    <  1:153750/251‑1 (MQ=255)
                                                                                     gcggcgAATTTGTCGCTCTGTTCGCTGCGTGAGGAAAGCACACTCAACGGCAACGCGCTCCACTGGCGCAGGCCGCGGATAACCTCAATCCCATAAACATCGGGCAGGCCGAGATCGAGAATAATCCAATCTGGCTTACAGGCTGCGGCTTCCAGCAAGCCGCGTTGCCGCGTTTCGGCCTCCAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAGTAGCCTGTTCATCTTcaa                                                                                                                             >  2:31690/1‑250 (MQ=255)
                                                                                       ggcgATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTcaacaa                                                                                                                          <  1:237127/251‑1 (MQ=255)
                                                                                             tttGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCCAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGa                                                                                                                     >  1:128310/1‑250 (MQ=255)
                                                                                              ttGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAAc                                                                                                                   >  1:150528/1‑251 (MQ=255)
                                                                                                        ctTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTc                                                                                                            <  1:100778/248‑1 (MQ=255)
                                                                                                        ctTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTc                                                                                                            >  2:100778/1‑248 (MQ=255)
                                                                                                          tCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGGCGCGGATAAACTCAATCCCAGCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTTCCGCTTCCAGCAAGCCGCGTTGCATCGTTTCGTCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGGAGAAAGTGACGAATAGCCTGTTCATCTGCAACAATCAGAACGTTTGTcac                                                                                                          <  1:98283/248‑1 (MQ=255)
                                                                                                           cGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGACTAGCCTGTTCATCTTCAACAATCAGCACGTTTGTCACatat                                                                                                      >  2:145257/1‑251 (MQ=255)
                                                                                                                                      gCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAAGTCTTCAAGTTCAgggg                                                                           <  1:145257/251‑1 (MQ=255)
                                                                                                                                       cACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCCCATATCCTCATGAAATTCTTAAAGTTCAGGGGc                                                                          >  2:224066/1‑251 (MQ=255)
                                                                                                                                            cgcTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCAGCGCCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCGGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGttt                                                                     <  1:142057/251‑1 (MQ=255)
                                                                                                                                                                 ggATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCACGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGATTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGa                                                 >  1:109423/1‑250 (MQ=255)
                                                                                                                                                                                      cagcGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCCGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCcggtc                           <  2:148618/248‑1 (MQ=255)
                                                                                                                                                                                      cATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCcggtc                           <  2:150528/251‑1 (MQ=255)
                                                                                                                                                                                           ggCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACTCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTGCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGTTACCACCTTCCGTTCGGTTg                      <  2:34431/251‑1 (MQ=255)
                                                                                                                                                                                             cAGGCCGAGATCGAGAATAATCAAATCTGGCTTACCGGTTTCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGa                     >  2:155148/1‑250 (MQ=255)
                                                                                                                                                                                              aGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAAc                   <  1:124233/251‑1 (MQ=255)
                                                                                                                                                                                                               aatCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGccccc  >  1:110538/1‑251 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
GCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCC  >  NC_000913/721380‑721837

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑