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

CAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAG  >  NC_000913/3980903‑3980993                                                  |                                           cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑g                                           <  1:213100/50‑1 (MQ=255)  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gg                                          >  1:2409779/1‑50 (MQ=255)  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gg                                          >  1:1924838/1‑50 (MQ=255)       tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                     >  1:1999752/1‑49 (MQ=255)             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                <  1:1406134/50‑1 (MQ=255)              tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                               >  1:168720/1‑50 (MQ=255)               gggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                              >  1:1570190/1‑50 (MQ=255)                        gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCCTTGGCGTCGGcc                     >  1:2093076/1‑50 (MQ=255)                        gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGcc                     >  1:1414195/1‑50 (MQ=255)                            gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                 <  1:1167234/50‑1 (MQ=255)                             aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                <  1:420092/50‑1 (MQ=255)                                tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg             >  1:2198354/1‑50 (MQ=255)                                 atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg            >  1:1094714/1‑50 (MQ=255)                                 atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg            >  1:1433581/1‑50 (MQ=255)                                  tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg           <  1:2625095/50‑1 (MQ=255)                                    gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg         >  1:1900508/1‑50 (MQ=255)                                     aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc        <  1:1028605/50‑1 (MQ=255)                                      aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc       >  1:2122264/1‑50 (MQ=255)                                       cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccg      >  1:622591/1‑50 (MQ=255)                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc     >  1:427108/1‑48 (MQ=39)                                            cGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg  >  1:2457077/1‑49 (MQ=39)                                                  |                                           CAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAG  >  NC_000913/3980903‑3980993

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