Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,949,114	C→A	G286W (GGG→TGG)	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,114	0	C	A	100.0%	163.7 / NA	48	G286W (GGG→TGG)	yifB	putative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (22/26);  total (22/26)

TTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAG  >  NC_000913/3948977‑3949243                                                                                                                                          |                                                                                                                                  ttttGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATg                                                                                                                                 >  2:305488/1‑140 (MQ=255)    tGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGaa                                                                                                                              <  2:415194/140‑1 (MQ=255)         gtgtGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTgcg                                                                                                                         <  1:469786/140‑1 (MQ=255)          tgtgCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTgcgc                                                                                                                        >  2:123999/1‑140 (MQ=255)          tgggCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTgcgc                                                                                                                        >  2:8408/4‑140 (MQ=255)              cGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCAcc                                                                                                                    >  2:336188/1‑140 (MQ=255)                     gATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTAcc                                                                                                             >  2:226043/1‑140 (MQ=255)                        ggATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATc                                                                                                          >  1:33601/1‑140 (MQ=255)                               cccGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTa                                                                                                   <  1:502276/140‑1 (MQ=255)                                     ttCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATg                                                                                             <  2:567570/140‑1 (MQ=255)                                     ttCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATg                                                                                             <  1:336188/140‑1 (MQ=255)                                       cAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCa                                                                                           <  1:86807/140‑1 (MQ=255)                                       cAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCa                                                                                           <  1:226043/140‑1 (MQ=255)                                              ctctCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGtgat                                                                                    >  2:664149/1‑140 (MQ=255)                                               tctcGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGtgatg                                                                                   >  1:107143/1‑140 (MQ=255)                                                  cGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAgg                                                                                >  2:225130/1‑140 (MQ=255)                                                    cAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTg                                                                              <  2:698361/140‑1 (MQ=255)                                                            ccAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAAc                                                                      >  1:139935/1‑140 (MQ=255)                                                                  tACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAgg                                                                                >  1:343351/1‑124 (MQ=255)                                                                  tACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAgg                                                                                <  2:343351/124‑1 (MQ=255)                                                                     gCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTgc                                                             >  2:124175/1‑140 (MQ=255)                                                                         ttCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTgccgcc                                                         >  1:439915/1‑140 (MQ=255)                                                                         ttCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTgccgcc                                                         >  1:561465/1‑140 (MQ=255)                                                                          tCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCa                                                        >  2:433421/1‑140 (MQ=255)                                                                          tCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCa                                                        <  2:33601/140‑1 (MQ=255)                                                                           cAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCAt                                                       <  1:433421/140‑1 (MQ=255)                                                                           cAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCAt                                                       <  2:139935/140‑1 (MQ=255)                                                                           cAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCAt                                                       <  1:666844/140‑1 (MQ=255)                                                                                    gTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGt                                              <  2:107143/140‑1 (MQ=255)                                                                                        cTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTa                                             >  2:332088/1‑137 (MQ=255)                                                                                        cTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTa                                             <  1:332088/137‑1 (MQ=255)                                                                                        cTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTg                                          >  2:216321/1‑140 (MQ=255)                                                                                           atcaAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGAtt                                       <  1:216321/140‑1 (MQ=255)                                                                                           atcaAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGAtt                                       <  1:455188/140‑1 (MQ=255)                                                                                                      gCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTAcc                            <  1:225130/140‑1 (MQ=255)                                                                                                        aCGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAg                          >  2:521809/1‑140 (MQ=255)                                                                                                         cGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGa                         <  2:265708/140‑1 (MQ=255)                                                                                                         cGCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGa                         <  1:521809/140‑1 (MQ=255)                                                                                                          gCCGTTATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGAc                        >  1:404039/1‑140 (MQ=255)                                                                                                              ttATGCGCCAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGACTTaa                    <  1:305488/140‑1 (MQ=255)                                                                                                                     ccAGCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGACTTAATATCgca             <  1:664149/140‑1 (MQ=255)                                                                                                                        gCGAATTTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGACTTAATATTgcagca          >  2:644017/1‑140 (MQ=255)                                                                                                                        gCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGACTTAATATCgcagca          >  1:552804/1‑140 (MQ=255)                                                                                                                        gCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGACTTAATATCgcagca          <  1:124175/140‑1 (MQ=255)                                                                                                                        gCGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGACTTAATATCgcagca          <  2:404039/140‑1 (MQ=255)                                                                                                                         cGAAATTTCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCAc         >  2:430640/1‑140 (MQ=255)                                                                                                                                tCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAg  <  1:8408/140‑1 (MQ=255)                                                                                                                                tCACCGGGCCATGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAg  <  1:123999/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                  TTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAG  >  NC_000913/3948977‑3949243

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