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%	42.3 / NA	11	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 . (7/4);  total (7/4)

AAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTT  >  NZ_CP009273/949963‑950103                                                                      |                                                                        aaaGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTaa                                                    >  1:324575/1‑90 (MQ=255)    gCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTa                                                     >  1:128511/1‑86 (MQ=255)    gCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTa                                                     <  2:128511/86‑1 (MQ=255)              tGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTAAATATTTGATTATc                                       >  2:916225/1‑90 (MQ=255)                                  tatCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTa                                                     <  1:989702/56‑1 (MQ=255)                                  tatCTCGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTa                                                     >  2:989702/1‑56 (MQ=255)                                      tcGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAAt               >  1:548493/1‑90 (MQ=255)                                       cGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACgg                                                         <  1:370908/47‑1 (MQ=255)                                       cGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACgg                                                         >  2:370908/1‑47 (MQ=255)                                       cGCCTGGCTTATAGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAAtt              >  1:304252/1‑90 (MQ=255)                                                   aGGCCCGTAACTCGCATGG‑TTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGtt  <  1:168790/90‑1 (MQ=255)                                                                      |                                                                        AAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTT  >  NZ_CP009273/949963‑950103

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