Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,684,806	T→C	intergenic (‑349/+209)	qseE ← / ← purL	two component system sensor histidine kinase QseE/GlrK/phosphoribosylformylglycinamidine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,684,806	0	T	C	100.0%	43.5 / NA	13	intergenic (‑349/+209)	qseE/purL	two component system sensor histidine kinase QseE/GlrK/phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/5);  total (8/5)

TGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAATCAGCCACTGTGTCGCTAAAAAGAGACAACTTAAGATAAACTTATTAGATAATATTAAAATCAATGAGTTAAGTGTCGCCAGAAAGCGACACGGCAAACCACCCATTGTCGTGATTTACAGACACAAAAA  >  NZ_CP009273/2684691‑2684935                                                                                                                    |                                                                                                                                  tGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAACCAGCCACTGTGTCGCTaaaaa                                                                                                              >  2:423600/1‑137 (MQ=255)   aGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAACCAGCCACTGTGTCGCTAAAAAgaga                                                                                                          <  1:139248/139‑1 (MQ=255)   aGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAACCAGCCACTGTGTCGCTAAAAAgaga                                                                                                          <  2:234411/139‑1 (MQ=255)                     tGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAACCAGCCACTGTGTCGCTAAAAAGAGACAACTTAAGATAAACtta                                                                                        >  2:233207/1‑139 (MQ=255)                        aGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAACCAGCCACTGTGTCGCTAAAAAGAGACAACTTAAGATAAACttat                                                                                       >  2:301787/1‑137 (MQ=255)                          tAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAACCAGCCACTGTGTCGCTAAAAAGAGACAACTTAAGATAAACttat                                                                                       >  2:313621/1‑135 (MQ=255)                                  ccAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAACCAGCCACTGTGTCGCTAAAAAGAGACAACTTAAGATAAACTTATTAGATAATAtt                                                                            >  1:299096/1‑138 (MQ=255)                                         aGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAACCAGCCACTGTGTCGCTAAAAAGAGACAACTTAAGATAAACTTATTAGATAATATTAAAATCaa                                                                    <  1:233207/139‑1 (MQ=255)                                                  tttGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAACCAGCCACTGTGTCGCTAAAAAGAGACAACTTAAGATAAACTTATTAGATAATATTAAAATCAATGAGTTAAg                                                           >  1:18670/1‑139 (MQ=255)                                                     ggTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAACCAGCCACTGTGTCGCTAAAAAGAGACAACTTAAGATAAACTTATTAGATAATATTAAAATCAATGAGTTAAgtgt                                                        >  2:113361/1‑139 (MQ=255)                                                                                            tctcCTCCGCTGTTTGAATGTAAACCAGCCACTGTGTCGCTAAAAAGAGACAACTTAAGATAAACTTATTAGATAATATTAAAATCAATGAGTTAAGTGTCGCCAGAAAGCGACACGGCAAACCACCCATTGTCGTGa                  <  1:113361/138‑1 (MQ=255)                                                                                                ctccGCTGTTTGAATGTAAACCAGCCACTGTGTCGCTAAAAAGAGACAACTTAAGATAAACTTATTAGATAATATTAAAATCAATGAGTTAAGTGTCGCCAGAAAGCGACACGGCAAACCACCCATTGTCGTGATTTa              >  1:1113/1‑138 (MQ=255)                                                                                                            aaTGTAAACCAGCCACTGTGTCGCTAAAAAGAGACAACTTAAGATAAACTTATTAGATAATATTAAAATCAATGAGTTAAGTGTCGCCAGAAAGCGACACGGCAAACCACCCATTGTCGTGATTTACAGACACaaaaa  <  2:255906/138‑1 (MQ=255)                                                                                                                    |                                                                                                                                  TGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAATCAGCCACTGTGTCGCTAAAAAGAGACAACTTAAGATAAACTTATTAGATAATATTAAAATCAATGAGTTAAGTGTCGCCAGAAAGCGACACGGCAAACCACCCATTGTCGTGATTTACAGACACAAAAA  >  NZ_CP009273/2684691‑2684935

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