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

GAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGC  >  NC_000913/3948662‑3948930                                                                                                                                           |                                                                                                                                    gAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gg                                                                                                                                   <  1:41619/140‑1 (MQ=255)     tGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggg                                                                                                                               <  2:582510/140‑1 (MQ=255)     tGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggg                                                                                                                               <  1:45973/140‑1 (MQ=255)     tGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggg                                                                                                                               <  1:679019/140‑1 (MQ=255)      ggCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggggg                                                                                                                              <  2:565048/140‑2 (MQ=255)       gCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCCCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggggt                                                                                                                             >  1:539410/1‑138 (MQ=255)            tttCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggg                                                                                                                               <  2:845436/133‑1 (MQ=255)            tttCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggg                                                                                                                               >  1:845436/1‑133 (MQ=255)             ttCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                                                                                                        <  1:166248/140‑1 (MQ=255)                  gACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCCCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                                                                                                                   <  2:539410/140‑1 (MQ=255)                     tATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                                                                                                                <  2:704834/140‑1 (MQ=255)                       tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCa                                                                                                              >  1:93991/1‑140 (MQ=255)                       tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCa                                                                                                              >  1:840649/1‑140 (MQ=255)                       tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCa                                                                                                              >  1:741288/1‑140 (MQ=255)                       tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggggg                                                                                                                              <  2:902336/123‑2 (MQ=255)                       tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggggg                                                                                                                              >  1:902336/1‑122 (MQ=255)                            ccagGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                                                                                                          >  2:825863/1‑139 (MQ=255)                            ccagGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                                                                                                          <  1:825863/139‑1 (MQ=255)                               ggCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGCGGAAATGGGATCTCCAGTgagaga                                                                                                      >  1:114831/1‑140 (MQ=255)                                 cATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc                                                                                                    >  2:809048/1‑140 (MQ=255)                                    tCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa                                                                                                 <  2:741288/140‑1 (MQ=255)                                     cAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg                                                                                                >  2:832605/1‑140 (MQ=255)                                         ttATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt                                                                                            <  1:860929/140‑1 (MQ=255)                                          tATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc                                                                                           >  1:672921/1‑140 (MQ=255)                                          tATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc                                                                                           >  1:467829/1‑140 (MQ=255)                                                gCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGaaa                                                                                     <  2:672921/140‑1 (MQ=255)                                                gCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGaaa                                                                                     >  1:301203/1‑140 (MQ=255)                                                gCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGaaa                                                                                     >  2:257869/1‑140 (MQ=255)                                                 ccGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAg                                                                                    <  2:467829/140‑1 (MQ=255)                                                            ggCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGagag                                                                         <  1:173034/140‑1 (MQ=255)                                                                 tctctGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGt                                                                    <  1:832605/140‑1 (MQ=255)                                                                   tctGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTg                                                                  >  1:415893/1‑140 (MQ=255)                                                                       gcggcCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGAt                                                              <  1:76520/140‑1 (MQ=255)                                                                       gcggcCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGAt                                                              <  2:415893/140‑1 (MQ=255)                                                                             aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAGGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTa                                                        <  1:171126/140‑1 (MQ=255)                                                                             aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAGGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTa                                                        >  2:171126/1‑140 (MQ=255)                                                                             aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTa                                                        >  2:235809/1‑140 (MQ=255)                                                                             aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACATa                                                        >  2:364246/1‑140 (MQ=255)                                                                               tatgcgTTGTTTAACGGTGGCGCTGCTTTCACCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAAt                                                      <  2:589941/137‑1 (MQ=255)                                                                                       gttTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTc                                              >  2:406055/1‑140 (MQ=255)                                                                                       gttTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTc                                              >  2:201238/1‑140 (MQ=255)                                                                                                tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCa                                     >  1:348050/1‑140 (MQ=255)                                                                                                  gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGc                                   >  2:77679/1‑140 (MQ=255)                                                                                                    gctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCgg                                 >  1:905998/1‑140 (MQ=255)                                                                                                    gctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGCCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCgg                                 >  1:678731/1‑140 (MQ=255)                                                                                                        ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTAt                             >  1:563709/1‑140 (MQ=255)                                                                                                           tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGt                          <  1:235809/140‑1 (MQ=255)                                                                                                           tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGt                          <  1:257869/140‑1 (MQ=255)                                                                                                           tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGt                          <  2:301203/140‑1 (MQ=255)                                                                                                               ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTcc                      <  2:840649/140‑1 (MQ=255)                                                                                                               ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTcc                      <  1:835814/140‑1 (MQ=255)                                                                                                                   cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGa                  <  1:809048/140‑1 (MQ=255)                                                                                                                       aCCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAAt              <  1:831686/140‑1 (MQ=255)                                                                                                                                cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAg     >  2:868634/1‑140 (MQ=255)                                                                                                                                   aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGc  <  1:868634/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                    GAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGC  >  NC_000913/3948662‑3948930

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