Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,003,264	C→T	E115K (GAA→AAA)	fliC ←	flagellar filament structural protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,003,264	0	C	T	100.0%	109.7 / NA	30	E115K (GAA→AAA)	fliC	flagellar filament structural protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/16);  total (14/16)

AAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATAC  >  NC_000913/2003039‑2003497                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                          aaGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGAc                                                                                                                                                                                                                  <  1:9049/251‑1 (MQ=255)         ccTGGAGCAGTAGTGATAACTGTAGCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGGGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATGTCAGCCAGACGGGATTTAATTTTGTCCTGGAGAGAAGACAGATCAGACTCAGAGtt                                                                                                                                                                                                          <  2:308308/250‑1 (MQ=255)               gCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTAcc                                                                                                                                                                                                    <  1:191374/251‑1 (MQ=255)                                                               ccAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTGCCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAgttgttgtt                                                                                                                                                    <  2:263376/251‑1 (MQ=255)                                                                               cAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAgcgcg                                                                                                                                    <  2:8732/251‑1 (MQ=255)                                                                                     gCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCTCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGTACAGCGCGCCTTc                                                                                                                               <  1:380973/250‑1 (MQ=255)                                                                                      cTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCt                                                                                                                                                                                         >  2:57600/1‑191 (MQ=255)                                                                                      cTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCt                                                                                                                                                                                         <  1:57600/191‑1 (MQ=255)                                                                                            aTCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCt                                                                                                                       >  2:306407/1‑251 (MQ=255)                                                                                                   gtgatagtCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAAc                                                                                                                >  2:78084/1‑251 (MQ=255)                                                                                                          tCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCCGACGGGATTTAATTTTATCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGAt                                                                                                          >  1:412587/1‑250 (MQ=255)                                                                                                          tCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGcaca                                                                                                                                                                     <  2:18127/191‑1 (MQ=255)                                                                                                          tCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGcaca                                                                                                                                                                     >  1:18127/1‑191 (MQ=255)                                                                                                                  tCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGTAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGtt                                                                                                 >  2:280281/1‑251 (MQ=255)                                                                                                                     tttGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGc                                                                                              >  1:185198/1‑251 (MQ=255)                                                                                                                               ccTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGATTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGgcc                                                                                    >  2:284883/1‑251 (MQ=255)                                                                                                                                        tCATGGAGCCATTTTTTGCCAGCACTTTCACGCCGTTGAACTGGGTCTTACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTGGATTTCGGACAGCGCTCCTGCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGt                                                                            <  1:237782/250‑1 (MQ=255)                                                                                                                                         cATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTc                                                                           >  2:99403/1‑250 (MQ=255)                                                                                                                                           tGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGCTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAgg                                                                        >  1:264062/1‑251 (MQ=255)                                                                                                                                              ggCCATTTTTTGCCAGCACGTTCCCGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTGTGTCCTGGATAGAAGACATATCAGACTCAGAGTTAGTTCCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCTGTTGTCTTCGCAATGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCt                                                                     <  1:17671/250‑1 (MQ=255)                                                                                                                                                ccATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCGTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCttt                                                                   <  2:185198/251‑1 (MQ=255)                                                                                                                                                 cATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTa                                                                  <  1:284883/251‑1 (MQ=255)                                                                                                                                                       tttccAGCCCGTTCACTCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGtt                                                             <  2:3660/247‑1 (MQ=255)                                                                                                                                                                          gTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGCTTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGAc                                                                                                                                          >  1:231705/1‑154 (MQ=255)                                                                                                                                                                          gTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGGGGCCTGTACCGTCAGTTCACGCACACGCGGTAAGTTGTTGTTGATTTCGGAc                                                                                                                                          <  2:231705/154‑1 (MQ=255)                                                                                                                                                                                       tGACCAGATACGAGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGAACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGACTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGc                            >  1:397794/1‑251 (MQ=255)                                                                                                                                                                                                  gcgGTCAATTTCATCCAGACGGGATTTAATGTTGTCCGGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCt                 <  2:381688/251‑1 (MQ=255)                                                                                                                                                                                                             cATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCt                                                                                                                       >  2:13472/1‑138 (MQ=255)                                                                                                                                                                                                             cATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCt                                                                                                                       <  1:13472/138‑1 (MQ=255)                                                                                                                                                                                                                 cAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATAc  >  1:153705/1‑251 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                          AAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATAC  >  NC_000913/2003039‑2003497

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