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%	93.8 / NA	28	E389D (GAA→GAC)	tar	methyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/15);  total (13/15)

AGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTC  >  NZ_CP009273/1965696‑1965868                                                                                       |                                                                                       aGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCa                                                                                     >  2:580959/1‑90 (MQ=255)            ctctTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCAcgcgc                                                                          >  2:918494/1‑90 (MQ=255)            ctctTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCAcgcgc                                                                          >  2:546573/1‑90 (MQ=255)               ttttGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGc                                                                       >  1:531601/1‑90 (MQ=255)                     cgccGGGGCGCTGCGACTGGCAAGATTCCGCCCTCCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAAc                                                                 <  2:305678/90‑1 (MQ=25)                        cTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACcgc                                                              >  1:256865/1‑90 (MQ=255)                                 gCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAg                                                     >  1:795118/1‑90 (MQ=255)                                   gccTGGCAAGATTCCGCTCTTCTCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAgcg                                                   <  1:201025/88‑1 (MQ=21)                                    aCTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAgcgc                                                  >  1:723110/1‑90 (MQ=255)                                     cTGGCAAGATTCCCCCCTCCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg                                                 <  2:119586/90‑1 (MQ=255)                                     cTGGCAAGATTACGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg                                                 <  1:498820/90‑1 (MQ=255)                                               tACGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAg                                       <  1:124254/90‑1 (MQ=255)                                                    aCTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTgg                                  <  1:723266/90‑1 (MQ=255)                                                      ttCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGaa                                >  2:468732/1‑90 (MQ=255)                                                       tCACCCGCCACCACGGCAAAACCACGGGCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAg                               >  1:480469/1‑90 (MQ=255)                                                       tCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAg                                               >  1:806149/1‑74 (MQ=255)                                                       tCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAg                                               <  2:806149/74‑1 (MQ=255)                                                           ccGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCaa                           <  2:531601/90‑1 (MQ=255)                                                               caccacGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAg                               <  2:642101/82‑1 (MQ=255)                                                               caccacGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAg                               >  1:642101/1‑82 (MQ=255)                                                               caccacGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATAcc                       >  1:50182/1‑90 (MQ=255)                                                                accacGGCAAAACCACGGCCCTGGTCCCCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCg                      <  1:920932/90‑1 (MQ=255)                                                                      gCAAAACCACGGCCCTGGTCGCCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATa                <  1:23785/90‑1 (MQ=255)                                                                      gCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATa                <  1:468732/90‑1 (MQ=255)                                                                        aaaaCCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAc              <  1:424066/90‑1 (MQ=255)                                                                        aaaaCCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAc              <  1:92515/90‑1 (MQ=255)                                                                        aaaaCCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAc              <  1:374473/90‑1 (MQ=255)                                                                                    cTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTc  >  1:897680/1‑90 (MQ=255)                                                                                       |                                                                                       AGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTC  >  NZ_CP009273/1965696‑1965868

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