Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,948,805	(G)6→7	coding (1165/1521 nt)	yifB ←	putative magnesium chelatase YifB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,948,799	1	.	G	100.0%	24.2 / NA	8	G391R (GGC→CGC)	yifB	putative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (1/7);  total (1/7)

AACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGT  >  NC_000913/3948752‑3948840                                                 |                                          aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑g                                          >  1:1080068/1‑49 (MQ=255)  aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑g                                          <  1:1083991/48‑1 (MQ=255)   cGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggg                                        <  1:1870082/49‑1 (MQ=255)    ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggg                                       <  1:1376904/49‑1 (MQ=255)    ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggg                                       >  1:324158/1‑49 (MQ=255)        gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggggt                                   <  1:305491/49‑3 (MQ=255)               tttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATg                             >  1:1757426/1‑49 (MQ=255)                          aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                  <  1:1744869/49‑1 (MQ=255)                          aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                  <  1:511616/49‑1 (MQ=255)                           cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTg                 <  1:1519999/49‑1 (MQ=255)                           cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTg                 <  1:659451/49‑1 (MQ=255)                           cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTg                 <  1:943893/49‑1 (MQ=255)                                   ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc          <  1:1536960/48‑1 (MQ=255)                                          aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt  <  1:343303/49‑1 (MQ=255)                                                 |                                          AACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGT  >  NC_000913/3948752‑3948840

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