Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,931,448	G→T	L265F (TTG→TTT)	rbsK →	ribokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,931,448	0	G	T	100.0%	46.9 / NA	13	L265F (TTG→TTT)	rbsK	ribokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (7/6);  total (7/6)

CAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAA  >  NZ_CP009273/3931365‑3931529                                                                                    |                                                                                  cAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGgaa                                                                             >  1:174426/1‑90 (MQ=255)   gcgcGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGgaaga                                                                           >  2:161645/1‑90 (MQ=255)                  ccGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCaga                                                            >  2:160169/1‑90 (MQ=255)                  ccGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCaga                                                            >  2:347788/1‑90 (MQ=255)                        gCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGAt                                                      <  1:22728/90‑1 (MQ=255)                        gCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGAt                                                      <  1:316839/90‑1 (MQ=255)                                                      tACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGc                        <  1:10889/90‑1 (MQ=255)                                                        ccTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGaaaa                                                                        >  1:99735/1‑40 (MQ=39)                                                        ccTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGaaaa                                                                        <  2:99735/40‑1 (MQ=39)                                                                  tGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaa            >  1:223049/1‑90 (MQ=255)                                                                  tGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaa            >  2:195289/1‑90 (MQ=255)                                                                            aCGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACaa  <  1:161645/90‑1 (MQ=255)                                                                            aCGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACaa  <  1:245029/90‑1 (MQ=255)                                                                                    |                                                                                  CAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAA  >  NZ_CP009273/3931365‑3931529

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