Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,437,859	T→C	intergenic (‑8/‑91)	pdxB ← / → flk	erythronate‑4‑phosphate dehydrogenase/putative flagella assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,437,859	0	T	C	100.0%	108.8 / NA	32	intergenic (‑8/‑91)	pdxB/flk	erythronate‑4‑phosphate dehydrogenase/putative flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/23);  total (9/23)

AGCGCATCCGCGTCTGCCAGTTGAGCGACGGGGATTGGACGCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGTTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCGGCCAGGCATACTTGCGAAGATTTCAGGTATAAGGATACGTAATGATACAACCTATTTCCGGCC  >  NC_000913/2437730‑2437971                                                                                                                                  |                                                                                                                 aGCGCATCCGCGTCTGCCAGTTGAGCGACGGGGATTGGACGCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTAc                                                                                                               <  1:37554/133‑1 (MQ=255)     cATCCGCGTCTGCCAGTTGAGCGACGGGGATTGGACGCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGTATATTTTCATCAACAAAGATTTTCAAGTTTCTGCTACCTGt                                                                                                           >  2:195035/1‑133 (MQ=255)      aTCCGCGTCTGCCAGTTGAGCGACGGGGATTGGACGCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTa                                                                                                          <  1:34114/133‑1 (MQ=255)      aTCCGCGTCTGCCAGTTGAGCGACGGGGATTGGACGCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTa                                                                                                          <  2:423496/133‑1 (MQ=255)       tCCGCGTCTGCCAGTTGAGCGACGGGGATTGGACGCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTAt                                                                                                         >  1:384585/1‑133 (MQ=255)       tCCGCGTCTGCCAGTTGAGCGACGGGGATTGGACGCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCAAATTTTCATTAACGAGGTTTTTCACGATTGTCCTACCTGTAt                                                                                                         >  2:169438/1‑133 (MQ=255)           cgTCTGCCCGTTGAGCGACGGGCATTGGACTCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATgaga                                                                                                     <  1:149057/133‑1 (MQ=255)                gCCAGTTGAGCGACGGGGATTGGACGCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATgaga                                                                                                     >  1:128610/1‑128 (MQ=255)                gCCAGTTGAGCGACGGGGATTGGACGCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATgaga                                                                                                     <  2:128610/128‑1 (MQ=255)                 ccAGTTGAGCGACGGAGATTGGGCGCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGagag                                                                                               <  1:178792/133‑1 (MQ=255)                     ttGAGCGACGGGGATTGGACGCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTaa                                                                                           <  2:90808/133‑1 (MQ=255)                              ggggATTGGACGCCCGGGAACCACGGTCACGTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAg                                                                                  <  1:195035/133‑1 (MQ=255)                              ggggATTGGAAGCCCTGGAACCGCCGTCTCCTGACCCAAACCGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAg                                                                                  <  2:184195/133‑1 (MQ=255)                                ggATTGGACGTCCGGGAACCGCGGTCACTTCACCCAAACGGGTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAgtg                                                                                <  2:187678/133‑1 (MQ=255)                                ggATTGGACGCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAgtg                                                                                <  1:332263/133‑1 (MQ=255)                                           ccGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAgtg                                                                                >  2:346198/1‑122 (MQ=255)                                           ccGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAgtg                                                                                <  1:346198/122‑1 (MQ=255)                                               gAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTcgcg                                                                 <  2:32830/133‑1 (MQ=255)                                               gAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTcgcg                                                                 <  1:97507/133‑1 (MQ=255)                                                aaCCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCgg                                                                <  2:284611/133‑1 (MQ=255)                                                aaCCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCgg                                                                <  2:294205/133‑1 (MQ=255)                                                 aCCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCGGc                                                               <  1:90913/133‑1 (MQ=255)                                                                     aCGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCGGCCAGGCATACTTGCGAAGAtt                                           <  1:382336/133‑1 (MQ=255)                                                                       ggCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCGGCCAGGCATACTTGCGAAGATTTc                                         <  1:313161/133‑1 (MQ=255)                                                                         cTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCGGCCAGGCATACTTGCGAAGATTTCAg                                       >  1:316702/1‑133 (MQ=255)                                                                         cTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCCCGCGGCCAGGCGTATTTGCGAGGATTTCAg                                       >  1:248836/1‑133 (MQ=255)                                                                           aaaTAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCGGCCAGGCATACTTGCGAAGATTTCAGGt                                     <  1:511044/133‑1 (MQ=255)                                                                              tAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCGGCCAGGCATACTTGCGAAGATTTCAGGtata                                  <  2:384585/133‑1 (MQ=255)                                                                              tAAGTCGCGGCCATATGGCATATTCTCACCAAAAAGGATTTCCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCGGCCAGGCATACTTGCGAAGATTTCAGGtata                                  <  2:130489/133‑1 (MQ=255)                                                                                 gTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCGGCCTGGCATACATGCTAATATTTCATGTATAAgg                               >  1:134270/1‑133 (MQ=255)                                                                                                         atcaACAAGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCGGCCAGGCATACTTGCGAAGATTTCAGGTATAAGGATACGTAATGATACAACCTATTTc       <  2:316702/133‑1 (MQ=255)                                                                                                              caaGGATTTTCACGTTTGTGCTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCGGCCAGGCATACTTGCGAAGATGTCAGGTATTAGGGTACGTAATGATACAACCCATTTCCGGcc  >  2:206692/1‑133 (MQ=255)                                                                                                                                  |                                                                                                                 AGCGCATCCGCGTCTGCCAGTTGAGCGACGGGGATTGGACGCCCGGGAACCGCGGTCACCTCACCCAAACGGCTAAATAAGTCGCGGGCATAAGGCATATTTTCATCAACAAGGATTTTCACGTTTGTGTTACCTGTATGAGACGAGAGTTAACCGGACAAGTGTGCCATAATCTCGCGGCCAGGCATACTTGCGAAGATTTCAGGTATAAGGATACGTAATGATACAACCTATTTCCGGCC  >  NC_000913/2437730‑2437971

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