Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,965,782	T→G	E389D (GAA→GAC)	tar ←	methyl‑accepting chemotaxis protein II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,965,782	0	T	G	100.0%	33.8 / NA	12	E389D (GAA→GAC)	tar	methyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (5/7);  total (5/7)

GAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAC  >  NZ_CP009273/1965695‑1965856                                                                                        |                                                                           gAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTc                                                                          <  2:321386/90‑1 (MQ=255)  aGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCa                                                                         >  2:410862/1‑90 (MQ=255)  aGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCa                                                                         >  2:859302/1‑90 (MQ=255)                       gccTGGGCGCTGGGACTGGCAAGATTACGCACTTTACCCCCCACCACGGCAAAACCACGGCCCTGGGCACCCCCACGCGCGGCCTCAAcc                                                    >  1:1371/1‑90 (MQ=21)                                      cgggCAAGATTCCGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg                                     <  1:862936/88‑1 (MQ=255)                                      cTGGCAAGATTCCGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg                                     <  2:1371/90‑1 (MQ=255)                                      cTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg                                     <  1:859302/90‑1 (MQ=255)                                        ggCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAg                                   >  2:23587/1‑90 (MQ=255)                                                     aCTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTgg                      <  1:545790/90‑1 (MQ=255)                                                                 accacGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCg          <  2:470001/90‑1 (MQ=255)                                                                  ccacGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGGGGGTATTAGTCTGGGAGGCAATACCGt         >  1:691949/1‑90 (MQ=25)                                                                         aaaaCCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAc  <  2:402923/90‑1 (MQ=255)                                                                                        |                                                                           GAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAC  >  NZ_CP009273/1965695‑1965856

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