Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,949,783	(T)7→6	coding (187/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,949,777	0	T	.	100.0%	127.5 / NA	30	coding (193/1521 nt)	yifB	putative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (20/10);  total (20/10)

GCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTG  >  NC_000913/3949702‑3949851                                                                            |                                                                           gCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGAtttttt                                                                       >  1:772838/1‑81 (MQ=255)   aaCGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGAttttttcg                                                                     >  2:772835/1‑79 (MQ=255)        aaTGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGa                                                               <  1:772807/81‑1 (MQ=255)            gCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGatat                                                            <  1:772808/80‑1 (MQ=255)             cGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATAtt                                                           >  2:772845/1‑80 (MQ=255)                tAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCata                                                       >  1:772850/1‑81 (MQ=255)                   gTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATcc                                                    >  2:772840/1‑81 (MQ=255)                    tAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCg                                                   >  2:772856/1‑81 (MQ=255)                     aaaTCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGc                                                  >  1:772854/1‑81 (MQ=255)                      aaTCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGc                                                  >  1:772855/1‑80 (MQ=255)                      aaTCATATCGTCCCCCTTCTTTTGGCAGATAAGCTGGAGCCAGATTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCt                                                 <  1:772809/81‑1 (MQ=255)                               gTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAAt                                        >  2:772848/1‑81 (MQ=255)                               gTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAAt                                        >  2:772857/1‑81 (MQ=255)                                tCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATg                                       <  1:772810/81‑1 (MQ=255)                                  ccccTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGgc                                     <  2:772811/81‑1 (MQ=255)                                   cccTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGgcg                                    <  2:772812/81‑1 (MQ=255)                                     cttcttTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCt                                  >  2:772858/1‑81 (MQ=255)                                        cttTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTgc                                <  1:772813/80‑1 (MQ=255)                                          tttGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCa                             >  1:772846/1‑81 (MQ=255)                                              gCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCAcgcg                         <  1:772814/81‑1 (MQ=255)                                               cAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGa                        <  1:772815/81‑1 (MQ=255)                                                aGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGAt                       >  2:772859/1‑81 (MQ=255)                                                 gATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATc                      >  2:772853/1‑81 (MQ=255)                                                     aGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATcgc                    >  1:772847/1‑79 (MQ=255)                                                        tGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGct               >  1:772852/1‑81 (MQ=255)                                                          gAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGcttc             >  1:772863/1‑81 (MQ=255)                                                           aGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGcttct            <  2:772816/81‑1 (MQ=255)                                                               aGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTa        >  2:772861/1‑81 (MQ=255)                                                                 gTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTAcc      >  1:772868/1‑81 (MQ=255)                                                                    gATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCgtt   >  1:772866/1‑81 (MQ=255)                                                                    gATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCgtt   >  1:772849/1‑81 (MQ=255)                                                                     aTGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCgttg  >  1:772837/1‑81 (MQ=255)                                                                            |                                                                           GCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTG  >  NC_000913/3949702‑3949851

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