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

GATGGATTTCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCG  >  NC_000913/983488‑983746                                                                                                                                        |                                                                                                                            gATGGATTTCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGtaataa                                                                                                                          <  1:162021/139‑1 (MQ=255)                    ggTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       <  1:46452/139‑1 (MQ=255)                       aTCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                                 <  2:558638/126‑1 (MQ=255)                       aTCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                                 >  1:558638/1‑126 (MQ=255)                           ccttAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATg                                                                                                <  1:687887/137‑1 (MQ=255)                             ttAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                              <  1:492383/139‑1 (MQ=255)                                                       aTCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                    <  1:831026/139‑1 (MQ=255)                                                                cgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                    <  1:567487/130‑1 (MQ=255)                                                                cgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                    >  2:567487/1‑130 (MQ=255)                                                                 gcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                          <  1:630948/139‑1 (MQ=255)                                                                 cagggtggtGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCAGGCCAGCACGGCACTGGTAATAAGTGCTTAACGAAACAGAGGCAGTATGATCTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                          <  1:534133/136‑1 (MQ=255)                                                                  cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       <  2:527893/93‑1 (MQ=255)                                                                  cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       >  1:527893/1‑93 (MQ=255)                                                                  cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAggg                                                                                                        >  2:716634/1‑92 (MQ=255)                                                                  cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAggg                                                                                                        <  1:716634/92‑1 (MQ=255)                                                                      tggtgTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGtt                                                     <  1:188335/139‑1 (MQ=255)                                                                               atatTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       <  1:476774/80‑1 (MQ=255)                                                                               atatTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       >  2:476774/1‑80 (MQ=255)                                                                               atatTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                          <  1:181218/125‑1 (MQ=255)                                                                               atatTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                          >  2:181218/1‑125 (MQ=255)                                                                               atatTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCt                                            >  1:296044/1‑139 (MQ=255)                                                                               atatTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATCCCGGTCACCGTATTCTCCCAGAACTGTTCATTaactg                                            >  2:534133/1‑136 (MQ=255)                                                                                   tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTaa                                                                                                                  <  2:41592/65‑1 (MQ=255)                                                                                   tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTaa                                                                                                                  >  1:41592/1‑65 (MQ=255)                                                                                   tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       >  1:592998/1‑76 (MQ=255)                                                                                   tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       <  2:592998/76‑1 (MQ=255)                                                                                   tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTc                                                                           <  1:751356/104‑1 (MQ=255)                                                                                   tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTc                                                                           >  2:751356/1‑104 (MQ=255)                                                                                   tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgtg                                        >  1:789809/1‑139 (MQ=255)                                                                                   tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgtg                                        >  2:261557/1‑139 (MQ=255)                                                                                        aCGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATg                                                                                                <  1:829755/78‑1 (MQ=255)                                                                                        aCGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATg                                                                                                >  2:829755/1‑78 (MQ=255)                                                                                            agtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                                 >  1:812258/1‑57 (MQ=255)                                                                                            agtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                                 <  2:812258/57‑1 (MQ=255)                                                                                               aaCTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCg                            >  2:137348/1‑139 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       <  2:226954/57‑1 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       >  1:226954/1‑57 (MQ=255)                                                                                                      tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGcc                     >  1:476593/1‑139 (MQ=255)                                                                                                           atTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGt                           >  2:375042/1‑128 (MQ=255)                                                                                                           atTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGt                           <  1:375042/128‑1 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       >  2:389525/1‑42 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       <  1:389525/42‑1 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       >  2:743760/1‑42 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       >  1:430048/1‑42 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       <  1:743760/42‑1 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                       <  2:430048/42‑1 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                          <  1:793961/87‑1 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGaa                                                          >  2:793961/1‑87 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCa                  <  1:834567/127‑1 (MQ=255)                                                                                                                     ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCa                  >  2:834567/1‑127 (MQ=255)                                                                                                                         aGCACGGCACTGGTAATAATTGCTTAACGAAATAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTggcg  >  2:131991/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                            GATGGATTTCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCG  >  NC_000913/983488‑983746

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