Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	282,361	C→T	G28G (GGC→GGT)	yagE →	2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	282,361	0	C	T	100.0%	210.2 / NA	54	G28G (GGC→GGT)	yagE	2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/26);  total (28/26)

ATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCG  >  NC_000913/282213‑282508                                                                                                                                                     |                                                                                                                                                    aTGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTAc                                                                                                                                                   >  2:746666/1‑151 (MQ=255)   gATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTAcc                                                                                                                                                  <  2:522887/150‑1 (MQ=255)       gTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACcgccgc                                                                                                                                             >  1:538354/1‑151 (MQ=255)                   aacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcgacga                                                                                                                                 <  1:582362/151‑1 (MQ=255)                   aacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcgacga                                                                                                                                 <  2:450445/151‑1 (MQ=255)                   aacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcgacga                                                                                                                                 <  1:758077/151‑1 (MQ=255)                   aacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcgacg                                                                                                                                  >  1:355705/1‑150 (MQ=255)                    acgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGAt                                                                                                                                >  1:364381/1‑151 (MQ=255)                      gCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCTCCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCt                                                                                                                              >  2:199504/1‑151 (MQ=255)                      gCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCt                                                                                                                              >  1:460925/1‑151 (MQ=255)                           tCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATc                                                                                                                          >  2:75044/1‑150 (MQ=255)                                    tAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGc                                                                                                                 <  2:34325/150‑1 (MQ=255)                                     aGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCg                                                                                                                <  2:508550/150‑1 (MQ=255)                                       aaCAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTg                                                                                                             <  1:315494/151‑1 (MQ=255)                                         cAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGAc                                                                                                           <  1:734266/151‑1 (MQ=255)                                          aaaaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACg                                                                                                          >  1:395370/1‑151 (MQ=255)                                          aaaaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACg                                                                                                          >  1:73134/1‑151 (MQ=255)                                              tGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCt                                                                                                      >  1:614285/1‑151 (MQ=255)                                                  tCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGtt                                                                                                   <  1:180489/150‑1 (MQ=255)                                                   cAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGttct                                                                                                 >  1:733978/1‑151 (MQ=255)                                                       aaGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTcc                                                                                              <  2:679502/150‑1 (MQ=255)                                                        aGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTg                                                                                            >  2:283144/1‑151 (MQ=255)                                                        aGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTg                                                                                            >  1:517544/1‑151 (MQ=255)                                                          gagaTCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTggg                                                                                          >  2:708584/1‑151 (MQ=255)                                                           agaTCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGc                                                                                         <  2:409043/151‑1 (MQ=255)                                                                 cATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCg                                                                                     <  1:670237/149‑1 (MQ=255)                                                                 cATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCg                                                                                     <  1:498495/149‑1 (MQ=255)                                                                   tGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTgg                                                                                 >  1:375154/1‑151 (MQ=255)                                                                              gcgTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTccc                                                                      <  1:70207/151‑1 (MQ=255)                                                                               cgTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCACCGGTGGCGAGTTCTccc                                                                      >  1:765653/1‑150 (MQ=255)                                                                                   gttCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTc                                                                 >  1:144897/1‑151 (MQ=255)                                                                                    ttCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCg                                                                >  2:441985/1‑151 (MQ=255)                                                                                     tCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCgg                                                               >  1:369343/1‑151 (MQ=255)                                                                                     tCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCg                                                                <  1:313908/150‑1 (MQ=255)                                                                                           gAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCg                                                          >  1:688934/1‑150 (MQ=255)                                                                                                 ttCCCCCTATCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGTTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAagag                                                     <  2:162443/149‑1 (MQ=255)                                                                                                       cTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGc                                             >  1:515125/1‑151 (MQ=255)                                                                                                         gTCTCCACCATTTTTACCGCTGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGc                                             <  1:5950/149‑1 (MQ=255)                                                                                                         gTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCa                                           <  1:486766/151‑1 (MQ=255)                                                                                                         gTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCa                                           <  1:288313/151‑1 (MQ=255)                                                                                                          tctcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCa                                           <  2:437081/150‑1 (MQ=255)                                                                                                             ccgccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCAtt                                         <  2:698505/146‑1 (MQ=255)                                                                                                             ccaccaTTTTTCCCGCCGAGGGCTAGCTCGATAAGCCGGGTCCCGCCGCGCTGATCGACGATCTGATCAAAGCAGTCGTGGACGGCCTGTCCTTCCTGTGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTg                                        <  2:702568/150‑1 (MQ=255)                                                                                                                ccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGTGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCg                                    >  2:96236/1‑151 (MQ=255)                                                                                                                ccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCg                                    >  2:352507/1‑151 (MQ=255)                                                                                                                ccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCg                                    <  2:726567/151‑1 (MQ=255)                                                                                                                   tttttACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCtt                                 >  1:691594/1‑151 (MQ=255)                                                                                                                   tttttACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCtt                                 >  2:257561/1‑151 (MQ=255)                                                                                                                       tACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGc                              <  2:667787/150‑1 (MQ=255)                                                                                                                           gccgACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCg                         >  2:50842/1‑151 (MQ=255)                                                                                                                            ccgACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGa                        <  1:244499/151‑1 (MQ=255)                                                                                                                                               gCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTg     <  2:279757/151‑1 (MQ=255)                                                                                                                                                  gggTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCg  >  2:492799/1‑151 (MQ=255)                                                                                                                                                  gggTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCg  >  2:734304/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCG  >  NC_000913/282213‑282508

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