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%	37.6 / NA	10	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 . (5/5);  total (5/5)

TACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATT  >  NZ_CP009273/949956‑950091                                                                             |                                                            tACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATAc                                               <  1:210867/90‑1 (MQ=255) tACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATAc                                               <  1:402553/90‑1 (MQ=255)      aaaaaGAATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAATTCGCATGG‑‑TTTATGCAAATACGGAGTa                                         <  2:88885/88‑1 (MQ=255)       cAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTa                                         <  1:52409/90‑1 (MQ=255)                                    gATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAaaataaat            >  1:106228/1‑90 (MQ=255)                                    gATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAaaataaat            >  1:107427/1‑90 (MQ=255)                                        atatCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTa        >  1:4978/1‑90 (MQ=255)                                              cGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTa                                         >  1:364105/1‑51 (MQ=255)                                              cGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTa                                         <  2:364105/51‑1 (MQ=255)                                              cGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAAtt  >  1:258443/1‑90 (MQ=255)                                                                             |                                                            TACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATT  >  NZ_CP009273/949956‑950091

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