Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,346,086	A→G	I414T (ATC→ACC)	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,346,086	0	A	G	100.0%	127.9 / NA	36	I414T (ATC→ACC)	fumB	anaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/17);  total (19/17)

CTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCG  >  NC_000913/4345957‑4346214                                                                                                                                  |                                                                                                                                 ctgcCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCg                                                                                                                         >  1:795643/1‑139 (MQ=255)  tgcCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCgg                                                                                                                        <  2:795643/139‑1 (MQ=255)       ccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGtt                                                                                                                   >  2:794769/1‑139 (MQ=255)          tGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGGGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTcat                                                                                                                >  1:431183/1‑137 (MQ=255)                                  gAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCa                                                                                        >  2:568085/1‑139 (MQ=255)                                  gAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCa                                                                                        >  2:560931/1‑139 (MQ=255)                                       cATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTg                                                                                   >  2:50312/1‑139 (MQ=255)                                         tACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                                 <  1:568085/139‑1 (MQ=255)                                          aCGGCCTGCGGTGGTTGGGCCAAGTTAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCTGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCg                                                                                >  1:685358/1‑139 (MQ=255)                                           cGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgt                                                                               <  2:628227/139‑1 (MQ=255)                                            ggCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgtg                                                                              >  1:93677/1‑139 (MQ=255)                                            ggCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgtg                                                                              <  1:794769/139‑1 (MQ=255)                                                tGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCa                                                                          <  1:560931/139‑1 (MQ=255)                                                      ggtTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATAtct                                                                    <  1:359206/139‑1 (MQ=255)                                                           ggCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                               <  2:716814/139‑1 (MQ=255)                                                                          aTGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                                           >  2:288603/1‑96 (MQ=255)                                                                          aTGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                                           <  1:288603/96‑1 (MQ=255)                                                                           tGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCg                                               >  1:136593/1‑139 (MQ=255)                                                                                    ggCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTATCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGAc                                      >  2:500842/1‑139 (MQ=255)                                                                                    ggCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                                           <  1:681549/86‑1 (MQ=255)                                                                                    ggCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                                           >  2:681549/1‑86 (MQ=255)                                                                                     gCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACa                                     >  2:47868/1‑139 (MQ=255)                                                                                       aGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCttt                                                                                          <  2:86099/84‑1 (MQ=255)                                                                                       aGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCttt                                                                                          >  1:86099/1‑84 (MQ=255)                                                                                            ttttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAg                                          >  2:68630/1‑127 (MQ=255)                                                                                            ttttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAg                                          <  1:68630/127‑1 (MQ=255)                                                                                                        ccGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGa                                                                                                                       <  1:189471/38‑1 (MQ=255)                                                                                                        ccGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGa                                                                                                                       >  2:189471/1‑38 (MQ=255)                                                                                                        ccGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTa                  <  1:47868/139‑1 (MQ=255)                                                                                                                 agATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAg         <  2:136593/139‑1 (MQ=255)                                                                                                                   aTCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCt       <  1:171234/139‑1 (MQ=255)                                                                                                                     cGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                              <  2:774103/114‑1 (MQ=255)                                                                                                                     cGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                              >  1:774103/1‑114 (MQ=255)                                                                                                                        gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCg  <  2:93677/139‑1 (MQ=255)                                                                                                                        gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCg  >  2:101942/1‑139 (MQ=255)                                                                                                                        gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCg  >  1:527869/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                 CTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCG  >  NC_000913/4345957‑4346214

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