Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,021,189	T→G	A34A (GCT→GCG)	fliN →	flagellar motor switching and energizing component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,021,189	0	T	G	100.0%	161.9 / NA	45	A34A (GCT→GCG)	fliN	flagellar motor switching and energizing component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (28/17);  total (28/17)

ATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGC  >  NC_000913/2021042‑2021334                                                                                                                                                    |                                                                                                                                                  aTTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGcc                                                                                                                                                <  2:847705/151‑1 (MQ=255)  tttGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCg                                                                                                                                               >  1:205770/1‑151 (MQ=255)     gATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCg                                                                                                                                               >  1:882671/1‑148 (MQ=255)      aTTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGAc                                                                                                                                           >  2:516354/1‑151 (MQ=255)      aTTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGAc                                                                                                                                           >  2:812643/1‑151 (MQ=255)        tAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACgg                                                                                                                                         >  1:324553/1‑151 (MQ=255)         aaCCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACgg                                                                                                                                         >  1:665902/1‑150 (MQ=255)           cccGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGAGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGgtg                                                                                                                                       >  1:1028/1‑150 (MQ=255)             cGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTc                                                                                                                                    <  1:889548/151‑1 (MQ=255)               aTTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTcc                                                                                                                                   >  2:359199/1‑150 (MQ=255)               aTTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCa                                                                                                                                  <  2:719868/151‑1 (MQ=255)                      aTTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAAttt                                                                                                                           <  2:578471/151‑1 (MQ=255)                                 gAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGgtggtg                                                                                                                 >  1:270762/1‑150 (MQ=255)                                  aGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGAGgtggtggt                                                                                                               >  2:619131/1‑149 (MQ=255)                                         gCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCg                                                                                                        <  1:167012/151‑1 (MQ=255)                                             aaaTGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAAc                                                                                                    >  2:658966/1‑151 (MQ=255)                                              aaTGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACg                                                                                                   >  1:496737/1‑151 (MQ=255)                                               aTGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGt                                                                                                  <  1:249853/151‑1 (MQ=255)                                               aTGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGt                                                                                                  >  2:41084/1‑151 (MQ=255)                                                tGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGtt                                                                                                 >  2:843510/1‑151 (MQ=255)                                                 gAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGtt                                                                                                 <  1:524269/150‑1 (MQ=255)                                                    tGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGc                                                                                               >  1:715386/1‑149 (MQ=255)                                                     gACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAg                                                                                             <  2:712046/150‑1 (MQ=255)                                                       cATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGa                                                                                           <  1:439046/150‑1 (MQ=255)                                                            ataatCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGa                                                                                     <  1:115122/151‑1 (MQ=255)                                                               atCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCt                                                                                  <  2:258842/151‑1 (MQ=255)                                                                    gCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTa                                                                             <  2:124676/151‑1 (MQ=255)                                                                      cGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATg                                                                           >  2:376928/1‑151 (MQ=255)                                                                       gATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATgg                                                                          >  1:702554/1‑151 (MQ=255)                                                                            caacaaCGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATAtt                                                                     <  1:876383/151‑1 (MQ=255)                                                                             aacaacGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTc                                                                    >  1:585565/1‑151 (MQ=255)                                                                                aacGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCgg                                                                 >  2:737090/1‑151 (MQ=255)                                                                                   ggCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTc                                                               >  2:875163/1‑150 (MQ=255)                                                                                        aaTGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTg                                                         <  2:676666/151‑1 (MQ=255)                                                                                               gATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCg                                                  >  2:905099/1‑151 (MQ=255)                                                                                                    gtgGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCg                                                  <  2:916407/146‑1 (MQ=255)                                                                                                    gtgGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTg                                             >  1:44588/1‑151 (MQ=255)                                                                                                     tgGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTgg                                            >  1:696004/1‑151 (MQ=255)                                                                                                     tgGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTg                                             <  2:257840/150‑1 (MQ=255)                                                                                                                          cAAAAATCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCa                       <  1:432512/151‑1 (MQ=255)                                                                                                                             aaaTCAACCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCaa                      >  2:877956/1‑149 (MQ=255)                                                                                                                                   aCCAGCAGCAAAAGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGt              >  1:359237/1‑151 (MQ=255)                                                                                                                                            aaaaGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTg      >  1:1358/1‑150 (MQ=255)                                                                                                                                            aaaaGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTg      >  2:606976/1‑150 (MQ=255)                                                                                                                                               aGCGCGGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGc  >  1:294568/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                  ATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTGGGCTGAAGCGTTGAGCGAACAAAAATCAACCAGCAGCAAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGC  >  NC_000913/2021042‑2021334

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