Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	817,246	G→C	R68P (CGC→CCC)	moaA →	GTP 3',8'‑cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	817,246	0	G	C	100.0%	29.2 / NA	12	R68P (CGC→CCC)	moaA	GTP 3',8'‑cyclase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (7/5);  total (7/5)

CTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGCCGTCTTTACGCCGCGACTTTACCGATATCATCGCCGCTGTGCGGGAAAACGACGCTATCCGCCAGATTGCGGTCACAACCAATGGTTACCGTCTGGAACGCGATGTGGCGAGCTGGCGCGATGCGGGACTTACTGGCATTAACGTCAGTGTCGACAGTCTGGACGCCCGCCAGTTTCACGCTATTACCGGGCAGGATAAATT  >  NC_000913/817118‑817465                                                                                                                                 |                                                                                                                                                                                                                            cTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGCTGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCCCCTGACAGGAGGAGAGCCGTCTTTACGCCGCGACTTTCCCGATATCATCGCCGCTGTGCGGGAAAacgac                                                                                                                                                       >  1:172736/1‑199 (MQ=255)            ccTACTGCCTGCCGGATGTCTACAAACCGAGCGGCGTCACCAATAAAGGCGTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCCCCTGACAGGAGGAGAGGCGTCTTTACGCCGCTACTTTACCGATATCATCTCCGCTGTGCGGGAAAACGACGCTATCCGCCAGATTGCGGTCACAATCAATGGTTACCGTCTGGAACGCGATGTGGCGAGCTgg                                                                                        <  1:228198/251‑1 (MQ=255)            ccTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAAGTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCCCCTGACAGGAGGAGAGCCGTCTTGATGCCGCGACTTTACCGATATCACCGCCGCTGTGCGGGAAAACGACGCTATCCGCCAGATGGCGGTCACACCCAATGGTTACCGTCTGGAACGCGATGTGGCGAGCTgg                                                                                        <  1:325111/251‑1 (MQ=255)                  ncctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCCCCTGACAGGAGGAGAGCCGTCTTTACGCCGCGACTTTACCGATATCATCGCCGCTGTGCGGGAAAACGACGCTATCCGCCAGATTGCGGTCACAACCAATGGTTACCGTCTGGAACGCGATGTGGCGAGCTGGCGCGAt                                                                                  >  1:104766/2‑251 (MQ=255)                         ggATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTGCCGTCGCTGAAATTCGCCGGGTTACGCACGCCTTCGCCAGACTGGGCACCGAAGAAGTGACCCTGAGAGGAGGAGAACCGTCTTTACGCCGCGACTTTACCGATATCATCGCCGCTGTGCGGGAAAACGACGCtcgcaga                                                                                                                                              >  1:7655/1‑178 (MQ=255)                             ggCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCCCCTGACAGGAGGAGAGCCGTCTTTACGCCGCGACTTTACCGATATCATCGCCGCTGTGCGGGAAAACGACGCTATCCGCCAGATTGCGGTCACAACCAATGGTTACCGTCTGGAACGCGATGTGGCGAGCTGGCGCGATGCGGGa                                                                            >  1:133803/1‑246 (MQ=255)                                        gAGCGGCGTCACCAATAAAGGCTTTCTTGCCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCCCCTGACAGGAGGAGAGCCGTCTTTACGCCGCGACTTTACCGTTATCATCGCCGCTGTGCGGGAAAac                                                                                                                                                          <  1:21258/157‑1 (MQ=255)                                        gAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGNNCCTGAcn                                                                                                                                                                                                                      >  2:21258/1‑96 (MQ=255)                                                                                                  gcCTTCGCCAGACTGGGCACCGAAAAAGGGCCCCTGACAGGAGGAGAGCCGTCTTTACGCCGCGGCTTTACCGATATCATCGCCGCTGTGCGGGAAAACGACGCTATCCGCCAGATTGCGGTCACAACCAATGGTTACCGTCTGGAACGCGATGTGGCGAGCTGGCGCGATGCGGGACTTACTGGCATTAACGTCAGTGTCGACAGTCTGGACGCCCGCCAGTTTCACGCTATTACCGGGCAGGATAAAtt  >  1:201089/1‑251 (MQ=255)                                                                                                                       gAAAAAGTGCCCCTGACAGGAGGAGAGCCGTCTTTACGCCGCGACTTTACCGATATCATCGCCGCTGTGCGGGAAAac                                                                                                                                                          <  2:18252/78‑1 (MQ=255)                                                                                                                       gAAAAAGGGCCCCTGACAGGAGGAGAGCCGTCTTTACGCCGCGACTTTACCGATATCATCGCCGCTGTGCGGGAAAac                                                                                                                                                          >  1:18252/1‑78 (MQ=255)                                                                                                                            aGTGCCCCTGACAGGAGGAGAGCCGTCTTTACGCCGCGACTTTACCGATATCATCGCCGCTGTGCg                                                                                                                                                                 <  1:7087/66‑1 (MQ=255)                                                                                                                            aGTGCCCCTGACAGGAGGAGAGCCGTCTTTACGCCGCGACTTTACCGATATCATCGCCGCTGTGCg                                                                                                                                                                 >  2:7087/1‑66 (MQ=255)                                                                                                                                 |                                                                                                                                                                                                                            CTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGCCGTCTTTACGCCGCGACTTTACCGATATCATCGCCGCTGTGCGGGAAAACGACGCTATCCGCCAGATTGCGGTCACAACCAATGGTTACCGTCTGGAACGCGATGTGGCGAGCTGGCGCGATGCGGGACTTACTGGCATTAACGTCAGTGTCGACAGTCTGGACGCCCGCCAGTTTCACGCTATTACCGGGCAGGATAAATT  >  NC_000913/817118‑817465

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