Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	983,623	G→A	*183* (TAG→TAA)	ycbK →	M15A protease‑related family periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	983,623	0	G	A	100.0%	254.0 / NA	77	*183* (TAG→TAA)	ycbK	M15A protease‑related family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/37);  total (40/37)

GGATTTCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGC  >  NC_000913/983491‑983753                                                                                                                                     |                                                                                                                                   ggATTTCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTg                                                                                                                              >  1:385296/1‑139 (MQ=255)   aTTTCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCt                                                                                                                            <  1:501171/139‑1 (MQ=255)               aaGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgg                                                                                                                <  1:372477/139‑1 (MQ=255)                 ggTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  2:899037/139‑1 (MQ=255)                     tCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  1:580012/1‑135 (MQ=255)                     tCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  2:580012/135‑1 (MQ=255)                            aaGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACt                                                                                                   >  2:990377/1‑139 (MQ=255)                             aGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTaa                                                                                                                         >  2:1168509/1‑116 (MQ=255)                             aGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTaa                                                                                                                         <  1:1168509/116‑1 (MQ=255)                             aGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                                  <  2:1554313/139‑1 (MQ=255)                                  tatTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGta                                                                                             >  2:984782/1‑139 (MQ=255)                                   atTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  2:1996140/1‑121 (MQ=255)                                   atTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  1:1996140/121‑1 (MQ=255)                                            gCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGAGCAGTATGAACTATCGTATTATTCCGGt                                                                                   >  2:1430669/1‑139 (MQ=255)                                              cgcgTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGAGCAGTATGAACTATCGTATTATTCCGGTCa                                                                                 <  1:1430669/139‑1 (MQ=255)                                                    aTCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                           <  2:385296/139‑1 (MQ=255)                                                       tATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  1:14603/101‑1 (MQ=255)                                                       tATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  2:14603/1‑101 (MQ=255)                                                         tGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTcc                                                                      >  1:445926/1‑139 (MQ=255)                                                               cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  2:1703624/1‑93 (MQ=255)                                                               cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  1:1703624/93‑1 (MQ=255)                                                                   tggtgTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGtt                                                            <  2:445926/139‑1 (MQ=255)                                                                     gtgtAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTaataat                                                                                                                                <  1:1846648/69‑1 (MQ=255)                                                                     gtgtAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTaataat                                                                                                                                >  2:1846648/1‑69 (MQ=255)                                                                          ggATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACt                                                               <  1:112730/129‑1 (MQ=255)                                                                          ggATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACt                                                               >  2:112730/1‑129 (MQ=255)                                                                            atatTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  2:1646676/1‑80 (MQ=255)                                                                            atatTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  1:1646676/80‑1 (MQ=255)                                                                             tattatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAggg                                                                                                               >  1:508030/1‑78 (MQ=255)                                                                             tattatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAggg                                                                                                               <  2:508030/78‑1 (MQ=255)                                                                                tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  2:1810632/76‑1 (MQ=255)                                                                                tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  1:1810632/1‑76 (MQ=255)                                                                                tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                     >  1:1695319/1‑85 (MQ=255)                                                                                tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                     <  2:1695319/85‑1 (MQ=255)                                                                                tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                                  >  2:831553/1‑88 (MQ=255)                                                                                tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                                  <  1:831553/88‑1 (MQ=255)                                                                                  tCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGaaa                                                                                                                    >  2:788938/1‑68 (MQ=255)                                                                                  tCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGaaa                                                                                                                    <  1:788938/68‑1 (MQ=255)                                                                                  tCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                           <  1:1401260/109‑1 (MQ=255)                                                                                  tCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                           >  2:1401260/1‑109 (MQ=255)                                                                                    cACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGaa                                           >  2:1760658/1‑139 (MQ=255)                                                                                       gtagtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGTAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACaa                                        >  2:1506111/1‑139 (MQ=255)                                                                                       gtagtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgt                                                >  1:94110/1‑131 (MQ=255)                                                                                       gtagtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgt                                                <  2:94110/131‑1 (MQ=255)                                                                                                 tgtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTgg                              >  2:399438/1‑139 (MQ=255)                                                                                                  gtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAggg                                                                                                               <  1:1995695/57‑1 (MQ=255)                                                                                                  gtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAggg                                                                                                               >  2:1995695/1‑57 (MQ=255)                                                                                                   tgCATATTGATACCGGGCCAGCACGGCACTGGTTATAATTGCTTAACGAAACAggg                                                                                                               >  1:480948/1‑56 (MQ=255)                                                                                                   tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCtt                                                                                                                           <  2:1895456/44‑1 (MQ=255)                                                                                                   tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCtt                                                                                                                           >  1:1895456/1‑44 (MQ=255)                                                                                                   tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  2:1166752/57‑1 (MQ=255)                                                                                                   tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  1:1166752/1‑57 (MQ=255)                                                                                                   tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAggg                                                                                                               <  2:480948/56‑1 (MQ=255)                                                                                                      atatTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGTAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCa                         <  1:1506111/139‑1 (MQ=255)                                                                                                      atatTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCa                         <  1:990377/139‑1 (MQ=255)                                                                                                         ttGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTg                      <  1:984782/139‑1 (MQ=255)                                                                                                             tACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTAGTGTGAACAAACCCGTCTGGCCGCACTGGTCg                  >  1:1392563/1‑139 (MQ=255)                                                                                                               ccGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGAt                <  2:1777289/139‑1 (MQ=255)                                                                                                               ccGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCCGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGAt                >  2:2028645/1‑139 (MQ=255)                                                                                                                 gggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATcc              >  1:1587659/1‑139 (MQ=255)                                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAAc                                                                                                                   >  1:270797/1‑37 (MQ=255)                                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAAc                                                                                                                   <  2:270797/37‑1 (MQ=255)                                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  2:253466/1‑42 (MQ=255)                                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  1:253466/42‑1 (MQ=255)                                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                     <  1:1318874/51‑1 (MQ=255)                                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                     >  2:1318874/1‑51 (MQ=255)                                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTc                                                                                  >  2:1480052/1‑70 (MQ=255)                                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTc                                                                                  <  1:1480052/70‑1 (MQ=255)                                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                                 <  1:873327/87‑1 (MQ=255)                                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                                 >  2:873327/1‑87 (MQ=255)                                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCt             >  1:1979569/1‑139 (MQ=255)                                                                                                                        cACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTggcggc       >  2:1336809/1‑139 (MQ=255)                                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                                  >  1:27949/1‑43 (MQ=255)                                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                                  <  2:27949/43‑1 (MQ=255)                                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACt                                                               <  2:1841987/78‑1 (MQ=255)                                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACt                                                               >  1:1841987/1‑78 (MQ=255)                                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGt                                  <  2:1278153/107‑1 (MQ=255)                                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGt                                  >  1:1278153/1‑107 (MQ=255)                                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGc  >  2:562339/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                   GGATTTCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGC  >  NC_000913/983491‑983753

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