Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,139,721	Δ1 bp	coding (2501/2502 nt)	frwA ←	putative PTS enzyme: Hpr, enzyme I and II components

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,139,718	0	T	.	100.0%	209.8 / NA	44	intergenic (‑10/+2)	fsaB/frwA	fructose‑6‑phosphate aldolase 2/putative PTS enzyme: Hpr, enzyme I and II components
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (24/20);  total (24/20)

GCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTT  >  NC_000913/4139570‑4139863                                                                                                                                                     |                                                                                                                                                  gCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCttt                                                                                                                                                 >  1:414970/1‑151 (MQ=255)  cAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCttta                                                                                                                                                <  1:667527/151‑2 (MQ=255)  cAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCttt                                                                                                                                                 >  1:464630/1‑150 (MQ=255)   aGACGCGGCAGCACTTCCCATATGGACTCCTTGTTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGCCGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCtttac                                                                                                                                               >  2:203439/1‑149 (MQ=255)      cgcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTAcagt                                                                                                                                           <  1:620680/151‑1 (MQ=255)           cagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTAcagttcc                                                                                                                                        >  1:478437/1‑149 (MQ=255)            agcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTAcagttccagt                                                                                                                                     >  2:432128/1‑151 (MQ=255)            agcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTAcagttccagt                                                                                                                                     >  1:28620/1‑151 (MQ=255)             gcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTAcagttccagtt                                                                                                                                    >  1:654410/1‑151 (MQ=255)               aCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTAcagttccagttc                                                                                                                                   >  2:114874/1‑150 (MQ=255)               aCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCa                                                                                                                                  <  1:284996/151‑1 (MQ=255)                  tCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCAt                                                                                                                                 <  1:573695/149‑1 (MQ=255)                    ccATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTg                                                                                                                             <  1:612561/151‑1 (MQ=255)                      atatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTgca                                                                                                                           <  2:508166/151‑1 (MQ=255)                      atatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTgc                                                                                                                            >  1:636439/1‑150 (MQ=255)                      atatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTATCAAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTgc                                                                                                                            >  1:550769/1‑150 (MQ=255)                         tGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTgcagc                                                                                                                         >  2:458598/1‑150 (MQ=255)                          ggACCCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGcagcag                                                                                                                       <  2:430393/151‑1 (MQ=255)                              tCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTg                                                                                                                   >  1:27780/1‑151 (MQ=255)                                cTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTgg                                                                                                                  >  2:399987/1‑150 (MQ=255)                                  tGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCTGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGTCTGGCg                                                                                                                <  2:51293/150‑1 (MQ=255)                                           cGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCgg                                                                                                      <  2:613405/151‑1 (MQ=255)                                             aTAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCa                                                                                                    <  2:602921/151‑1 (MQ=255)                                                         ttAGTTGTAACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGAGtt                                                                                         >  2:655534/1‑150 (MQ=255)                                                               gTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAg                                                                                   <  2:471189/150‑1 (MQ=255)                                                                    cccGGGAATGGGGAAAATGCGGGCCAGGCCGTTGGCTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGtt                                                                             <  2:379497/151‑1 (MQ=255)                                                                    cccGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGtt                                                                             <  1:162211/151‑1 (MQ=255)                                                                     ccGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGtt                                                                             >  1:89356/1‑150 (MQ=255)                                                                     ccGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTa                                                                            <  1:207445/151‑1 (MQ=255)                                                                        gCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACgg                                                                         <  1:512841/151‑1 (MQ=255)                                                                          aaTGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGa                                                                        >  2:399855/1‑150 (MQ=255)                                                                           aTGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAAt                                                                      >  1:322860/1‑151 (MQ=255)                                                                              gggAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAAttct                                                                   >  1:520168/1‑151 (MQ=255)                                                                                          gCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCatgatg                                                       <  2:425239/151‑1 (MQ=255)                                                                                            cAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATgcg                                                     >  2:401797/1‑151 (MQ=255)                                                                                             aGACGTTCGACTTCTGCGACGCTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATgc                                                      >  2:26899/1‑149 (MQ=255)                                                                                                         tCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGgcgc                                        >  2:618219/1‑151 (MQ=255)                                                                                                            gCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCga                                      <  2:654389/150‑1 (MQ=255)                                                                                                              gACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGaa                                   <  1:222297/151‑1 (MQ=255)                                                                                                                   tcgCGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATgc                              >  1:116937/3‑151 (MQ=255)                                                                                                                      cGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATgcgc                            >  2:445636/1‑150 (MQ=255)                                                                                                                      cGGTGTCCAGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATgcgc                            >  2:224309/1‑150 (MQ=255)                                                                                                                              aGATACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTg                    >  1:137965/1‑150 (MQ=255)                                                                                                                                 tACAGTTCCATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCg                >  1:42921/1‑151 (MQ=255)                                                                                                                                        ccATGATGTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGc           >  2:402056/1‑149 (MQ=255)                                                                                                                                            gatgTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCgt      <  2:647288/150‑1 (MQ=255)                                                                                                                                             atgTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCgtgt    <  2:101715/151‑1 (MQ=255)                                                                                                                                               gTTTCC‑TTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGtgttt  <  1:249680/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTT  >  NC_000913/4139570‑4139863

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