Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	282,311	C→T	P12S (CCC→TCC)	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,311	0	C	T	100.0%	192.5 / NA	58	P12S (CCC→TCC)	yagE	CP4‑6 prophage; putative 2‑keto‑3‑deoxygluconate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/21);  total (37/21)

TCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCC  >  NC_000913/282180‑282440                                                                                                                                    |                                                                                                                                  tccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGtc                                                                                                                           >  2:946387/1‑140 (MQ=255) tccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGtc                                                                                                                           >  1:155667/1‑140 (MQ=255) tccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGtc                                                                                                                           >  2:944399/1‑140 (MQ=255) tccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGtc                                                                                                                           >  1:24651/1‑140 (MQ=255) tccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGt                                                                                                                            <  2:323318/139‑1 (MQ=255) tccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGt                                                                                                                            >  1:323318/1‑139 (MQ=255)    gATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTcc                                                                                                                        <  2:24651/140‑1 (MQ=255)     aTTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTcca                                                                                                                       >  2:425231/1‑140 (MQ=255)                       tGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTGACCGCCGACg                                                                                                     >  2:853400/1‑140 (MQ=255)                         cGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGc                                                                                                   <  2:194952/140‑1 (MQ=255)                         cGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGc                                                                                                   >  1:568223/1‑140 (MQ=255)                              tAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCt                                                                                              >  2:433083/1‑140 (MQ=255)                                gTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCg                                                                                            <  1:433083/140‑1 (MQ=255)                                  aTGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGAt                                                                                          <  2:362085/140‑1 (MQ=255)                                   tGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAg                                                                                                >  2:313046/1‑133 (MQ=255)                                   tGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAg                                                                                                <  1:313046/133‑1 (MQ=255)                                      tCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGc                                                                                      >  2:112810/1‑140 (MQ=255)                                      tCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGc                                                                                      >  1:277383/1‑140 (MQ=255)                                      tCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGc                                                                                      >  1:867499/1‑140 (MQ=255)                                         ttAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCgg                                                                                   >  2:785956/1‑140 (MQ=255)                                         ttAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCgg                                                                                   >  2:304985/1‑140 (MQ=255)                                         ttAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCgg                                                                                   >  1:198526/1‑140 (MQ=255)                                           aaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGc                                                                                 >  1:876167/1‑140 (MQ=255)                                                 acgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCAccgcc                                                                           >  2:119791/1‑140 (MQ=255)                                                   gaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACcgccgc                                                                         <  2:198526/140‑1 (MQ=255)                                                      cgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCCCCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCt                                                                      <  1:944399/140‑1 (MQ=255)                                                      cgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCt                                                                      <  2:568223/140‑1 (MQ=255)                                                      cgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCt                                                                      <  1:830915/140‑1 (MQ=255)                                                        cTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCt                                                                      <  2:626683/138‑1 (MQ=255)                                                        cTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCt                                                                      >  1:626683/1‑138 (MQ=255)                                                           ttCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATcg                                                                 <  2:867499/140‑1 (MQ=255)                                                                     tAGAACAAAATGATTCAGCAAGGAGATCTCATGTCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGAt                                                       >  1:137587/1‑140 (MQ=255)                                                                     tAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGAt                                                       >  1:647044/1‑140 (MQ=255)                                                                     tAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGAt                                                       >  2:265181/1‑140 (MQ=255)                                                                          cAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATcg                                                                 >  2:193820/1‑125 (MQ=255)                                                                          cAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATcg                                                                 <  1:193820/125‑1 (MQ=255)                                                                          cAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAg                                                  >  2:778559/1‑140 (MQ=255)                                                                            aaaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCa                                                <  1:853400/140‑1 (MQ=255)                                                                             aaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAg                                               >  1:421984/1‑140 (MQ=255)                                                                                   tCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTg                                         <  1:304985/140‑1 (MQ=255)                                                                                    cAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGa                                        >  1:199727/1‑140 (MQ=255)                                                                                     aGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGAc                                       >  2:101352/1‑140 (MQ=255)                                                                                       cAAGGAGATCTCATGTCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACgg                                     <  2:137587/140‑1 (MQ=255)                                                                                       cAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACgg                                     <  1:286178/140‑1 (MQ=255)                                                                                            agaTCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGt                                <  1:112810/140‑1 (MQ=255)                                                                                                 tcATGCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGttcttc                           >  2:22686/1‑140 (MQ=255)                                                                                                     gCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTgg                       >  1:645715/1‑140 (MQ=255)                                                                                                     gCCGCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTgg                       >  1:175133/1‑140 (MQ=255)                                                                                                        gCAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCa                    >  1:267627/1‑140 (MQ=255)                                                                                                         cAGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAg                   >  2:396546/1‑140 (MQ=255)                                                                                                          aGTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGc                  >  1:336905/1‑140 (MQ=255)                                                                                                           gTCCGCGTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAACCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCg                 >  1:103650/1‑140 (MQ=255)                                                                                                              cgcgTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCg                 <  2:641228/137‑1 (MQ=255)                                                                                                              cgcgTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCg                 >  1:641228/1‑137 (MQ=255)                                                                                                               gcgTTGTTCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTgg             <  1:785956/140‑1 (MQ=255)                                                                                                                   tgttCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAg         <  2:645715/140‑1 (MQ=255)                                                                                                                   tgttCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAg         <  2:267627/140‑1 (MQ=255)                                                                                                                      tCACGGGAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTtc      >  2:441194/1‑140 (MQ=255)                                                                                                                          gggAATCATTTCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTccc  >  2:864495/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                  TCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCC  >  NC_000913/282180‑282440

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