Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,912,293	G→A	Y453Y (TAC→TAT)	relA ←	GDP/GTP pyrophosphokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,912,293	0	G	A	100.0%	32.4 / NA	13	Y453Y (TAC→TAT)	relA	GDP/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/8);  total (5/8)

CCAGCTCGTCGTCAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTCACGGTCCTGTTTACGGAACCAGGCGTGAATTTTCGAACGCCCACGGCTGGTTGTGACGTAACCGAGGTTTGGGTTTAACCAGTCACGGCTGGGGTTCGGCTGTTTCTGGGTGATAATTTCAATCTGGTCGCCCATCTGCAGCTGGTAGGTGAACGGCACAATGCGCCCGCCAATTTTTGCCCCGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCC  >  NC_000913/2912103‑2912391                                                                                                                                                                                               |                                                                                                   ccAGCTCGTCGTCAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTCACGGTCCTGTTGACGGAACCAGGCGTGAATTTTCGAACGCCCACGGCTGTTTGTGACGTAACCGATGTTTGGGTTTAACCAGTCATGGCTGTGGTTCTGCTGTTTCTGGTTGATAATGTCAATCTGGTCGCCCATCTGCAGCTGATAGGTGAACGGcac                                                                                      <  1:379789/205‑1 (MQ=255)        gtcgtcAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTCACGGTCCTGTTTACGGAACCAGGCGTGAATTTTCGAACGCCCACGGCTGGTTGTGACGTAACCGAGGTTTGGGTTTAACCAGTCACGGCTGGGGTTCGGCTTTTTCTGGGTGATAATTTCAATCTGGTCGCCCATCTGCAGCTGATAGGTGAACGGCACAATGCGCCCGCCAATTTTTGCCCCGATGCAGCGGTGTCCGACATCACTGTGGa                                 <  1:39757/251‑1 (MQ=255)          cgtcAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTCACGGTCCTGTTTACGGAACCAGGCGTGAATTTTCGAACGCCCACGGCTGGTTGTGACGTAACCGAGGTTTGGGTTTAACCAGTCACGGCTGGGGTTCGGCTGTTTCTGGGTGATAATTTCAATCTGGTCGCCCATCTGCAGCTGATAGGTGAACGGCACAATGCGCCCGCCAATTTTTGCCCCGATGCAGCGGTGTCCGACATCACt                                      >  1:369513/1‑244 (MQ=255)                     tGCCGCCCAGCCAGAATGTTTTTGTCACGGTCCTGTTTACGGAACCAGGCGTGAATTTTCGAACGCCCACGGCTGGTTGTGACGTAACCGAGGTTTGGGTTTAACCAGTCACGGCTGGGGTTCGGCTGTTTCTGGGTGATAATTTCAATCTGGTCGCCCATCTGCAGCTGATAGGTGAACGGCACAATgc                                                                                 <  1:33577/190‑1 (MQ=255)                                  gAATGTTTTTGTCACGGTCCTGTTTACGGAACCAGGCGTGAATTTTCGAACGCCCACGGCTGGTTGTGACGTAACCGAGGTTTGGGTTTAACCAGTCACGGCTGGGGTTCGGCTGTTTCTGGGTGATAATTTCAATCTGGTCGCCCATCTGCAGCTGATAGGTGAACGGCACAAtn                                                                                  <  1:306763/176‑2 (MQ=255)                                                       gTTTACGGAACCAGGCGTGAATTTTCGAACTCCCCCGGCTGGTTGTGCCGTAACCGAGGTTTGGGTTTAACCAGTCACGGCTGGGGTTCGGCTGTTTCTTGGGGATAATTTCAATCTGGTCGCCCCTCTGCAGCTGATAGGTGa                                                                                             >  1:22915/1‑144 (MQ=255)                                                              gAACCAGGCGTGAATTTTCGAACGCCCACGGCTGGTTGTGACGTAGCCGAGGTTTGTGTTTAACCGGTCACGGCTGGGGTTCGGCTGTTTCTGGGTTATAATTTCAATCTGGTCGCCCATCTGCAGCTGATAGGTGAGCGGCACAATGCGCCCGCCAATTTTTGCCCCGATGCAgn                                                      <  1:90562/176‑2 (MQ=255)                                                                        tGAATTTTCGAACGCCCACGGCTGGTTGTGACGTAACAGAGGTTTGGGTTTAACCAGTCCCGGCTGGGGTTCGGATGCTTCTGGGTGATAATTTCAATCTGGTCGCCCATCTGCAGCTGATAGGTGAACGGCACAATgcgc                                                                               >  1:203008/1‑141 (MQ=255)                                                                                                           nccGAGGTTTGGGTTTAACCAGTCACGGCTGGGGTTCGGCTGTTTCTGGGTGATAATTTCAATCTGGTCGCCCATCTGCAGCTGATAGGTGAACGGCACAATGcgcc                                                                              >  1:304279/2‑107 (MQ=255)                                                                                                           aCCGAGGTTTGGGTTTAACCAGTCACGGCTGGGGTTCGGCTGTTTCTGGGTGATAATTTCAATCTGGTCGCCCATCTGCAGCTGATAGGTGAACGGCACAATGcgcc                                                                              <  2:304279/107‑1 (MQ=255)                                                                                                                                                   tGTTTNNGGNNNATNNTTTCAATCTGGTCGCCCATCTGCNGCTGATAGGTGAACGGCACAATGCGCCCGNCNANTTTTNNCCCNATGCAGCGGTNTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGCTcc  >  2:379551/1‑143 (MQ=255)                                                                                                                                                   ggTTTCTGGGTGATAATGTCAATTTGGTCGCCCATCTGCAGCTGATAGGTGAACGGCACAATGGGCCCTGCAATTTTTGCCCCGATGCAGCGGTGTCCGACATCGCTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATcc  <  1:379551/142‑1 (MQ=255)                                                                                                                                                      tnCTGGGNGATANTTTCAATCTGGTCGCNCATCTGCAGCTGATAGGTGNNCNGCACANTgcgc                                                                               <  2:203008/63‑1 (MQ=38)                                                                                                                                                                                               |                                                                                                   CCAGCTCGTCGTCAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTCACGGTCCTGTTTACGGAACCAGGCGTGAATTTTCGAACGCCCACGGCTGGTTGTGACGTAACCGAGGTTTGGGTTTAACCAGTCACGGCTGGGGTTCGGCTGTTTCTGGGTGATAATTTCAATCTGGTCGCCCATCTGCAGCTGGTAGGTGAACGGCACAATGCGCCCGCCAATTTTTGCCCCGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCC  >  NC_000913/2912103‑2912391

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