Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,421,451	A→G	K20R (AAG→AGG)	ulaB →	L‑ascorbate specific PTS enzyme IIB component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,421,451	0	A	G	100.0%	31.6 / NA	12	K20R (AAG→AGG)	ulaB	L‑ascorbate specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (6/6);  total (6/6)

TCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGGCTGTGTGTGGCAACGGACAAGGCAGTTCCATGATCATGAAGATGAAAGTGGACCAGTTTTTAACCCAATCAAACATTGACCATACGGTAAACAGCTGCGCGGTTGGCGAGTACAAAAGCGAGTTGAGTGGCGCGGATATCATCATCGCTTCTACGCACATTGCGGGCGAAATCACCGTGACCGGCAACAAATACGTGGTTGGCGTGCGCAACATGCTCTCTCCTGCCGACTTTGGCCCGAAACTGCTGGAAGTGATCAAAGAGC  >  NC_000913/4421310‑4421675                                                                                                                                              |                                                                                                                                                                                                                                 tcttcGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGGCTGTGTGTGGCAACGGACAAGGCAGTTCCATGATCATGAGGATGAAAGTGGACCAGTTTTTAACCCAATCAAACATTGACCATACGGTAAACAGCTGCGCGGTTGGc                                                                                                                                                               <  1:10708/209‑1 (MQ=255)                               aagaTGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGGCTGTGTGTGGCAACGGACAAGGCAGTTCCATGATCATGAGGATGAAAGTGGCCCAGTTTTTAACCCAATCAAACATTGACCATACGGTAAACAGCTGCGCGGTTGGCGAGTACAAAAGCg                                                                                                                                                  >  1:188676/1‑192 (MQ=255)                                    gCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGGCTGTGTGTGGCAACGGACAAGGCAGTTCCATGATCAGGAGGATGAAAGTGGACCAGTTTTTAACCCAATCAAACATTGACCATACGGTAAACAGCTGCGCGGTTGGCGAGTACAAAAGCGAGTTGAGTGGCGCGGATATCATCATCGCTTCTACGCACATTGCGGGCGAAATCACCGTGACCgn                                                                                  <  1:83147/251‑2 (MQ=255)                                                    naGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGGCTGTGTGTGGCAACGGACAAGGCAGTTCCATGATCATGAGGATGAAAGTGGACCAGTTTTTAACCCAATCAAACATTGACCATACGGTAAACAGCTGCGCGGTTGGCGAGTACAAAAGCGAGTTGAGTGGcgcg                                                                                                                                    >  1:93902/2‑185 (MQ=255)                                                                                                     ggCTGTGTGTGGCAACGGACAAGGCAGTGCCATGATCATTAGGATGAAAGTGGACCAGTTTTTAACCCAATCAAACAGTGACCATACGGTAAACAGCTGCGCGGTTGGCGAGTACAAAAGCGAGTTGAGTGGCGCGGATATCATCATCGCTTCTACGTACATTGCGGGCGAAATCACCTTGACAGGCAACAAATACTTGGTTGGCGTGCGCAACATGc                                                  <  1:220359/218‑1 (MQ=255)                                                                                                         gtgtgtgGCAACGGACAAGGCAGTTCCATGATCATGAGGATGAAAGTGTACCAGTTTTTAACCCAATCAAACATTGACCATACGGTAAACAGCTGCGCGGTTGGCGAGTACAAAAGCGAGTTGAGTGGCGCGGATATCATCATCGCTTCTACGCACATTGCGGGCGAAATCACCGTGACCGGCAACAAATACGTGGTTGGCGTGCGCAACATGCTCTCTCCTGCCGACTTTGGCCCGAGACTGCTGGAAg              <  1:88638/250‑1 (MQ=255)                                                                                                                  ggcGGACAAGGCAGTTACATGATCATGAGGATGAAAGTGGACCAGTTTTTAACCCAATCAAACATTGACCATACGGTAAACAGCTGCGCGGTTGGCGAGTACAAAAGCGAGTTGAGTGGCGCGGATATCATCATCGCTTCTACGCACATTGCGGGCGAAATCACCGTGACCGGc                                                                                 >  1:268795/3‑174 (MQ=255)                                                                                                                    nggACAAGGCAGTTCCATGATCATGAGGATGAAAGTGGACCAGTTTTTAACCCAATCAAACATTGACCATACGGTAAACAGCTGCGCGGTTGGCGAGTACAAAAGCGAGTTGAGTGGCGCGGATATCATCATCGCTTCTACGCACATTGCGGGCGAAATCACCGTGACCGGCAACAAATACGTGGTTGGCGTGCGCAACATGCTCTCTCCTGCCGACTTTGGCCCGAAACTGCTGGAAGTGATCAAAGAGc  >  1:202556/2‑251 (MQ=255)                                                                                                                     ggACAAGGCAGTTCCATGATCATGAGGATGAAAGTGGACCAGTTTTTAACCCAATCAAACATTGACCATACGGTAAACAGCTGCGCGGTTGGCGAGTACAAAAGCGAGTTGAGTGGCGCGGATATCATCATCGCTTCTACGcac                                                                                                            <  1:111997/144‑1 (MQ=255)                                                                                                                     ggACAAGGCAGTTCCATGATCATGAGGATGAAAGTGGACCAGTTTTTAACCCAATCAAACATTGACCATACGGTAAACAGCTGCGCGGTTGGCGAGTACAAAAGCGAGTTGAGTGGCGCGGATATCATCATCGCTTCTACGCACATTGCGGGCGAAATCACCGTGACCGGCAACAAATACGTGGTTGGCGTGCGCAACATGCTCTCTCCTg                                         >  1:89073/1‑211 (MQ=255)                                                                                                                               gTTCCNTGATCNTNNGGGTGAAAGTGGACCCGTTNTTAACCCAATCAAACNTTGANCATACNGTna                                                                                                                                                                                >  2:142035/1‑66 (MQ=255)                                                                                                                               gTTCCATGATCATGAGTATGAAAGTGGACCAGTTTTTAACCCAATCAAACATTGACCATACGGTAAACAGCTGCGCGGTTGTCGAGTACAAAAGCGATTTGAGTGGCGCGGAGATCATCATCGCTTCTACGCACATTGCGGGCGAAATCACCGGGACCGGCAACAAATACGTGGTTGGCGTGCGCAACATGc                                                  <  1:142035/192‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                                                                                                 TCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGGCTGTGTGTGGCAACGGACAAGGCAGTTCCATGATCATGAAGATGAAAGTGGACCAGTTTTTAACCCAATCAAACATTGACCATACGGTAAACAGCTGCGCGGTTGGCGAGTACAAAAGCGAGTTGAGTGGCGCGGATATCATCATCGCTTCTACGCACATTGCGGGCGAAATCACCGTGACCGGCAACAAATACGTGGTTGGCGTGCGCAACATGCTCTCTCCTGCCGACTTTGGCCCGAAACTGCTGGAAGTGATCAAAGAGC  >  NC_000913/4421310‑4421675

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