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

CAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAG  >  NC_000913/3980903‑3980993                                                  |                                           cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                             >  1:1037250/1‑49 (MQ=255) cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                             <  1:2054662/49‑1 (MQ=255)  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑g                                           >  1:164056/1‑49 (MQ=255)   gagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gg                                          >  1:1317365/1‑49 (MQ=255)             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                 >  1:294821/1‑49 (MQ=39)                 gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                             <  1:204159/49‑1 (MQ=39)                  cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                            >  1:697337/1‑49 (MQ=255)                         cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGcc                     >  1:721905/1‑49 (MQ=39)                          tCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                    >  1:1405179/1‑49 (MQ=39)                          tCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                    >  1:1592210/1‑49 (MQ=39)                             aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                 <  1:1837685/49‑1 (MQ=39)                             aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                 <  1:659899/49‑1 (MQ=39)                              cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                <  1:588957/49‑1 (MQ=39)                                       cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc       <  1:1673255/49‑1 (MQ=39)                                       cTCAACGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc       <  1:657064/49‑1 (MQ=38)                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc    <  1:1728678/49‑1 (MQ=39)                                             gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg  >  1:1255440/1‑48 (MQ=39)                                                  |                                           CAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAG  >  NC_000913/3980903‑3980993

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