Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,346,121	C→T	E402E (GAG→GAA)	fumB ←	fumarase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,346,121	0	C	T	100.0%	15.8 / NA	8	E402E (GAG→GAA)	fumB	fumarase B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (4/4);  total (4/4)

TACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTC  >  NC_000913/4345934‑4346243                                                                                                                                                                                            |                                                                                                                           tACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGTTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGTCGGTGAGCGACAAACg                                                              >  1:213728/1‑250 (MQ=255)                       gctgcCGCCGTGGGATTGCCGCAGATNCACGTANNAGNCCCTACGGCCTGCGGCGGTNGGNNNANNNGNCCCTGNTGNNTAACCGGCANGGGGTTTCGCCGGACCCGCGTAGTCGATCGGGTGATCTTTGTTGTACTGCGGGNGTTCTTTACNNGCGTCAATCAGTTCTTTCAGCNTGGNNTNNGNAatnccn                                                                                                 >  2:57632/1‑189 (MQ=255)                                       tGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACATGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCCGTTCTTTCAGCTTGGCGTGTGCAATATCTAGGCCCACGCTAATGGTGTCGGTGAGCGACAAACGAGTGGATCCCGGGTATTGCGAAAGCTGGGCGAGGAtctc                       >  1:224455/1‑251 (MQ=255)                                                             ccATACGGCCTGCGGTGGTTGGGGCAAGTGAACCTGAGGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGTTCTTTCAGCTTGGCGTGGGCAATATCTCGGCCCACGATAATGGTGTCGGTGAGCGACAAACGAGTGGATTCTGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCTGTTAAGGTc  <  1:57632/250‑1 (MQ=255)                                                                                                  tGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGTTCTTTCAGCTTGGCGTGTGCa                                                                                                       <  2:154543/112‑1 (MQ=255)                                                                                                  nggATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGTTCTTTCAGCTTGGCGTGTGCa                                                                                                       >  1:154543/2‑112 (MQ=255)                                                                                                                            gACCCGCGTAGTAGATCGTGTGATCGTTGATGTACTGCGGAATTTCTTTACCGGCGTCAATCAGTTCTTTCAGCTTGGCGTGTGCAATATCTCgg                                                                                              <  1:120480/95‑1 (MQ=255)                                                                                                                                                            tACTGCGTAAGTTCTTTACCGGCTTCAATCATTTCTTTCAGCTTGGCGTTTGCAAGATCTCGGCCCACGATAATCGTGTCGGTGAGCGACATACGAGTGCATACCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAg     <  1:114021/152‑1 (MQ=255)                                                                                                                                                                                            |                                                                                                                           TACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTC  >  NC_000913/4345934‑4346243

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