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.T93.3% 86.4 / ‑0.2 30coding (895/2685 nt)kdpDfused sensory histidine kinase in two‑component regulatory system with KdpE: signal sensing protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (13/15);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CGAGAATAATCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCG  >  NC_000913/723388‑723644
                                                                                                                                     |                                                                                                                             
cGAGAATAATCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑tttttt                                                                                                                        >  1:462224/1‑139 (MQ=255)
                               ggCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTc                                                                                          <  2:462224/139‑1 (MQ=255)
                                 cATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTTCCGGTAAGCGGTgcagggca                                                                                                  <  1:1178216/130‑1 (MQ=255)
                                 cATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTTCCGGTAAGCGGTgcagggca                                                                                                  >  2:1178216/1‑130 (MQ=255)
                                   tAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACa                                                                                      >  2:728113/1‑139 (MQ=255)
                                                cTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGcc                                                                         >  1:610732/1‑139 (MQ=255)
                                                cTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑tt                                                                                                                            >  1:93245/1‑87 (MQ=255)
                                                cTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑tt                                                                                                                            <  2:93245/87‑1 (MQ=255)
                                                     tctTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACg                                                                    <  1:567418/139‑1 (MQ=255)
                                                           gCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTAccc                                                              <  2:1359464/139‑1 (MQ=255)
                                                               gATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGtt                                                                                            >  2:945875/1‑105 (MQ=255)
                                                               gATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGtt                                                                                            <  1:945875/105‑1 (MQ=255)
                                                                          gtTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATa                                                                                    >  2:401133/1‑102 (MQ=255)
                                                                          gtTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATa                                                                                    <  1:401133/102‑1 (MQ=255)
                                                                                   tcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAg                                                                                                              <  2:1011196/67‑1 (MQ=255)
                                                                                   tcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAg                                                                                                              >  1:1011196/1‑67 (MQ=255)
                                                                                   tcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCt                                      >  2:603659/1‑139 (MQ=255)
                                                                                   tcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCATCATACACCGCGTGCCAGACg                                                                    <  1:515712/108‑1 (MQ=255)
                                                                                   tcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCATCATACACCGCGTGCCAGACg                                                                    >  2:515712/1‑108 (MQ=255)
                                                                                     cgcgCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACACCAGCCGCGCCGCTgc                                    <  1:323902/139‑1 (MQ=255)
                                                                                                      ccAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGcc                                                                         <  1:490244/85‑1 (MQ=255)
                                                                                                      ccAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGcc                                                                         >  2:490244/1‑85 (MQ=255)
                                                                                                        aGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTgcagggca                                                                                                  <  1:1135294/58‑1 (MQ=255)
                                                                                                        aGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTgcagggca                                                                                                  >  2:1135294/1‑58 (MQ=255)
                                                                                                        aGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAgccg                                             <  1:655198/111‑1 (MQ=255)
                                                                                                        aGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAgccg                                             >  2:655198/1‑111 (MQ=255)
                                                                                                                    cTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTg                                                                                                         >  2:518766/1‑39 (MQ=37)
                                                                                                                    cTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTg                                                                                                         <  1:518766/39‑1 (MQ=37)
                                                                                                                    cTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                              >  2:292104/1‑50 (MQ=255)
                                                                                                                    cTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                              <  1:292104/50‑1 (MQ=255)
                                                                                                                     tGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGc    >  2:544004/1‑139 (MQ=255)
                                                                                                                     tGAGAATGGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGc    <  1:961883/139‑1 (MQ=255)
                                                                                                                       agaATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCg  <  1:603659/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                             
CGAGAATAATCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCG  >  NC_000913/723388‑723644

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