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%	227.8 / NA	70	*183* (TAG→TAA)	ycbK	M15A protease‑related family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/32);  total (38/32)

GATGGATTTCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGC  >  NC_000913/983488‑983753                                                                                                                                        |                                                                                                                                   gATGGATTTCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGtaataa                                                                                                                                 <  1:762610/139‑1 (MQ=255)         tCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTaa                                                                                                                         >  1:1553194/1‑139 (MQ=255)                 gAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAg                                                                                                                 >  1:739782/1‑139 (MQ=255)                    ggTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  1:427558/139‑1 (MQ=255)                    ggTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  1:154736/139‑1 (MQ=255)                       aTCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAg                                                                                                           >  2:803052/1‑139 (MQ=255)                         cgcgTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTa                                                                                                         >  2:904952/1‑139 (MQ=255)                             ttAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                     <  1:834011/139‑1 (MQ=255)                                aGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                                  <  2:1553194/139‑1 (MQ=255)                                                ccGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  2:16771/111‑1 (MQ=255)                                                ccGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  1:16771/1‑111 (MQ=255)                                                ccGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTc                                                                                  <  1:904952/139‑1 (MQ=255)                                                      tATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCa                                                                            <  2:698528/139‑1 (MQ=255)                                                       aTCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                           <  1:875144/139‑1 (MQ=255)                                                         cTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATtc                                                                         <  1:803052/139‑1 (MQ=255)                                                             gcgcgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  2:177723/98‑1 (MQ=255)                                                             gcgcgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  1:177723/1‑98 (MQ=255)                                                                cgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGa                                                                  >  1:1262032/1‑139 (MQ=255)                                                                 gcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                                 <  2:1262032/139‑1 (MQ=255)                                                                  cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTcc                                                                      >  1:1413023/1‑133 (MQ=255)                                                                  cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTcc                                                                      <  2:1413023/133‑1 (MQ=255)                                                                    ggtggtGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTg                                                              <  1:861307/139‑1 (MQ=255)                                                                      tggtgTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  2:313962/1‑89 (MQ=255)                                                                      tggtgTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  1:313962/89‑1 (MQ=255)                                                                         tgtAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTg                                                              >  1:416742/1‑134 (MQ=255)                                                                         tgtAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTg                                                              <  2:416742/134‑1 (MQ=255)                                                                           tAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAAc                                                                                                                   <  1:1560469/79‑1 (MQ=255)                                                                           tAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAAc                                                                                                                   >  2:1560469/1‑79 (MQ=255)                                                                               atatTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCAc                                                                                >  2:1264580/1‑110 (MQ=255)                                                                               atatTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCAc                                                                                <  1:1264580/110‑1 (MQ=255)                                                                                   tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  1:623785/76‑1 (MQ=255)                                                                                   tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  2:623785/1‑76 (MQ=255)                                                                                   tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAtt                                                                          >  1:155594/1‑112 (MQ=255)                                                                                   tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAtt                                                                          <  2:155594/112‑1 (MQ=255)                                                                                   tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgtg                                               >  2:880095/1‑139 (MQ=255)                                                                                        aCGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAAc                                          >  2:1345715/1‑139 (MQ=255)                                                                                        aCGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAAc                                          >  1:1256017/1‑139 (MQ=255)                                                                                           tagtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACaaa                                       <  1:1345715/139‑1 (MQ=255)                                                                                                aCTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGt                                  >  2:10055/1‑139 (MQ=255)                                                                                                    tgtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                           >  2:760129/1‑94 (MQ=255)                                                                                                    tgtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                           <  1:760129/94‑1 (MQ=255)                                                                                                     gtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATg                                                                                                       >  2:643191/1‑65 (MQ=255)                                                                                                     gtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATg                                                                                                       <  1:643191/65‑1 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTaa                                                                                                                         <  2:975295/46‑1 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTaa                                                                                                                         >  1:975295/1‑46 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  2:1481453/57‑1 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  1:1481453/1‑57 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGTGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgt                                                <  2:897588/119‑1 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATg                                                                                                       <  2:1436804/64‑1 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATg                                                                                                       >  1:1436804/1‑64 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgt                                                <  2:1284405/119‑1 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgt                                                >  1:1284405/1‑119 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgt                                                >  1:897588/1‑119 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGcc                            >  2:1408176/1‑139 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAggg                                                                                                               <  2:303812/41‑1 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAggg                                                                                                               >  1:303812/1‑41 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                     >  2:1228383/1‑51 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                     <  1:1228383/51‑1 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                                 >  1:1349666/1‑87 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                                 <  2:1349666/87‑1 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCt             >  2:1287548/1‑139 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCt             >  1:901926/1‑139 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCt             >  2:609890/1‑139 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCt             >  1:1027058/1‑139 (MQ=255)                                                                                                                      gCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              <  2:249505/41‑1 (MQ=255)                                                                                                                      gCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                              >  1:249505/1‑41 (MQ=255)                                                                                                                       ccAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                     <  2:1280452/49‑1 (MQ=255)                                                                                                                       ccAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                     >  1:1280452/1‑49 (MQ=255)                                                                                                                           cACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTggcggc       >  2:825928/1‑139 (MQ=255)                                                                                                                                cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGc  >  2:1221682/1‑139 (MQ=255)                                                                                                                                cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGc  >  1:1067686/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                   GATGGATTTCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGC  >  NC_000913/983488‑983753

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