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

CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980902‑3980995                                                   |                                             ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                               <  1:7119313/50‑1 (MQ=255) ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                               <  1:7255485/50‑1 (MQ=255) ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                               >  1:3017404/1‑50 (MQ=255)  cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑g                                             >  1:1155407/1‑50 (MQ=255)  cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑g                                             <  1:6291600/50‑1 (MQ=255)  cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑g                                             <  1:4580645/50‑1 (MQ=255)  cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑g                                             >  1:3910298/1‑50 (MQ=255)   agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gg                                            >  1:2219133/1‑50 (MQ=255)    gagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggg                                           >  1:1745427/1‑50 (MQ=255)      gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggg                                         <  1:6782073/50‑1 (MQ=255)       cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gggggg                                        <  1:3725001/50‑1 (MQ=255)       cTACAGCGTGAGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gggggg                                        <  1:7005547/50‑1 (MQ=255)        tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                       >  1:4598488/1‑49 (MQ=255)        tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                       >  1:1601463/1‑49 (MQ=255)        tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                       >  1:6280267/1‑49 (MQ=255)        tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                       >  1:362180/1‑49 (MQ=255)         aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                       >  1:401538/1‑48 (MQ=255)         aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                       >  1:2576093/1‑48 (MQ=255)         aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg                                       >  1:8085221/1‑48 (MQ=255)          cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gggggggc                                      <  1:7234850/49‑3 (MQ=255)             cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCTCTGGGGGGGCACCa                                   >  1:89725/1‑50 (MQ=255)             cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                   >  1:8607381/1‑50 (MQ=255)             cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                   <  1:4374773/50‑1 (MQ=255)             cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                   >  1:8364640/1‑50 (MQ=255)              gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  >  1:179078/1‑50 (MQ=255)              gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  <  1:115801/50‑1 (MQ=255)              gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  >  1:5514547/1‑50 (MQ=255)              gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  <  1:7605824/50‑1 (MQ=255)               tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                                 <  1:2343691/50‑1 (MQ=255)               tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                                 <  1:4423123/50‑1 (MQ=255)               tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                                 <  1:7881332/50‑1 (MQ=255)                gggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                >  1:8023510/1‑50 (MQ=255)                gggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                >  1:5155583/1‑50 (MQ=255)                gggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                >  1:7738751/1‑50 (MQ=255)                 ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                               >  1:312414/1‑50 (MQ=255)                 ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                               <  1:7724210/50‑1 (MQ=255)                 ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                <  1:48202/49‑1 (MQ=39)                  gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                              <  1:571743/50‑1 (MQ=255)                  gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                              >  1:8088723/1‑50 (MQ=255)                  gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                              >  1:4994580/1‑50 (MQ=255)                  gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                              >  1:5580792/1‑50 (MQ=255)                   cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                             >  1:7100719/1‑50 (MQ=255)                   cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                             >  1:8272731/1‑50 (MQ=255)                    tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                            >  1:490980/1‑50 (MQ=255)                    tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                            >  1:3799174/1‑50 (MQ=255)                    tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                            >  1:7346250/1‑50 (MQ=255)                     ggAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           >  1:4869429/1‑50 (MQ=255)                     ggAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           >  1:8220499/1‑50 (MQ=255)                      gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                            <  1:4746939/48‑1 (MQ=255)                      gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           <  1:507560/49‑1 (MQ=255)                      gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           <  1:559388/49‑1 (MQ=255)                      gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           <  1:221139/49‑1 (MQ=255)                      gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           <  1:2654687/49‑1 (MQ=255)                      gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           <  1:3306889/49‑1 (MQ=255)                      gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                          >  1:7786874/1‑50 (MQ=255)                         gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGcc                       <  1:1317454/50‑1 (MQ=255)                          cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                      <  1:1138561/50‑1 (MQ=255)                          cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                      <  1:8516341/50‑1 (MQ=255)                          cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                      <  1:2375028/50‑1 (MQ=255)                           tCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTg                     >  1:45428/1‑50 (MQ=255)                           tCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTg                     >  1:6882846/1‑50 (MQ=255)                            cGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                    >  1:5427712/1‑50 (MQ=255)                             gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTg                     <  1:664077/48‑1 (MQ=39)                             gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   <  1:8217874/50‑1 (MQ=255)                             gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                    <  1:1920382/49‑1 (MQ=39)                              aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                  <  1:6783048/50‑1 (MQ=255)                              aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                  <  1:2885389/50‑1 (MQ=255)                              aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   >  1:8175172/1‑49 (MQ=39)                               cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                 <  1:7738056/50‑1 (MQ=255)                                 tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg               >  1:3504433/1‑50 (MQ=255)                                 tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg               <  1:4084552/50‑1 (MQ=255)                                   tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg             >  1:3773176/1‑50 (MQ=255)                                   tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg             >  1:3802705/1‑50 (MQ=255)                                   tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg             >  1:4377602/1‑50 (MQ=255)                                     gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg           <  1:6144180/50‑1 (MQ=255)                                     gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg           <  1:7384062/50‑1 (MQ=255)                                     gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg            <  1:1371067/49‑1 (MQ=39)                                     gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg            <  1:4303565/49‑1 (MQ=39)                                      aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc          >  1:7737832/1‑50 (MQ=255)                                       aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         <  1:5518959/50‑1 (MQ=255)                                       aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         <  1:2796162/50‑1 (MQ=255)                                       aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         >  1:6419067/1‑50 (MQ=255)                                          cATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:6169620/1‑49 (MQ=39)                                           aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc      <  1:801912/49‑1 (MQ=39)                                           aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:1558001/1‑48 (MQ=39)                                           aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:995623/1‑48 (MQ=39)                                           aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa     <  1:2999385/50‑1 (MQ=255)                                              gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt   <  1:5760778/49‑1 (MQ=39)                                              gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa  >  1:1898431/1‑50 (MQ=255)                                                   |                                             CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980902‑3980995

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