Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	282,806	C→T	R177C (CGC→TGC)	yagE →	CP4‑6 prophage; putative 2‑keto‑3‑deoxygluconate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	282,806	0	C	T	100.0%	189.0 / NA	55	R177C (CGC→TGC)	yagE	CP4‑6 prophage; putative 2‑keto‑3‑deoxygluconate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/31);  total (24/31)

CACGCTGCCGGTGATGCTCTATAACTTCCCGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGC  >  NC_000913/282670‑282933                                                                                                                                         |                                                                                                                                cACGCTGCCGGTGATGCTCTATAACTTCCCGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTgca                                                                                                                              >  1:612773/1‑140 (MQ=255)   cGCTGCCGGTGATGCTCTATAACTTCCCGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTgcagc                                                                                                                            <  1:54007/140‑1 (MQ=255)   cGCTGCCGGTGATGCTCTATAACTTCCCGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTgcagc                                                                                                                            <  1:638350/140‑1 (MQ=255)   cGCTGCCGGTGATGCTCTATAACTTCCCGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTgcagc                                                                                                                            <  2:325786/140‑1 (MQ=255)   cGCTGCCGGTGATGCTCTATAACTTCCCGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTgcagc                                                                                                                            <  1:330809/140‑1 (MQ=255)        ccgacgATGCTCTATAACTTCCCGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGAt                                                                                                                       >  2:407922/6‑140 (MQ=255)                  tctATAACTTCCCGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTc                                                                                                             >  1:287901/1‑140 (MQ=255)                      taACTTCCCGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAg                                                                                                         >  1:1023375/1‑140 (MQ=255)                      taACTTCCCGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAg                                                                                                         >  1:351121/1‑140 (MQ=255)                         cTTCCCGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTg                                                                                                      <  2:1144833/140‑1 (MQ=255)                           tCCCGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGcc                                                                                                    <  2:351121/140‑1 (MQ=255)                            cccGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGccc                                                                                                   >  2:11145/1‑140 (MQ=255)                              cGGCGCTGACCGGGCAGGATCTGACTTCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCAt                                                                                                 <  2:287901/140‑1 (MQ=255)                              cGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATGCCGTCAAAGGTGCCCAt                                                                                                 >  2:1089992/1‑140 (MQ=255)                              cGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGt                                                                                                              <  2:153176/127‑1 (MQ=255)                              cGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGt                                                                                                              >  1:153176/1‑127 (MQ=255)                              cGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCAt                                                                                                 <  2:189871/140‑1 (MQ=255)                              cGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCAt                                                                                                 >  2:607557/1‑140 (MQ=255)                              cGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCAt                                                                                                 >  1:638006/1‑140 (MQ=255)                              cGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCAt                                                                                                 <  1:406786/140‑1 (MQ=255)                                gcgcTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATcc                                                                                               <  2:638006/140‑1 (MQ=255)                                 cgcTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCg                                                                                              <  2:990556/140‑1 (MQ=255)                                   cTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTg                                                                                                      <  2:54618/130‑1 (MQ=255)                                   cTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTg                                                                                                      >  1:54618/1‑130 (MQ=255)                                     gACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACt                                                                                          <  1:407922/140‑1 (MQ=255)                                      aCCGGGCAGGATCTGACTCCGGCGCTGGTGAAAGCCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCa                                                                                            <  2:505328/137‑1 (MQ=255)                                      aCCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCa                                                                                            >  1:505328/1‑137 (MQ=255)                                          ggCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCa                                                                                            <  2:943020/133‑1 (MQ=255)                                          ggCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCa                                                                                            >  1:943020/1‑133 (MQ=255)                                                aTCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTg                                                                                  <  1:173445/137‑1 (MQ=255)                                                aTCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTg                                                                                  >  2:173445/1‑137 (MQ=255)                                                 tCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGctct                                                                              <  2:170863/140‑1 (MQ=255)                                                     aCTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCgg                                                                          >  2:1086018/1‑140 (MQ=255)                                                      cTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGc                                                                         <  1:584239/140‑1 (MQ=255)                                                        ccGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTa                                                                       <  1:1089992/140‑1 (MQ=255)                                                         cGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTac                                                                      >  1:312232/1‑140 (MQ=255)                                                         cGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTa                                                                       >  2:1138904/1‑139 (MQ=255)                                                         cGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTa                                                                       <  1:1138904/139‑1 (MQ=255)                                                             gcTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTacgacg                                                                  <  2:312232/140‑1 (MQ=255)                                                               tGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGat                                                                <  2:612773/140‑1 (MQ=255)                                                                 gTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTC‑CCGTGCTCTGCGGCTACGACGatcat                                                             <  1:11145/140‑1 (MQ=255)                                                                  tGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGatcat                                                             <  1:607557/140‑1 (MQ=255)                                                                           tCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAAt                                                    >  2:774481/1‑140 (MQ=255)                                                                           tCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAAt                                                    >  2:705497/1‑140 (MQ=255)                                                                                   tCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGatcat                                                             >  1:305062/1‑123 (MQ=255)                                                                                   tCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGatcat                                                             <  2:305062/123‑1 (MQ=255)                                                                                   tCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCctgct                                            >  1:1215857/1‑140 (MQ=255)                                                                                               attatCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGa                                <  1:1086018/140‑1 (MQ=255)                                                                                                     ggCATCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGc                          >  2:397115/1‑140 (MQ=255)                                                                                                        aTCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGat                       >  2:649892/1‑140 (MQ=255)                                                                                                        aTCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGat                       >  1:455020/1‑140 (MQ=255)                                                                                                         tCAAAGACACCATCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGata                      >  1:597741/1‑140 (MQ=255)                                                                                                                    aTCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCgg           <  1:97569/140‑1 (MQ=255)                                                                                                                    aTCGACTCCGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCgg           <  1:649892/140‑1 (MQ=255)                                                                                                                            cGTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTg   <  2:1215857/140‑1 (MQ=255)                                                                                                                             gTCGCCCACCTGTGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGc  <  1:397115/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                CACGCTGCCGGTGATGCTCTATAACTTCCCGGCGCTGACCGGGCAGGATCTGACTCCGGCGCTGGTGAAAACCCTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGC  >  NC_000913/282670‑282933

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