Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,020,172	T→C	V29A (GTT→GCT)	fliM →	flagellar motor switch protein FliM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,020,172	0	T	C	100.0%	58.2 / NA	18	V29A (GTT→GCT)	fliM	flagellar motor switch protein FliM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/7);  total (11/7)

ATGTCACCGACGTGCTGTATACCGCTTTTATTCTGCGATAACGACATGGGCGATAGTATTCTTTCTCAAGCTGAAATTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGTTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGTTTGCAGGCGCTGGAAATCATTAATGAGCGCTTTGCCCGCCATTTTCGTATGGGG  >  NC_000913/2020042‑2020290                                                                                                                                   |                                                                                                                       aTGTCACCGACGTGCTGTATACCGCTTTTATTCTGCGATAACGACATGGGCGATAGTATTCTTTCTCAAGCTGAAATTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTa                                                                                                                      >  1:405404/1‑133 (MQ=255)                tGTATACCGCTTTTATTCTGCGATAACGACATGGGCGATAGTATTCTTTCTCAAGCTGAAATTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACa                                                                                                       <  1:368533/133‑1 (MQ=255)                tGTATACCGCTTTTATTCTGCGATAACGACATGGGCGATAGTATTCTTTCTCAAGCTGAAATTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACa                                                                                                       <  1:448412/133‑1 (MQ=255)                         cTTTTATTCTGCGATAACGACATGGGCGATAGTATTCTTTCTCAAGCTGAAATTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGt                                                                                              >  2:327009/1‑133 (MQ=255)                                    cGATAACGACATGGGCGATAGTATTCTTTCTCAAGCTGAAATTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCCGTGCTAGCGGCGTAAGCGACATTCGTCCGTCCGATCCGAAt                                                                                   >  2:200772/1‑133 (MQ=255)                                                                tctcAAGCTGAAATTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGtt                                                       >  1:172353/1‑133 (MQ=255)                                                                tctcAAGCTGAAATTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGtt                                                       >  1:317579/1‑133 (MQ=255)                                                                tctcAAGCTGAAATTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGtt                                                       >  1:410948/1‑133 (MQ=255)                                                                 ctcAAGCTGAAATTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGttt                                                      >  1:326961/1‑133 (MQ=255)                                                                   cAAGCTGAAATTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGTTTGc                                                    >  2:505397/1‑133 (MQ=255)                                                                       cTGAAATTCATGCGCTGTAGAAGGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGTTTGCAGgc                                                <  1:242063/133‑1 (MQ=255)                                                                         gAAATTGATGCGCTGCTGCATGGTGACAGCCAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGTTTGCAGgcgc                                              <  1:200772/133‑1 (MQ=255)                                                                           aaTTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGTTTGCAGGCGCTg                                            <  1:505397/133‑1 (MQ=255)                                                                           aaTTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGTTTGCAGGCGCTg                                            <  2:326961/133‑1 (MQ=255)                                                                             ttGATGCGCTATTGACTGGTGACATCGAAGTCTAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGTTTGCAGGCGCTGGa                                          <  2:172353/133‑1 (MQ=255)                                                                                                  aCAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCTAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGTTTGCAGGCGCTGGAAATCATTAATGAGCGCTTTGc                     >  2:432975/1‑133 (MQ=255)                                                                                                  aCAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGTTTGCAGGCGCTGGAAATCATTAATGAGCGCTTTGc                     >  1:407956/1‑133 (MQ=255)                                                                                                  aCAGCGAAGTCAAAGACGAACCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGTTTGCAGGCGCTGGAAAGCATTAATGAGCGCTTTGc                     >  1:176663/1‑133 (MQ=255)                                                                                                                    acccGACAGCCAGTGCTAGCGGCGAAAGTGAAATTCGCCCGTACGATCCAAATACCCAACGGGGGGTTGTGCGCGAACGTTTGCAGGCGCTGGAAATCATTAATGAGCGTTTTGCCCGCCATTTTCGTATggg   <  2:405404/131‑1 (MQ=255)                                                                                                                    aaCCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGTTTGCAGGCGCGGAAAATCATTAATGGGCGTTTTGCCCGCCATTTTCGTATggg   <  2:410948/133‑1 (MQ=255)                                                                                                                     aCCGACAGCCAGTGCTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGTTTGCAGGCGCTGGAAATCATTAATGAGCGCTTTGCCCGCCATTTTCGTATgggg  >  2:365370/1‑133 (MQ=255)                                                                                                                                   |                                                                                                                       ATGTCACCGACGTGCTGTATACCGCTTTTATTCTGCGATAACGACATGGGCGATAGTATTCTTTCTCAAGCTGAAATTGATGCGCTGTTGAATGGTGACAGCGAAGTCAAAGACGAACCGACAGCCAGTGTTAGCGGCGAAAGTGACATTCGTCCGTACGATCCGAATACCCAACGACGGGTTGTGCGCGAACGTTTGCAGGCGCTGGAAATCATTAATGAGCGCTTTGCCCGCCATTTTCGTATGGGG  >  NC_000913/2020042‑2020290

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