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.T96.4% 169.4 / ‑0.6 55coding (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/2);  new base T (24/29);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GCGACCGAGAATAATCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAA  >  NC_000913/723383‑723652
                                                                                                                                           |                                                                                                                                     
gCGACCGAGAATAATCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑t                                                                                                                                     >  2:73225/1‑139 (MQ=255)
         aataatCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑ttttttttt                                                                                                                             <  2:347351/138‑2 (MQ=255)
         aataatCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑ttttttttt                                                                                                                             >  1:347351/1‑137 (MQ=255)
                ttGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAg                                                                                                                      <  2:114658/139‑1 (MQ=255)
                ttGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAg                                                                                                                      <  1:24011/139‑1 (MQ=255)
                   ccGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCgg                                                                                                                   >  2:557580/1‑139 (MQ=255)
                     gggaTTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATAAAAAAGTGTTGCCGTCCCCGCCCCCATTTCCTGCGCCAGACGTAAGGCGTTGAGAATTGCCCGACTT‑‑TTTTTTTTCCGGTAAGCGGTg                                                                                                                 <  2:1450829/137‑1 (MQ=255)
                     gagaTTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTg                                                                                                                 <  1:733512/139‑1 (MQ=255)
                     gagaTTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTg                                                                                                                 <  1:1367918/139‑1 (MQ=255)
                           aTGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTgcagggc                                                                                                           <  1:474936/139‑1 (MQ=255)
                             gTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTgcagggcag                                                                                                         <  1:557580/139‑1 (MQ=255)
                                 aCGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTg                                                                                                                 >  2:305713/1‑127 (MQ=255)
                                 aCGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTg                                                                                                                 <  1:305713/127‑1 (MQ=255)
                                  cGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGtt                                                                                                    >  2:1321092/1‑139 (MQ=255)
                                  cGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGtt                                                                                                    <  2:369122/139‑1 (MQ=255)
                                  cGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGtt                                                                                                    <  1:73225/139‑1 (MQ=255)
                                                     cTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGcc                                                                                 <  1:1321092/139‑1 (MQ=255)
                                                        tctcTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGa                                                                              <  2:637440/139‑1 (MQ=255)
                                                          tctTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACg                                                                            >  1:608495/1‑139 (MQ=255)
                                                                gCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTAccc                                                                      >  2:1092557/1‑139 (MQ=255)
                                                                 cTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAg                                                                                                                      <  2:485557/90‑1 (MQ=255)
                                                                 cTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAg                                                                                                                      >  1:485557/1‑90 (MQ=255)
                                                                              tgtTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAg                                                        >  1:1334565/1‑139 (MQ=255)
                                                                              tgtTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAg                                                        >  1:763837/1‑139 (MQ=255)
                                                                              tgtTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAg                                                        >  2:447214/1‑139 (MQ=255)
                                                                              tgtTGCCGTCTCCGCGCCCAGTTCATGCGCCAGACGTAAGGCGCTGAGAATGGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCTAGCCGTGACGCCAg                                                        >  2:899189/1‑139 (MQ=255)
                                                                                    cGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAgccgcgc                                                  >  2:62839/1‑139 (MQ=255)
                                                                                      tctcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGCCGTTTTTTTTTTCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAgccgcgccg                                                <  2:763837/139‑1 (MQ=255)
                                                                                        tcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAg                                                                                                                      <  1:303903/67‑1 (MQ=255)
                                                                                        tcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAg                                                                                                                      >  2:303903/1‑67 (MQ=255)
                                                                                        tcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                                      >  1:761843/1‑83 (MQ=255)
                                                                                        tcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAgggg                                                                                                      <  2:761843/83‑1 (MQ=255)
                                                                                        tcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGttt                                                                                                   <  1:1090298/86‑1 (MQ=255)
                                                                                        tcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGttt                                                                                                   >  2:1090298/1‑86 (MQ=255)
                                                                                           gcgcCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTgcg                                           <  1:1092557/139‑1 (MQ=255)
                                                                                           gcgcCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTgcg                                           <  2:487296/139‑1 (MQ=255)
                                                                                                   ttCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTgcagggca                                                                                                          >  2:1386758/1‑68 (MQ=255)
                                                                                                   ttCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTgcagggca                                                                                                          <  1:1386758/68‑1 (MQ=255)
                                                                                                           ccAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGCTTTTCGc                           <  2:608495/139‑1 (MQ=255)
                                                                                                             aGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTg                         >  1:171791/1‑139 (MQ=255)
                                                                                                              gACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTgcagggcagg                                                                                                        <  2:199173/59‑1 (MQ=255)
                                                                                                              gACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTgcagggcagg                                                                                                        >  1:199173/1‑59 (MQ=255)
                                                                                                                cGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCg                      <  1:62839/139‑1 (MQ=255)
                                                                                                                   aaGGCGCTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGgtg                   >  2:1074521/1‑139 (MQ=255)
                                                                                                                         cTGAGAATTGCCCGACGTTTTTTTTTTCCAGTAAGCGGTgcagggcag                                                                                                         <  2:943466/48‑1 (MQ=25)
                                                                                                                         cTGAGAATTGCCCGACGTTTTTTTTTTCCAGTAAGCGGTgcagggcag                                                                                                         >  1:943466/1‑48 (MQ=37)
                                                                                                                         cTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTg                                                                                                                 <  2:1457877/39‑1 (MQ=37)
                                                                                                                         cTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTg                                                                                                                 >  1:1457877/1‑39 (MQ=37)
                                                                                                                         cTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATgg             >  2:975075/1‑139 (MQ=255)
                                                                                                                         cTGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATgg             >  2:551658/1‑139 (MQ=255)
                                                                                                                          tGAGAATTGCCCGACAT‑‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGc            <  1:143871/139‑1 (MQ=255)
                                                                                                                          tGAGAATTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGc            <  2:1193593/139‑1 (MQ=255)
                                                                                                                               aTTGCCCGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATg       <  1:1074521/139‑1 (MQ=255)
                                                                                                                                   cccGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGtt                                                                                                    <  1:1095578/42‑1 (MQ=255)
                                                                                                                                   cccGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGtt                                                                                                    >  2:1095578/1‑42 (MQ=38)
                                                                                                                                   cccGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGa         >  1:1105426/1‑133 (MQ=255)
                                                                                                                                   cccGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGa         <  2:1105426/133‑1 (MQ=255)
                                                                                                                                    ccGACGTTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCCGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATaa  <  1:899189/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                     
GCGACCGAGAATAATCTTGCCGAGATTATGTTCACGGGCATAACGCACTACCGCTTTCTCTTCCGCTGGATCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACG‑‑TTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAA  >  NC_000913/723383‑723652

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