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%	120.6 / NA	27	coding (193/1521 nt)	yifB	putative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (18/9);  total (18/9)

CAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTT  >  NC_000913/3949703‑3949850                                                                           |                                                                          cAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGAttttttc                                                                     >  1:722614/1‑80 (MQ=255)  aaCGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGAttttttcg                                                                    >  2:722606/1‑79 (MQ=255)   aCGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGAttttttcgc                                                                   >  2:722620/1‑78 (MQ=255)        aTGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGa                                                              >  2:722604/1‑80 (MQ=255)           gCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATAtt                                                          <  1:722581/81‑1 (MQ=255)           gCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATAtt                                                          <  1:722580/81‑1 (MQ=255)             gATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCa                                                        >  2:722617/1‑81 (MQ=255)              aTAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCat                                                       <  2:722582/81‑1 (MQ=255)                aGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCatat                                                     >  2:722616/1‑81 (MQ=255)                 ggTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATc                                                    >  2:722613/1‑81 (MQ=255)                           atCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGa                                           >  1:722643/1‑80 (MQ=255)                            tCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATa                                         <  1:722583/81‑1 (MQ=255)                                 ccccTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATgg                                     <  1:722584/80‑1 (MQ=255)                                   ccTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGgcgc                                  >  1:722627/1‑81 (MQ=255)                                     ttcttTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTg                                >  2:722615/1‑81 (MQ=255)                                      tcttTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTgc                               >  2:722625/1‑81 (MQ=255)                                          ttGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCAc                           >  1:722626/1‑81 (MQ=255)                                          ttGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCAc                           >  2:722619/1‑81 (MQ=255)                                           tGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCAcg                          >  2:722629/1‑81 (MQ=255)                                             gCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCAcgc                         <  2:722585/80‑1 (MQ=255)                                                 aTCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATcg                    <  1:722586/81‑1 (MQ=255)                                                    aGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATcgc                   >  2:722622/1‑79 (MQ=255)                                                     gCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAg                <  1:722587/81‑1 (MQ=255)                                                     gCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAg                >  2:722637/1‑81 (MQ=255)                                                      cTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTTATAATGGCGCTGCGCACGCGATCGCGAGc               >  1:722632/1‑81 (MQ=255)                                                         gAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGcttc            <  2:722588/81‑1 (MQ=255)                                                             cAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCttt         >  2:722642/1‑80 (MQ=255)                                                              aGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTa       >  1:722623/1‑81 (MQ=255)                                                                 ttGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCg    >  2:722618/1‑81 (MQ=255)                                                                   gATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCgtt  >  2:722634/1‑81 (MQ=255)                                                                           |                                                                          CAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTT  >  NC_000913/3949703‑3949850

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: ‑