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

TAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTC  >  NC_000913/3948751‑3948841                                                  |                                           tAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                             >  1:31738/1‑49 (MQ=255)  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                             <  1:1148185/48‑1 (MQ=255)  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑g                                           >  1:1340689/1‑49 (MQ=255)   aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑g                                           <  1:2262901/48‑1 (MQ=255)   aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑g                                           <  1:1754747/48‑1 (MQ=255)       tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggg                                      >  1:994219/1‑49 (MQ=255)         gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggggt                                    <  1:1099794/49‑3 (MQ=255)         gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggggt                                    <  1:2420384/49‑3 (MQ=255)         gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggggt                                    >  1:2275689/1‑47 (MQ=255)               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                <  1:2768849/48‑1 (MQ=255)               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTa                                 <  1:546131/47‑1 (MQ=255)               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                               >  1:1450919/1‑49 (MQ=255)               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                               <  1:1732033/49‑1 (MQ=255)                tttCTCCCGGCACTACCGTTTTACTCAAATTGCCGGGGGGGTGGTAATg                              >  1:30101/1‑49 (MQ=38)                 ttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                             >  1:2276681/1‑49 (MQ=255)                  tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                             >  1:1554054/1‑48 (MQ=255)                  tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATg                              <  1:2447237/47‑1 (MQ=255)                   ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATggg                            <  1:2299496/48‑1 (MQ=255)                    tcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                          >  1:2652210/1‑49 (MQ=255)                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                           <  1:1476929/46‑1 (MQ=255)                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtct                        <  1:1743555/49‑1 (MQ=255)                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtc                         <  1:2457771/48‑1 (MQ=255)                       cGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                       <  1:1495691/49‑1 (MQ=255)                        ggCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                       >  1:57092/1‑48 (MQ=255)                        ggCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                       >  1:1981181/1‑48 (MQ=255)                         gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                       >  1:2101514/1‑47 (MQ=255)                         gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCa                     <  1:2670209/49‑1 (MQ=255)                          cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                    >  1:495251/1‑49 (MQ=255)                           aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                   >  1:2436546/1‑49 (MQ=255)                            cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAGGGGAtctcaa                     >  1:1157824/1‑44 (MQ=255)                             tACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                   >  1:566749/1‑47 (MQ=255)                             tACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                   >  1:1737802/1‑47 (MQ=255)                              aCCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                   >  1:262415/1‑46 (MQ=255)                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                 >  1:107981/1‑47 (MQ=255)                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                 <  1:763007/47‑1 (MQ=255)                                cGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                 >  1:2099595/1‑46 (MQ=255)                                 gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgagaga             <  1:435436/49‑1 (MQ=255)                                 gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgagaga             <  1:532219/49‑1 (MQ=255)                                 gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                 >  1:1759958/1‑45 (MQ=255)                                  ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                 >  1:993118/1‑44 (MQ=255)                                  ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGAt            <  1:2217151/49‑1 (MQ=255)                                    ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc           <  1:1795621/48‑1 (MQ=255)                                    ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc           <  1:1516026/48‑1 (MQ=255)                                    ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc           <  1:378390/48‑1 (MQ=255)                                     tACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg          <  1:588310/48‑1 (MQ=255)                                      aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGa         <  1:1598826/48‑1 (MQ=255)                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa        <  1:1377647/48‑1 (MQ=255)                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa        <  1:1872399/48‑1 (MQ=255)                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg       <  1:460937/49‑1 (MQ=255)                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg       >  1:2777541/1‑49 (MQ=255)                                        tCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc      >  1:1720017/1‑49 (MQ=255)                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc      <  1:2029206/47‑1 (MQ=255)                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCgg    >  1:1800970/1‑49 (MQ=255)                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg     <  1:1709475/48‑1 (MQ=255)                                           aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt   >  1:984135/1‑49 (MQ=255)                                           aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt   <  1:1727562/49‑1 (MQ=255)                                            aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc  <  1:2733570/49‑1 (MQ=255)                                                  |                                           TAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTC  >  NC_000913/3948751‑3948841

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