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

CAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAG  >  NC_000913/3980903‑3980993                                                  |                                           cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                             <  1:277635/49‑1 (MQ=255) cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                             <  1:1639496/49‑1 (MQ=255) cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                             >  1:1510066/1‑49 (MQ=255)  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑g                                           >  1:1055537/1‑49 (MQ=255)  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑g                                           >  1:1704901/1‑49 (MQ=255)  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑g                                           >  1:393858/1‑49 (MQ=255)   gagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gg                                          >  1:801693/1‑49 (MQ=255)    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggg                                         <  1:875507/49‑1 (MQ=255)    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggg                                         >  1:1613107/1‑49 (MQ=255)     gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gggg                                        >  1:473154/1‑49 (MQ=255)       tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gggggg                                      >  1:2763026/1‑49 (MQ=255)        aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                     >  1:1176891/1‑48 (MQ=255)         cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gggggggc                                    <  1:1602837/49‑3 (MQ=255)         cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gggggggc                                    <  1:904098/49‑3 (MQ=255)                ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                               <  1:150951/48‑1 (MQ=39)                ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                              <  1:878239/49‑1 (MQ=39)                ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                              >  1:2059299/1‑49 (MQ=39)                ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                              <  1:2440917/49‑1 (MQ=39)                 gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                             >  1:935060/1‑49 (MQ=39)                 gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                             <  1:2726015/49‑1 (MQ=39)                  cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                            >  1:2121839/1‑49 (MQ=255)                  cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                            <  1:193383/49‑1 (MQ=255)                   tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                           >  1:1041010/1‑49 (MQ=255)                      aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                        >  1:2250481/1‑49 (MQ=39)                      aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                        <  1:2623024/49‑1 (MQ=39)                       aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                       >  1:1194713/1‑49 (MQ=39)                        gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGc                      >  1:2093524/1‑49 (MQ=39)                         cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGcc                     <  1:2300088/49‑1 (MQ=39)                          tCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                    >  1:2205431/1‑49 (MQ=39)                            gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                  <  1:2407004/49‑1 (MQ=39)                            gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                  >  1:1535071/1‑49 (MQ=39)                             aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                 >  1:140564/1‑49 (MQ=39)                                 atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg             >  1:598065/1‑49 (MQ=39)                                 atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg             >  1:1677654/1‑49 (MQ=39)                                  tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg             <  1:503130/48‑1 (MQ=39)                                    gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg          <  1:1569774/49‑1 (MQ=39)                                    gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg          <  1:835418/49‑1 (MQ=39)                                    gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg          >  1:2591269/1‑49 (MQ=39)                                     aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg         <  1:805465/49‑1 (MQ=39)                                     aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg         >  1:1068452/1‑49 (MQ=39)                                       cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc       >  1:1770521/1‑49 (MQ=39)                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc    <  1:420391/49‑1 (MQ=39)                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc     >  1:439584/1‑48 (MQ=39)                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc     >  1:1380477/1‑48 (MQ=39)                                           tCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa   <  1:868630/49‑1 (MQ=39)                                            cGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg  >  1:2725441/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 < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑