Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,015,409	A→G	E177G (GAA→GGA)	fliG →	flagellar motor switch protein FliG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,015,409	0	A	G	100.0%	125.1 / NA	36	E177G (GAA→GGA)	fliG	flagellar motor switch protein FliG
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/15);  total (21/15)

ACGTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGAAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACGGCAGCCGAAAT  >  NC_000913/2015332‑2015484                                                                              |                                                                            aCGTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGt                                                                          <  2:358317/81‑1 (MQ=255)   gTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTAc                                                                        <  1:358318/81‑1 (MQ=255)   gTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTAc                                                                        >  2:358349/1‑81 (MQ=255)   gTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTAc                                                                        >  1:358360/1‑81 (MQ=255)    tCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACt                                                                       <  2:358319/81‑1 (MQ=255)       gcgcCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGaa                                                                    >  2:358357/1‑81 (MQ=255)          ccACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATgg                                                                 >  1:358356/1‑81 (MQ=255)          ccACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATgg                                                                 >  1:358359/1‑81 (MQ=255)               acgTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGc                                                            >  2:358379/1‑81 (MQ=255)                 gTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTc                                                          <  2:358320/81‑1 (MQ=255)                   gATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGa                                                        >  2:358369/1‑81 (MQ=255)                   gATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGa                                                        >  1:358362/1‑81 (MQ=255)                     tGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACg                                                      >  1:358370/1‑81 (MQ=255)                      gTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACgg                                                     >  2:358364/1‑81 (MQ=255)                       ttGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGt                                                    <  2:358321/81‑1 (MQ=255)                         gCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTc                                                   <  1:358322/80‑1 (MQ=255)                            tATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGaa                                               >  2:358373/1‑81 (MQ=255)                              tCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAAt                                              <  1:358323/80‑1 (MQ=255)                               cGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAAtct                                            >  1:358365/1‑81 (MQ=255)                                gCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAAtct                                            >  2:358366/1‑80 (MQ=255)                                  cACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAAtctc                                           >  1:358354/1‑79 (MQ=255)                                        tGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAgcgc                                     <  1:358324/79‑1 (MQ=255)                                         ggcggcGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCgcag                                   <  1:358325/80‑1 (MQ=255)                                                    gccagccGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGggcgg                       <  2:358326/81‑1 (MQ=255)                                                       agccGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGt                    <  2:358327/81‑1 (MQ=255)                                                           gcgcTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTgaga                <  2:358329/81‑1 (MQ=255)                                                           gcgcTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTgag                 >  1:358375/1‑80 (MQ=255)                                                           gcgcTGGCGGAGCTGACCGGAGTACTGAATGACTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTgaga                <  2:358328/81‑1 (MQ=255)                                                            cgcTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGaa               >  2:358381/1‑81 (MQ=255)                                                             gcTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAAc              <  1:358330/81‑1 (MQ=255)                                                               tGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACgg            >  1:358378/1‑81 (MQ=255)                                                               tGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACgg            >  2:358374/1‑81 (MQ=255)                                                                     aGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACGGCAGcc       >  1:358367/1‑80 (MQ=255)                                                                      gCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACGGCAGCCg      >  1:358376/1‑80 (MQ=255)                                                                         gACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACGGCAGCCGAAAt  >  2:358372/1‑81 (MQ=255)                                                                          aCCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACGGCAGCCGaaa   <  1:358331/79‑1 (MQ=255)                                                                              |                                                                            ACGTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGAAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACGGCAGCCGAAAT  >  NC_000913/2015332‑2015484

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