Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,484,897 (C)6→7 coding (524/3594 nt) evgS → hybrid sensory histidine kinase in two‑component regulatory system with EvgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,484,8911.C98.2% 181.6 / ‑3.7 54coding (518/3594 nt)evgShybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (25/28);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑A‑CCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAG  >  NC_000913/2484754‑2485023
                                                                                                                                           |                                                                                                                                    
ttCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑A‑c                                                                                                                                     >  2:84680/1‑139 (MQ=255)
ttCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑A‑c                                                                                                                                     <  1:305667/139‑1 (MQ=255)
     cgccTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑A‑cccccc                                                                                                                                >  1:624525/1‑139 (MQ=255)
                 atGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTaa                                                                                                                     <  1:100156/139‑1 (MQ=255)
                  tGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAAAGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGATGTAAt                                                                                                                    >  1:479946/1‑139 (MQ=255)
                    aCATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑A‑c                                                                                                                                     >  2:382676/1‑119 (MQ=255)
                    aCATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑A‑c                                                                                                                                     <  1:382676/119‑1 (MQ=255)
                       ttGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATtcatc                                                                                                               >  2:704882/1‑139 (MQ=255)
                       ttGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATtcatc                                                                                                               <  2:57782/139‑1 (MQ=255)
                                       aTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                               <  1:668572/139‑1 (MQ=255)
                                       aTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                               <  2:8815/139‑1 (MQ=255)
                                        ttGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGca                                                                                              >  2:542708/1‑139 (MQ=255)
                                         tGATAATTTCCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATTCGACCGCGTACCTCACCAAACCCCGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                             <  2:498306/139‑1 (MQ=255)
                                         tGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                             >  2:715331/1‑139 (MQ=255)
                                         tGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                             >  1:151312/1‑139 (MQ=255)
                                                tACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAt                                                                                      >  1:745922/1‑139 (MQ=255)
                                                  ccTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAtct                                                                                    >  2:90048/1‑139 (MQ=255)
                                                     ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAtct                                                                                    <  2:769757/136‑1 (MQ=255)
                                                     ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAtct                                                                                    >  1:769757/1‑136 (MQ=255)
                                                        cGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                              <  1:406287/139‑1 (MQ=255)
                                                         ggCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTAc                                                                             <  2:124725/139‑1 (MQ=255)
                                                           cgcTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑A‑c                                                                                                                                     >  1:416327/1‑80 (MQ=255)
                                                           cgcTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑A‑c                                                                                                                                     <  2:416327/80‑1 (MQ=255)
                                                                   accaccCTTCACGACTCAATGCGACCGCTTACCTCACCAAATCCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTatat                                                                   <  2:125192/139‑1 (MQ=255)
                                                                   accaccCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTatat                                                                   <  2:745922/139‑1 (MQ=255)
                                                                      accCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAg                                                                <  1:84680/139‑1 (MQ=255)
                                                                          ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                               <  1:180877/104‑1 (MQ=255)
                                                                          ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                               >  2:180877/1‑104 (MQ=255)
                                                                           tCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCAtt                                                           >  2:247802/1‑139 (MQ=255)
                                                                           tCAAGAATCAATGCGACCGCTTACCTCACGAAAAACAGTAAATATTGCTAGGGTAGGAAATTAA‑CGCACCAGACGGGGTAAGTCATAAATCATTTCCAAAAGAAACAATTATATCTTTTACAAATTAATATCAGGAAtt                                                           >  2:222176/1‑139 (MQ=255)
                                                                                    aaTGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGt                                                  >  2:411117/1‑139 (MQ=255)
                                                                                           ccGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGtct                                                <  2:143544/134‑1 (MQ=255)
                                                                                           ccGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGtct                                                >  1:143544/1‑134 (MQ=255)
                                                                                              cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                             >  2:502278/1‑86 (MQ=255)
                                                                                              cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                             <  1:502278/86‑1 (MQ=255)
                                                                                              cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTggg                                        >  2:717494/1‑139 (MQ=255)
                                                                                                tACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGca                                      <  1:411117/139‑1 (MQ=255)
                                                                                                       ccAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGAt                               <  1:715331/139‑1 (MQ=255)
                                                                                                        cAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGAtt                              >  1:52693/1‑139 (MQ=255)
                                                                                                         aaaaCCAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTa                             <  1:704882/139‑1 (MQ=255)
                                                                                                             ccAGTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACttt                         >  1:178363/1‑139 (MQ=255)
                                                                                                                gTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATtcat                                                                                                                >  1:186132/1‑49 (MQ=255)
                                                                                                                gTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATtcat                                                                                                                <  2:186132/49‑1 (MQ=255)
                                                                                                                gTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGa                                <  1:754950/129‑1 (MQ=255)
                                                                                                                gTAAATATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGa                                >  2:754950/1‑129 (MQ=255)
                                                                                                                  aaaTATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                              <  2:472518/81‑1 (MQ=255)
                                                                                                                  aaaTATTGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                              >  1:472518/1‑81 (MQ=255)
                                                                                                                        tGCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTaa              <  2:624525/139‑1 (MQ=255)
                                                                                                                         gCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAc             >  2:424140/1‑139 (MQ=255)
                                                                                                                         gCTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAc             >  1:451395/1‑139 (MQ=255)
                                                                                                                          cTCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGAATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAca            >  2:408243/1‑139 (MQ=255)
                                                                                                                           tCGGGTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAcat           >  2:40203/1‑139 (MQ=255)
                                                                                                                               gTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTatat                                                                   <  1:395466/79‑1 (MQ=255)
                                                                                                                               gTAGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGAAACAATTATCTCTTTTACAAATTTatat                                                                   >  2:395466/1‑79 (MQ=255)
                                                                                                                                 aGCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTAc     <  1:429851/139‑1 (MQ=255)
                                                                                                                                  gCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTAcc    <  2:479946/139‑1 (MQ=255)
                                                                                                                                  gCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTAcc    <  2:433371/139‑1 (MQ=255)
                                                                                                                                  gCAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTAcc    <  1:90048/139‑1 (MQ=255)
                                                                                                                                   cAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAt                                                                                      >  2:629439/1‑56 (MQ=255)
                                                                                                                                   cAAATT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAt                                                                                      <  1:629439/56‑1 (MQ=255)
                                                                                                                                    aaaTT‑ACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACcag  <  2:75864/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                    
TTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATT‑A‑CCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAG  >  NC_000913/2484754‑2485023

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