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
                                                                                                                            |                                                                                                                            
tCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAg                                                                                                               <  2:456929/139‑1 (MQ=255)
     ccGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTg                                                                                                          >  2:86881/1‑139 (MQ=255)
            tATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTgcagggca                                                                                                   <  1:493013/139‑1 (MQ=255)
             aTGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTgcagggcag                                                                                                  <  1:58416/139‑1 (MQ=255)
              tGTTCACGGGC‑‑AACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                               <  1:454151/139‑1 (MQ=255)
                ttCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                               <  1:86881/139‑1 (MQ=255)
                  cACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGtt                                                                                             <  1:149110/139‑1 (MQ=255)
                   aCGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGttt                                                                                            <  2:96206/139‑1 (MQ=255)
                   aCGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGttt                                                                                            >  1:487277/1‑139 (MQ=255)
                    cGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                               <  2:688226/136‑1 (MQ=255)
                    cGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                               >  1:688226/1‑136 (MQ=255)
                     gggCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAg                                                                                                               >  1:1060/1‑118 (MQ=255)
                     gggCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAg                                                                                                               <  2:1060/118‑1 (MQ=255)
                                aCTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACcgc                                                                               <  2:687235/139‑1 (MQ=255)
                                  tACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGt                                                                             >  1:580087/1‑139 (MQ=255)
                                            tctTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCt                                                                   >  1:446828/1‑139 (MQ=255)
                                             ctTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTa                                                                  <  2:487277/139‑1 (MQ=255)
                                              ttCTGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGATGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTAc                                                                 >  1:143537/1‑139 (MQ=255)
                                                        tCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGATGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGGCGCGa                                                       <  2:143537/139‑1 (MQ=255)
                                                                   tGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAg                                                                                                               <  2:415132/72‑1 (MQ=255)
                                                                   tGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAg                                                                                                               >  1:415132/1‑72 (MQ=255)
                                                                    gCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGttttttt                                                                                                                        >  1:67445/1‑63 (MQ=255)
                                                                    gCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGttttttt                                                                                                                        <  2:67445/63‑1 (MQ=255)
                                                                          tcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTgc                                     >  1:110510/1‑139 (MQ=255)
                                                                                 ccAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGAcc                              <  2:435971/139‑1 (MQ=255)
                                                                                    gTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGt                           >  2:320299/1‑139 (MQ=255)
                                                                                           cgcCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGc                    <  2:446828/139‑1 (MQ=255)
                                                                                           cgcCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGc                    >  1:56094/1‑139 (MQ=255)
                                                                                            gcCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCt                   <  1:320299/139‑1 (MQ=255)
                                                                                               aGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                               >  1:47617/1‑60 (MQ=255)
                                                                                               aGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                               <  2:47617/60‑1 (MQ=255)
                                                                                                gACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAg                                                                                                               >  1:134544/1‑43 (MQ=255)
                                                                                                gACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAg                                                                                                               <  2:134544/43‑1 (MQ=255)
                                                                                                             gagaATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGa  <  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: