Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,016,917	G→A	D122N (GAC→AAC)	fliI →	flagellar export ATPase FliI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,016,917	0	G	A	100.0%	24.0 / NA	11	D122N (GAC→AAC)	fliI	flagellar export ATPase FliI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (6/5);  total (6/5)

GCTCGGCGCAACCTGTGTCATTGAGCGCCAGAACGGCAGCGAAACGCACGAAGTAGAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGG  >  NC_000913/2016694‑2017132                                                                                                                                                                                                                                |                                                                                                                                                                                                                        gCTCGGCACAACCTGCGTCATTGAGCGCCAGAACGGCCGCGAAACGCACGAAGTAGAAAGCGAAGTCGGTGGCTTTACCGGTCAGCGGCTGCTTTTCATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCACGGCGCGCGTGTTTATGCCTAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGAACGGCAGCGGTAAAc                                                                                                                                                                                                          >  1:316689/1‑239 (MQ=255)         naaCCTGTGTCATTGAGCGCCAGAACGGCAGCGAAACGCACGAAGTAGAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCCAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGAACGGCAGCGGTAAACCGCTCTATGGCCTGCCCTcc                                                                                                                                                                                      >  1:192233/2‑251 (MQ=255)                          agcCAGACCGGCAGCGAAACGCACGAAGTAGAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGAACGGCAGCGGTAAACAGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAAc                                                                                                                                                                     >  1:86560/2‑251 (MQ=255)                                                      tAGAAAGCGAAGTAGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTAGAAGGTGTCATGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCAGGCGTTATTAGGTCGCGTTCTGAACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGCAACCGGTGCGCTGATTAc                                                                                                                                                      >  1:269134/1‑238 (MQ=255)                                                                                           gTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGNACGNNNNNggnnnnnnnnnn                                                                                                                                                                                                     >  2:199613/1‑144 (MQ=255)                                                                                                 aaTGCCGCTGGCGGAAGTCGAAGGTGTCCTGACCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGATGCAAAGCGGCAAGCAGTTGCCGCTAGGTCCGGCGTTATTAGGTCGCGTTCTGAACGGCAGCGGTAACCCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCACCGTACACCGCTTGAACATGTGCTGGACACCGGCg                                                                                              >  1:284421/1‑251 (MQ=255)                                                                                                   tGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGTGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGAACGGCAGCGGTAAACCg                                                                                                                                                                                                        <  1:212853/143‑1 (MQ=255)                                                                                                                                        gcgTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGAACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTgggn                                                       <  1:199613/251‑2 (MQ=255)                                                                                                                                                               nggCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGAACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGAtt                                                                                                                                                        >  1:203696/2‑131 (MQ=255)                                                                                                                                                               cGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGAACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGAtt                                                                                                                                                        <  2:203696/131‑1 (MQ=255)                                                                                                                                                                         gCTGCAAAGCGTCAATCAGTGGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGATCGGCAGCGTTAAACCGCTCGATGGCCTGCCCGCCCCCGAAGCGACGGAAATCGGTGCGCTGATTTCCCCGCCATTTAACCCGTTGCAACg                                                                                                                              <  1:290783/147‑1 (MQ=255)                                                                                                                                                                                              gCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGAACGGCAGCGGTAAACCGCTCGATGGCCGGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTgggg  <  1:270147/250‑1 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                        GCTCGGCGCAACCTGTGTCATTGAGCGCCAGAACGGCAGCGAAACGCACGAAGTAGAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGG  >  NC_000913/2016694‑2017132

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