Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,136,783	T→C	V170V (GTT→GTC)	flgI →	putative flagellar basal body protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,136,783	0	T	C	100.0%	153.5 / NA	41	V170V (GTT→GTC)	flgI	putative flagellar basal body protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/24);  total (17/24)

CGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCGCGTGCGTGGATATGGCAGCGCCAC  >  NC_000913/1136647‑1136914                                                                                                                                         |                                                                                                                                    cGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTATCCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCAt                                                                                                                                   <  1:537381/139‑1 (MQ=255) cGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCAt                                                                                                                                   >  2:548420/1‑139 (MQ=255)  gCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATACTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCAtt                                                                                                                                  <  1:899587/139‑1 (MQ=255)               gACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGc                                                                                                                     <  1:1200511/139‑1 (MQ=255)                aCAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGcc                                                                                                                    <  2:1147891/139‑1 (MQ=255)                  aGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCcca                                                                                                                  <  2:732344/139‑1 (MQ=255)                    tCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCccagc                                                                                                                <  1:613522/139‑1 (MQ=255)                               cgcTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGt                                                                                                     <  1:106139/139‑1 (MQ=255)                                  tGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGG‑CGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCggg                                                                                                 >  2:1498075/1‑139 (MQ=255)                                   ggCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCggg                                                                                                 >  1:736949/1‑139 (MQ=255)                                    gcgcAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGa                                                                                                <  1:1500547/139‑1 (MQ=255)                                     cgcAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGaa                                                                                               <  1:548420/139‑1 (MQ=255)                                     cgcAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGG‑CGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAAt                                                                                              <  1:1498075/139‑1 (MQ=255)                                           gCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCt                                                                                         <  2:468094/139‑1 (MQ=255)                                                atTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAAtt                                                                                    <  1:614631/139‑1 (MQ=255)                                                         ggcggcGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTa                                                                           >  2:496470/1‑139 (MQ=255)                                                           cggcgCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCcca                                                                                                                  <  2:680295/98‑1 (MQ=255)                                                           cggcgCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCcca                                                                                                                  >  1:680295/1‑98 (MQ=255)                                                           cggcgCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAac                                                                         >  1:104088/1‑139 (MQ=255)                                                               gcAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAacgacg                                                                     <  1:651756/139‑1 (MQ=255)                                                                                             cAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTg                                       <  1:496470/139‑1 (MQ=255)                                                                                                   aaCCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCa                                 >  2:8098/1‑139 (MQ=255)                                                                                                           gAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAAtt                                                                                    >  2:330962/1‑80 (MQ=255)                                                                                                           gAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAAtt                                                                                    <  1:330962/80‑1 (MQ=255)                                                                                                              cGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCgcgtg                      >  2:554008/1‑139 (MQ=255)                                                                                                               ggTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCgcgtgc                     <  2:1135377/139‑1 (MQ=255)                                                                                                                        aTCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTg                                                                                  >  1:410803/1‑69 (MQ=255)                                                                                                                        aTCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTg                                                                                  <  2:410803/69‑1 (MQ=255)                                                                                                                           aCCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGa                                                                                                >  2:1076831/1‑52 (MQ=255)                                                                                                                           aCCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGa                                                                                                <  1:1076831/52‑1 (MQ=255)                                                                                                                             cAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGgcg                                                    <  2:86100/94‑1 (MQ=255)                                                                                                                             cAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGgcg                                                    >  1:86100/1‑94 (MQ=255)                                                                                                                                  gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTg                                       <  1:321358/102‑1 (MQ=255)                                                                                                                                  gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTg                                       >  2:321358/1‑102 (MQ=255)                                                                                                                                  gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCgcgtgc                     <  1:1483411/120‑1 (MQ=255)                                                                                                                                  gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCgcgtgc                     >  2:699091/1‑120 (MQ=255)                                                                                                                                  gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCgcgtgc                     >  2:1483411/1‑120 (MQ=255)                                                                                                                                  gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCgcgtgc                     <  1:699091/120‑1 (MQ=255)                                                                                                                                  gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCGCGTGCGTGGATATGGCAGCGCCAc  >  1:13685/1‑139 (MQ=255)                                                                                                                                  gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAA‑TTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCgcgtgc                     >  2:659026/1‑119 (MQ=255)                                                                                                                                  gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAA‑TTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCgcgtgc                     <  1:659026/119‑1 (MQ=255)                                                                                                                                         |                                                                                                                                    CGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCGCGTGCGTGGATATGGCAGCGCCAC  >  NC_000913/1136647‑1136914

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