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

GAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980905‑3980995                                                |                                             gagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggg                                           >  1:1494560/1‑50 (MQ=255)  agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gggg                                          >  1:824101/1‑50 (MQ=255)  agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gggg                                          >  1:3860372/1‑50 (MQ=255)    cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gggggg                                        >  1:414209/1‑50 (MQ=255)     tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                       >  1:2670260/1‑49 (MQ=255)           gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  >  1:143499/1‑50 (MQ=255)            tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                                 <  1:2709904/50‑1 (MQ=255)             gggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                >  1:2692648/1‑50 (MQ=255)              ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                               >  1:1334613/1‑50 (MQ=255)                cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                             >  1:91259/1‑50 (MQ=255)                cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                             >  1:3170460/1‑50 (MQ=255)                cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                             >  1:1160171/1‑50 (MQ=255)                   gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                            <  1:2021502/48‑1 (MQ=255)                   gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                            <  1:1619488/48‑1 (MQ=255)                   gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                          >  1:1423550/1‑50 (MQ=255)                    aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                         >  1:1476739/1‑50 (MQ=255)                    aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                         >  1:3899322/1‑50 (MQ=255)                      gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGcc                       >  1:2161783/1‑50 (MQ=255)                       cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                      >  1:862753/1‑50 (MQ=255)                       cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                      <  1:3233551/50‑1 (MQ=255)                          gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTg                     <  1:2839491/48‑1 (MQ=39)                           aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                  <  1:3717830/50‑1 (MQ=255)                           aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   >  1:1334480/1‑49 (MQ=39)                            cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                 <  1:1568956/50‑1 (MQ=255)                             atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                <  1:2569391/50‑1 (MQ=255)                               atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg              >  1:3126906/1‑50 (MQ=255)                               atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg              >  1:1082720/1‑50 (MQ=255)                               atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg              >  1:1587509/1‑50 (MQ=255)                                tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg             >  1:2186908/1‑50 (MQ=255)                                   aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc          >  1:3757603/1‑50 (MQ=255)                                   aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc          >  1:2879989/1‑50 (MQ=255)                                       cATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:1475861/1‑49 (MQ=39)                                        aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc      <  1:1810821/49‑1 (MQ=39)                                        aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa     <  1:3186986/50‑1 (MQ=255)                                           gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt   <  1:3480671/49‑1 (MQ=39)                                           gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa  >  1:2466585/1‑50 (MQ=255)                                                |                                             GAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980905‑3980995

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: ‑