Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,013,820	G→A	A198T (GCC→ACC)	fliF →	flagellar basal‑body MS‑ring and collar protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,013,820	0	G	A	100.0%	151.2 / NA	45	A198T (GCC→ACC)	fliF	flagellar basal‑body MS‑ring and collar protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/25);  total (20/25)

GGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCGCCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGGCGCGATCTTAATGACGCTCAGTTGAAATATGCCAGCGATGTCGA  >  NC_000913/2013690‑2013944                                                                                                                                   |                                                                                                                             ggCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTg                                                                                                                     >  2:170799/1‑140 (MQ=255) ggCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTg                                                                                                                     >  1:510467/1‑140 (MQ=255)              ccGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCggg                                                                                                        <  2:398048/140‑1 (MQ=255)                 tCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGaaa                                                                                                     >  2:269174/1‑140 (MQ=255)                  cTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAAc                                                                                                    >  2:350282/1‑140 (MQ=255)                       ttCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCAc                                                                                               <  1:2658/140‑1 (MQ=255)                        tCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACg                                                                                              >  2:13520/1‑140 (MQ=255)                           tCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCtg                                                                                           <  1:593615/140‑1 (MQ=255)                               tGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCtggtgg                                                                                       >  1:923138/1‑140 (MQ=255)                               tGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCtggtgg                                                                                       >  1:770669/1‑140 (MQ=255)                               tGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCtgg                                                                                          >  2:688080/1‑137 (MQ=255)                               tGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCtgg                                                                                          <  1:688080/137‑1 (MQ=255)                                 aaCAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCtg                                                                                           <  2:669437/134‑1 (MQ=255)                                 aaCAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCtg                                                                                           >  1:669437/1‑134 (MQ=255)                                 aaCAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGAt                                                                                     <  1:247394/140‑1 (MQ=255)                                 aaCAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGAt                                                                                     <  1:13520/140‑1 (MQ=255)                                      aaaTCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAggg                                                                                <  1:350282/140‑1 (MQ=255)                                      aaaTCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAggg                                                                                <  1:170799/140‑1 (MQ=255)                                        aTCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCg                                                                              <  2:697725/140‑1 (MQ=255)                                           cccTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGAc                                                                           >  2:570268/1‑140 (MQ=255)                                                cTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGAt                                                                                     >  1:927730/1‑125 (MQ=255)                                                cTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGAt                                                                                     <  2:927730/125‑1 (MQ=255)                                                   cATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTa                                                                   >  1:300412/1‑140 (MQ=255)                                                       ggTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAAccc                                                               <  2:300412/140‑1 (MQ=255)                                                         tGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAg                                                             <  1:570268/140‑1 (MQ=255)                                                          gACGGTAAATCTGTTACCCGGCCGCGCACTCGATAAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGt                                                            <  1:687519/140‑1 (MQ=255)                                                             ggTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAAcc                                                                >  1:374969/1‑133 (MQ=255)                                                             ggTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAAcc                                                                <  2:374969/133‑1 (MQ=255)                                                             ggTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCa                                                         >  1:163540/1‑140 (MQ=255)                                                                       ttACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCgg                                               >  1:476262/1‑140 (MQ=255)                                                                           ccGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGgcgc                                           >  2:790681/1‑140 (MQ=255)                                                                              gCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGTAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGgc                                             >  2:529886/1‑135 (MQ=255)                                                                              gCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGgc                                             <  1:529886/135‑1 (MQ=255)                                                                               ccGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGgc                                             <  2:557808/134‑1 (MQ=255)                                                                               ccGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGgc                                             >  1:557808/1‑134 (MQ=255)                                                                               ccGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGGCGCGATc                                       >  1:270157/1‑140 (MQ=255)                                                                                  cgcACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGgcg                                            >  2:379281/1‑132 (MQ=255)                                                                                  cgcACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGgcg                                            <  1:379281/132‑1 (MQ=255)                                                                                        cGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGGCGCGATCTTAATGACg                              <  2:476262/140‑1 (MQ=255)                                                                                        cGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGGCGCGATCTTAATGACg                              <  2:262518/140‑1 (MQ=255)                                                                                        cGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGGCGCGATCTTAATGACg                              <  2:923138/140‑1 (MQ=255)                                                                                          aTGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGGCGCGATCTTAATGACGCt                            <  1:790681/140‑1 (MQ=255)                                                                                                             cATTGTGCATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGGCGCGATCTTAATGACGCTCAGTTGAAATATGCCAGCg         <  2:270157/140‑1 (MQ=255)                                                                                                                    cATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGGCGCGATCTTAATGACGCTCAGTTGAAATATGCCAGCGATGTCGa  <  2:485181/140‑1 (MQ=255)                                                                                                                    cATCTGGTTTCCAGCACCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGGCGCGATCTTAATGACGCTCAGTTGAAATATGCCAGCGATGTCGa  <  1:269174/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                             GGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAGCGCCGTTGCTGGTCTGCCGCCGGGAAACGTCACGCTGGTGGATCAGGGCGGACATCTGTTAACCCAGTCCAATACCAGCGGGCGCGATCTTAATGACGCTCAGTTGAAATATGCCAGCGATGTCGA  >  NC_000913/2013690‑2013944

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