Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,948,805	(G)6→7	coding (1165/1521 nt)	yifB ←	magnesium chelatase family protein and putative transcriptional regulator

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%	145.0 / NA	41	G391R (GGC→CGC)	yifB	magnesium chelatase family protein and putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (22/19);  total (22/19)

TTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAA  >  NC_000913/3948665‑3948917                                                                                                                                        |                                                                                                                       ttGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggg                                                                                                                    >  2:967363/1‑139 (MQ=255) ttGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggg                                                                                                                    >  2:404345/1‑139 (MQ=255) ttGGCGTATTTCCCGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggg                                                                                                                    <  1:1617047/139‑1 (MQ=255)   ggCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggg                                                                                                                  <  2:12629/139‑1 (MQ=255)    gCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggggg                                                                                                                 <  1:967363/139‑2 (MQ=255)    gCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggggg                                                                                                                 <  1:20769/139‑2 (MQ=255)               gACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                                                                                                       <  2:1358254/139‑1 (MQ=255)                   aTCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                                                                                                   <  1:971310/139‑1 (MQ=255)                    tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTcc                                                                                                  >  2:1331266/1‑139 (MQ=255)                    tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTcc                                                                                                  >  1:158788/1‑139 (MQ=255)                    tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggggg                                                                                                                 >  1:46544/1‑122 (MQ=255)                    tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggggg                                                                                                                 <  1:680892/123‑2 (MQ=255)                    tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggggg                                                                                                                 <  2:46544/123‑2 (MQ=255)                    tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggggg                                                                                                                 >  2:680892/1‑122 (MQ=255)                            ggCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgagag                                                                                          >  2:359963/1‑139 (MQ=255)                             gCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgagaga                                                                                         >  1:1229957/1‑139 (MQ=255)                                   aGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg                                                                                   <  2:268544/139‑1 (MQ=255)                                   aGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg                                                                                   <  2:1229957/139‑1 (MQ=255)                                       tATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                                                                                                       <  1:1630855/115‑1 (MQ=255)                                       tATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                                                                                                       >  2:1630855/1‑115 (MQ=255)                                             gCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGaa                                                                         <  1:404345/139‑1 (MQ=255)                                               cGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAg                                                                       >  1:286625/1‑139 (MQ=255)                                               cGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gg                                                                                                                      <  2:1727718/91‑1 (MQ=255)                                               cGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gg                                                                                                                      >  1:1727718/1‑91 (MQ=255)                                                           cgcTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                                                                                                         >  1:758791/1‑77 (MQ=255)                                                           cgcTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                                                                                                         <  2:758791/77‑1 (MQ=255)                                                           cgcTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTg                                                                                              >  1:661401/1‑104 (MQ=255)                                                           cgcTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTg                                                                                              <  2:661401/104‑1 (MQ=255)                                                                    gcggcCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTgaga                                                  <  2:158788/139‑1 (MQ=255)                                                                          aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGTC‑ggggggg                                                                                                                 <  2:1473839/69‑2 (MQ=255)                                                                          aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGTC‑ggggggg                                                                                                                 >  1:1473839/1‑68 (MQ=255)                                                                          aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                                                                                           <  1:1110944/76‑1 (MQ=255)                                                                          aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                                                                                           >  2:1110944/1‑76 (MQ=255)                                                                          aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc                                                                              >  2:379116/1‑105 (MQ=255)                                                                          aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGaa                                                                         >  1:1582322/1‑110 (MQ=255)                                                                          aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGaa                                                                         <  2:1582322/110‑1 (MQ=255)                                                                          aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAGGGGATCTCCAGTGAGAGATCGAAGCGGTc                                                                              <  1:379116/105‑1 (MQ=255)                                                                          aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggggg                                                                                                                 >  1:1087965/1‑68 (MQ=255)                                                                          aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggggg                                                                                                                 <  2:1087965/69‑2 (MQ=255)                                                                          aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggg                                                                                                                  <  1:372808/68‑1 (MQ=255)                                                                          aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggggg                                                                                                                  >  2:372808/1‑68 (MQ=255)                                                                              cgcgTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg                                                                                   >  1:655066/1‑96 (MQ=255)                                                                              cgcgTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg                                                                                   <  2:655066/96‑1 (MQ=255)                                                                               gcgTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGt                                       <  1:1331266/139‑1 (MQ=255)                                                                                  ttgttTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                                                                                                       >  2:1586593/1‑72 (MQ=255)                                                                                  ttgttTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                                                                                                       <  1:1586593/72‑1 (MQ=255)                                                                                         aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATggg                                                                                                        >  2:735966/1‑64 (MQ=255)                                                                                         aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATggg                                                                                                        <  1:735966/64‑1 (MQ=255)                                                                                             tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGc                         >  1:1722285/1‑139 (MQ=255)                                                                                               gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                                                                                                   <  2:945429/63‑1 (MQ=255)                                                                                               gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                                                                                                   >  1:945429/1‑63 (MQ=255)                                                                                               gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gg                                                                                                                      <  1:1427339/43‑1 (MQ=255)                                                                                               gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gg                                                                                                                      >  2:1427339/1‑43 (MQ=255)                                                                                                      tttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTAt                >  1:1341606/1‑139 (MQ=255)                                                                                                               gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGAt                                                 <  2:322170/97‑1 (MQ=255)                                                                                                               gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGAt                                                 >  1:322170/1‑97 (MQ=255)                                                                                                               gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCt       >  2:485718/1‑139 (MQ=255)                                                                                                                    aCCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                                                                                                <  1:267164/45‑1 (MQ=255)                                                                                                                    aCCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                                                                                                >  2:267164/1‑45 (MQ=255)                                                                                                                    aCCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATaa  >  2:1453989/1‑139 (MQ=255)                                                                                                                    aCCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATaa  >  1:986434/1‑139 (MQ=255)                                                                                                                                  aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGagag                                                            <  2:587247/67‑1 (MQ=255)                                                                                                                                  aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGagag                                                            >  1:587247/1‑67 (MQ=255)                                                                                                                                        |                                                                                                                       TTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAA  >  NC_000913/3948665‑3948917

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