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

AGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGCGG  >  NC_000913/983519‑983755                                                                                                         |                                                                                                                                     aGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                                    >  1:1280202/1‑139 (MQ=255)           gCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                                          <  1:794913/107‑1 (MQ=255)           gCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                                          >  2:794913/1‑107 (MQ=255)                gCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGt                                                                                     >  1:644781/1‑139 (MQ=255)                        aTCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                             <  2:644781/139‑1 (MQ=255)                                 cgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGa                                                                    >  2:1416698/1‑139 (MQ=255)                                   cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTaa                                                                                                                           <  1:413565/82‑1 (MQ=255)                                   cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTaa                                                                                                                           >  2:413565/1‑82 (MQ=255)                                   cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                <  2:1408023/93‑1 (MQ=255)                                   cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                >  1:1408023/1‑93 (MQ=255)                                   cAGGCGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAAc                                                                  >  1:115346/1‑139 (MQ=255)                                     ggtggtGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                       <  2:1710387/100‑1 (MQ=255)                                     ggtggtGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                       >  1:1710387/1‑100 (MQ=255)                                         gtgtAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCa                                                            <  2:1280202/139‑1 (MQ=255)                                                   ttatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgt                                                  >  2:1321964/1‑139 (MQ=255)                                                    tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                >  2:1521166/1‑76 (MQ=255)                                                    tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                <  1:1521166/76‑1 (MQ=255)                                                    tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCa                                                            <  1:1307587/128‑1 (MQ=255)                                                    tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCa                                                            >  2:1307587/1‑128 (MQ=255)                                                    tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgtg                                                 >  2:1024205/1‑139 (MQ=255)                                                    tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgtg                                                 >  1:1415714/1‑139 (MQ=255)                                                                      gtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                                          <  2:737525/48‑1 (MQ=255)                                                                      gtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                                          >  1:737525/1‑48 (MQ=255)                                                                       tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                <  1:1293940/57‑1 (MQ=255)                                                                       tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                >  2:1293940/1‑57 (MQ=255)                                                                       tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                             >  1:287418/1‑92 (MQ=255)                                                                       tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                             <  2:287418/92‑1 (MQ=255)                                                                       tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAAt                                                       >  1:1660742/1‑114 (MQ=255)                                                                       tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAAt                                                       <  2:1660742/114‑1 (MQ=255)                                                                       tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGcc                              >  2:203666/1‑139 (MQ=255)                                                                            atTGATACCGGGCCAGCACGGCACTGGTAATAATTACTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACt                         >  1:581023/1‑139 (MQ=255)                                                                                 tACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                             >  1:1309634/1‑82 (MQ=255)                                                                                 tACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                             <  2:1309634/82‑1 (MQ=255)                                                                                  aCCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgg                                                                                                                  <  1:630371/44‑1 (MQ=255)                                                                                  aCCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgg                                                                                                                  >  2:630371/1‑44 (MQ=255)                                                                                   ccGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGAt                  <  1:203666/139‑1 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGaaa                                                                                                                      <  1:49776/36‑1 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGaaa                                                                                                                      >  2:49776/1‑36 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                >  2:404275/1‑42 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                <  1:404275/42‑1 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                <  1:1602578/42‑1 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                >  2:1299839/1‑42 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                >  2:1602578/1‑42 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                <  1:1299839/42‑1 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                       <  1:1041388/51‑1 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                       >  2:1041388/1‑51 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                             <  2:1658896/77‑1 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                             >  1:1658896/1‑77 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgt                                                  >  2:377250/1‑104 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgt                                                  <  1:377250/104‑1 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGt                      <  1:450868/132‑1 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGt                      >  2:450868/1‑132 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCt               >  2:186777/1‑139 (MQ=255)                                                                                      ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCt               >  1:490782/1‑139 (MQ=255)                                                                                               ggCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGAt      <  1:1321964/139‑1 (MQ=255)                                                                                                 cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                             >  2:1639521/1‑66 (MQ=255)                                                                                                 cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                             <  1:1639521/66‑1 (MQ=255)                                                                                                 cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAAt                                                       >  2:1317533/1‑88 (MQ=255)                                                                                                 cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAAt                                                       <  1:1317533/88‑1 (MQ=255)                                                                                                 cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCg                    >  2:1127760/1‑123 (MQ=255)                                                                                                 cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCg                    <  1:1127760/123‑1 (MQ=255)                                                                                                 cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGc    >  2:934580/1‑139 (MQ=255)                                                                                                   cTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATAATTCCGGTCACCGCAATCTCCCAGAACTGTTCATTAATCAGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGCgg  >  1:1103086/1‑139 (MQ=255)                                                                                                         |                                                                                                                                     AGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGCGG  >  NC_000913/983519‑983755

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