Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,980,232	C→T	R135H (CGC→CAC)	ygeA ←	Asp/Glu_racemase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,980,232	0	C	T	100.0%	44.4 / NA	14	R135H (CGC→CAC)	ygeA	Asp/Glu_racemase family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/5);  total (9/5)

CGCGTGACGCTTCGGTAAATTGCCCCAGACACAGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGCGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCG  >  NC_000913/2980110‑2980367                                                                                                                           |                                                                                                                                        cgcgTGACGCTTCGGTAAATTGCCCCAGACACAGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATaaaaa                                                                                                                                 <  1:639339/131‑1 (MQ=255) cgcgTGACGCTTCGGTAAATTGCCCCAGACACAGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATaaaaa                                                                                                                                 >  2:639339/1‑131 (MQ=255)                                              taatCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGc                                                                            <  2:71894/139‑1 (MQ=255)                                                                   gTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGt                                                       <  2:163324/139‑1 (MQ=255)                                                                              tCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCaa                                            >  1:318374/1‑139 (MQ=255)                                                                               cAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAAt                                           >  1:20156/1‑139 (MQ=255)                                                                                  gAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAAtgtg                                        >  1:305716/1‑139 (MQ=255)                                                                                             ttGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCa                             >  1:289445/1‑139 (MQ=255)                                                                                                aTGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAAc                                                                                                              >  1:366847/1‑55 (MQ=255)                                                                                                aTGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAAc                                                                                                              <  2:366847/55‑1 (MQ=255)                                                                                                  ggAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTg                        >  1:154/1‑139 (MQ=255)                                                                                                  ggAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTg                        <  2:318374/139‑1 (MQ=255)                                                                                                                  aGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACTTGACTCAATgg        >  2:45246/1‑139 (MQ=255)                                                                                                                  aGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATgg        >  1:538891/1‑139 (MQ=255)                                                                                                                        ccgtGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGAGACAAGGGAGTCCGTCGGGCCCCGGTAATTTCTAGAACGGAGGCATACGCATTGTGTAAACAAGGCAGAGTGCAATATGGCTCAATTGCATCCg  <  1:13194/139‑1 (MQ=255)                                                                                                                           |                                                                                                                                        CGCGTGACGCTTCGGTAAATTGCCCCAGACACAGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGCGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCG  >  NC_000913/2980110‑2980367

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