Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,133,140	T→G	S189R (AGT→AGG)	flgE →	flagellar hook protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,133,140	0	T	G	100.0%	200.3 / NA	56	S189R (AGT→AGG)	flgE	flagellar hook protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (33/23);  total (33/23)

AATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGT  >  NC_000913/1132994‑1133283                                                                                                                                                   |                                                                                                                                                aaTACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGAcaggcagg                                                                                                                                             >  1:501947/1‑151 (MQ=255)  aTACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAggg                                                                                                                                            >  1:389034/1‑151 (MQ=255)    aCCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGATCGGTGACTGTTTTCGGCAGGCAggg                                                                                                                                            >  2:592862/1‑149 (MQ=255)       cTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATg                                                                                                                                       >  2:802115/1‑151 (MQ=255)        tGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGc                                                                                                                                      >  1:706905/1‑151 (MQ=255)          aTGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTc                                                                                                                                    >  2:586801/1‑151 (MQ=255)            ggCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTcat                                                                                                                                  <  1:571172/151‑1 (MQ=255)             gCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTcatg                                                                                                                                 <  1:344226/151‑1 (MQ=255)             gCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTcatg                                                                                                                                 <  1:403203/151‑1 (MQ=255)                   aaaaCTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTcatgacatg                                                                                                                            >  2:444590/1‑150 (MQ=255)                         accaccACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCt                                                                                                                     >  2:782386/1‑151 (MQ=255)                         accaccACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCt                                                                                                                     <  1:452557/151‑1 (MQ=255)                          ccaccaCGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTa                                                                                                                    <  1:367114/151‑1 (MQ=255)                           caccacGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTAc                                                                                                                   >  1:681368/1‑151 (MQ=255)                           caccacGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTAc                                                                                                                   >  2:287111/1‑151 (MQ=255)                               acGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTtg                                                                                                               <  1:901494/151‑1 (MQ=255)                               acGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTtg                                                                                                               >  2:413095/1‑151 (MQ=255)                                      cGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGTAGa                                                                                                         >  2:590439/1‑150 (MQ=255)                                       gATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGAcc                                                                                                       >  1:464696/1‑151 (MQ=255)                                         tGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCgg                                                                                                     >  1:388130/1‑151 (MQ=255)                                            aGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCgggg                                                                                                   >  2:934577/1‑150 (MQ=255)                                             gATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGat                                                                                                 >  2:99831/1‑151 (MQ=255)                                                       aaTTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGc                                                                                       >  1:405155/1‑151 (MQ=255)                                                          tCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAg                                                                                     <  1:145651/150‑1 (MQ=255)                                                           ccAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGt                                                                                   >  1:527314/1‑151 (MQ=255)                                                             aGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCt                                                                                 <  1:270194/151‑1 (MQ=255)                                                               tGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTa                                                                                >  1:755782/1‑150 (MQ=255)                                                                gATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTaca                                                                              <  2:162737/151‑1 (MQ=255)                                                                  tCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACAcc                                                                            >  1:732823/1‑151 (MQ=255)                                                                    cGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACAccc                                                                           >  2:161745/1‑150 (MQ=255)                                                                     gCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACAccc                                                                           >  1:123751/1‑149 (MQ=255)                                                                       ttCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAgg                                                                        >  1:710250/1‑150 (MQ=255)                                                                           tACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAg                                                                    >  1:55807/1‑150 (MQ=255)                                                                            aCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAg                                                                    <  2:706925/149‑1 (MQ=255)                                                                                ttACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTg                                                               >  1:545468/1‑150 (MQ=255)                                                                                    gCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCa                                                          <  2:716312/151‑1 (MQ=255)                                                                                     ccATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCaa                                                         <  1:169128/151‑1 (MQ=255)                                                                                     ccATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCa                                                          >  1:806236/1‑150 (MQ=255)                                                                                      cATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCaa                                                         >  2:96310/1‑150 (MQ=255)                                                                                            gcgcCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAt                                                  <  2:816028/151‑1 (MQ=255)                                                                                            gcgcCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGAGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAt                                                  <  1:486575/151‑1 (MQ=255)                                                                                             cgcCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACTCCCAGGATAGCAGTGATCCAAACAGCAt                                                  <  1:70730/150‑1 (MQ=255)                                                                                             cgcCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAtt                                                 <  1:533187/151‑1 (MQ=255)                                                                                                 aGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGa                                             <  2:457669/151‑1 (MQ=255)                                                                                                  gCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGaa                                            <  1:381235/151‑1 (MQ=255)                                                                                                      tGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGAc                                         >  1:515157/1‑150 (MQ=255)                                                                                                          gATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGa                                    <  2:444472/151‑1 (MQ=255)                                                                                                              gCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGacaac                                >  2:879769/1‑151 (MQ=255)                                                                                                                tataACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAacac                              >  1:614425/1‑151 (MQ=255)                                                                                                                 ataACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGacaaca                               <  1:628999/149‑1 (MQ=255)                                                                                                                  taACAAAAAAGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTg                            >  1:484129/1‑151 (MQ=255)                                                                                                                        aaaaGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAAttt                      >  1:730963/1‑151 (MQ=255)                                                                                                                         aaaGGTTCGGTGACTGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTa                     <  2:389048/151‑1 (MQ=255)                                                                                                                                       tGTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGcac        >  2:255028/1‑150 (MQ=255)                                                                                                                                        gTTTTCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACAt      <  2:484733/151‑1 (MQ=255)                                                                                                                                          tttCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACAtt     <  1:600842/150‑1 (MQ=255)                                                                                                                                            tCGACAGGCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGt  >  2:499330/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                AATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGT  >  NC_000913/1132994‑1133283

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