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

TACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980909‑3980995                                            |                                             tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                       >  1:2301088/1‑49 (MQ=255)   cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                       <  1:1749948/48‑2 (MQ=255)      cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                   >  1:703082/1‑50 (MQ=255)      cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                   <  1:4432162/50‑1 (MQ=255)      cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                   <  1:3568635/50‑1 (MQ=255)       gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  >  1:1656811/1‑50 (MQ=255)        tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                                 <  1:523250/50‑1 (MQ=255)        tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                                 <  1:2202462/50‑1 (MQ=255)          ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                <  1:1463239/49‑1 (MQ=39)              ggAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           >  1:2957559/1‑50 (MQ=255)              ggAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           >  1:2297785/1‑50 (MQ=255)               gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                            <  1:3383546/48‑1 (MQ=255)               gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                          <  1:1387849/50‑1 (MQ=255)               gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                          >  1:3612507/1‑50 (MQ=255)               gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                          >  1:4047756/1‑50 (MQ=255)                 aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGc                        <  1:4136506/50‑1 (MQ=255)                 aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGc                        <  1:2011016/50‑1 (MQ=255)                   cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                      <  1:767044/50‑1 (MQ=255)                      gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTg                     <  1:681457/48‑1 (MQ=39)                      gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   <  1:2615796/50‑1 (MQ=255)                      gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   >  1:3939572/1‑50 (MQ=255)                       aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                  <  1:2910470/50‑1 (MQ=255)                        cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                 <  1:2035970/50‑1 (MQ=255)                              gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg            <  1:3446393/49‑1 (MQ=39)                                aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         >  1:3270890/1‑50 (MQ=255)                                aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         <  1:2834523/50‑1 (MQ=255)                                aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         >  1:1850471/1‑50 (MQ=255)                                   cATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:2573173/1‑49 (MQ=39)                                    aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa     <  1:2456330/50‑1 (MQ=255)                                    aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa     <  1:1384605/50‑1 (MQ=255)                                       gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt   <  1:2046193/49‑1 (MQ=39)                                       gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa  >  1:2153964/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 < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑