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

GGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGCGGAAA  >  NC_000913/983507‑983758                                                                                                                     |                                                                                                                                        ggTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACa                                                                                                                       >  2:171125/1‑135 (MQ=255) ggTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACa                                                                                                                       <  1:171125/135‑1 (MQ=255) ggTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                   <  2:793315/139‑1 (MQ=255)                       gCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                   >  2:905002/1‑117 (MQ=255)                       gCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                   <  1:905002/117‑1 (MQ=255)                                             cgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCa                                                                                                                 <  2:18348/97‑1 (MQ=255)                                             cgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCa                                                                                                                 >  1:18348/1‑97 (MQ=255)                                             cgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATTAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGa                                                                       >  1:652778/1‑139 (MQ=255)                                             cgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCgg                                                                                         >  1:1205358/1‑121 (MQ=255)                                             cgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCgg                                                                                         <  2:1205358/121‑1 (MQ=255)                                             cgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGa                                                                       >  2:336157/1‑139 (MQ=255)                                               cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                                             >  1:881080/1‑83 (MQ=255)                                               cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGta                                                                                                  >  1:1373654/1‑110 (MQ=255)                                               cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGta                                                                                                  <  2:1373654/110‑1 (MQ=255)                                               cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                                <  1:580852/128‑1 (MQ=255)                                               cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                                >  2:580852/1‑128 (MQ=255)                                               cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGAGCCAGCACGGCACGGGTAATAATTGCTTAAc                                                                                                                             <  2:881080/83‑1 (MQ=255)                                                  gtggtgTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGt                                                                  >  1:1551298/1‑139 (MQ=255)                                                               ttatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgt                                                     <  1:336157/139‑1 (MQ=255)                                                                tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCa                                                               >  2:762125/1‑128 (MQ=255)                                                                tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCa                                                               <  1:762125/128‑1 (MQ=255)                                                                tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgtg                                                    >  1:396818/1‑139 (MQ=255)                                                                       gtagtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACaa                                             >  1:303322/1‑139 (MQ=255)                                                                       gtagtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACaa                                             >  2:1687870/1‑139 (MQ=255)                                                                                 tgtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAAc                                                                                                                        >  1:1553220/1‑54 (MQ=255)                                                                                 tgtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAAc                                                                                                                        <  2:1553220/54‑1 (MQ=255)                                                                                   tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                   >  2:1290292/1‑57 (MQ=255)                                                                                   tgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                   <  1:1290292/57‑1 (MQ=255)                                                                                      atatTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATTAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCa                              <  2:652778/139‑1 (MQ=255)                                                                                            aTACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGaaa                                                                                                                         >  2:1534408/1‑42 (MQ=255)                                                                                            aTACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGaaa                                                                                                                         <  1:1534408/42‑1 (MQ=255)                                                                                            aTACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTc                        <  2:396818/139‑1 (MQ=255)                                                                                               ccGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGAt                     <  2:1551298/139‑1 (MQ=255)                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                   >  2:1659366/1‑42 (MQ=255)                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                   <  1:1659366/42‑1 (MQ=255)                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                                <  1:85014/77‑1 (MQ=255)                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                                >  2:85014/1‑77 (MQ=255)                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAAc                                                                     <  1:1479343/88‑1 (MQ=255)                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAAc                                                                     >  2:1479343/1‑88 (MQ=255)                                                                                                  ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCt                  >  1:587455/1‑139 (MQ=255)                                                                                                   gCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                                       >  2:1478565/1‑53 (MQ=255)                                                                                                   gCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                                       <  1:1478565/53‑1 (MQ=255)                                                                                                    ccAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGc                                  >  1:876627/1‑121 (MQ=255)                                                                                                    ccAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGc                                  <  2:876627/121‑1 (MQ=255)                                                                                                      aGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                                >  1:100951/1‑73 (MQ=255)                                                                                                      aGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                                <  2:100951/73‑1 (MQ=255)                                                                                                        cACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTTTCCCAGAACTGTTCATTAATCTGGtgtg                                                  <  2:1005984/101‑1 (MQ=255)                                                                                                        cACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTTTCCCAGAACTGTTCATTAATCTGGtgtg                                                  >  1:1005984/1‑101 (MQ=255)                                                                                                             cACTGGTAATAATTGCTTAACGAAACTGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGAt         <  2:23151/137‑1 (MQ=255)                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACTGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGc       >  2:166632/1‑139 (MQ=255)                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgt                                                     >  2:526280/1‑93 (MQ=255)                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgt                                                     <  1:526280/93‑1 (MQ=255)                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTggcgg             <  2:1342356/133‑1 (MQ=255)                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTggcgg             >  1:1342356/1‑133 (MQ=255)                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGAt         >  1:23151/1‑137 (MQ=255)                                                                                                             cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGc       >  2:872041/1‑139 (MQ=255)                                                                                                                  gtaataATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGCGGaaa  >  1:1205504/1‑139 (MQ=255)                                                                                                                     |                                                                                                                                        GGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGCGGAAA  >  NC_000913/983507‑983758

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