Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,042,663	A→G	R516R (CGT→CGC)	etk ←	tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,042,663	0	A	G	100.0%	195.3 / NA	56	R516R (CGT→CGC)	etk	tyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/27);  total (29/27)

TTTTTGATCGGACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGAT  >  NC_000913/1042525‑1042803                                                                                                                                           |                                                                                                                                             tttttGATCGGACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTcc                                                                                                                                  >  2:465769/1‑151 (MQ=255)      gATCGGACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCAcagc                                                                                                                             <  2:644727/151‑1 (MQ=255)          ggACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGaa                                                                                                                         <  2:766269/151‑1 (MQ=255)                    aTCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGt                                                                                                               >  1:1635/1‑151 (MQ=255)                       aCCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTAt                                                                                                            >  1:75180/1‑151 (MQ=255)                       aCCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTAt                                                                                                            >  2:727198/1‑151 (MQ=255)                               cAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGcc                                                                                                    <  1:544504/151‑1 (MQ=255)                                agagTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCa                                                                                                   <  2:736110/151‑1 (MQ=255)                                   gTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAgg                                                                                                 <  1:766439/150‑1 (MQ=255)                                     tGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGaa                                                                                               <  1:1582/150‑1 (MQ=255)                                     tGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGa                                                                                                >  1:16328/1‑149 (MQ=255)                                      gAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAg                                                                                              <  1:457292/150‑1 (MQ=255)                                        aCTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAgggg                                                                                           <  1:338073/151‑1 (MQ=255)                                           gACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGata                                                                                        <  2:35087/151‑1 (MQ=255)                                            aCAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGatat                                                                                       <  2:142156/151‑1 (MQ=255)                                                aCGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATAttttt                                                                                   >  1:28229/1‑151 (MQ=255)                                                           cTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATgg                                                                         <  1:408151/150‑1 (MQ=255)                                                           cTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGgc                                                                        >  2:765799/1‑151 (MQ=255)                                                           cTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGgc                                                                        >  2:570923/1‑151 (MQ=255)                                                           cTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTAAGATGgc                                                                        >  2:42155/1‑151 (MQ=255)                                                             gTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGgcgc                                                                      >  2:384066/1‑151 (MQ=255)                                                                   cGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGctgctg                                                                 >  1:750589/1‑150 (MQ=255)                                                                     tCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGAtt                                                              >  2:86511/1‑151 (MQ=255)                                                                      cGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGAtt                                                              <  1:385639/150‑1 (MQ=255)                                                                         aCCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGa                                                           >  1:85533/1‑150 (MQ=255)                                                                          ccGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGaaa                                                         <  1:476677/151‑1 (MQ=255)                                                                          ccGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGaa                                                          >  2:697857/1‑150 (MQ=255)                                                                          ccGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGaa                                                          >  1:130115/1‑150 (MQ=255)                                                                                 tCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAtt                                                  >  2:617500/1‑151 (MQ=255)                                                                                   aTTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAtttt                                                <  1:224130/151‑1 (MQ=255)                                                                                        aaTGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTa                                           <  1:517781/151‑1 (MQ=255)                                                                                          tGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACg                                         >  2:624262/1‑151 (MQ=255)                                                                                           gTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACg                                         <  1:79775/150‑1 (MQ=255)                                                                                             tATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCa                                       <  1:202209/150‑1 (MQ=255)                                                                                               ttCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGAc                                    >  2:229684/1‑151 (MQ=255)                                                                                                    cGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTg                               <  2:450001/151‑1 (MQ=255)                                                                                                     gTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGc                              >  2:579108/1‑151 (MQ=255)                                                                                                      tctcCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAtttt                                                >  1:417282/1‑132 (MQ=255)                                                                                                      tctcCATCAGAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAtttt                                                <  2:417333/132‑1 (MQ=255)                                                                                                         ccATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGGTGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCg                             <  1:755352/148‑1 (MQ=255)                                                                                                         ccATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGttt                          <  1:670001/151‑1 (MQ=255)                                                                                                           atcatAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTAt                        >  2:309479/1‑151 (MQ=255)                                                                                                              atAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCa                     <  2:8024/151‑1 (MQ=255)                                                                                                                aGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGc                   >  1:560479/1‑151 (MQ=255)                                                                                                                    aaaTGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACt               <  1:737895/151‑1 (MQ=255)                                                                                                                     aaTGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTc              >  2:142185/1‑151 (MQ=255)                                                                                                                          aGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAg                                      >  2:117765/1‑122 (MQ=255)                                                                                                                          aGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAg                                      <  1:117733/122‑1 (MQ=255)                                                                                                                          aGACTGGATCGAAACGCGCGAACGGCTTCCACCGCAGAATCCGCCGGGTTATCTACCGACAGGAAGGGGCTATTTTTAGAACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGATGGGTTTATCCAGCGACTCGGAc          >  1:125319/1‑150 (MQ=255)                                                                                                                            aCTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACAtt       >  2:31015/1‑151 (MQ=255)                                                                                                                             cTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACAtt       <  1:399801/150‑1 (MQ=255)                                                                                                                               ggTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTgg     <  2:530264/150‑1 (MQ=255)                                                                                                                               ggTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTgg     <  2:363579/150‑1 (MQ=255)                                                                                                                                gTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTggg    >  2:316884/1‑150 (MQ=255)                                                                                                                                gTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTggg    >  1:758074/1‑150 (MQ=255)                                                                                                                                  tCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGAt  >  2:352413/1‑150 (MQ=255)                                                                                                                                           |                                                                                                                                             TTTTTGATCGGACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGAT  >  NC_000913/1042525‑1042803

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