Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	283,014	T→C	L246P (CTG→CCG)	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	283,014	0	T	C	100.0%	228.5 / NA	66	L246P (CTG→CCG)	yagE	2‑keto‑3‑deoxy gluconate (KDG) aldolase, CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/33);  total (33/33)

TACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTCGCCGCTGGACGAGCCGCGCAAGGCGCA  >  NC_000913/282883‑283149                                                                                                                                    |                                                                                                                                        tACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAAtt                                                                                                                                  >  2:157934/1‑139 (MQ=255)   cccTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTcc                                                                                                                                <  2:1367691/139‑1 (MQ=255)    ccTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCg                                                                                                                               >  2:266589/1‑139 (MQ=255)      tgctgcTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCa                                                                                                                             <  1:1146749/139‑1 (MQ=255)      tgctgcTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCa                                                                                                                             >  2:1484143/1‑139 (MQ=255)       gctgctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAg                                                                                                                            <  1:266589/139‑1 (MQ=255)        ctgctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGa                                                                                                                           >  2:1814291/1‑139 (MQ=255)           ctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGt                                                                                                                        <  2:1086515/139‑1 (MQ=255)           ctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGt                                                                                                                        <  1:157934/139‑1 (MQ=255)           ctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGt                                                                                                                        >  1:1434161/1‑139 (MQ=255)             cggcggCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTAt                                                                                                                      >  1:1460904/1‑139 (MQ=255)               gcggcgACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCa                                                                                                                    <  1:1312610/139‑1 (MQ=255)                          gCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCa                                                                                                                             >  2:212718/1‑119 (MQ=255)                          gCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCa                                                                                                                             <  1:212718/119‑1 (MQ=255)                                   gCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGcc                                                                                                        >  1:1312690/1‑131 (MQ=255)                                   gCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGcc                                                                                                        <  2:1312690/131‑1 (MQ=255)                                    cGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                                               <  1:716524/139‑1 (MQ=255)                                    cGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                                               <  1:1510902/139‑1 (MQ=255)                                    cGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                                               <  2:726001/139‑1 (MQ=255)                                       gCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGt                                                                                            <  1:1814291/139‑1 (MQ=255)                                        cGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCaaa                                                                                                                                    >  2:242967/1‑98 (MQ=255)                                        cGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCaaa                                                                                                                                    <  1:242967/98‑1 (MQ=255)                                        cGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                                               >  1:1742728/1‑135 (MQ=255)                                        cGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                                               <  2:1742728/135‑1 (MQ=255)                                           cAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGAtt                                                                                        <  1:1177792/139‑1 (MQ=255)                                                ttGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAaga                                                                                   <  1:1370793/139‑1 (MQ=255)                                                   ccccGCAGGTGTCGGTGACTCTTCTGAAAGCCTGGCGCGACGGAGACGTGGCGAAAGCGGCCGGGTATGATCAGACCTGCCCGCACCTTCCGCAGATGTATCAGCTGGGTACGCCGGTTGTTAACGAGATTACAGAGGc                                                                                >  1:363410/1‑139 (MQ=255)                                                      cGCAGTCGTCGGTAAATCTTCTGAAATCCTGGGTCGACGGCGACGTGGCGAACGCGGCCGGGTATCATCAGTCCTTGCCGAAAATCCAGCAGACGTATCAGCTGGATACGCAGTTTGTGAACGTGATTAAAGAGGCGAt                                                                             <  1:1677140/139‑1 (MQ=255)                                                        cAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTcc                                                                                                                                >  1:1866071/1‑86 (MQ=255)                                                        cAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTcc                                                                                                                                <  2:1866071/86‑1 (MQ=255)                                                        cAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                                               <  1:1866557/119‑1 (MQ=255)                                                        cAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                                               >  2:1866557/1‑119 (MQ=255)                                                           gtgtCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCa                                                                                                                                      <  2:1251560/77‑1 (MQ=255)                                                           gtgtCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCa                                                                                                                                      >  1:1251560/1‑77 (MQ=255)                                                                ggTGAATCTTCTGATAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCATATTCCGCAGATGTATCAg                                                                                                                   >  2:1141527/1‑91 (MQ=255)                                                                ggTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAg                                                                                                                   <  1:1141527/91‑1 (MQ=255)                                                                  tGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCgg                                                                 <  2:1434161/139‑1 (MQ=255)                                                                         tctGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCt                                                          >  2:1699290/1‑139 (MQ=255)                                                                           tGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTTTCATCAGACCTTGCCGCAAATTcc                                                                                                                                <  1:1275698/67‑1 (MQ=255)                                                                           tGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTcc                                                                                                                                >  2:1275698/1‑67 (MQ=255)                                                                              aaGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGtctc                                                     >  1:1633864/1‑139 (MQ=255)                                                                                gCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCa                                                                                                                             >  2:894018/1‑65 (MQ=255)                                                                                gCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCa                                                                                                                             <  1:894018/65‑1 (MQ=255)                                                                                gCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGcc                                                                                                        >  1:1845324/1‑86 (MQ=255)                                                                                gCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGcc                                                                                                        <  2:1845324/86‑1 (MQ=255)                                                                                    ggCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCa                                                                                                                             <  1:1173374/61‑1 (MQ=255)                                                                                    ggCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCa                                                                                                                             >  2:1173374/1‑61 (MQ=255)                                                                                    ggCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCcacgca                                               >  2:1087431/1‑139 (MQ=255)                                                                                    ggCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCcacgca                                               >  2:656738/1‑139 (MQ=255)                                                                                           cGGGGAGGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCGGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGtgctg                                        >  1:823925/1‑139 (MQ=255)                                                                                           cGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGtgctg                                        >  1:969211/1‑139 (MQ=255)                                                                                           cGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGtgctg                                        >  1:485655/1‑139 (MQ=255)                                                                                                     gCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGGTGTATCAGCTGGATACGGCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCcctcct                              <  1:1209190/139‑1 (MQ=255)                                                                                                     gCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCt                              <  1:1484143/139‑1 (MQ=255)                                                                                                         aaGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCt                                                          >  2:1716016/1‑107 (MQ=255)                                                                                                         aaGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCt                                                          <  1:1716016/107‑1 (MQ=255)                                                                                                              gCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAg                                                                                                                            >  2:183932/1‑36 (MQ=255)                                                                                                              gCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAg                                                                                                                            <  1:183932/36‑1 (MQ=255)                                                                                                              gCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                                               <  1:1389218/65‑1 (MQ=255)                                                                                                              gCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                                               >  2:1389218/1‑65 (MQ=255)                                                                                                                            gACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTCGCCGCTGGACGAGCCGCGCAAg       >  2:32612/1‑139 (MQ=255)                                                                                                                            gACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTCGCCGCTGGACGAGCCGCGCAAg       >  1:881619/1‑139 (MQ=255)                                                                                                                             aCCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTc                             <  1:1727681/116‑1 (MQ=255)                                                                                                                             aCCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTc                             >  2:1727681/1‑116 (MQ=255)                                                                                                                              ccTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTCGCCGCTGGACGAGCCGCGCAAGgc     >  2:595537/1‑139 (MQ=255)                                                                                                                               cTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCTACGCACGTGCTGCCGCCCGCCTCGCCGCTGGACGAGCCGCGCAAGgcg    <  2:823925/139‑1 (MQ=255)                                                                                                                                 tgccgcAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTCGCCGCTGGACGAGCCGCGCAAGGCGCa  >  1:819090/1‑139 (MQ=255)                                                                                                                                 tgccgcAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTCGCCGCTGGACGAGCCGCGCAAGGCGCa  >  1:1018710/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                        TACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTCGCCGCTGGACGAGCCGCGCAAGGCGCA  >  NC_000913/282883‑283149

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