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

TTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGA  >  NC_000913/3948750‑3948843                                                   |                                             ttAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                               >  1:2024086/1‑50 (MQ=255)   aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gg                                            >  1:1821854/1‑50 (MQ=255)      ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggg                                         >  1:359677/1‑50 (MQ=255)                ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                  <  1:1859366/48‑1 (MQ=255)                ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATg                                >  1:1224141/1‑50 (MQ=255)                ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATg                                <  1:1835700/50‑1 (MQ=255)                   tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATggg                              <  1:1160960/49‑1 (MQ=255)                   tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                               <  1:1029878/48‑1 (MQ=255)                       ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                         <  1:1412868/50‑1 (MQ=255)                       ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                         <  1:2344382/50‑1 (MQ=255)                         ggCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGTATCTCCa                       <  1:2464925/50‑1 (MQ=255)                         ggCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                         >  1:1384242/1‑48 (MQ=255)                           cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     <  1:2065729/50‑1 (MQ=255)                           cACTACCGTGTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     <  1:221430/50‑1 (MQ=255)                            aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  1:1703559/1‑49 (MQ=255)                            aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTg                    <  1:377855/50‑1 (MQ=255)                             cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  1:2518658/1‑48 (MQ=255)                              tACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  1:1785777/1‑47 (MQ=255)                               aCCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   >  1:2452279/1‑48 (MQ=255)                                 cGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgagaga               <  1:2571145/50‑1 (MQ=255)                                  gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGAt              <  1:908656/50‑1 (MQ=255)                                   ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   >  1:2575205/1‑44 (MQ=255)                                   ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc             <  1:508362/50‑1 (MQ=255)                                    tttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc             >  1:2156615/1‑49 (MQ=255)                                     ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg            <  1:1769574/49‑1 (MQ=255)                                      tACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg            >  1:1241008/1‑48 (MQ=255)                                       aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa          <  1:2203097/49‑1 (MQ=255)                                          cAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCgg      >  1:1550395/1‑50 (MQ=255)                                           aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCgg      <  1:774562/49‑1 (MQ=255)                                           aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt     <  1:1427840/50‑1 (MQ=255)                                            aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc    <  1:989508/50‑1 (MQ=255)                                             aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc    >  1:2461953/1‑49 (MQ=255)                                              aTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCga  <  1:1454046/50‑1 (MQ=255)                                                   |                                             TTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGA  >  NC_000913/3948750‑3948843

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