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.C93.8% 156.6 / 3.2 48coding (518/3594 nt)evgShybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (23/22);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

TTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑CCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCA  >  NC_000913/2484754‑2485022
                                                                                                                                           |                                                                                                                                    
ttCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑c                                                                                                                                    <  1:474815/139‑1 (MQ=255)
ttCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTCACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTa                                                                                                                                       <  2:718900/138‑1 (MQ=255)
ttCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTCACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTa                                                                                                                                       >  1:718900/1‑138 (MQ=255)
  cgccgcCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑ccc                                                                                                                                  >  1:632443/1‑139 (MQ=255)
    agacagctTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCCCCCTTCACGACTCAATGCGACCGCATACCTCGCCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑ccccc                                                                                                                                <  1:798304/133‑1 (MQ=255)
                    aCATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑c                                                                                                                                    >  2:1210026/1‑119 (MQ=255)
                    aCATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑c                                                                                                                                    <  1:1210026/119‑1 (MQ=255)
                     cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtca                                                                                                                >  2:639887/1‑139 (MQ=255)
                     cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑ccccccca                                                                                                                             >  1:477208/1‑123 (MQ=255)
                     cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑ccccccca                                                                                                                             <  2:477208/125‑3 (MQ=255)
                     cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑c                                                                                                                                    <  1:1019973/118‑1 (MQ=255)
                     cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑c                                                                                                                                    >  2:1019973/1‑118 (MQ=255)
                     cATTGCTGCAACCAAACCATTGATAATTACATTTCCGCCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTTCTCGGGTACAAAATTAA‑‑CCCCCCAGACTAGGtctgtct                                                                                                                >  2:798304/1‑133 (MQ=255)
                             cAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCAt                                                                                                        >  2:111302/1‑139 (MQ=255)
                               accaaaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCAttt                                                                                                      >  1:614757/1‑139 (MQ=255)
                                        ttGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGca                                                                                             >  2:195990/1‑139 (MQ=255)
                                          gATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaac                                                                                           >  1:219024/1‑139 (MQ=255)
                                                      tCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtcat                                                                                                               <  1:338203/107‑1 (MQ=255)
                                                      tCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtcat                                                                                                               >  2:338203/1‑107 (MQ=255)
                                                       ccGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑ccccc                                                                                                                                <  1:131073/88‑1 (MQ=255)
                                                       ccGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑ccccc                                                                                                                                >  2:131073/1‑88 (MQ=255)
                                                          gcgcTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACa                                                                           <  2:219024/139‑1 (MQ=255)
                                                             cTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAAt                                                                        >  1:331940/1‑139 (MQ=255)
                                                                gTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTta                                                                     >  1:34007/1‑139 (MQ=255)
                                                                  aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTtata                                                                   <  2:34007/139‑1 (MQ=255)
                                                                   accagcCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTatat                                                                  <  1:821813/139‑1 (MQ=255)
                                                                     caccCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑cccccccaga                                                                                                                           <  2:167567/79‑5 (MQ=255)
                                                                     caccCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑cccccccaga                                                                                                                           >  1:167567/1‑75 (MQ=255)
                                                                        ccTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtcat                                                                                                               >  2:928684/1‑89 (MQ=255)
                                                                        ccTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtcat                                                                                                               <  1:928684/89‑1 (MQ=255)
                                                                          ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCAt                                                           <  1:111302/139‑1 (MQ=255)
                                                                               gACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGc                                                       >  2:500576/1‑139 (MQ=255)
                                                                               gACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                      <  2:778726/139‑1 (MQ=255)
                                                                                      tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGtct                                               >  2:86489/1‑139 (MQ=255)
                                                                                      tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGtct                                               >  1:477797/1‑139 (MQ=255)
                                                                                            cGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTg                                         >  2:622506/1‑139 (MQ=255)
                                                                                               ttACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGc                                      <  2:477797/139‑1 (MQ=255)
                                                                                                tACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGca                                     <  1:528350/139‑1 (MQ=255)
                                                                                                      aCCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACg                                                                                                                          <  2:355612/48‑1 (MQ=255)
                                                                                                      aCCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACg                                                                                                                          >  1:355612/1‑48 (MQ=255)
                                                                                                        cAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑cccccccaga                                                                                                                           >  1:405269/1‑40 (MQ=255)
                                                                                                        cAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑cccccccaga                                                                                                                           <  2:405269/44‑5 (MQ=255)
                                                                                                          aaaCCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTAc                           <  2:331940/139‑1 (MQ=255)
                                                                                                            aCCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                      <  1:1258805/110‑1 (MQ=255)
                                                                                                            aCCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                      >  2:1258805/1‑110 (MQ=255)
                                                                                                                gTAAATATTGCTCGGGTAGCAAATTACCCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAt                      <  1:500576/139‑1 (MQ=255)
                                                                                                                gTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGcaca                                   <  1:355317/125‑1 (MQ=255)
                                                                                                                gTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGcaca                                   >  2:355317/1‑125 (MQ=255)
                                                                                                                 tAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAg                                            >  1:1530683/1‑115 (MQ=255)
                                                                                                                 tAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACAAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAg                                            <  2:1530683/115‑1 (MQ=255)
                                                                                                                      atTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGgtagt               <  2:1301319/139‑1 (MQ=255)
                                                                                                                          cTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                             <  1:679054/73‑1 (MQ=255)
                                                                                                                          cTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                             >  2:679054/1‑73 (MQ=255)
                                                                                                                               gTAGCAAATTACCCCCCCAGACGAGGTAATTCatcaatca                                                                                                         <  2:399255/40‑1 (MQ=255)
                                                                                                                               gTAGCAAATTACCCCCCCAGACGAGGTAATTCatcaatca                                                                                                         >  1:399255/1‑40 (MQ=255)
                                                                                                                               gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGca                                     <  2:900815/108‑1 (MQ=255)
                                                                                                                               gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGca                                     >  1:900815/1‑108 (MQ=255)
                                                                                                                               gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGa                               >  2:1375807/1‑114 (MQ=255)
                                                                                                                               gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGa                               <  1:1375807/114‑1 (MQ=255)
                                                                                                                                 aGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTAc    <  1:639887/139‑1 (MQ=255)
                                                                                                                                   cAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGca                                     >  2:623325/1‑104 (MQ=255)
                                                                                                                                   cAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGca                                     <  1:623325/104‑1 (MQ=255)
                                                                                                                                   cAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTAcaa  >  2:250942/1‑137 (MQ=255)
                                                                                                                                   cAAATTA‑‑CCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                              <  1:1511509/46‑1 (MQ=255)
                                                                                                                                   cAAATTA‑‑CCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                              >  2:1511509/1‑46 (MQ=255)
                                                                                                                                           |                                                                                                                                    
TTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑CCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCA  >  NC_000913/2484754‑2485022

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: