Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,092,339	T→A	I10F (ATC→TTC)	frvB ←	putative PTS enzyme IIBC component FrvB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,092,339	0	T	A	100.0%	114.2 / NA	34	I10F (ATC→TTC)	frvB	putative PTS enzyme IIBC component FrvB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (16/18);  total (16/18)

TCCCGGCAAACCGCGCGCGATCATCACCGCTCAGGCCACGCCCGGTAGCGAGAATGACGTAATCGGCAGCGGCAATCTCTTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGATCGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGC  >  NC_000913/4092123‑4092577                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               tcCCGGCAAACCGCGCGCGATCATCACCGCTCAGGCCACGCCCGGTAGCGAGAATGACGTAATCGGCAGCGGCAATCTCTTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGcc                                                                                                                                                                                                                 <  1:165155/248‑1 (MQ=255)          aCCGCGCGCGATCATCACCGCTCAGGCCACGCCCGGTAGCGAGAATGACGTAATCGGCAGCGGCAATCTCTTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGa                                                                                                                                                                                                     <  1:226937/251‑1 (MQ=255)            cgcgcgcgATCATCACCGCTCAGGCCACGCCCGGTAGCGAGAATGACGTAATCGGCAGCGGCAATCTCTTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGc                                                                                                                                                                                                   >  1:234995/1‑251 (MQ=255)            cgcgcgcgATCATCACCGCTCAGGCCACGCCCGGTAGCGAGAATGACGTAATCGGCAGCGGCAATCTCTTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGc                                                                                                                                                                                                   <  1:370904/251‑1 (MQ=255)                         caCCGCTCAGGCCACGCCCGGTAGCGAGAATGACGTAATCGGCAGCGGCAATCTCTTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGAt                                                                                                                                                                                       <  1:109397/250‑1 (MQ=255)                                 aGGCCACGCCCGGTAGCGAGAATGACGTAATCGGCAGCGGCAATCTCTTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGGTCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCCTTTGAGCTGAGGGTTTGATTTAACAAgg                                                                                                                                                                              >  2:62245/1‑251 (MQ=255)                                  ggCCACGCCCGGTAGCGAGAATGACGTAATCGGCAGCGGCAATCTCTTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGGTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAAAAAGGc                                                                                                                                                                             >  2:115483/1‑251 (MQ=255)                                           cGGTAGCGAGAATGACGTAATCGGCAGCGGCAATCTCTTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTg                                                                                                                                                                    >  1:164937/1‑251 (MQ=255)                                               aGCGAGAATGACGTAATCGGCAGCGGCAATCTCTTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTc                                                                                                                                                                 <  2:143470/250‑1 (MQ=255)                                                          cGTAATCGGCAGCGGCAATCTCTTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACGTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTg                                                                                                                                                     <  2:164937/251‑1 (MQ=255)                                                          cGTAATCGGCAGCGGCAATCTCTTCGCTGGATAAGCGGTTTGCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCGTTCGGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTg                                                                                                                                                     <  2:360903/251‑1 (MQ=255)                                                                  gCAGCGGCAATCTCGTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTGATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGTCGGGGCAGGTGGTGAACGCGACAATACGTAAGTATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCTTGTTTCAGtt                                                                                                                                             <  2:343815/251‑1 (MQ=255)                                                                        acAATCTCTTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCGTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGaa                                                                                                                                       <  2:149806/250‑1 (MQ=255)                                                                               ctTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATaa                                                                                                                                    <  2:333807/247‑1 (MQ=255)                                                                                      ggATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACTGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTtgt                                                                                                                           <  2:308951/249‑1 (MQ=255)                                                                                      ggATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTtgt                                                                                                                           <  1:246603/249‑1 (MQ=255)                                                                                            gCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAAtt                                                                                                                   <  1:374747/251‑1 (MQ=255)                                                                                                   tCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAATTTTTCGac                                                                                                            >  2:84777/1‑251 (MQ=255)                                                                                                          cACTGGACACTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGCCAATACGTAAGGATGACTCCCTCGCCCCTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAACCCCTGGTCGGGATCGCCCTGTTGCGGTTGATGAATAAATTCCTTGTGATTAATTTTGCGaccaactgt                                                                                                     >  2:217911/1‑244 (MQ=255)                                                                                                           aCTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGGTGGGGAAAGCGAAAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGTGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAATTTTGCGACACAGtgtg                                                                                                    >  1:218481/1‑251 (MQ=255)                                                                                                                    ttGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGGTTGATTTAACAAGGCAAGCACCTGGTCGGGATCGCCCTGTTGCCGTTGATGAATAAATTCCTTGTGAATAATTTTGCGACACAGTGTGCCGATGAtt                                                                                           >  1:184010/1‑251 (MQ=255)                                                                                                                      gAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCTATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAATTTTGCGACACCGTGTGCCGATGAtttt                                                                                         >  2:254587/1‑251 (MQ=255)                                                                                                                        gTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTTAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGGTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAATTTTGCGACACAGTGTGCCGATGAttttt                                                                                        >  2:266789/1‑249 (MQ=255)                                                                                                                           tCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGGAGGTGTGAGCGAGCCCGGCGGGGCAGTTGGTGAACTCGACAATACGTAAGGATGACTCCATCGCCACTTCCTGATGAGCTGAGTGTTTGATTGAACAAGGCAAGCACCTGGGCTGTAGCGCCCTGTTGCAGTTGATGAATTAATGCCTTGTGAATAATTTTGCGACACAGTGTGCCGATGATTTTgccc                                                                                     <  1:28024/250‑4 (MQ=255)                                                                                                                                  ttATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCct                                                                              <  1:115917/250‑1 (MQ=255)                                                                                                                                                       gggCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGGGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCCCCTGGTCGGTATCGCCCTGTTGCCGTTGATGAATAAATTCCTTGTGAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCCCGcc                                                        >  2:340073/1‑251 (MQ=255)                                                                                                                                                                           tCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAATTGTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAAtt                                      <  2:179929/249‑1 (MQ=255)                                                                                                                                                                                 accaTGTCGGTGTGAGCGATCCAGGCGGGGGAGTTGGGGAACGCGCCAATACGTAAGGAGGACTCCATCGCCACTTCCTTATGCGGTGAGGGTTTGATTTAACAAGGCAAGCACCTGGGCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATACTTTTGCGACACAGTGTGCCGATGATTTTGACCTGATTCTTTTCGCCGCGTTGCGGCACGCCGAGGCAGATCCAGCAATTGACggatt                              >  2:87021/1‑247 (MQ=255)                                                                                                                                                                                   caTGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGAACGCGACAATACGTAAGGGTGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCg                            >  2:188870/1‑251 (MQ=255)                                                                                                                                                                                                gCGATCCCGGCGGGGAAGTTGGTGAACGCGACAATACGTAAGGATAACTCCATCGCCACTTCCTTATGAGCTGAGGGGTTGATTTAACAAGGAAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGGGAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCtt               >  1:131343/1‑251 (MQ=255)                                                                                                                                                                                                        ggcggGGCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTc        >  1:147504/1‑250 (MQ=255)                                                                                                                                                                                                           ggggCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCCTCGCCAATTCCTTATGAGCTGAGGGATTGATTTAACAATGCAAGCACATTGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCactc    >  1:38943/1‑248 (MQ=255)                                                                                                                                                                                                             ggCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGc  <  1:62245/251‑1 (MQ=255)                                                                                                                                                                                                             ggCAGTTGGTGAACGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGc  <  1:340099/251‑1 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               TCCCGGCAAACCGCGCGCGATCATCACCGCTCAGGCCACGCCCGGTAGCGAGAATGACGTAATCGGCAGCGGCAATCTCTTCGCTGGATAAGCGGTTTTCAACGCCACTGGACCCTTGAGTTTCCACTTTTATGGTATGACCGAGAGAACGGGCTTTCTGTTCCAGGGCTTCCGCCACCATGTAGGTGTGAGCGATCCCGGCGGGGCAGTTGGTGATCGCGACAATACGTAAGGATGACTCCATCGCCACTTCCTTATGAGCTGAGGGTTTGATTTAACAAGGCAAGCACCTGGTCGGTATCGCCCTGTTGCAGTTGATGAATAAATTCCTTGTGAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGC  >  NC_000913/4092123‑4092577

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