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%	111.6 / NA	29	L265F (TTG→TTT)	rbsK	ribokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/14);  total (15/14)

AGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAAC  >  NZ_CP009273/3931366‑3931530                                                                                   |                                                                                   aGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGGGCGTTAATCACGGCATTTCTGgaag                                                                             >  1:611244/1‑90 (MQ=255)    gcgTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGaaa                                                                          <  2:361832/90‑1 (MQ=255)            ggATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTg                                                                  >  1:65675/1‑90 (MQ=255)                 ccGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCaga                                                             >  1:22191/1‑90 (MQ=255)                       gCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGAt                                                       <  2:747037/90‑1 (MQ=255)                        cAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGAtt                                                      <  1:581947/90‑1 (MQ=255)                                    aCCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGga                                                                               >  1:692619/1‑53 (MQ=255)                                    aCCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGga                                                                               <  2:692619/53‑1 (MQ=255)                                    aCCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAt                                          >  1:24921/1‑90 (MQ=255)                                        ttGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAtgctg                                      >  1:216402/1‑90 (MQ=255)                                        ttGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAtgctg                                      >  2:630603/1‑90 (MQ=255)                                        ttGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAtgctg                                      >  1:306925/1‑90 (MQ=255)                                        ttGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAtgctg                                      >  2:404729/1‑90 (MQ=255)                                        ttGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAtgctg                                      >  2:129532/1‑90 (MQ=255)                                                    aTACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTg                          <  2:611244/90‑1 (MQ=255)                                                          ttAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTaa                    <  2:357531/90‑1 (MQ=255)                                                          ttAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTaa                    >  2:109512/1‑90 (MQ=255)                                                          ttAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTaa                    <  2:100457/90‑1 (MQ=255)                                                                 tGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaa             >  2:250137/1‑90 (MQ=255)                                                                 tGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaa             >  2:73117/1‑90 (MQ=255)                                                                 tGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaa             >  1:302619/1‑90 (MQ=255)                                                                  gCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaaa            <  2:22191/90‑1 (MQ=255)                                                                   cGTTAATCGCGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAg           <  2:218809/90‑1 (MQ=255)                                                                         tCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGcac     <  2:70238/90‑1 (MQ=255)                                                                         tCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGcac     <  1:161530/90‑1 (MQ=255)                                                                         tCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGcac     <  2:773390/90‑1 (MQ=255)                                                                           aCGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACaa   <  2:252579/90‑1 (MQ=255)                                                                           aCGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACaa   <  2:16035/90‑1 (MQ=255)                                                                            cGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAAc  >  2:563921/1‑90 (MQ=255)                                                                                   |                                                                                   AGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAAC  >  NZ_CP009273/3931366‑3931530

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