Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,015,254	G→T	E125D (GAG→GAT)	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,254	0	G	T	100.0%	121.8 / NA	37	E125D (GAG→GAT)	fliG	flagellar motor switch protein FliG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (19/18);  total (19/18)

TATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGT  >  NC_000913/2015118‑2015379                                                                                                                                         |                                                                                                                              tATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCAt                                                                                                                            >  1:476568/1‑140 (MQ=255)         aaaGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAAt                                                                                                                    >  1:638284/1‑140 (MQ=255)           aGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTa                                                                                                                  >  2:795210/1‑140 (MQ=255)           aGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTa                                                                                                                  >  2:280377/1‑140 (MQ=255)           aGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTa                                                                                                                  >  1:433007/1‑140 (MQ=255)            gCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTAt                                                                                                                 <  1:373492/140‑1 (MQ=255)            gCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTAt                                                                                                                 <  2:67132/140‑1 (MQ=255)             ctctGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATc                                                                                                                <  2:370757/140‑1 (MQ=255)                 ggctGAAGAACGTGCCACCAGCCTGCTGGAAAATAATCTAGAAATTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGcca                                                                                                            <  2:600063/137‑1 (MQ=255)                            gTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTg                                                                                                 <  1:68074/140‑1 (MQ=255)                            gTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTg                                                                                                 >  1:94354/1‑140 (MQ=255)                             tGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGc                                                                                                >  2:330194/1‑140 (MQ=255)                                cgccAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATc                                                                                             >  2:539678/1‑140 (MQ=255)                                                   tttCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAgccgc                                                                          <  1:539678/139‑1 (MQ=255)                                                    ttCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAgccgc                                                                          >  1:73500/1‑139 (MQ=255)                                                    ttCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAgccgc                                                                          <  2:73500/139‑1 (MQ=255)                                                     tcttGAATCTGGCGATACCGCCAGCGGTATTGAAACGCACAACTTTATGGGTCAAAAGTGCGCCGCCGATCTGATTCGCGATGATTATCGGCAAATTATCCCCCCCATTCTGGTGCATCTGAAGCGCGCCTAAGccgccg                                                                        <  2:288005/140‑1 (MQ=255)                                                     tctcGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGccgccg                                                                        <  1:795210/140‑1 (MQ=255)                                                       tcGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGat                                                                      >  1:748328/1‑140 (MQ=255)                                                       tcGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGat                                                                      >  1:606350/1‑140 (MQ=255)                                                          aaaCTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATAtt                                                                   <  2:491500/140‑1 (MQ=255)                                                           aaCTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTc                                                                  <  2:433007/140‑1 (MQ=255)                                                                gcgATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTTGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCg                                                             <  1:512273/140‑1 (MQ=255)                                                                     cccGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCgttgtt                                                        <  2:13509/139‑1 (MQ=255)                                                                                     aaaCGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTgcgc                                        <  1:176757/140‑1 (MQ=255)                                                                                      aaCGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGcc                                       <  1:474184/140‑1 (MQ=255)                                                                                       aCGCTCAACTTTATGGAGCCCCAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCa                                      <  2:53595/140‑1 (MQ=255)                                                                                       aCGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCa                                      <  2:564760/140‑1 (MQ=255)                                                                                          cTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCacga                                   >  2:827705/1‑140 (MQ=255)                                                                                            cAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCacgacg                                 <  1:86712/140‑1 (MQ=255)                                                                                               cTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCACGACGTGa                              >  2:834883/1‑140 (MQ=255)                                                                                                    tGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCACGACGTGATGTTg                         >  1:411888/1‑140 (MQ=255)                                                                                                         ccACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCACGACGTGATGTTGCGTAt                    >  1:174437/1‑140 (MQ=255)                                                                                                         ccACAGAGCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCACGACGTGATGTTGCGTAt                    >  2:887543/1‑140 (MQ=255)                                                                                                              gagCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCACGACGTGATGTTGCGTATCGCCa               >  1:440439/1‑140 (MQ=255)                                                                                                              gagCGCCGCCGATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCACGACGTGATGTTGCGTATCGCCa               <  2:891392/140‑1 (MQ=255)                                                                                                                    cgccgATCTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTg         >  2:268576/1‑140 (MQ=255)                                                                                                                           cTGATTCGCGATGATCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCACAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGt  >  1:649960/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                              TATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGT  >  NC_000913/2015118‑2015379

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