Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,228,724	A→G	T216A (ACT→GCT)	ycgM →	putative isomerase/hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,228,724	0	A	G	100.0%	140.4 / NA	40	T216A (ACT→GCT)	ycgM	putative isomerase/hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/22);  total (18/22)

AAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAACTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGGCAAAGTAA  >  NC_000913/1228610‑1228858                                                                                                                   |                                                                                                                                       aaGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTgc                                                                                                                >  2:1166667/1‑139 (MQ=255)                      ggTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATAc                                                                                                                       <  1:798638/111‑1 (MQ=255)                      ggTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATAc                                                                                                                       >  2:798638/1‑111 (MQ=255)                       gTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAAc                                                                                          >  1:513723/1‑139 (MQ=255)                          aCGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTg                                                                                       <  2:17949/139‑1 (MQ=255)                                 gtTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCt                                                                                >  2:624782/1‑139 (MQ=255)                                  ttGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTcgcg                                                                                                                                  <  1:1122962/88‑1 (MQ=255)                                  ttGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTcgcg                                                                                                                                  >  2:1122962/1‑88 (MQ=255)                                  ttGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCtg                                                                               <  2:513723/139‑1 (MQ=255)                                    gACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCtgtg                                                                             <  2:1039809/139‑1 (MQ=255)                                    gACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCtgtg                                                                             <  1:493125/139‑1 (MQ=255)                                        ggCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTAGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAg                                                                         <  1:750986/139‑1 (MQ=255)                                          cACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                       >  2:753911/1‑139 (MQ=255)                                          cACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                       <  1:753911/139‑1 (MQ=255)                                          cACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                       <  1:624782/139‑1 (MQ=255)                                          cACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                       <  1:459773/139‑1 (MQ=255)                                               cTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                                  >  2:301968/1‑123 (MQ=255)                                               cTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                                  <  1:301968/123‑1 (MQ=255)                                                      cGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGt                                                           <  1:1166667/139‑1 (MQ=255)                                                       gTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                       >  2:286689/1‑126 (MQ=255)                                                       gTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                       <  1:286689/126‑1 (MQ=255)                                                         cGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCg                                                        >  2:1064012/1‑139 (MQ=255)                                                                      aGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                       >  2:345709/1‑111 (MQ=255)                                                                      aGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                       <  1:345709/111‑1 (MQ=255)                                                                           tgatgaGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGc                                                                                      >  1:822696/1‑91 (MQ=255)                                                                           tgatgaGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGc                                                                                      <  2:822696/91‑1 (MQ=255)                                                                                  cTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCtg                               >  1:296114/1‑139 (MQ=255)                                                                                     aCAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                                  <  1:631312/85‑1 (MQ=255)                                                                                     aCAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                                  >  2:631312/1‑85 (MQ=255)                                                                                         tCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCg                        <  1:790893/139‑1 (MQ=255)                                                                                         tCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCg                        <  1:1129033/139‑1 (MQ=255)                                                                                          cACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAacac                                  <  2:1041903/128‑1 (MQ=255)                                                                                          cACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAacac                                  >  1:1041903/1‑128 (MQ=255)                                                                                            cTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATg                     >  2:1055202/1‑139 (MQ=255)                                                                                               tCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCtgat                             >  2:215363/1‑128 (MQ=255)                                                                                               tCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCtgat                             <  1:215363/128‑1 (MQ=255)                                                                                                     ggCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCt            >  2:980785/1‑139 (MQ=255)                                                                                                             ttGACAGCTCGCGTTTTGTAATACTTTTTGCCGGCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                       <  1:1099551/72‑1 (MQ=255)                                                                                                             ttGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                       >  2:1099551/1‑72 (MQ=255)                                                                                                               gACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGGCAAAGTaa  <  1:1055202/139‑1 (MQ=255)                                                                                                                   |                                                                                                                                       AAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAACTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGGCAAAGTAA  >  NC_000913/1228610‑1228858

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