Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 723,528 (T)8→7 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,5210T.93.8% 117.7 / 3.6 32coding (894/2685 nt)kdpDfused sensory histidine kinase in two‑component regulatory system with KdpE: signal sensing protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (12/18);  total (13/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

TCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGA  >  NC_000913/723397‑723645
                                                                                                                            |                                                                                                                            
tCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAg                                                                                                               <  2:456929/139‑1 (MQ=255)
     ccGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTg                                                                                                          >  2:86881/1‑139 (MQ=255)
            tATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTgcagggca                                                                                                   <  1:493013/139‑1 (MQ=255)
             aTGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTgcagggcag                                                                                                  <  1:58416/139‑1 (MQ=255)
              tGTTCACGGGC‑‑AACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                               <  1:454151/139‑1 (MQ=255)
                ttCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                               <  1:86881/139‑1 (MQ=255)
                  cACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGtt                                                                                             <  1:149110/139‑1 (MQ=255)
                   aCGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGttt                                                                                            <  2:96206/139‑1 (MQ=255)
                   aCGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGttt                                                                                            >  1:487277/1‑139 (MQ=255)
                    cGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                               <  2:688226/136‑1 (MQ=255)
                    cGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                               >  1:688226/1‑136 (MQ=255)
                     gggCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAg                                                                                                               >  1:1060/1‑118 (MQ=255)
                     gggCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAg                                                                                                               <  2:1060/118‑1 (MQ=255)
                                aCTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACcgc                                                                               <  2:687235/139‑1 (MQ=255)
                                  tACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGt                                                                             >  1:580087/1‑139 (MQ=255)
                                            tctTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCt                                                                   >  1:446828/1‑139 (MQ=255)
                                             ctTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTa                                                                  <  2:487277/139‑1 (MQ=255)
                                              ttCTGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGATGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTAc                                                                 >  1:143537/1‑139 (MQ=255)
                                                        tCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGATGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGGCGCGa                                                       <  2:143537/139‑1 (MQ=255)
                                                                   tGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAg                                                                                                               <  2:415132/72‑1 (MQ=255)
                                                                   tGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAg                                                                                                               >  1:415132/1‑72 (MQ=255)
                                                                    gCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGttttttt                                                                                                                        >  1:67445/1‑63 (MQ=255)
                                                                    gCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGttttttt                                                                                                                        <  2:67445/63‑1 (MQ=255)
                                                                          tcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTgc                                     >  1:110510/1‑139 (MQ=255)
                                                                                 ccAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGAcc                              <  2:435971/139‑1 (MQ=255)
                                                                                    gTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGt                           >  2:320299/1‑139 (MQ=255)
                                                                                           cgcCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGc                    <  2:446828/139‑1 (MQ=255)
                                                                                           cgcCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGc                    >  1:56094/1‑139 (MQ=255)
                                                                                            gcCAGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCt                   <  1:320299/139‑1 (MQ=255)
                                                                                               aGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                               >  1:47617/1‑60 (MQ=255)
                                                                                               aGACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                               <  2:47617/60‑1 (MQ=255)
                                                                                                gACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAg                                                                                                               >  1:134544/1‑43 (MQ=255)
                                                                                                gACGTAAGGCGCTGAGAATTGCCCGACG‑TTTTTTTCCGGTAAg                                                                                                               <  2:134544/43‑1 (MQ=255)
                                                                                                             gagaATTGCCCGACG‑TTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGa  <  2:580087/139‑1 (MQ=255)
                                                                                                                            |                                                                                                                            
TCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGA  >  NC_000913/723397‑723645

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