Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 717,981 C→T 18.6% L630L (CTG→CTA 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
*CP009273717,9810CT18.6% 31.3 / 2.9 21L630L (CTG→CTAkdpDfused sensory histidine kinase in two‑component regulatory system with KdpE: signal sensing protein
Reads supporting (aligned to +/- strand):  ref base C (8/9);  new base T (2/2);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

AATCATGCGAAAGCGCCGCCAGCAGGGCGTTGCGGATCTGTTCACGTTCGCTTGCCATCCGCGCCTGTTCTTCGCTGGCGGTTAGCGTCAGCCGCTCAAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGC  >  CP009273/717848‑718126
                                                                                                                                     |                                                                                                                                                 
aaTCATGCGAAAGCGCCGCCAGCAGGGCGTTGCGGATCTGTTCACGTTCGCTTGCCATCCGCGCCTGTTCTTCGCTGGCGGTTAGCGTCAGCCGCTCAAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTcc                                                                                                                                    >  1:190815/1‑149 (MQ=255)
  tCATGCGAAAGCGCCGCCAGCAGGGCGTTGCGGATCTGTTCACGTTCGCTTGCCATCCGCGCCTGTTCTTCGCTGGCGGTTAGCGTCAGCCGCTCAAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCgg                                                                                                                                  >  2:99121/1‑149 (MQ=255)
    agtcGAAAGGGCCGCCAGCAGGGCGTTGCGGATCTGTTCACGTTCGCTTGCCATCCGCGCCTGTTCTTCGCTGGCGGTTAGCGTCAGCCGCTCAAGTGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCGGTTCCGGGa                                                                                                                                <  2:113580/146‑1 (MQ=255)
               ccgccAGCAGGGCGTTGCGGATCTGTTCACGTTCGCTTGCCATCCGCGCCTGTTCTTCGCTGGCGGTTAGCGTCAGCCGCTCAAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTgg                                                                                                                     >  1:135111/1‑149 (MQ=255)
                  ccAGCAGGGCGTTGCGGATCTGTTCACGTTCGCTTGCCATCCGCGCCTGTTCTTCGCTGGCGGTTAGCGTCAGCCGCTCAAGGGCATTGGCGTCTAACAGTGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGgcgc                                                                                                                  <  1:88125/149‑1 (MQ=255)
                                cGGAGCTGTTCACGTTCGCTTGCCATCCGCGCCTGTTCTTCGCTGGCGGTTAGCGTCAGCCGCTCAAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCGGTTCCGGGATCATCAACTGGCGCAGGTGCCCCGGCTc                                                                                                    >  2:211430/1‑149 (MQ=255)
                                 ggATCTGTTCACGTTCGCTTGCCATCCGCGCCTGTTCTTCGCTGGCGGTTAGCGTCAGCCGCTCAAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCGACTGGCGCAGCGTCCCCGGCTcc                                                                                                   >  2:24435/1‑149 (MQ=255)
                                           aCGGTCGATTGCCATCCGCGCCGGTTTTTCGCTGGCGGTTAGCGTCAGCCGCTCAAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTTCGGGATCATCAACTGGCGCAGATTCCCCGGCTccaccaccacca                                                                                         <  1:109622/149‑1 (MQ=255)
                                             gTTCGCTTGCCATCCGCGCCTGTTCTTCGCTGGCGGTTAGCGTCAGCCGCTCAAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGc                                                                                       <  1:120566/149‑1 (MQ=255)
                                                           cgcgCCTGTTCTTCGCTGGCGGTTAGCGTCAGCCGCTCAAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTtctc                                                                         <  2:135111/149‑1 (MQ=255)
                                                                  gTTCTTCGCTGGCGGTTAGCGTCAGCCGCTCCAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCtt                                                                  >  1:245380/1‑149 (MQ=255)
                                                                      ttcGCTGGCGGTTAGCGTCAGCCGCTCAAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGTAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTa                                                              >  1:224333/1‑149 (MQ=255)
                                                                      ttcGCTGGCGGTTAGCGTCAGCCGCTCAAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGTAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTa                                                              >  1:224244/1‑149 (MQ=255)
                                                                                              cTCAAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTacac                                      <  2:190815/149‑1 (MQ=255)
                                                                                                 aaGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAACATCTGGTACGGTACACCgg                                   >  1:238632/1‑149 (MQ=255)
                                                                                                  aGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCGCCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCggg                                  <  1:28665/149‑1 (MQ=255)
                                                                                                            cGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCg                        <  1:99121/149‑1 (MQ=255)
                                                                                                               cTAACAGCGTAAACGTCTCCAGTAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGTCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTg                     <  1:320062/149‑1 (MQ=255)
                                                                                                               cTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTg                     >  1:234260/1‑149 (MQ=255)
                                                                                                                tAACAGAGGAAACATCTTCCGTCGGCGCCGCTGTCCCGGGGTCACCACCTCGCGGAGGTCCCTCGTGCCCACCACCACCGGCCCGGTGGTTTCCTCGCCGCTTTTTAGCGGCAAAAGCTGGTACGGTACGCCGGGTAACGTGTCGGTgc                    <  2:127248/149‑1 (MQ=255)
                                                                                                                                tcCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCgcaggcag    <  1:211430/149‑1 (MQ=255)
                                                                                                                                  ccgcagGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGTCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCaggcaggc  <  1:24435/147‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                 
AATCATGCGAAAGCGCCGCCAGCAGGGCGTTGCGGATCTGTTCACGTTCGCTTGCCATCCGCGCCTGTTCTTCGCTGGCGGTTAGCGTCAGCCGCTCAAGGGCATTGGCGACTAACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGC  >  CP009273/717848‑718126

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