Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,345,727	A→G	F534L (TTC→CTC)	fumB ←	anaerobic class I fumarate hydratase (fumarase B)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,345,727	0	A	G	100.0%	171.1 / NA	45	F534L (TTC→CTC)	fumB	anaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/19);  total (26/19)

TTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTAC  >  NC_000913/4345592‑4345856                                                                                                                                        |                                                                                                                                  ttCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTc                                                                                                                                >  1:782110/1‑139 (MQ=255)           aGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGt                                                                                                                              <  2:1369633/131‑1 (MQ=255)           aGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGt                                                                                                                              >  1:1369633/1‑131 (MQ=255)           aGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACTTTTg                                                                                                                      >  2:541192/1‑139 (MQ=255)           aGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTg                                                                                                                      >  2:1381666/1‑139 (MQ=255)               tGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCAt                                                                                                                  <  2:387244/139‑1 (MQ=255)                gAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATc                                                                                                                 <  1:1381666/139‑1 (MQ=255)                   aTTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGAc                                                                                                              <  1:1262103/139‑1 (MQ=255)                    ttACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGAcc                                                                                                             <  1:431473/139‑1 (MQ=255)                            gCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATaaa                                                                                                     >  1:802863/1‑139 (MQ=255)                                                 aTCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCggg                                                                                              <  2:244505/125‑1 (MQ=255)                                                 aTCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCggg                                                                                              >  1:244505/1‑125 (MQ=255)                                                        ggCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTc                                                                         >  2:877654/1‑139 (MQ=255)                                                           agcaTGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCa                                                                                                                   >  2:1067646/1‑94 (MQ=255)                                                           agcaTGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCa                                                                                                                   <  1:1067646/94‑1 (MQ=255)                                                           agcaTGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAg                                                                      >  1:447681/1‑139 (MQ=255)                                                           agcaTGCTGCCAGGCGCTCGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAg                                                                      >  2:824439/1‑139 (MQ=255)                                                                       ggCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCa                                                            >  1:692003/1‑137 (MQ=255)                                                                       ggCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCa                                                            <  2:692003/137‑1 (MQ=255)                                                                           cTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGa                                                                                                               >  2:392956/1‑82 (MQ=255)                                                                           cTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGa                                                                                                               <  1:392956/82‑1 (MQ=255)                                                                             gggCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGc                                                    >  1:975285/1‑139 (MQ=255)                                                                              ggCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCt                                                   >  1:53884/1‑139 (MQ=255)                                                                                          tACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGAt                                                                    >  1:358543/1‑110 (MQ=255)                                                                                          tACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGAt                                                                    <  2:358543/110‑1 (MQ=255)                                                                                                 tGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGc                                                                 >  2:1221660/1‑106 (MQ=255)                                                                                                 tGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGc                                                                 <  1:1221660/106‑1 (MQ=255)                                                                                                      ttCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATa                                                                   >  1:1119245/1‑99 (MQ=255)                                                                                                      ttCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATa                                                                   <  2:1119245/99‑1 (MQ=255)                                                                                                      ttCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAg                                                                  <  1:1404488/100‑1 (MQ=255)                                                                                                      ttCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAg                                                                  >  2:1404488/1‑100 (MQ=255)                                                                                                      ttCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGAt                           >  1:1191026/1‑139 (MQ=255)                                                                                                       tCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAg                                                                  >  2:635874/1‑99 (MQ=255)                                                                                                       tCGCGCACTGTTTGTTGACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAg                                                                  <  1:635874/99‑1 (MQ=255)                                                                                                                    gttgACGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTctgc             <  2:53884/139‑1 (MQ=255)                                                                                                                        aCGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTgcg         <  2:447681/139‑1 (MQ=255)                                                                                                                         cGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACg                                     >  2:938118/1‑110 (MQ=255)                                                                                                                         cGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACg                                     <  1:938118/110‑1 (MQ=255)                                                                                                                         cGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTgcgc        >  2:1456867/1‑139 (MQ=255)                                                                                                                         cGATTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTgcgc        >  1:13865/1‑139 (MQ=255)                                                                                                                           aTTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCa      >  1:948163/1‑139 (MQ=255)                                                                                                                           aTTTGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCa      >  1:1050512/1‑139 (MQ=255)                                                                                                                             ttGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAg                                                                  <  2:285537/77‑1 (MQ=255)                                                                                                                             ttGCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAg                                                                  >  1:285537/1‑77 (MQ=255)                                                                                                                               gCTGGAAGAGGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGTTTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTAc  <  2:1050512/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  TTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTAC  >  NC_000913/4345592‑4345856

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