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%	168.7 / NA	48	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 (23/25);  total (23/25)

CAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGC  >  NC_000913/282880‑283119                                                                                                                                       |                                                                                                          cAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCaa                                                                                                       >  1:755680/1‑139 (MQ=255) cAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCaa                                                                                                       >  2:474564/1‑139 (MQ=255)    tACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAAtt                                                                                                    >  2:899600/1‑139 (MQ=255)    tACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAAtt                                                                                                    >  2:1144794/1‑139 (MQ=255)      cccTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTcc                                                                                                  <  1:899600/139‑1 (MQ=255)      cccTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTcc                                                                                                  <  2:1149952/139‑1 (MQ=255)       ccTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCg                                                                                                 <  2:532738/139‑1 (MQ=255)       ccTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCg                                                                                                 <  1:474564/139‑1 (MQ=255)         tgctgcTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGTATCTTCTGAAAGCCTGGCGCGACGCGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAACTTCCGCa                                                                                               <  2:745860/139‑1 (MQ=255)         tgctgcTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCa                                                                                               <  1:330993/139‑1 (MQ=255)          gctgctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAg                                                                                              <  2:807931/139‑1 (MQ=255)           ctgctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGa                                                                                             >  1:541041/1‑139 (MQ=255)              ctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGt                                                                                          <  1:1144794/139‑1 (MQ=255)              ctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGt                                                                                          >  2:159993/1‑139 (MQ=255)                 ggcggcGACGGGGCTATATCGGCGAGCGGCAACATTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCTACGGGACCGTGGCGAAACCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATc                                                                                       >  1:460597/1‑139 (MQ=255)                    ggcgACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGc                                                                                    <  2:627224/139‑1 (MQ=255)                    ggcgACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGc                                                                                    <  1:943381/139‑1 (MQ=255)                             gCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGc                                                                           <  2:54617/139‑1 (MQ=255)                              cGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGcc                                                                          <  1:159993/139‑1 (MQ=255)                                   tCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTtg                                                                     <  2:665399/139‑1 (MQ=255)                                       cGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                 <  2:82369/139‑1 (MQ=255)                                           cGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTg                                                             >  2:536450/1‑139 (MQ=255)                                           cGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTg                                                             >  2:729646/1‑139 (MQ=255)                                                aCTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTaa                                                        >  2:827136/1‑139 (MQ=255)                                                          gCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATc                                              <  2:384165/139‑1 (MQ=255)                                                           cAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGa                                                                  <  2:594257/118‑1 (MQ=255)                                                           cAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGa                                                                  >  1:594257/1‑118 (MQ=255)                                                             ggTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTg                                           <  2:541041/139‑1 (MQ=255)                                                                     tGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                 <  1:1074979/109‑1 (MQ=255)                                                                     tGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                 >  2:1074979/1‑109 (MQ=255)                                                                         tcttctGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGgt                                  >  2:975350/1‑136 (MQ=255)                                                                         tcttctGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGgt                                  <  1:975350/136‑1 (MQ=255)                                                                                    ccTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATa                                                                              >  2:1199199/1‑81 (MQ=255)                                                                                    ccTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATa                                                                              <  1:1199199/81‑1 (MQ=255)                                                                                     cTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                 <  1:104583/93‑1 (MQ=255)                                                                                     cTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                 >  2:104583/1‑93 (MQ=255)                                                                                       ggCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCa                                                                                               >  2:492358/1‑61 (MQ=255)                                                                                       ggCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCa                                                                                               <  1:492358/61‑1 (MQ=255)                                                                                        gcgcGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCg                                                                                                 >  2:943576/1‑58 (MQ=255)                                                                                        gcgcGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCg                                                                                                 <  1:943576/58‑1 (MQ=255)                                                                                                   aCGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTggcggc     >  1:1083006/1‑137 (MQ=255)                                                                                                   aCGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTgccgcc     >  1:890106/1‑139 (MQ=255)                                                                                                      tGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGTGGTCGTCCTGTCTCCACGCACGTGCTGccgcccgc  >  2:314751/1‑139 (MQ=255)                                                                                                      tGGCGAAAGCGGCCGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGccgcccgc  >  2:1185879/1‑139 (MQ=255)                                                                                                                  ccGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATc                                              >  1:509433/1‑83 (MQ=255)                                                                                                                  ccGGGTATCATCAGACCTTGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATc                                              <  2:509433/83‑1 (MQ=255)                                                                                                                                   ttGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTtg                                                                     >  1:373622/1‑43 (MQ=255)                                                                                                                                   ttGCCGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTtg                                                                     <  2:373622/43‑1 (MQ=255)                                                                                                                                    tgccgcAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa                                                                 <  1:719387/46‑1 (MQ=255)                                                                                                                                       |                                                                                                          CAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGC  >  NC_000913/282880‑283119

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: ‑