Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	635,213	A→G	K186K (AAA→AAG)	SAUSA300_RS02990 →	acetyl‑CoA acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	635,213	0	A	G	100.0%	43.5 / NA	17	K186K (AAA→AAG)	SAUSA300_RS02990	acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/1);  total (16/1)

GAAAATGTGGCCAAGATGTATGATGTTTCAAGAGAATTACAAGATGAATTTGCTTATCGAAGTCATCAATTGACAGCGGAAAATGTAAAAAATGGAAATATTTCTCAGGAAATATTACCTATAACCGTTAAAGGAGAAATATTCAACACTGATGAAAGTCTAAAATCACAT  >  NC_007793/635124‑635294                                                                                          |                                                                                  gAAAATGTGGCCAAGATGTATGATGTTTCAAGAGAATTACAAGATGAATTTGCTTATCGAAGTCATCAATTGACAGCGGAAAATGTAAAGAATGGAAata                                                                         >  1:3641692/1‑100 (MQ=255)       ttgGCCAAGATGTATGATGTTTCAAGAGAATTACAAGATGAATTTGCTTATCGAAGTCATCAATTGACAGCGGAAAATGTAAAGAATGGAAATATTtctc                                                                   >  1:2799078/2‑100 (MQ=255)       gtgGCCAAGATGTATGATGTTTCAAGAGAATTACAAGATGAATTTGCTTATCGAAGTCATCAATTGACAGCGGAAAATGTAAAGAATGGAAagatcggaa                                                                   >  1:2893456/1‑95 (MQ=255)       gtgGCCAAGATGTATGATGTTTCAAGAGAATTACAAGATGAATTTGCTTATCGAAGTCATCAATTGACAGCGGAAAATGTAAAGAATGGAAATATTtctc                                                                   >  1:281402/1‑100 (MQ=255)       ctgGCCAAGATGTATGATGTTTCAAGAGAATTACAAGATGAATTTGCTTATCGAAGTCATCAATTGACAGCGGAAAATGTAAAGAATGGAAATATTtctc                                                                   >  1:2514803/2‑100 (MQ=255)         tgCCAAGATGTATGATGTTTCAAGAGAATTACAAGATGAATTTGCTTATCGAAGTCATCAATTGACAGCGGAAAATGTAAAGAATGGAAATATTTCTCAg                                                                 >  1:5077582/2‑100 (MQ=255)         ggCCAAGATGTATGATGTTTCAAGAGAATTACAAGATGAATTTGCTTATCGAAGTCATCAATTGACAGCGGAAAATGTAAAGAATGGAAATATTTCTCAg                                                                 >  1:1291276/1‑100 (MQ=255)         ggCCAAGATGTATGATGTTTCAAGAGAATTACAAGATGAATTTGCTTATCGAAGTCATCAATTGACAGCGGAAAATGTAAAGAATGGAAATATTTCTCAg                                                                 >  1:4954431/1‑100 (MQ=255)         ggCCAAGATGTATGATGTTTCAAGAGAATTACAAGATGAATTTGCTTATCGAAGTCATCAATTGACAGCGGAAAATGTAAAGAATGGAAATATTTCTCAg                                                                 >  1:3488645/1‑100 (MQ=255)                  gTATGATGTTTCAAGAGAATTACAAGATGAATTTGCTTATCGAAGTCATCAATTGACAGCGGAAAATGTAAAGAATGGAAATATTTCTCAGGAAATATTa                                                        >  1:2112044/1‑100 (MQ=255)                                               caTTTGCTTATCGAAGTCATCAATTGACAGCGGAAAATGTAAAGAATGGAAATATTTCTCAGGAAATATTACCTATAACCGTTAAAGTAGAAATATTCaa                           >  1:2460848/2‑100 (MQ=255)                                                                   caATTGACAGCGGAAAATGTAAAGAATGGAAATATTTCTCAGGAAATATTACCTATAACCGTTAAAGGAGAAATATTCAACACTGATGAAAGTCTAAAAt       >  1:261617/1‑100 (MQ=255)                                                                       tGACAGCGGAAAATGTAAAGAATGGAAATATTTCTCAGGAAATATTACCTATAACCGTTAAAGGAGAAATATTCAACACTGATGAAAGTCTAAAATcaca   >  1:1604366/1‑100 (MQ=255)                                                                        gACAGCGGAAAATGTAAAGAATGGAAATATTTCTCAGGAAATATTACCTATAACCGTTAAAGGAGAAATATTCAACACTGATGAAAGTCTAAAATCACat  >  1:4875351/1‑100 (MQ=255)                                                                        gACAGCGGAAAATGTAAAGAATGGAAATATTTCTCAGGAAATATTACCTATAACCGTTAAAGGAGAAATATTCAACACTGATGAAAGTCTAAAATCACat  >  1:1773831/1‑100 (MQ=255)                                                                        gACAGCGGAAAATGTAAAGAATGGAAATATTTCTCAGGAAATATTACCTATAACCGTTAAAGGAGAAATATTCAACACTGATGAAAGTCTAAAATCACat  >  1:1750502/1‑100 (MQ=255)                                                                        gACAGCGGAAAATGTAAAGAATGGAAATATTTCTCAGGAAATATTACCTATAACCGTTAAAGGAGAAATATTCAACACTGATGAAAGTCTAAAATCACat  <  1:5267123/100‑1 (MQ=255)                                                                                          |                                                                                  GAAAATGTGGCCAAGATGTATGATGTTTCAAGAGAATTACAAGATGAATTTGCTTATCGAAGTCATCAATTGACAGCGGAAAATGTAAAAAATGGAAATATTTCTCAGGAAATATTACCTATAACCGTTAAAGGAGAAATATTCAACACTGATGAAAGTCTAAAATCACAT  >  NC_007793/635124‑635294

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