Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,438,548	T→C	V200A (GTG→GCG)	flk →	putative flagella assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,438,548	0	T	C	100.0%	208.6 / NA	60	V200A (GTG→GCG)	flk	putative flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/30);  total (30/30)

GCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGTGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTAGATGC  >  NC_000913/2438419‑2438683                                                                                                                                  |                                                                                                                                        gcgAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                                >  2:869967/1‑139 (MQ=255)  cgAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                               <  2:715283/139‑1 (MQ=255)  cgAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                               >  2:634801/1‑139 (MQ=255)  cgAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                               >  1:304235/1‑139 (MQ=255)     ggCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAg                                                                                                                            >  1:40129/1‑139 (MQ=255)                 tCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCg                                                                                                                >  2:1471261/1‑139 (MQ=255)                  ccTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                               <  1:634801/139‑1 (MQ=255)                    tttATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                                             <  1:801579/139‑1 (MQ=255)                    tttATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                                             <  1:400526/139‑1 (MQ=255)                              gcggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                                             >  2:28096/1‑129 (MQ=255)                              gcggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                                             <  1:28096/129‑1 (MQ=255)                                ggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCAt                                                                                                 >  1:121714/1‑139 (MQ=255)                                        aTGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                                         <  2:889503/139‑1 (MQ=255)                                        aTGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                                         <  2:1166670/139‑1 (MQ=255)                                            gcTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAAt                                                                                                                       >  1:1028129/1‑105 (MQ=255)                                            gcTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAAt                                                                                                                       <  2:1028129/105‑1 (MQ=255)                                              tAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCt                                                                                   <  1:1471261/139‑1 (MQ=255)                                                 aaCAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACa                                                                                >  1:1298000/1‑139 (MQ=255)                                                     gCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGgcg                                                                            >  1:255240/1‑139 (MQ=255)                                                        ggTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGt                                                                         >  1:1321945/1‑139 (MQ=255)                                                          tGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCa                                                                       <  2:1068236/139‑1 (MQ=255)                                                            aCCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                               >  2:395847/1‑81 (MQ=255)                                                            aCCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                               <  1:395847/81‑1 (MQ=255)                                                                      cggcggCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCa                                                           >  2:1360928/1‑139 (MQ=255)                                                                           gcAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                               >  1:799253/1‑82 (MQ=255)                                                                           gcAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                               <  2:799253/82‑1 (MQ=255)                                                                           gcAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCaaa                                                      >  1:1195651/1‑139 (MQ=255)                                                                                  gggAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGcc                                               <  1:869967/139‑1 (MQ=255)                                                                                   ggAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCg                                              <  2:304235/139‑1 (MQ=255)                                                                                       cccAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGc                                          <  2:121714/139‑1 (MQ=255)                                                                                        ccAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCt                                         <  2:40129/139‑1 (MQ=255)                                                                                          aGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAAcc                                                                                                   >  2:758902/1‑79 (MQ=255)                                                                                          aGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAAcc                                                                                                   <  1:758902/79‑1 (MQ=255)                                                                                             aaaCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAAt                                    <  2:552855/139‑1 (MQ=255)                                                                                             aaaCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAAt                                    <  2:1195651/139‑1 (MQ=255)                                                                                               aCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                                <  2:119553/45‑1 (MQ=255)                                                                                               aCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                                >  1:119553/1‑45 (MQ=255)                                                                                                   aTCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAAt                                                    >  2:1012653/1‑117 (MQ=255)                                                                                                   aTCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAAt                                                    <  1:1012653/117‑1 (MQ=255)                                                                                                   aTCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGAc                                            <  2:1428920/125‑1 (MQ=255)                                                                                                   aTCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGAc                                            >  1:1428920/1‑125 (MQ=255)                                                                                                         cAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGa                        <  2:255240/139‑1 (MQ=255)                                                                                                           aTCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                                             >  1:379651/1‑52 (MQ=255)                                                                                                           aTCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                                             <  2:379651/52‑1 (MQ=255)                                                                                                             cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTg                                                                                              <  1:927544/65‑1 (MQ=255)                                                                                                             cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTg                                                                                              >  2:927544/1‑65 (MQ=255)                                                                                                             cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCt                    >  2:264550/1‑139 (MQ=255)                                                                                                                   tGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAaaca              <  2:1321945/139‑1 (MQ=255)                                                                                                                       aCTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCt          >  1:370558/1‑139 (MQ=255)                                                                                                                        cTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                                         >  2:668565/1‑59 (MQ=255)                                                                                                                        cTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                                         <  1:668565/59‑1 (MQ=255)                                                                                                                          ttCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAgg                                                                              >  2:523385/1‑68 (MQ=255)                                                                                                                          ttCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAgg                                                                              <  1:523385/68‑1 (MQ=255)                                                                                                                          ttCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGaa                                     <  1:1345677/109‑1 (MQ=255)                                                                                                                          ttCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGaa                                     >  2:1345677/1‑109 (MQ=255)                                                                                                                          ttCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGa                                      >  1:1285127/1‑108 (MQ=255)                                                                                                                          ttCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGa                                      <  2:1285127/108‑1 (MQ=255)                                                                                                                          ttCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTa       >  2:292695/1‑139 (MQ=255)                                                                                                                             cGGGGCGAAAGATGGCGAGTTTATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTAGAt    >  1:638928/1‑139 (MQ=255)                                                                                                                               gggCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTAGATGc  >  2:29992/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                        GCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGTGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTAGATGC  >  NC_000913/2438419‑2438683

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