Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	282,993	C→T	A239V (GCC→GTC)	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,993	0	C	T	100.0%	135.8 / NA	36	A239V (GCC→GTC)	yagE	2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/20);  total (16/20)

GCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTCGCCGCTGGACG  >  NC_000913/282858‑283133                                                                                                                                        |                                                                                                                                             gCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAg                                                                                                                                 <  2:188527/149‑1 (MQ=255) gCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGAc                                                                                                                               >  2:564956/1‑151 (MQ=255) gCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGAc                                                                                                                               <  1:556197/151‑1 (MQ=255)  cGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGAcc                                                                                                                              <  1:355744/151‑1 (MQ=255)   ggCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCt                                                                                                                             >  2:472226/1‑151 (MQ=255)          acgaTCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGTACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCa                                                                                                                      <  1:410220/151‑1 (MQ=255)                         aTACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGt                                                                                                        >  2:50860/1‑150 (MQ=255)                                  tgctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTgg                                                                                               >  2:493259/1‑150 (MQ=255)                                   gctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAACCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGAt                                                                                             <  2:448739/151‑1 (MQ=255)                                           gcgACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGtt                                                                                     <  1:153126/151‑1 (MQ=255)                                                  ggCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAAc                                                                              >  2:550681/1‑151 (MQ=255)                                                   gCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACg                                                                             >  1:432995/1‑151 (MQ=255)                                                     gATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTg                                                                           <  1:20181/151‑1 (MQ=255)                                                     gATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTg                                                                           <  2:318935/151‑1 (MQ=255)                                                      atatCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGa                                                                          <  1:371235/151‑1 (MQ=255)                                                          cGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTa                                                                       >  1:76279/1‑150 (MQ=255)                                                              gAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAagag                                                                  <  1:366846/151‑1 (MQ=255)                                                              gAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAag                                                                    >  1:395038/1‑149 (MQ=255)                                                               aGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTGTCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAgg                                                                 <  2:38052/151‑1 (MQ=255)                                                                   gCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGAt                                                             <  2:422835/151‑1 (MQ=255)                                                                         ttGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGATCGTGATTAAAGAGGCGATCGTGct                                                       >  2:353213/1‑151 (MQ=255)                                                                                 cAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGgtc                                               <  2:550099/151‑1 (MQ=255)                                                                                   ggTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGgtcgt                                             <  2:212673/151‑1 (MQ=255)                                                                                      gtCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTcc                                           <  1:50861/150‑1 (MQ=255)                                                                                      gtCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCt                                          >  2:129896/1‑151 (MQ=255)                                                                                      gtCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCt                                          >  1:283695/1‑151 (MQ=255)                                                                                      gtCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCt                                          >  2:113459/1‑151 (MQ=255)                                                                                                        aGCCTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGtgctg                        >  1:425171/1‑151 (MQ=255)                                                                                                          ccTGGCGCGACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTgcc                      <  2:515915/151‑1 (MQ=255)                                                                                                                gcgACGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTTATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGccgcccgc                >  2:39391/1‑151 (MQ=255)                                                                                                                   aCGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTc             <  2:43806/151‑1 (MQ=255)                                                                                                                   aCGGGGACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTc             <  2:247140/151‑1 (MQ=255)                                                                                                                       ggACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTcgccg         <  2:616264/151‑1 (MQ=255)                                                                                                                        gACGTGGCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTcgccgc        >  2:563197/1‑151 (MQ=255)                                                                                                                              gCGAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTCGCCGCTGGACg  <  1:169253/151‑1 (MQ=255)                                                                                                                                gAAAGCGGTCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTCGCCGCTGGACg  >  1:260866/1‑149 (MQ=255)                                                                                                                                        |                                                                                                                                             GCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGAGGCGATCGTGCTCTGCGGTCGTCCTGTCTCCACGCACGTGCTGCCGCCCGCCTCGCCGCTGGACG  >  NC_000913/282858‑283133

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