Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,019,460	G→A	V359I (GTC→ATC)	fliK →	flagellar hook‑length control protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,019,460	0	G	A	100.0%	140.6 / NA	42	V359I (GTC→ATC)	fliK	flagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/27);  total (15/27)

AGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGA  >  NC_000913/2019316‑2019598                                                                                                                                                 |                                                                                                                                           aGTGGCATTCAGGTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACACCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCAtctc                                                                                                                                        <  1:91409/149‑1 (MQ=255)  gTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATctct                                                                                                                                       <  1:901002/149‑1 (MQ=255)   tGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCGTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCtt                                                                                                                                      <  2:702263/149‑1 (MQ=255)            gTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGc                                                                                                                             >  2:293235/1‑149 (MQ=255)               aGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCgtg                                                                                                                          <  2:601503/149‑1 (MQ=255)                 ggCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGTGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTa                                                                                                                        >  2:610576/1‑149 (MQ=255)                    aaaGCAACATCAGTGGCGAAAGCTTTAGTGGGCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACa                                                                                                                     <  1:904993/149‑1 (MQ=255)                     aaGCAACATAAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAg                                                                                                                    <  1:464348/149‑1 (MQ=255)                      aGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCACTTTACAAGGGCGTGTAACAgg                                                                                                                   >  1:448424/1‑149 (MQ=255)                         aaCATCAGTGGCGAAAGCTTTAGTGGTCAGAAGCAGGCCGCTTCCCAGCAACAGAAAAGCCAACGCAAAGCAAACCATGAACCTCTGGCGGGGGAAGCCGACGATACGCTTCCGCTTCCCAAATCTTTACAAGGGCGTGTAACAGGCaa                                                                                                                >  1:49847/1‑149 (MQ=255)                              cAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAgcg                                                                                                           >  1:417014/1‑149 (MQ=255)                                  ggCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGt                                                                                                       <  2:717538/149‑1 (MQ=255)                                   gCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                                                                                      <  2:417014/149‑1 (MQ=255)                                              aGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGc                                                                                           <  2:549135/149‑1 (MQ=255)                                                 ggTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTa                                                                                        <  2:339279/149‑1 (MQ=255)                                                  gTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTaa                                                                                       <  1:293235/149‑1 (MQ=255)                                                    gagcagcatGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACg                                                                                     <  2:844318/148‑1 (MQ=255)                                                    cagcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACg                                                                                     <  2:387536/149‑1 (MQ=255)                                                        agcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCAGGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCag                                                                                 <  2:55731/149‑1 (MQ=255)                                                         gcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCGCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTGGATATTTTCGCCTAACGTCcga                                                                                >  2:1033967/1‑149 (MQ=255)                                                             gCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGCCGACGCTACGCTTCCGGTTCCCATCTATTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGAGAGAGGTa                                                                            >  2:1069142/1‑149 (MQ=255)                                                             gCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTa                                                                            <  1:610576/149‑1 (MQ=255)                                                             gCCGCTTCCCAGCAACACCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCCGGGGAAGACGACGAAACGCATCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCTATGAAAAATTCGCCTAAAGTCAGAGGTa                                                                            >  2:111414/1‑149 (MQ=255)                                                              ccGCTTCCCAGCACCAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAg                                                                           <  2:1162243/149‑1 (MQ=255)                                                               cGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGc                                                                          <  2:988624/149‑1 (MQ=255)                                                                 cTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACTCTTCCGGTTCCCATCTCTTTACAATGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCAc                                                                        <  2:162425/149‑1 (MQ=255)                                                                        ccaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATc                                                                 <  2:997427/148‑1 (MQ=255)                                                                          aacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCcg                                                               >  2:651826/1‑149 (MQ=255)                                                                                   gCCAACGCACAGCAACCCAGTACCCTCTCGCGGGCGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTc                                                      <  2:65372/149‑1 (MQ=255)                                                                                     cAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTccc                                                    <  2:846329/149‑1 (MQ=255)                                                                                         gCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCt                                                >  2:826747/1‑149 (MQ=255)                                                                                             aGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTttgt                                            >  2:981675/1‑149 (MQ=255)                                                                                              tcAAACCATGAACGTCTGGCGGGGGAAGACGACGATTCGCTTCCGGTTCACATCTCTTTTCAAGGGCGTGTAAGAGGCAACAGCTGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCAttgtt                                           <  1:801708/148‑1 (MQ=255)                                                                                              gCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGGCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTttgtt                                           <  2:1092074/149‑1 (MQ=255)                                                                                                                gCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGAt                         >  2:363797/1‑149 (MQ=255)                                                                                                                   ggggAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAAt                      <  1:981675/149‑1 (MQ=255)                                                                                                                    gggAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAAtt                     <  1:651826/149‑1 (MQ=255)                                                                                                                       acgacgagGATAGGGATCCGGTTCCCATCTCTGTCCAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTATCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGc                  <  2:194707/147‑1 (MQ=255)                                                                                                                               gaTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAga          >  1:794418/1‑148 (MQ=255)                                                                                                                               gaTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGc          >  2:189428/1‑149 (MQ=255)                                                                                                                               gaTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGc          >  2:602853/1‑149 (MQ=255)                                                                                                                                       tCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGa  <  1:764661/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                           AGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGA  >  NC_000913/2019316‑2019598

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑