Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,948,805	(G)6→7	coding (1165/1521 nt)	yifB ←	putative magnesium chelatase YifB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,948,799	1	.	G	100.0%	186.9 / NA	47	G391R (GGC→CGC)	yifB	putative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (15/32);  total (15/32)

TAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGA  >  NC_000913/3948751‑3948843                                                  |                                             tAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                               <  1:1694767/49‑1 (MQ=255)  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gg                                            >  1:4198165/1‑50 (MQ=255)  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gg                                            >  1:1377002/1‑50 (MQ=255)     ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggg                                         >  1:2599487/1‑50 (MQ=255)     ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggg                                         >  1:1653137/1‑50 (MQ=255)       tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggggg                                       >  1:1053008/1‑49 (MQ=255)        ggCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggggt                                      >  1:2857633/1‑48 (MQ=255)            ctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTa                                   <  1:3553783/50‑1 (MQ=255)             tgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                  >  1:4165456/1‑50 (MQ=255)              gctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                 >  1:1710545/1‑50 (MQ=255)               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                 <  1:2007030/49‑1 (MQ=255)               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATg                                <  1:1267980/50‑1 (MQ=255)               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATg                                <  1:1098517/50‑1 (MQ=255)                  tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                               <  1:4023437/48‑1 (MQ=255)                   ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                             <  1:267059/49‑1 (MQ=255)                   ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                             <  1:1201792/49‑1 (MQ=255)                   ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                            >  1:1703698/1‑50 (MQ=255)                   ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                            >  1:2555671/1‑50 (MQ=255)                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                         <  1:4380456/50‑1 (MQ=255)                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                         <  1:2361491/50‑1 (MQ=255)                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtc                           <  1:3288433/48‑1 (MQ=255)                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtc                           <  1:1232498/48‑1 (MQ=255)                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtc                           <  1:2174758/48‑1 (MQ=255)                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                            <  1:1370526/47‑1 (MQ=255)                         gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                      >  1:2748912/1‑50 (MQ=255)                         gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                      >  1:4246925/1‑50 (MQ=255)                         gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                      >  1:905066/1‑50 (MQ=255)                           aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTg                    <  1:2446382/50‑1 (MQ=255)                            cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   <  1:3516779/50‑1 (MQ=255)                            cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   <  1:1155103/50‑1 (MQ=255)                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgaga                 <  1:3787627/49‑1 (MQ=255)                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgaga                 <  1:1285429/49‑1 (MQ=255)                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgag                  <  1:2985819/48‑1 (MQ=255)                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgag                  <  1:4365938/48‑1 (MQ=255)                                cGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgagaga               <  1:3150801/50‑1 (MQ=255)                                cGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgagaga               <  1:104862/50‑1 (MQ=255)                                cGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   >  1:1947707/1‑46 (MQ=255)                                 gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGAt              <  1:1490667/50‑1 (MQ=255)                                  ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc             <  1:2273938/50‑1 (MQ=255)                                   tttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc             >  1:690290/1‑49 (MQ=255)                                   tttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc             >  1:2714170/1‑49 (MQ=255)                                    ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg            >  1:3624358/1‑49 (MQ=255)                                     tACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg            >  1:4204319/1‑48 (MQ=255)                                     tACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGa           <  1:3351700/49‑1 (MQ=255)                                      aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa          <  1:4108043/49‑1 (MQ=255)                                      aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa          <  1:4208353/49‑1 (MQ=255)                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc        <  1:1710322/50‑1 (MQ=255)                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg       <  1:466105/48‑1 (MQ=255)                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg       <  1:265357/48‑1 (MQ=255)                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg       <  1:1306751/48‑1 (MQ=255)                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt     >  1:1868302/1‑50 (MQ=255)                                            aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc    >  1:1958047/1‑49 (MQ=255)                                             aTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc    >  1:1343298/1‑48 (MQ=255)                                             aTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCga  <  1:341989/50‑1 (MQ=255)                                                  |                                             TAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGA  >  NC_000913/3948751‑3948843

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