Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	950,036	(T)5→4	intergenic (‑114/+292)	focA ← / ← ycaO	formate transporter FocA/30S ribosomal protein S12 methylthiotransferase accessory factor YcaO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	950,032	0	T	.	100.0%	32.7 / NA	9	intergenic (‑110/+296)	focA/ycaO	formate transporter FocA/30S ribosomal protein S12 methylthiotransferase accessory factor YcaO
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (4/5);  total (4/5)

ACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTTGAA  >  NZ_CP009273/949950‑950106                                                                                   |                                                                           acgaGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGttttatgc                                                                     <  2:528660/90‑5 (MQ=255)      aTACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATa                                                               <  1:491502/90‑1 (MQ=255)       tACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATAc                                                              <  2:1553/90‑1 (MQ=255)       tACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATAc                                                              <  2:417622/90‑1 (MQ=255)               aaGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTaaa                                                      >  1:392726/1‑90 (MQ=255)                         gATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTa                                                        <  1:341777/78‑1 (MQ=255)                         gATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTa                                                        >  2:341777/1‑78 (MQ=255)                              gaAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATAc                                                              <  1:472285/67‑1 (MQ=255)                              gaAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATAc                                                              >  2:472285/1‑67 (MQ=255)                                                                   cccGTAACTCGCATGG‑TTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTTGaa  >  1:630171/1‑90 (MQ=255)                                                                                   |                                                                           ACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTTGAA  >  NZ_CP009273/949950‑950106

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