Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,133,535	A→G	N321S (AAC→AGC)	flgE →	flagellar hook protein FlgE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,133,535	0	A	G	100.0%	62.0 / NA	20	N321S (AAC→AGC)	flgE	flagellar hook protein FlgE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/10);  total (10/10)

AGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAG  >  NC_000913/1133412‑1133657                                                                                                                            |                                                                                                                           aGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTCCTGGGGCAGATTGTACGGGCGAACTGTGCCGGCAACGAAgg                                                                                                                   >  2:184913/1‑133 (MQ=255)                  ggCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCg                                                                                                  >  2:70142/1‑133 (MQ=255)                    cGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGAc                                                                                                <  2:259197/133‑1 (MQ=255)                                ttATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACgtctggtct                                                                                    >  2:314082/1‑133 (MQ=255)                                  atcaaatcaaTGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGc                                                                                  >  1:323214/1‑133 (MQ=255)                                      aatcaaTGATGACGGTACGGTTGTCGGCTACTGTTCCAACGTACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACg                                                                              <  2:244392/133‑1 (MQ=255)                                       atcaaTGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTTGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGc                                                                             <  1:12169/133‑1 (MQ=255)                                                   ggtacggtTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCtggcg                                                                 <  2:430234/133‑1 (MQ=255)                                                        ggtTGTCGTCAACTACTCGAACGAACAAACGCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCtggcgtggcg                                                            <  1:241477/133‑1 (MQ=255)                                                        ggtTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCtggcgtggcg                                                            <  2:299114/133‑1 (MQ=255)                                                          ttGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCt                                                          >  1:462617/1‑133 (MQ=255)                                                            gTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGt                                                        >  1:512348/1‑133 (MQ=255)                                                                gCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTggg                                                    <  1:404154/133‑1 (MQ=255)                                                                                aCAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCcgggacggga                                    <  1:70142/133‑1 (MQ=255)                                                                                        aaCTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTgg                            >  2:465710/1‑133 (MQ=255)                                                                                             ctgGGGCAGATTGTACTGGCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCAccc                       >  2:559/1‑133 (MQ=255)                                                                                                               gCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGTGTGGTGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACgg             <  1:410141/125‑1 (MQ=255)                                                                                                               gCGAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACgg             >  2:410141/1‑125 (MQ=255)                                                                                                                cGAACTTTGCGAGCAACTAAGGTCTGGCATCCGAAGCCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGCCCCTGTTGGGGACCGCCGGGACGGGCACCTTTGGCACCCTGACCAACGGTGCGCTGGa    >  2:270072/1‑133 (MQ=255)                                                                                                                 gAACTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGaa   <  1:314082/133‑1 (MQ=255)                                                                                                                  aaCTTTGCCAGCAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAg  >  1:69110/1‑133 (MQ=255)                                                                                                                            |                                                                                                                           AGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAG  >  NC_000913/1133412‑1133657

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