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%	139.3 / NA	39	T216A (ACT→GCT)	ycgM	putative isomerase/hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/23);  total (16/23)

CGCTGATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAACTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGGCAAAGTA  >  NC_000913/1228590‑1228857                                                                                                                                       |                                                                                                                                      cGCTGATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTcg                                                                                                                                   <  2:560443/139‑1 (MQ=255)   cTGATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTcgcg                                                                                                                                 <  2:1158543/139‑1 (MQ=255)                   gCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACttttt                                                                                                                 >  1:547805/1‑139 (MQ=255)                     aaGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTgc                                                                                                               <  1:37470/139‑1 (MQ=255)                                                 ttGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCAt                                                                                   <  1:720208/139‑1 (MQ=255)                                                   gtgtTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                                 <  2:509793/139‑1 (MQ=255)                                                    tgtTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                                 <  2:1275602/138‑1 (MQ=255)                                                    tgtTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                                 >  1:1275602/1‑138 (MQ=255)                                                    tgtTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGc                                                                                >  1:28212/1‑139 (MQ=255)                                                              cACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                      <  2:1307621/139‑1 (MQ=255)                                                              cACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                      <  2:547805/139‑1 (MQ=255)                                                              cACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                      >  1:357372/1‑139 (MQ=255)                                                                        ggCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTg                                                                                      <  1:259539/113‑1 (MQ=255)                                                                        ggCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTg                                                                                      >  2:259539/1‑113 (MQ=255)                                                                          cGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGt                                                          <  2:357372/139‑1 (MQ=255)                                                                               gCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGtt                                                     >  1:109044/1‑139 (MQ=255)                                                                                    ttacAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTcgcg                                                                                                                                 >  1:710921/4‑58 (MQ=255)                                                                                    ttacAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTcgcg                                                                                                                                 <  2:710921/55‑1 (MQ=255)                                                                                          aGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                      >  2:256297/1‑111 (MQ=255)                                                                                          aGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                      <  1:256297/111‑1 (MQ=255)                                                                                              gTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAgcg                                                                                                   >  1:1001051/1‑78 (MQ=255)                                                                                              gTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAgcg                                                                                                   <  2:1001051/78‑1 (MQ=255)                                                                                              gTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCgg                                      <  2:540601/139‑1 (MQ=255)                                                                                               tgatgaGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                                 <  1:462724/95‑1 (MQ=255)                                                                                               tgatgaGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                                 >  2:462724/1‑95 (MQ=255)                                                                                                         aCAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCaaa                                                                                           <  1:43582/75‑1 (MQ=255)                                                                                                         aCAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCaaa                                                                                           >  2:43582/1‑75 (MQ=255)                                                                                                         aCAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGaa                                    <  1:755736/130‑1 (MQ=255)                                                                                                         aCAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGaa                                    >  2:755736/1‑130 (MQ=255)                                                                                                         aCAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCtgatg                           >  2:566349/1‑139 (MQ=255)                                                                                                         aCAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCtgatg                           <  1:113450/139‑1 (MQ=255)                                                                                                               aCTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGAt                     >  1:1059767/1‑139 (MQ=255)                                                                                                                cTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCgg                                      <  2:399762/121‑1 (MQ=255)                                                                                                                cTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCgg                                      >  1:399762/1‑121 (MQ=255)                                                                                                                           cATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTgg         <  2:28212/139‑1 (MQ=255)                                                                                                                           cATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTgg         <  1:566349/139‑1 (MQ=255)                                                                                                                               cTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                                 >  2:187552/1‑63 (MQ=255)                                                                                                                               cTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                                 <  1:187552/63‑1 (MQ=255)                                                                                                                                  tGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGAGACAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGGCAAAGTa  <  1:234258/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      CGCTGATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAACTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGGCAAAGTA  >  NC_000913/1228590‑1228857

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