Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,012,918	A→G	D351G (GAC→GGC)	fliK →	flagellar hook‑length control protein FliK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,012,918	0	A	G	100.0%	93.3 / NA	28	D351G (GAC→GGC)	fliK	flagellar hook‑length control protein FliK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/16);  total (12/16)

GAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGT  >  NZ_CP009273/2012833‑2013004                                                                                      |                                                                                       gAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAggcgac                                                                                    >  2:665422/1‑90 (MQ=255)      cTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATAc                                                                               <  2:1111425/90‑1 (MQ=255)      cTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATAc                                                                               <  1:677330/90‑1 (MQ=255)        ttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGc                                                                             <  2:581953/90‑1 (MQ=255)        ttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGc                                                                             >  2:264744/1‑90 (MQ=255)              gTCAGCAGCAGGCCGCTTCCCAGCAACAGCGAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCgg                                                                       >  1:803973/1‑90 (MQ=255)                cagcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGtt                                                                     >  1:634934/1‑90 (MQ=255)                  gcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTcc                                                                   <  1:842994/90‑1 (MQ=255)                      cagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtc                                                               <  2:935678/90‑1 (MQ=255)                      cagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtc                                                               <  2:199806/90‑1 (MQ=255)                      cagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtc                                                               <  2:1049923/90‑1 (MQ=255)                      cagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtc                                                               <  1:238113/90‑1 (MQ=255)                        ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc                                                             <  1:1020931/90‑1 (MQ=255)                        ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc                                                             >  1:135010/1‑90 (MQ=255)                        ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc                                                             <  1:196669/90‑1 (MQ=255)                        ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc                                                             <  1:415859/90‑1 (MQ=255)                        ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc                                                             <  2:377860/90‑1 (MQ=255)                               tCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACaa                                                      <  1:709933/90‑1 (MQ=255)                                       acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCgtgt                                              <  1:218574/90‑1 (MQ=255)                                       acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTCACAAGGGCgtgt                                              >  2:916672/1‑90 (MQ=255)                                          gcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAAc                                           >  1:498043/1‑90 (MQ=255)                                          gcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAAc                                           <  2:822648/90‑1 (MQ=255)                                                        caGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgcggcg                             >  1:760193/1‑90 (MQ=255)                                                                         ctGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTaa            >  2:284307/1‑90 (MQ=255)                                                                         ctGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTaa            >  1:395229/1‑90 (MQ=255)                                                                              gggggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCa       <  2:728543/90‑1 (MQ=255)                                                                                gggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCaga     >  2:139849/1‑90 (MQ=255)                                                                                gggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCaga     >  1:489300/1‑90 (MQ=255)                                                                                   aaGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGt  >  2:170440/1‑90 (MQ=255)                                                                                      |                                                                                       GAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGT  >  NZ_CP009273/2012833‑2013004

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