Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,347,002	C→T	G109S (GGC→AGC)	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,347,002	0	C	T	100.0%	84.3 / NA	26	G109S (GGC→AGC)	fumB	anaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/13);  total (13/13)

GCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGT  >  NC_000913/4346873‑4347130                                                                                                                                  |                                                                                                                                 gcgcCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTg                                                                                                                         >  2:998507/1‑139 (MQ=255)  cgccgcATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTgg                                                                                                                        <  2:9389/139‑1 (MQ=255)                    gcgcAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCt                                                                                                      >  2:95186/1‑139 (MQ=255)                          ggtATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCACGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcg                                                                                                >  2:742237/3‑139 (MQ=255)                           ttATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcgg                                                                                               <  2:902822/139‑1 (MQ=255)                                        aGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGgcag                                                                                                             <  2:97813/112‑1 (MQ=255)                                        aGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGgcag                                                                                                             >  1:97813/1‑112 (MQ=255)                                               ataGACGCCATTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAg                                                                           >  2:989770/1‑139 (MQ=255)                                                  gACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAAt                                                                        >  1:1154719/1‑139 (MQ=255)                                                    cGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGt                                                                                   >  1:409720/1‑126 (MQ=255)                                                    cGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGt                                                                                   <  2:409720/126‑1 (MQ=255)                                                       cTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAg                                                                   >  1:1171898/1‑139 (MQ=255)                                                       cTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAg                                                                   >  2:503687/1‑139 (MQ=255)                                                                 cGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTAc                                                         >  1:294231/1‑139 (MQ=255)                                                                               cgccgccgGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCt                                           <  2:309886/139‑1 (MQ=255)                                                                               cgccgccgGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCt                                           <  2:1171898/139‑1 (MQ=255)                                                                               cgccgccgGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCt                                           <  1:503687/139‑1 (MQ=255)                                                                               cgccgccgGTCCACACGCGCTGGCCTTTATTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCt                                           <  1:742237/139‑1 (MQ=255)                                                                                 ccgccgGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTgg                                         <  2:294231/139‑1 (MQ=255)                                                                                    ccgGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCAGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCtt                                      <  2:458407/139‑1 (MQ=255)                                                                                                 gcTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTgaag                         >  2:194902/1‑139 (MQ=255)                                                                                                  cTGGCCTTTTTTACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTgaaga                        >  1:586515/1‑139 (MQ=255)                                                                                                            ttACCGACGATGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCaa              <  2:52796/139‑1 (MQ=255)                                                                                                                 gacgaTGATCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGt         <  1:989770/139‑1 (MQ=255)                                                                                                                    gatgatCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTc      >  2:724503/1‑139 (MQ=255)                                                                                                                        atCGCGGTGCTGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGgt  <  1:95186/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                 GCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGT  >  NC_000913/4346873‑4347130

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