Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,345,751	G→A	L526L (CTG→TTG)	fumB ←	fumarase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,345,751	0	G	A	100.0%	33.2 / NA	13	L526L (CTG→TTG)	fumB	fumarase B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/8);  total (5/8)

GTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGC  >  NC_000913/4345526‑4345974                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                gtTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGGGGCCAGGGGGATCCCCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTTCTGAGTGCAGTTCGCTCACTGTTTGTTGACGATTTGCTGTAAGAAGTCGTTACCTTTGTCATCGACCAAGATAAACGCCGGGAAATCTTCTACt                                                                                                                                                                                                        <  1:42752/251‑1 (MQ=255)    aaGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTGACTTAGGGCAGTTCGCGCACTGTTTGTTGACGATTGGCGGGAAGAAGTCGTTACCTTTGTCATCGACCAAGATAAACGCCGGGAAATCTTCTACTTcn                                                                                                                                                                                                     <  1:62755/251‑2 (MQ=255)              ggCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCGCCAGCAAGCGTGAAATTACAATCGCAAAAGATGTGGCCAGGGGTATCACCTGGCAGCATGTTTCCAGGCGCGGGGCCGAAGAGTTGACTTAATGCAGTTCGCGCACTTTTTGTTGGCGATTGGCGGGAAGAAGCCGTTACCTTGGTCATCGGCCCAGATAAACGCCGGGAAATCTTCTACTTc                                                                                                                                                                                                      <  1:293503/240‑1 (MQ=255)                  ggAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAAGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAgn                                                                                                                                                                                       <  1:194191/250‑2 (MQ=255)                                                     gAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAAGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCg                                                                                                                                                                      <  1:68261/233‑1 (MQ=255)                                                                                                      ggTGGCCAGGGGGATCACCTGGCAGCATGATGCCAGGCGCTGGGCCGAACAGGATAATTAGTGCAGCTCGCGCCCTGTATGTTGACGATTTGATGGAAGAAGCCGTTATCTTTGTCATAGAACAAGATAAACGACGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCACGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTctgc                                                                                                                                   >  1:314594/1‑219 (MQ=255)                                                                                                       gTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTAGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAAGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCAGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCCGTACCGCCGCCGGACCGCCGATGCt                                                                                                   >  1:207425/1‑250 (MQ=255)                                                                                                               ggggATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAAGATAACCGCCGGGCAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAg                                                                                             >  1:149636/1‑248 (MQ=255)                                                                                                                        cTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTAGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAAGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTAcc                                                                                                                       <  1:330630/213‑1 (MQ=255)                                                                                                                                                                        ttCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTTCCTTTGTCATCGACCAAGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCGACCTAGGTAGAATCCGCCATTTTTATGACACTCGTCgg                                                             <  1:259735/223‑1 (MQ=255)                                                                                                                                                                                                          agTCTTTACCTTTTTCATCGACCAAGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATATCTTCCATCCCCAGCGCCGGATAAGCGACGCACTCCAGAGGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGTGATTgc  <  1:69998/248‑1 (MQ=255)                                                                                                                                                                                                          agTCGTTACCTTTGTCATCGACCAAGNTAACCGNNGGNAAATCTTCTACTTCGATTTNCCNNNTNNNTNCCATANCCNNCTCCGGATANGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCNCCGCTGCTANNGAGGTAGAAGACGCCGTGTTTNTGANNCNNGNCggnt                                                           >  2:69998/1‑189 (MQ=255)                                                                                                                                                                                                                             gACCAAGATAAACGCCGGGAAATCTTCTACTTCGATTTTCAAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATc                       >  1:215202/1‑208 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                GTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGC  >  NC_000913/4345526‑4345974

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