Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,684,019 C→T 10.5% R380W (CGG→TGG)  rstB → sensory histidine kinase RstB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,684,0190CT10.5% 120.7 / 5.9 51R380W (CGG→TGG) rstBsensory histidine kinase RstB
Reads supporting (aligned to +/- strand):  ref base C (21/23);  new base T (2/3);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

CCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATG  >  NC_000913/1683886‑1684152
                                                                                                                                      |                                                                                                                                     
ccATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCCTTGAACCTTTTGTACGCCTCGATCGCCGACGGGATCgtt                                                                                                                              >  2:250126/1‑141 (MQ=255)
 gaTTCAACGGTTGAACCCAGCCTGCTACTGTCGGGGAATAGAGCGACATGAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgagg                                                                                                                             <  1:1101814/142‑4 (MQ=255)
 cATTTAGCGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATCGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgagg                                                                                                                             <  2:506318/143‑4 (MQ=255)
    tCATCGGTTGAATCCAGCCTGCTACTGTCGGGGAATAGAGCGACAGTAATAGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgaggcga                                                                                                                          <  1:531647/143‑7 (MQ=255)
             gAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTGGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTg                                                                                                                 <  2:225730/143‑1 (MQ=255)
              aaaCCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGc                                                                                                                >  1:549904/1‑143 (MQ=255)
               aaCCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCggcggccacg                                                                                                               >  2:366898/1‑139 (MQ=255)
                 ccAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGTGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCggg                                                                                                             <  2:280311/143‑1 (MQ=255)
                    gCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctg                                                                                                          >  1:559566/1‑143 (MQ=255)
                       tGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGGTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctgggg                                                                                                       <  1:432253/143‑1 (MQ=255)
                       cgccaCTGTCGGGGAATAGAGCGACATTAATCGTCGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACACATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctgggg                                                                                                       <  1:544181/139‑1 (MQ=255)
                          tACTGTCGGGGAATAGAGCGACATTAATTGGTGAGGATGATGGCCCAGGGATTGCCCCAGAGAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggctg                                                                                                    <  1:96044/143‑1 (MQ=255)
                           aCTGTCGGGGAATAGAGCGAAATAAATTTGTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGggctggggctgg                                                                                                   <  2:442737/143‑1 (MQ=255)
                              gTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCaa                                                                                                <  1:111765/143‑1 (MQ=255)
                                  gggAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGt                                                                                            <  1:1049370/143‑1 (MQ=255)
                                   gaaaTAGGGCGGCGTTAATTGTTGAGGTTGATGCCCCAGGCATTGCCCCAGAAAACCGCGAACATATCTTCGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTc                                                                                           <  1:227078/141‑1 (MQ=255)
                                    gagtAGAGCGACAGTAACTGTGGGGGATGATGGCCCAGGGATTGCCCCAGAAATCCGCGAACATATCTTTGAACCTTTTGTTCGCCGCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTcc                                                                                          <  1:586274/140‑1 (MQ=255)
                                         gagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCActct                                                                                     <  1:1109842/143‑1 (MQ=255)
                                           gCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAAGTGTCCACTCtat                                                                                   >  1:539890/1‑143 (MQ=255)
                                                    atttctcaGGACGATTGCCCAGGGAGTGCCCCAGAAAACCGCTAACCTATCTTTGACCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGTTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACCGGc                                                                          <  2:987910/136‑1 (MQ=255)
                                                          gAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATggg                                                                    <  1:888409/143‑1 (MQ=255)
                                                          gAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCCCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTGTAGCACTGGCCGTggg                                                                    >  1:854184/1‑143 (MQ=255)
                                                                gatgGTTCAGGGATTGCCCCAGAAAACCGCGAGCAGATCTTTGAACCTGTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTAc                                                              <  1:531420/143‑1 (MQ=255)
                                                                gatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTAc                                                              >  1:840889/1‑143 (MQ=255)
                                                                 atgGCCCAGGGCTGGCGCCAGAAAACCGCGAACATATCTTTGAACCTTTGTTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACg                                                             <  1:591500/143‑1 (MQ=255)
                                                                   ggCCCACGGCTTGCCCGCGAACTCGTCGAACTTAGGTTTGCACCTTGTGTCTCGCTCGATCCCGGCCGGGATATCTCAACCGGAGGCTGCGGGCTGTGGTTGGTAATTGTCCATTCCATAGCACTGGCAATGGGCGGTACGGt                                                           <  2:1034010/143‑1 (MQ=255)
                                                                    gCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCCACCGGCGGCTGCGGGCTGGGGCTGGCAAATCTCCACGCTATAGCACTCGCAATGGGCGGTACCGtc                                                          >  1:248053/1‑142 (MQ=255)
                                                                               gCCCCCGGAAACCGCGAAGAGATCTGTGAACCTTTTGTGCGCCTCGATCCCAGCGGGGATCGTTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGGAATGGGCGGTACGGTTAATTGTGacac                                               <  2:563149/143‑1 (MQ=255)
                                                                                 cccAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCa                                             >  1:116456/1‑143 (MQ=255)
                                                                                  ccAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAg                                            <  1:113225/143‑1 (MQ=255)
                                                                                       aaaCCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGGTAGTTGTGACACCAGCGGAc                                       >  1:563285/1‑143 (MQ=255)
                                                                                         aCCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGGTGATTGTGGAGCCAGCGAAccg                                     >  1:698565/1‑141 (MQ=255)
                                                                                           cgcgAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTcgg                                   >  2:107325/1‑143 (MQ=255)
                                                                                             cgAACATATCTTTGAACCTTGTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGggtg                                 <  2:487795/143‑1 (MQ=255)
                                                                                               accATATCTTTGAACCTTTTGTCCGCCTCGACCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCAGCACGGTTAATTGTGACACCAGCGAACTGggtggt                               <  1:267298/141‑1 (MQ=255)
                                                                                                 cATATCTTTGATCCTTTTGTCCGCCTCGATCCCAGCCGGGATCGCTCAGCCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGc                             <  1:1101349/143‑1 (MQ=255)
                                                                                                 cATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGc                             >  2:1005445/1‑143 (MQ=255)
                                                                                                   tatCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAGTTGTCCACTCTATAGCACTGGCAATGGGCGGTAAGGATAATTGGGAGACCAGCGAACTGGGTGGGGccc                           >  2:66665/1‑143 (MQ=255)
                                                                                                   tatCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATCGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGccc                           >  1:109119/1‑143 (MQ=255)
                                                                                                   tatCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGccc                           >  1:972497/1‑143 (MQ=255)
                                                                                                       tttGAACCTTTTGTTCGCCTCGAGCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCGGGCAATTGTCCTCTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCtt                       <  2:776543/143‑1 (MQ=255)
                                                                                                          gAATCTTTAGTTCGCGACGATCCCAGCAGGGATCGGTCAACCGGCGGCGGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTtctc                    <  1:275802/143‑1 (MQ=255)
                                                                                                             ccTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCCATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGGTAATTGTGACACCCGCGAACTGGGTGGCGCCCGCTTCTCGtt                 >  1:393113/1‑143 (MQ=255)
                                                                                                             ccTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGtt                 <  1:1005445/143‑1 (MQ=255)
                                                                                                             ccTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGtt                 >  2:361693/1‑143 (MQ=255)
                                                                                                               ttttGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTa               >  2:586864/1‑143 (MQ=255)
                                                                                                               ttttGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTa               >  1:843491/1‑143 (MQ=255)
                                                                                                               ttttGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTa               <  2:840889/143‑1 (MQ=255)
                                                                                                                         gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGGGACGAGCACACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGt      >  1:515513/4‑143 (MQ=255)
                                                                                                                           gcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGGTTAGCTGGCCGTTa    >  2:460915/2‑143 (MQ=255)
                                                                                                                             gATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGGTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGGTTAGCGGGCCGtgctc  >  1:94876/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                     
CCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATG  >  NC_000913/1683886‑1684152

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