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%	69.6 / NA	21	R135H (CGC→CAC)	ygeA	Asp/Glu_racemase family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/9);  total (12/9)

TAAGCGCGTGACGCTTCGGTAAATTGCCCCAGACACAGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGCGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAAT  >  NC_000913/2980106‑2980359                                                                                                                               |                                                                                                                                tAAGCGCGTGACGCTTCGGTAAATTGCCCCAGACACAGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCt                                                                                                                     >  2:523189/1‑139 (MQ=255)    gcgcGTGACGCTTCGGTAAATTGCCCCAGACACAGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGtt                                                                                                                  <  2:307420/139‑1 (MQ=255)             gCTTCGGTAAATTGCCCCAGACACAGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGgtgt                                                                                                         <  1:44718/139‑1 (MQ=255)                           ccccAGACACAGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGtt                                                                                                                  <  2:458816/116‑1 (MQ=255)                           ccccAGACACAGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGtt                                                                                                                  >  1:458816/1‑116 (MQ=255)                                     aGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCacac                                                                                 >  1:254802/1‑139 (MQ=255)                                     aGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCacac                                                                                 >  2:207065/1‑139 (MQ=255)                                     aGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCacac                                                                                 >  1:183216/1‑139 (MQ=255)                                          ttcAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTc                                                                            >  1:561333/1‑139 (MQ=255)                                           tcAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCa                                                                           >  1:493927/1‑139 (MQ=255)                                           tcAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCa                                                                           >  1:29854/1‑139 (MQ=255)                                                aataatCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCg                                                                      <  2:29854/139‑1 (MQ=255)                                                   aatCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGcc                                                                   <  1:523189/139‑1 (MQ=255)                                                                           aTCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCa                                           <  2:183216/139‑1 (MQ=255)                                                                                cTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGc                                      <  2:254802/139‑1 (MQ=255)                                                                                 ttCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCa                                     <  2:493927/139‑1 (MQ=255)                                                                                 ttCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCa                                     >  2:119511/1‑139 (MQ=255)                                                                                                 ttGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCg                                                                      <  1:543329/90‑1 (MQ=255)                                                                                                 ttGATGGAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCg                                                                      >  2:543329/1‑90 (MQ=255)                                                                                                      ggAAAATTGTTCCGTCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTg                >  1:537073/1‑139 (MQ=255)                                                                                                                    tCAGCCGCCCGTGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAAt  >  2:587572/1‑139 (MQ=255)                                                                                                                               |                                                                                                                                TAAGCGCGTGACGCTTCGGTAAATTGCCCCAGACACAGTTCTTCAAAAATAATCTGATTAATTTTCGCCCGTTCATCCGCTTCAGGAATAAGACAATTGATGGAAAATTGTTCCGTCAGCCGCCCGCGATAAAAATCCTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAAT  >  NC_000913/2980106‑2980359

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