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

ATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGG  >  NC_000913/983501‑983744                                                                                                                           |                                                                                                                          atTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGaaa                                                                                                           <  2:1212036/139‑1 (MQ=255)  ttGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                               >  1:495213/1‑134 (MQ=255)  ttGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                               <  2:495213/134‑1 (MQ=255)          aTCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAg                                                                                                  >  1:861058/1‑139 (MQ=255)                                  gCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGt                                                                          >  2:631954/1‑139 (MQ=255)                                    cgcgTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                     <  1:446094/110‑1 (MQ=255)                                    cgcgTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                     >  2:446094/1‑110 (MQ=255)                                                    gcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                        <  2:861058/139‑1 (MQ=255)                                                     cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTAt                                                                                        <  1:353690/106‑1 (MQ=255)                                                     cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTAt                                                                                        >  2:353690/1‑106 (MQ=255)                                                     cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATtct                                                               >  2:819969/1‑131 (MQ=255)                                                     cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATtct                                                               <  1:819969/131‑1 (MQ=255)                                                           gtgtAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGTTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                  <  2:406226/122‑1 (MQ=255)                                                           gtgtAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGTTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                  >  1:406226/1‑122 (MQ=255)                                                                ggATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAAt                                            <  1:631954/139‑1 (MQ=255)                                                                  atatTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                        <  2:44893/125‑1 (MQ=255)                                                                  atatTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                        >  1:44893/1‑125 (MQ=255)                                                                      tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                               >  2:1051495/1‑66 (MQ=255)                                                                      tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                               <  1:1051495/66‑1 (MQ=255)                                                                      tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTcc                                                             >  2:591077/1‑116 (MQ=255)                                                                      tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTcc                                                             <  1:591077/116‑1 (MQ=255)                                                                      tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgtg                                      >  1:52649/1‑139 (MQ=255)                                                                      tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGGAATAATTGCTTAACGAAAc                                                                                                          <  1:955672/71‑1 (MQ=255)                                                                      tatCCACGTAGTAACATTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAAc                                                                                                          >  2:955672/1‑71 (MQ=255)                                                                        tCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                  >  2:221261/1‑109 (MQ=255)                                                                        tCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                  <  1:221261/109‑1 (MQ=255)                                                                        tCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGtgtg                                    >  1:213124/1‑139 (MQ=255)                                                                             gtagtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGgtaata                                                                                                                         >  2:785431/1‑49 (MQ=255)                                                                             gtagtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAggg                                                                                                      <  1:1150535/68‑1 (MQ=255)                                                                             gtagtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAggg                                                                                                      >  2:1150535/1‑68 (MQ=255)                                                                             gtagtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACaa                               >  1:1022313/1‑139 (MQ=255)                                                                               agtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAAt                                            <  2:853072/124‑1 (MQ=255)                                                                               agtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCCGTCACCGCATTCTCCCAGAACTGTTCATTAAt                                            >  1:853072/1‑124 (MQ=255)                                                                                 taACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCa                                                 >  1:1211077/1‑117 (MQ=255)                                                                                 taACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCa                                                 <  2:1211077/117‑1 (MQ=255)                                                                                         tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                  >  1:920796/1‑92 (MQ=255)                                                                                         tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                  <  2:920796/92‑1 (MQ=255)                                                                                         tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGcc                   >  1:510914/1‑139 (MQ=255)                                                                                         tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGcc                   >  1:721267/1‑139 (MQ=255)                                                                                              atTGATACCGGGCCAGCACGGCACTGGTAATACTTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACt              >  2:851571/1‑139 (MQ=255)                                                                                                  aTACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                  <  2:1114747/83‑1 (MQ=255)                                                                                                  aTACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                  >  1:1114747/1‑83 (MQ=255)                                                                                                  aTACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCCTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTc          <  2:1022313/139‑1 (MQ=255)                                                                                                   tACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                     <  1:1006288/47‑1 (MQ=255)                                                                                                   tACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                     >  2:1006288/1‑47 (MQ=255)                                                                                                        ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                     <  1:1158640/42‑1 (MQ=255)                                                                                                        ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                     >  2:1158640/1‑42 (MQ=255)                                                                                                        ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGt                                                                                                 <  2:77939/46‑1 (MQ=255)                                                                                                        ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGt                                                                                                 >  1:77939/1‑46 (MQ=255)                                                                                                        ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                  <  2:495402/77‑1 (MQ=255)                                                                                                        ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                  >  1:495402/1‑77 (MQ=255)                                                                                                        ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCt    >  1:1391124/1‑139 (MQ=255)                                                                                                        ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCt    >  2:359968/1‑139 (MQ=255)                                                                                                        agCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                         >  1:1288573/2‑54 (MQ=255)                                                                                                        agCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                         <  2:1288573/53‑1 (MQ=255)                                                                                                          ccAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                     <  2:543506/40‑1 (MQ=255)                                                                                                          ccAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                     >  1:543506/1‑40 (MQ=255)                                                                                                          ccAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTgg  >  1:1288959/1‑139 (MQ=255)                                                                                                                           |                                                                                                                          ATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGG  >  NC_000913/983501‑983744

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