Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	908,760	A→T	C179S (TGC→AGC)	ltaE ←	low‑specificity L‑threonine aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	908,760	0	A	T	100.0%	43.7 / NA	16	C179S (TGC→AGC)	ltaE	low‑specificity L‑threonine aldolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (9/7);  total (9/7)

CCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCATATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTC  >  NC_000913/908514‑908985                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                  ccGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCATATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCCGCGAATGGCACGTTTAATGTAATCTCGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCTGCCg                                                                                                                                                                                                                               >  1:442654/1‑251 (MQ=255)                                                         aGGGCATATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAAGCACAATAGGGCGTGATCTCTTGCAGTTCGCTGCCGTAAGCCGCCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTgcg                                                                                                                                                                       <  1:393865/251‑1 (MQ=255)                                                                  atCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTTCGCCAGCGAATGGCACGTTTAATGTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCTGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCggn                                                                                                                                                               <  1:110389/250‑2 (MQ=255)                                                                              gCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCTGCCGTAAGccac                                                                                                                                                                                                                       >  1:431147/1‑182 (MQ=255)                                                                                               gcgcATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCTGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTccc                                                                                                                                 >  1:163388/1‑251 (MQ=255)                                                                                                   ntCCCGCCACCTGTCATTTTCCGCCAGAGAATGGCACGTTTAATGTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCTGCCg                                                                                                                                                                                                                               >  1:110888/2‑153 (MQ=255)                                                                                                        gCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTAATCACGATTACCGACGAGTAATGAACAGACTGGCGTCCCAAGACCTTTCGAAAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCTGCCGTAAGCCACCACGGCATTCAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAac                                                                                                                         >  1:216562/1‑250 (MQ=255)                                                                                                                                                                                    ntCCCAAGACATTTAGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCTGCCGTAAGccacca                                                                                                                                                                                                                     >  1:130398/2‑82 (MQ=255)                                                                                                                                                                                    gTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCTGCCGTAAg                                                                                                                                                                                                                           >  2:278551/1‑76 (MQ=255)                                                                                                                                                                                    gTCCCAAGACCTTGCGACAGGCAAATGGTGAACGAATCCCAATATTGCGTGATCTCTTTCAGTTCGCTGCCGTaan                                                                                                                                                                                                                           <  1:278551/76‑2 (MQ=255)                                                                                                                                                                                      cacAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTGCGCTGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCg                                          <  1:37284/249‑1 (MQ=255)                                                                                                                                                                                                  cGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCTGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATcc                               >  1:226012/1‑250 (MQ=255)                                                                                                                                                                                                   gACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCTGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAACTGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATcc                               >  1:259775/1‑249 (MQ=255)                                                                                                                                                                                                              ggTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCTGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACGGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCAGCGCCACTTTATCCAGCGGTAGCGTGc                  <  1:160581/251‑1 (MQ=255)                                                                                                                                                                                                                       aTCACAATATTGCGTGATCTCTTTCAGTTCGCTGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCa                                                                                                                                          <  1:14454/122‑1 (MQ=255)                                                                                                                                                                                                                               atTGCGTGATCTCTTTCAGTTCGCTGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCTCTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTc  <  1:202683/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                  CCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCATATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTC  >  NC_000913/908514‑908985

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