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%	296.5 / NA	74	G391R (GGC→CGC)	yifB	putative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (36/38);  total (36/38)

TTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTC  >  NC_000913/3948750‑3948841                                                   |                                            ttAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                              >  1:5964144/1‑50 (MQ=255) ttAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                              >  1:16202/1‑50 (MQ=255)   aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gg                                          >  1:3877202/1‑50 (MQ=255)      ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑ggggg                                       <  1:3992517/50‑1 (MQ=255)      ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑ggggg                                       >  1:4290702/1‑50 (MQ=255)      ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑ggggg                                       <  1:5152756/50‑1 (MQ=255)       gTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gggggg                                      <  1:5982760/50‑1 (MQ=255)       gTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gggggg                                      >  1:8015765/1‑50 (MQ=255)        tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑ggggggg                                     >  1:8585802/1‑49 (MQ=255)         ggCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gggggggt                                    >  1:3435918/1‑48 (MQ=255)         ggCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gggggggt                                    >  1:6738550/1‑48 (MQ=255)          gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gggggggt                                    <  1:2931919/49‑3 (MQ=255)          gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gggggggt                                    <  1:1614181/49‑3 (MQ=255)              tgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTaa                                >  1:7107944/1‑50 (MQ=255)                ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAAt                               <  1:8541753/49‑1 (MQ=255)                ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAAt                               <  1:859708/49‑1 (MQ=255)                ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAAt                               <  1:4639575/49‑1 (MQ=255)                ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATg                              >  1:6048649/1‑50 (MQ=255)                ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATg                              <  1:101790/50‑1 (MQ=255)                 tttCTCCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATgg                             >  1:3172395/1‑50 (MQ=255)                 tttCTCCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATgg                             >  1:1294295/1‑50 (MQ=255)                 tttCTCCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATgg                             >  1:5095480/1‑50 (MQ=255)                  ttCTCCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATgg                             >  1:1234570/1‑49 (MQ=255)                  ttCTCCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATgg                             >  1:1229285/1‑49 (MQ=255)                   tctcCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATggg                            <  1:994709/49‑1 (MQ=255)                   tctcCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATggg                            <  1:4522257/49‑1 (MQ=255)                    ctcCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGa                           <  1:2301959/49‑1 (MQ=255)                    ctcCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGa                           <  1:6346203/49‑1 (MQ=255)                    ctcCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAt                          >  1:4359306/1‑50 (MQ=255)                     tcCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAtc                         >  1:7243423/1‑50 (MQ=255)                     tcCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAtc                         >  1:4892040/1‑50 (MQ=255)                     tcCCGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAtc                         >  1:5982092/1‑50 (MQ=255)                      cccGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAtct                        >  1:4306731/1‑50 (MQ=255)                       ccGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAtct                        <  1:4360654/49‑1 (MQ=255)                       ccGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAtc                         <  1:7766533/48‑1 (MQ=255)                       ccGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAt                          <  1:921557/47‑1 (MQ=255)                       ccGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAt                          <  1:3577456/47‑1 (MQ=255)                        cGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAtctc                       >  1:794734/1‑49 (MQ=255)                        cGGCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAtctc                       >  1:4937682/1‑49 (MQ=255)                         ggCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAtctc                       >  1:2228668/1‑48 (MQ=255)                         ggCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAtctc                       >  1:6615973/1‑48 (MQ=255)                         ggCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAtctc                       >  1:252541/1‑48 (MQ=255)                         ggCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGAtctc                       >  1:5641572/1‑48 (MQ=255)                         ggCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCa                     <  1:6360411/50‑1 (MQ=255)                          gCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAg                    >  1:7429399/1‑50 (MQ=255)                          gCACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAg                    <  1:4136679/50‑1 (MQ=255)                           cACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGt                   <  1:4500398/50‑1 (MQ=255)                           cACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGt                   <  1:824900/50‑1 (MQ=255)                           cACTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGt                   <  1:7132181/50‑1 (MQ=255)                            aCTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGt                   >  1:4747956/1‑49 (MQ=255)                            aCTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTg                  <  1:3364150/50‑1 (MQ=255)                             cTACCGTTTTACTCAAAATGCCGGGGGGGGTGGTAATGGGATCTCCAGt                   >  1:3822730/1‑49 (MQ=255)                             cTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGt                   >  1:2566825/1‑48 (MQ=255)                             cTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTga                 <  1:634424/50‑1 (MQ=255)                             cTACCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTga                 <  1:771810/50‑1 (MQ=255)                               aCCGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTga                 >  1:386147/1‑48 (MQ=255)                                ccGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTgagag              <  1:5979978/50‑1 (MQ=255)                                ccGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTgaga               <  1:6409172/49‑1 (MQ=255)                                ccGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTgag                <  1:4450402/48‑1 (MQ=255)                                 cGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTgagaga             <  1:4642088/50‑1 (MQ=255)                                 cGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTga                 >  1:7274098/1‑46 (MQ=255)                                 cGTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTga                 >  1:8053052/1‑46 (MQ=255)                                  gTTTTACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGAt            <  1:2576145/50‑1 (MQ=255)                                   ttttACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATc           <  1:6547412/50‑1 (MQ=255)                                   ttttACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATc           <  1:7340131/50‑1 (MQ=255)                                   ttttACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATc           <  1:3459758/50‑1 (MQ=255)                                    tttACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATc           >  1:5555743/1‑49 (MQ=255)                                     ttACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg          <  1:7745744/49‑1 (MQ=255)                                     ttACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg          <  1:7647864/49‑1 (MQ=255)                                      tACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg          >  1:3669463/1‑48 (MQ=255)                                      tACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg          >  1:8374863/1‑48 (MQ=255)                                      tACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg          >  1:4174372/1‑48 (MQ=255)                                      tACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg          >  1:6845676/1‑48 (MQ=255)                                      tACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg          >  1:4559427/1‑48 (MQ=255)                                      tACTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGa         <  1:6238251/49‑1 (MQ=255)                                        cTCAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg       <  1:3303744/49‑1 (MQ=255)                                        cTCAAAATGCCG‑GGGGGGTGGTAACGGGATCTCCAGTGAGAGATCGAAGc      >  1:6963267/1‑50 (MQ=255)                                          cAAAATGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCgg    >  1:4177578/1‑50 (MQ=255)                                           aaaaTGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt   <  1:7247235/50‑1 (MQ=255)                                           aaaaTGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt   <  1:1372256/50‑1 (MQ=255)                                           aaaaTGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt   <  1:3806240/50‑1 (MQ=255)                                           aaaaTGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt   <  1:8151093/50‑1 (MQ=255)                                            aaaTGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc  <  1:645861/50‑1 (MQ=255)                                            aaaTGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc  <  1:7388909/50‑1 (MQ=255)                                             aaTGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc  >  1:6808371/1‑49 (MQ=255)                                             aaTGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc  >  1:1020195/1‑49 (MQ=255)                                              aTGCCG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc  >  1:7000157/1‑48 (MQ=255)                                                   |                                            TTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTC  >  NC_000913/3948750‑3948841

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