Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,980,957	(G)6→7	coding (71/1386 nt)	yifK →	putative transporter YifK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,980,951	1	.	G	100.0%	45.4 / NA	13	L22R (CTG→CGG)	yifK	putative transporter YifK
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (7/6);  total (7/6)

CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980902‑3980995                                                   |                                             ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                               >  1:127317/1‑50 (MQ=255)  cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑g                                             <  1:742691/50‑1 (MQ=255)   agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gg                                            <  1:1217917/50‑1 (MQ=255)                 ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                <  1:1099626/49‑1 (MQ=39)                  gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                              <  1:1978880/50‑1 (MQ=255)                  gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                              <  1:979861/50‑1 (MQ=255)                      gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                            <  1:1155135/48‑1 (MQ=255)                        aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGc                        >  1:1399508/1‑50 (MQ=255)                        aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGc                        >  1:1969513/1‑50 (MQ=255)                             gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   <  1:1424669/50‑1 (MQ=255)                              aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                  <  1:1580002/50‑1 (MQ=255)                              aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   >  1:1094229/1‑49 (MQ=39)                                         tCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:1591917/1‑50 (MQ=255)                                         tCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:191465/1‑50 (MQ=255)                                             cGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg    >  1:2292332/1‑49 (MQ=39)                                              gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa  >  1:702967/1‑50 (MQ=255)                                                   |                                             CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980902‑3980995

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