Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,132,321	C→T	T156T (ACC→ACT)	flgD →	flagellar hook assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,132,321	0	C	T	100.0%	176.4 / NA	44	T156T (ACC→ACT)	flgD	flagellar hook assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/23);  total (21/23)

TGCGGGAACCGGCAGTGAAGAAGGGGCTGTGACCACGACCACGCCGTTTGGTGTTGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACCGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGG  >  NC_000913/1132177‑1132464                                                                                                                                                 |                                                                                                                                                tGCGGGAACCGGCAGTGAAGAAGGGGCTGTGACCACGACCACGCCGTTTGGTGTTGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGActgccg                                                                                                                                             <  1:20529/149‑1 (MQ=255)  gCGGGAACCGGCAGTGAAGAAGGGGCTGTGACCACGACCACGCCGTTTGGTGTTGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCgg                                                                                                                                            >  2:272405/1‑149 (MQ=255)     ggAACCGGCAGTGAAGAAGGGGCTGTGACCACGACCACGCCGTTTGGTGTTGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTc                                                                                                                                       >  1:46151/1‑151 (MQ=255)        aCCGGCAGTGAAGAAGGGGCTGTGACCACGACCACGCCGTTTGGTGTTGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTcaca                                                                                                                                    <  2:207337/151‑1 (MQ=255)           ggCAGTGAAGAAGGGGCTGTGACCACGACCACGCCGTTTGGTGTTGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGtt                                                                                                                                 <  1:577760/151‑1 (MQ=255)                      agGGGCTGTGACCACGACCACGCCGTTTGGTGTTGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGAc                                                                                                                      <  1:316604/151‑1 (MQ=255)                                accacgaccacgCCGTTTGGTGTTGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTg                                                                                                             >  1:590831/1‑150 (MQ=255)                                 ccacgaccacgCCGTTTGGTGTTGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGAc                                                                                                           >  1:440888/1‑151 (MQ=255)                                    cgaccacgCCGTTTGGTGTTGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTg                                                                                                         >  1:661777/1‑150 (MQ=255)                                           gCCGTTTGGTGTTGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGgcagc                                                                                                 <  1:571774/151‑1 (MQ=255)                                            ccGTTTGGTGTTGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGgcagca                                                                                                <  2:506185/151‑1 (MQ=255)                                                  ggTGTTGGGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTgcgc                                                                                          <  1:261006/151‑1 (MQ=255)                                                       tGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAAc                                                                                     >  2:389838/1‑151 (MQ=255)                                                        gAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACg                                                                                    <  2:353280/151‑1 (MQ=255)                                                        gAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACg                                                                                    <  2:398012/151‑1 (MQ=255)                                                        gAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCAATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAAc                                                                                     <  2:350261/150‑1 (MQ=255)                                                              cAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCtt                                                                              <  1:674539/151‑1 (MQ=255)                                                              cAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCtt                                                                              <  2:570745/151‑1 (MQ=255)                                                              cAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACATTGACTGATGGCAGCACTGCGCCGAACGGTTCtt                                                                              <  1:4275/151‑1 (MQ=255)                                                                  aGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACa                                                                           >  2:319690/1‑150 (MQ=255)                                                                   ggcggcAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAAt                                                                         <  1:164043/151‑1 (MQ=255)                                                                        cAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGc                                                                    <  1:590705/151‑1 (MQ=255)                                                                         aGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCg                                                                   <  2:265064/151‑1 (MQ=255)                                                                            cAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGAtt                                                                >  1:573633/1‑151 (MQ=255)                                                                              aaGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAg                                                              <  1:319314/151‑1 (MQ=255)                                                                               aGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAgc                                                             <  2:389012/151‑1 (MQ=255)                                                                               aGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAgc                                                             <  1:389849/151‑1 (MQ=255)                                                                                      gCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGt                                                       >  1:488004/1‑150 (MQ=255)                                                                                         accaTCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACg                                                   <  2:442183/151‑1 (MQ=255)                                                                                          ccaTCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACgg                                                  >  1:677878/1‑151 (MQ=255)                                                                                            aTCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACggtg                                                >  2:263331/1‑151 (MQ=255)                                                                                               aCCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACggtggt                                              >  1:28605/1‑150 (MQ=255)                                                                                                  gATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTacac                                          >  1:670142/1‑151 (MQ=255)                                                                                                    tAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAGCGGTGGTAcaca                                         >  2:153580/1‑150 (MQ=255)                                                                                                       aaaTGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTg                                     >  1:595847/1‑151 (MQ=255)                                                                                                         aTGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGt                                   <  1:171848/151‑1 (MQ=255)                                                                                                         aTGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGt                                   <  2:552316/151‑1 (MQ=255)                                                                                                                    gtgCGCACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGc                        >  2:458926/1‑151 (MQ=255)                                                                                                                         cACCATTGATATTGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAg                   <  1:594576/151‑1 (MQ=255)                                                                                                                                     tGGTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTg       >  2:176936/1‑151 (MQ=255)                                                                                                                                      ggTGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGc      <  2:421724/151‑1 (MQ=255)                                                                                                                                        tGAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGc      >  2:612635/1‑149 (MQ=255)                                                                                                                                         gAACTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAgg   >  2:335177/1‑151 (MQ=255)                                                                                                                                         aaaCTGACTGCCGGAGATCACAGTTTCACCAAGGACGGTACGTTGACTGATAGCAGCACTGCGCAGAACGGTTCTTACAATGTAGCGATTAGCGCCAGAAACGGAGGTACACAACTGGTTCTCCAGCCGCAGCATCGTGCTCTGGTGCAgg   >  2:545750/2‑151 (MQ=255)                                                                                                                                          aaCTGACTGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAggg  >  2:627343/1‑151 (MQ=255)                                                                                                                                           actgactgCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGt                                                       >  1:13015/1‑97 (MQ=255)                                                                                                                                           actgactgCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGt                                                       <  2:13015/97‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                TGCGGGAACCGGCAGTGAAGAAGGGGCTGTGACCACGACCACGCCGTTTGGTGTTGAGCTGCAACAGGCGGCAGACAAAGTTACGGCCACCATCACCGATAAAAATGGCGCGGTTGTGCGCACCATTGATATTGGTGAACTGACCGCCGGAGTTCACAGTTTCACCTGGGACGGTACGTTGACTGATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGG  >  NC_000913/1132177‑1132464

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