Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,562,963	(T)7→8	coding (641/861 nt)	yjiK ←	uncharacterized protein YjiK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,562,956	1	.	T	100.0%	185.0 / NA	53	N216K (AAT→AAA)	yjiK	uncharacterized protein YjiK
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (22/31);  total (22/31)

CGCAACCCCTTCGGCCTGTTTGATGTTATGCGACAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTA‑TTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACAGCCCGTTAACTTTGTAGACCTCAATCGGATTTT  >  NC_000913/4562821‑4563089                                                                                                                                         |                                                                                                                                      cGCAACCCCTTCGGCCTGTTTGATGTTATGCGACAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTA‑tttt                                                                                                                                   >  2:443166/1‑140 (MQ=255)           tCGGCCTGTTTGATGTTATGCGACAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGt                                                                                                                          >  2:88655/1‑140 (MQ=255)                          ttATGCGACAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCgg                                                                                                           <  1:382735/140‑1 (MQ=255)                           tATGCGACAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGa                                                                                                          <  1:710299/140‑1 (MQ=255)                             tGCGACAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTgcg                                                                                                               >  1:556254/1‑133 (MQ=255)                             tGCGACAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTgcg                                                                                                               <  2:556254/133‑1 (MQ=255)                             tGCGACAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGaca                                                                                                        >  1:117077/1‑140 (MQ=255)                             tGCGACAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGaca                                                                                                        >  2:23228/1‑140 (MQ=255)                             tGCGACAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGaca                                                                                                        >  2:489835/1‑140 (MQ=255)                             tGCGACAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGaca                                                                                                        >  2:587331/1‑140 (MQ=255)                             tGCGACAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGaca                                                                                                        >  2:720285/1‑140 (MQ=255)                                  cAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACAcatca                                                                                                   >  1:176947/1‑140 (MQ=255)                                  cAAACCTCAACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACAcatca                                                                                                   >  2:615225/1‑140 (MQ=255)                                   aaaCCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACAcatcat                                                                                                  <  2:838072/140‑1 (MQ=255)                                   aaaCCTCGACTGCCTTTGGTTAGCGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACAcatcat                                                                                                  <  1:691287/140‑1 (MQ=255)                                            cTGTCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTaa                                                                                         <  1:5264/140‑1 (MQ=255)                                                 tttGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAAt                                                                                       >  2:1024556/1‑137 (MQ=255)                                                 tttGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAAt                                                                                       <  1:1024556/137‑1 (MQ=255)                                                 tttGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGc                                                                                    >  1:290714/1‑140 (MQ=255)                                                         gTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAg                                                                            >  2:733093/1‑140 (MQ=255)                                                          tGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGt                                                                           >  2:267016/1‑140 (MQ=255)                                                          tGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGt                                                                           >  2:756937/1‑140 (MQ=255)                                                                 tCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGt                                                                    <  2:228680/140‑1 (MQ=255)                                                                 tCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGt                                                                    <  1:387136/140‑1 (MQ=255)                                                                   aCCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCt                                                                  <  1:529248/140‑1 (MQ=255)                                                                      gATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTg                                                               <  1:267016/140‑1 (MQ=255)                                                                      gATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTg                                                               <  1:733093/140‑1 (MQ=255)                                                                      gATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTg                                                               <  1:756937/140‑1 (MQ=255)                                                                       aTAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTg                                                               >  2:95248/1‑139 (MQ=255)                                                                       aTAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTg                                                               <  1:95248/139‑1 (MQ=255)                                                                        tAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCt                                                             <  2:482799/140‑1 (MQ=255)                                                                         aaCCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTg                                                            <  1:489835/140‑1 (MQ=255)                                                                         aaCCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTg                                                            <  2:117077/140‑1 (MQ=255)                                                                         aaCCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTg                                                            <  1:23228/140‑1 (MQ=255)                                                                              cGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGAtgtg                                                       >  1:172462/1‑140 (MQ=255)                                                                                      aGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCAt                                               <  1:383039/140‑1 (MQ=255)                                                                                      aGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCAt                                               <  1:587331/140‑1 (MQ=255)                                                                                         gTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGc                                            >  1:413998/1‑140 (MQ=255)                                                                                                 cTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAAc                                    <  1:88655/140‑1 (MQ=255)                                                                                                 cTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAAc                                    <  2:176947/140‑1 (MQ=255)                                                                                                  tGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTaa                                                    <  2:164666/123‑1 (MQ=255)                                                                                                  tGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTaa                                                    >  1:164666/1‑123 (MQ=255)                                                                                                  tGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACa                                   >  1:924245/1‑140 (MQ=255)                                                                                                   ggAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACAg                                  <  1:44237/140‑1 (MQ=255)                                                                                                        gcgcGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACAGCGCGt                             <  2:172462/140‑1 (MQ=255)                                                                                                        gcgcGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACAGCCCGt                             >  2:1021017/1‑140 (MQ=255)                                                                                                        gcgcGAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACAGCCCGt                             <  2:413998/140‑1 (MQ=255)                                                                                                         cgcgAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACAGCCCGtt                            <  1:542718/140‑1 (MQ=255)                                                                                                           cgAGATTCATGAGACAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACAGCCCGTTaa                          <  1:323931/140‑1 (MQ=255)                                                                                                                     gagaCAGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACAGCCCGTTAACTTTGTAGAc                <  1:443166/140‑1 (MQ=255)                                                                                                                          aGCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACAGCCCGTTAACTTTGTAGACCTCaa           >  1:587075/1‑140 (MQ=255)                                                                                                                           gCACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACAGCCCGTTAACTTTGTAGACCTCAAt          <  2:402502/140‑1 (MQ=255)                                                                                                                            cACCAGAAGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACAGCCCGTTAACTTTGTAGACCTCAATc         >  1:946816/1‑140 (MQ=255)                                                                                                                                   aGCGTATTTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACAGCCCGTTAACTTTGTAGACCTCAATCGGAtttt  <  2:946816/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      CGCAACCCCTTCGGCCTGTTTGATGTTATGCGACAAACCTCGACTGCCTTTGGTTAGTGACATTTCACCGATAACCTCGCCAACCAGGGTTACTTCCTGGAGTGCGCGAGATTCATGAGACAGCACCAGAAGCGTA‑TTTTTTTGCTGGTTAAATTCTGCGCCGGACACATCATCTAAAGTAAATTGCCGTTGTAGTGCCTCGTCTTTGCTGATGTGTAATTCATTGCTACTTAACAGCCCGTTAACTTTGTAGACCTCAATCGGATTTT  >  NC_000913/4562821‑4563089

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: ‑