Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	268,138	C→T	V275M (GTG→ATG)	ykfA ←	CP4‑6 prophage; putative GTP‑binding protein YkfA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	268,138	0	C	T	100.0%	185.3 / NA	56	V275M (GTG→ATG)	ykfA	CP4‑6 prophage; putative GTP‑binding protein YkfA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/28);  total (28/28)

CAGACTCCTCATAACGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCACCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGC  >  NC_000913/268009‑268269                                                                                                                                  |                                                                                                                                    cAGACTCCTCATAACGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTgcgc                                                                                                                           <  2:356687/140‑1 (MQ=255) cAGACTCCTCATAACGGGAAAACAAATAAAAGGGAACCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTgcgc                                                                                                                           <  1:941953/140‑1 (MQ=255)  aGACTCCTCATAACGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCa                                                                                                                          <  2:440313/140‑1 (MQ=255)  aGACTCCTCATAACGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCa                                                                                                                          <  1:459091/140‑1 (MQ=255)   gACTCCTCATAACGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCAc                                                                                                                         >  2:905881/1‑140 (MQ=255)   gACTCCTCATAACGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCAc                                                                                                                         >  2:772545/1‑140 (MQ=255)   gACTCCTCATAACGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCAc                                                                                                                         >  1:712045/1‑140 (MQ=255)    aCTCCTCATAACGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACg                                                                                                                        >  1:1003045/1‑140 (MQ=255)    aCTCCTCATAACGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACg                                                                                                                        <  2:953512/140‑1 (MQ=255)    aCTCCTCATAACGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACg                                                                                                                        >  1:933344/1‑140 (MQ=255)     ctcctcATAACGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACgg                                                                                                                       >  2:690913/1‑140 (MQ=255)        ctcATAACGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCt                                                                                                                    >  2:501233/1‑140 (MQ=255)             aaCGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATg                                                                                                               <  2:590982/140‑1 (MQ=255)                gggAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACg                                                                                                            >  2:966315/1‑140 (MQ=255)                  gAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGc                                                                                                          >  2:318912/1‑140 (MQ=255)                  gAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGc                                                                                                          >  2:153556/1‑140 (MQ=255)                   aaaaCAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGcc                                                                                                         <  2:753727/140‑1 (MQ=255)                         aaTAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCAcc                                                                                                   <  1:368041/140‑1 (MQ=255)                          aTAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCg                                                                                                  >  1:341232/1‑140 (MQ=255)                          aTAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCg                                                                                                  <  1:432403/140‑1 (MQ=255)                                gggATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCgg                                                                                            >  2:852219/1‑140 (MQ=255)                                    tCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTAt                                                                                        >  1:964782/1‑140 (MQ=255)                                     ccccGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATc                                                                                       <  2:523263/140‑1 (MQ=255)                                     ccccGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATc                                                                                       <  2:781201/140‑1 (MQ=255)                                        cGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAgg                                                                                    >  2:547946/1‑140 (MQ=255)                                          tCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTAACAgggg                                                                                  >  1:681354/1‑140 (MQ=255)                                           cGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGa                                                                                 <  2:341232/140‑1 (MQ=255)                                             cATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGaaa                                                                               <  1:1029414/140‑1 (MQ=255)                                                  gCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCAc                                                                          <  1:966315/140‑1 (MQ=255)                                                    gACAGGGGCAGGGAAGTAGCAGGGATGGGTTAAATACTCAGAAGAAGAATTCCCGGACGGCGCGGGCCACTGCCCCCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTg                                                                        <  1:466309/140‑1 (MQ=255)                                                    gACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTg                                                                        <  1:590107/140‑1 (MQ=255)                                                       gggggCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATg                                                                     <  2:722413/139‑1 (MQ=255)                                                          ggCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTa                                                                  <  2:561241/140‑1 (MQ=255)                                                                  aGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTAc                                                          <  2:933344/140‑1 (MQ=255)                                                                    tAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACgg                                                        >  1:846288/1‑140 (MQ=255)                                                                    tAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGGACgg                                                        >  2:683611/1‑140 (MQ=255)                                                                         gggATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCg                                                   >  1:252734/1‑140 (MQ=255)                                                                         gggATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCg                                                   >  1:273548/1‑140 (MQ=255)                                                                          ggATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGa                                                  >  2:522843/1‑140 (MQ=255)                                                                                         tCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAAt                                   <  2:720899/140‑1 (MQ=255)                                                                                         tCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAAt                                   >  2:137771/1‑140 (MQ=255)                                                                                         tCAGAAGAAGAAATCCCAGACGGCACGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAAt                                   >  2:770987/1‑140 (MQ=255)                                                                                           agaagaagaAATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCg                                 >  2:44546/1‑140 (MQ=255)                                                                                              agaagaaATCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCGCTa                              >  1:545476/1‑140 (MQ=255)                                                                                                   aaaTCCCAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGc                         <  1:852219/140‑1 (MQ=255)                                                                                                       cccAGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCCATCGAGTACGGCCCCGACGGCTTCACCAAAATCGCAACGAGCCTGt                     >  1:640689/1‑140 (MQ=255)                                                                                                          aGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTcc                  <  2:712045/140‑1 (MQ=255)                                                                                                          aGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTcc                  <  2:363520/140‑1 (MQ=255)                                                                                                          aGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTcc                  <  2:273548/140‑1 (MQ=255)                                                                                                          aGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTcc                  <  2:846288/140‑1 (MQ=255)                                                                                                          aGACGGCGCGGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTcc                  <  2:252734/140‑1 (MQ=255)                                                                                                                cgcgGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGGAAAGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCCAAATCGCTACGCGCCTGTTCCCGGACt            >  1:891879/1‑140 (MQ=25)                                                                                                                cgcgGGCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACt            >  1:227645/1‑140 (MQ=255)                                                                                                                     gCCACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACcgt       <  1:1002077/140‑1 (MQ=255)                                                                                                                      ccACTGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACcgtc      <  1:935210/140‑1 (MQ=255)                                                                                                                          tGACACCATCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGc  >  2:316346/1‑140 (MQ=255)                                                                                                                                  |                                                                                                                                    CAGACTCCTCATAACGGGAAAACAAATAAAAGGGATCCCCGTCGCATCGGCGACAGGGGCAGGGAAGTAACAGGGATGGGTTAAATACTCAGAAGAAGAAATCCCAGACGGCGCGGGCCACTGACACCACCGTGGTGCGCACGGCCTGAATGACGGCCCGCACCGGGGCGGGTATCAGGGGAAAGGCACTGATGCTATCGAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGC  >  NC_000913/268009‑268269

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