Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,237,358	C→T	S154S (AGC→AGT)	yjbG →	capsule biosynthesis GfcC family protein YjbG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,237,358	0	C	T	100.0%	205.8 / NA	58	S154S (AGC→AGT)	yjbG	capsule biosynthesis GfcC family protein YjbG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/24);  total (34/24)

GGGCAGGCAAAAAATCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGCCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGTTTACCCGGACGGACGCACGCAAAAAGCGC  >  NC_000913/4237226‑4237488                                                                                                                                      |                                                                                                                                   gggCAGGCAAAAAATCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTg                                                                                                                             >  1:1168126/1‑140 (MQ=255)      ggaaaacaaCCATTCTTGCACGCGATATCCCCCGCGTTGCCGAGGGCGGTAG‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAag                                                                                                                        <  2:325162/137‑1 (MQ=255)        aaaaaaaTCAATCTTGACCCCGATACCGCCCGCGTTGCCGAACGCGGTAC‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAagca                                                                                                                      <  1:934080/139‑1 (MQ=255)         aaaaaaTCAATCTTGACCCCGATATCGCCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAagcag                                                                                                                     <  2:761248/140‑1 (MQ=255)            aaaaCAATCTTGCACCCGGTAACCCCCCGGTTCCCCAACGGGGTAACCCCGCCGTTGC‑GGGCAACTACACGCTGTGGGTCGGACCACCACCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCa                                                                                                                  <  1:1317404/140‑1 (MQ=255)             aatcaatcTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCAt                                                                                                                 <  1:876012/140‑1 (MQ=255)                 aatcTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAagcag                                                                                                                     >  2:206693/1‑132 (MQ=255)                 aatcTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAagcag                                                                                                                     <  1:206693/132‑1 (MQ=255)                 aatcTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCAc                                                                                                             >  1:950752/1‑140 (MQ=255)                  atcTTGATCCCGATATCGCCCGCGTTGCCGAGCGCGGTAC‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACt                                                                                                            <  2:879485/140‑1 (MQ=255)                             gATATCGTCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGa                                                                                                 >  1:1447982/1‑140 (MQ=255)                               tatCGTCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACg                                                                                               >  2:1099011/1‑140 (MQ=255)                               tatCGTCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACg                                                                                               >  2:574809/1‑140 (MQ=255)                                  cgcccGCTTTGCCAACCGCGGTAC‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTgg                                                                                            <  1:897603/137‑1 (MQ=255)                                    tCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGa                                                                                                 <  1:540383/133‑1 (MQ=255)                                    tCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGa                                                                                                 >  2:540383/1‑133 (MQ=255)                                    tCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCg                                                                                          >  2:334363/1‑140 (MQ=255)                                    tCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCg                                                                                          <  1:1061204/140‑1 (MQ=255)                                    tCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCg                                                                                          >  2:63885/1‑140 (MQ=255)                                    tCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCg                                                                                          >  1:11048/1‑140 (MQ=255)                                    tCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCg                                                                                          >  2:1082998/1‑140 (MQ=255)                                    tCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCg                                                                                          >  2:1353694/1‑140 (MQ=255)                                    tCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCg                                                                                          >  2:1098370/1‑140 (MQ=255)                                    tCCGCGTTGCCGAACGCGCTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGCCGGACCACCGCCGTCCACGGGCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCg                                                                                          >  2:667557/1‑140 (MQ=255)                                    cccGCGTTGCCGACCGCGGTAC‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCg                                                                                          <  2:110366/139‑1 (MQ=255)                                     ccGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGa                                                                                         >  1:450925/1‑140 (MQ=255)                                         gTTGCTGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTa                                                                                     <  1:1034612/140‑1 (MQ=255)                                         gTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCTCCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTa                                                                                     <  2:950752/140‑1 (MQ=255)                                                       tAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCaga                                                                       >  1:375311/1‑140 (MQ=255)                                                               ccgTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCt                                                                <  2:1447982/140‑1 (MQ=255)                                                                 gTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCa                                                              >  1:195824/1‑140 (MQ=255)                                                                 gTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCa                                                              <  1:1098370/140‑1 (MQ=255)                                                                  ttGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGAcc                                                                          <  1:620307/127‑1 (MQ=255)                                                                  ttGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGAcc                                                                          >  2:620307/1‑127 (MQ=255)                                                                         gCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCg                                                      <  1:63885/140‑1 (MQ=255)                                                                             cTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATc                                                  <  1:645308/140‑1 (MQ=255)                                                                             cTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATc                                                  >  1:451108/1‑140 (MQ=255)                                                                             cTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATc                                                  >  1:377881/1‑140 (MQ=255)                                                                              tACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCg                                                 <  1:1082998/140‑1 (MQ=255)                                                                              tACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCg                                                 <  2:195824/140‑1 (MQ=255)                                                                                       tgGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTAcgc                                        >  1:312366/1‑140 (MQ=255)                                                                                       tgGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTAcgc                                        >  2:944974/1‑140 (MQ=255)                                                                                        gggTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTAcgcg                                       >  1:1046809/1‑140 (MQ=255)                                                                                          gTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCgtg                                     <  1:1099011/140‑1 (MQ=255)                                                                                           tCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCgtgg                                    <  2:1168126/140‑1 (MQ=255)                                                                                           tCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCgtgg                                    <  2:375311/140‑1 (MQ=255)                                                                                                 caccGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGtt                              >  2:202318/1‑140 (MQ=255)                                                                                                 caccGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGAGGtt                              >  2:396710/1‑140 (MQ=255)                                                                                                   ccgccgTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGTTTa                            >  1:1342180/1‑140 (MQ=255)                                                                                                   ccgccgTCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGTTTa                            >  2:956917/1‑140 (MQ=255)                                                                                                         tCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGTTTACCcgga                      >  2:1060007/1‑140 (MQ=255)                                                                                                         tCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGTTTACCcgga                      >  1:514303/1‑140 (MQ=255)                                                                                                         tCCACGGTCACGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGTTTACCcgga                      >  1:1092403/1‑140 (MQ=255)                                                                                                                   cGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGTTTACCCGGACGGacgcacg            >  1:988163/1‑140 (MQ=255)                                                                                                                   cGTTGTTCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGTTTACCCGGACGGacgcacg            >  1:411678/1‑140 (MQ=255)                                                                                                                        ttCGGGCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGTTTACCCGGACGGACGCACGCaaaa       >  1:1170733/1‑140 (MQ=255)                                                                                                                            ggCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGTTTACCCGGACGGACGCACGCAAAAAgcg   >  1:311671/1‑140 (MQ=255)                                                                                                                             gCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGTTTACCCGGACGGACGCACGCAAAAAgcgc  >  2:523673/1‑140 (MQ=255)                                                                                                                             gCTTATCAGTCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGTTTACCCGGACGGACGCACGCAAAAAgcgc  <  1:944974/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                   GGGCAGGCAAAAAATCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAA‑CCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGCCGTCCTGGCAAGCAGCCATTCACTCCAGGTCGCGACGTGGCGAGCTATCTCTCTGACCAGAGCCTGCTCAGCGGTGCGGATCGCAGCTACGCGTGGGTGGTTTACCCGGACGGACGCACGCAAAAAGCGC  >  NC_000913/4237226‑4237488

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