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%	51.5 / NA	14	L22R (CTG→CGG)	yifK	putative transporter YifK
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (8/6);  total (8/6)

TACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980909‑3980995                                            |                                             tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                       >  1:913657/1‑49 (MQ=255)  aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                       >  1:2068700/1‑48 (MQ=255)      cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                   >  1:556118/1‑50 (MQ=255)       gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  <  1:1539381/50‑1 (MQ=255)       gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  >  1:891178/1‑50 (MQ=255)              ggAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           >  1:978449/1‑50 (MQ=255)               gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                          <  1:40354/50‑1 (MQ=255)                  gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGcc                       >  1:736074/1‑50 (MQ=255)                   cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                      <  1:241284/50‑1 (MQ=255)                    tCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTg                     >  1:90228/1‑50 (MQ=255)                     cGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                    >  1:1319377/1‑50 (MQ=255)                            tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg             <  1:994906/50‑1 (MQ=255)                                 cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccg        <  1:1706975/50‑1 (MQ=255)                                    aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:1100622/1‑48 (MQ=39)                                    aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa     <  1:2065126/50‑1 (MQ=255)                                       gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa  >  1:1012332/1‑50 (MQ=255)                                            |                                             TACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980909‑3980995

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