Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 723,528 (T)8→9 coding (887/2685 nt) kdpD ← fused sensory histidine kinase in two‑component regulatory system with KdpE: signal sensing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913723,5201.T97.4% 117.3 / ‑3.4 38coding (895/2685 nt)kdpDfused sensory histidine kinase in two‑component regulatory system with KdpE: signal sensing protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (17/20);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

CCGAGAATAATCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCA  >  NC_000913/723387‑723649
                                                                                                                                      |                                                                                                                                  
ccGAGAATAATCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑ttttt                                                                                                                              >  2:451870/1‑139 (MQ=255)
            ttGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAg                                                                                                                   <  1:352794/139‑1 (MQ=255)
            ttGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAg                                                                                                                   <  2:628926/139‑1 (MQ=255)
             tGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGc                                                                                                                  >  1:508420/1‑139 (MQ=255)
               ccGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCgg                                                                                                                >  1:149408/1‑139 (MQ=255)
                 gagaTTATGTTTACGGGCATAACGCACGACCGCTTTCTCTTCCGCTGTATCAGAAAGTGTTGCGGGCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATGGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTg                                                                                                              <  2:613926/139‑1 (MQ=255)
                 gagaTTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTg                                                                                                              <  2:477524/139‑1 (MQ=255)
                  agaTTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTTCCGGTAAGCGGTg                                                                                                              <  2:576024/139‑1 (MQ=255)
                            cACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                                   <  1:478024/139‑1 (MQ=255)
                            cACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                                   <  1:447643/139‑1 (MQ=255)
                            cACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                                   <  1:752491/139‑1 (MQ=255)
                             aCGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGtt                                                                                                 <  2:296426/139‑1 (MQ=255)
                              cGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGtt                                                                                                 >  2:503431/1‑139 (MQ=255)
                              cGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGCTCCTGCGCCAGTCGTAAGGCGCCGAGTATTGCCCGCCGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGtt                                                                                                 >  2:921174/1‑139 (MQ=255)
                                  cATAACGCACTACCGCTTTCTCTTCCACTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCaa                                                                                             >  1:881581/1‑139 (MQ=255)
                                    tAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGCGTTTCAACa                                                                                           >  1:675023/1‑139 (MQ=255)
                                                 cTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGcc                                                                              <  2:932245/139‑1 (MQ=255)
                                                      tctTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACg                                                                         >  1:77788/1‑139 (MQ=255)
                                                                           gtTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTTCGGTAAGCGGTgcagggca                                                                                                       >  2:577320/1‑88 (MQ=255)
                                                                           gtTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTgcagggca                                                                                                       <  1:577320/88‑1 (MQ=255)
                                                                            ttGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAg                                                                                                                   <  2:570740/75‑1 (MQ=255)
                                                                            ttGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAg                                                                                                                   >  1:570740/1‑75 (MQ=255)
                                                                            ttGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTgcagggcag                                                                                                      >  2:951519/1‑88 (MQ=255)
                                                                            ttGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTgcagggcag                                                                                                      <  1:951519/88‑1 (MQ=255)
                                                                            ttGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTgcagggca                                                                                                       >  2:522435/1‑87 (MQ=255)
                                                                            ttGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTgcagggca                                                                                                       <  1:522435/87‑1 (MQ=255)
                                                                              gCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAgccgc                                                 >  1:512709/1‑139 (MQ=255)
                                                                               ccGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAgccgcg                                                >  2:723365/1‑139 (MQ=255)
                                                                                      cgcgCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTgc                                         >  2:48895/1‑139 (MQ=255)
                                                                                                      gcCAGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCg                         <  2:512709/139‑1 (MQ=255)
                                                                                                         aGACGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTg                      >  2:696755/1‑139 (MQ=255)
                                                                                                           aCGTAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGcc                    <  2:508420/139‑1 (MQ=255)
                                                                                                              tAAGGCGCTGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGgt                 <  1:503431/139‑1 (MQ=255)
                                                                                                                      tGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGAc                                   <  2:386743/113‑1 (MQ=255)
                                                                                                                      tGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGAc                                   >  1:386743/1‑113 (MQ=255)
                                                                                                                      tGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGc         <  1:723365/139‑1 (MQ=255)
                                                                                                                      tGAGAATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGc         <  1:696755/139‑1 (MQ=255)
                                                                                                                        agaATTGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCg       <  1:451870/139‑1 (MQ=255)
                                                                                                                             tGCCCGACGT‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCa  >  2:1052356/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                  
CCGAGAATAATCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCA  >  NC_000913/723387‑723649

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑