Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	289,154	A→C	I3M (ATT→ATG)	yagI ←	CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	289,154	0	A	C	100.0%	203.1 / NA	56	I3M (ATT→ATG)	yagI	CP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (29/27);  total (29/27)

AGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGC  >  NC_000913/289009‑289292                                                                                                                                                  |                                                                                                                                           aGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCgg                                                                                                                                       >  2:261524/1‑151 (MQ=255) aGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCgg                                                                                                                                       >  2:409248/1‑151 (MQ=255)  gCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGc                                                                                                                                      <  2:597413/151‑1 (MQ=255)     gTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCAt                                                                                                                                    >  2:692095/1‑150 (MQ=255)       gAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATaa                                                                                                                                  <  2:775605/150‑1 (MQ=255)              cAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCa                                                                                                                          <  2:605803/151‑1 (MQ=255)                  gTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTc                                                                                                                      >  1:770993/1‑151 (MQ=255)                  gTTTTTAGCAGCGAGTGGAGGGAACTCATGCTCAGCCCCCTCCGTTTGCTGATGTCGGTGATCTTAAGCTCGGTTGCCTGCTCGTTTAACAGGTCTGGGTTCTGCAACGCACGTTCAACAGACTGAATCCACGGCAGAATGCCTGCATGTc                                                                                                                      >  1:329532/1‑151 (MQ=255)                     tttAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCa                                                                                                                    <  2:418805/150‑1 (MQ=255)                      ttAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCAc                                                                                                                   <  2:518835/150‑1 (MQ=255)                      ttAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCAc                                                                                                                   >  1:681509/1‑150 (MQ=255)                           agcGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTgg                                                                                                              >  2:213946/1‑150 (MQ=255)                           agcGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGa                                                                                                             <  1:406619/151‑1 (MQ=255)                                  ggAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAg                                                                                                       <  2:902030/150‑1 (MQ=255)                                  ggAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAg                                                                                                       <  1:286128/150‑1 (MQ=255)                                   gAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGc                                                                                                     <  2:373541/151‑1 (MQ=255)                                     gggTACTCTTGCTCAGCCCCATCAGTTTGCTGAAGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCg                                                                                                    >  2:923623/1‑150 (MQ=255)                                         aCTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAAc                                                                                               >  2:887897/1‑151 (MQ=255)                                               gCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCt                                                                                         <  1:24050/151‑1 (MQ=255)                                                  cAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCttt                                                                                       <  1:209559/150‑1 (MQ=255)                                                    gCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCttttt                                                                                     <  1:777626/150‑1 (MQ=255)                                                       ccATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCggg                                                                                 <  1:814179/151‑1 (MQ=255)                                                          tCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATg                                                                              >  1:848084/1‑151 (MQ=255)                                                          tCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATg                                                                              >  2:929252/1‑151 (MQ=255)                                                           cAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGa                                                                             <  2:133619/151‑1 (MQ=255)                                                            aGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGaa                                                                            <  1:874177/151‑1 (MQ=255)                                                                 gCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAAt                                                                        <  1:553613/150‑1 (MQ=255)                                                                     aTGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATg                                                                   <  2:257724/151‑1 (MQ=255)                                                                          ggTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATg                                                              <  1:494290/151‑1 (MQ=255)                                                                           gTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATgg                                                             >  2:288826/1‑151 (MQ=255)                                                                            tGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATgg                                                             >  2:932726/1‑150 (MQ=255)                                                                            tGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATgg                                                             >  2:206036/1‑150 (MQ=255)                                                                            tGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATgg                                                             <  1:614881/150‑1 (MQ=255)                                                                             gATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATgg                                                             >  1:602172/1‑149 (MQ=255)                                                                                  tAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCg                                                        <  1:759256/149‑1 (MQ=255)                                                                                     gCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGtt                                                   >  2:271338/1‑151 (MQ=255)                                                                                       tCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTtct                                                 >  2:488076/1‑151 (MQ=255)                                                                                            ggCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTtctccctc                                            <  2:337193/151‑1 (MQ=255)                                                                                            ggCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTtctccctc                                            <  2:169307/151‑1 (MQ=255)                                                                                              ccTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTtctccctctc                                          <  2:329365/151‑1 (MQ=255)                                                                                                   tCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCtgt                                     >  1:165551/1‑151 (MQ=255)                                                                                                     gTTGAACAGGTCGAGGATCCGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTgg                                   <  2:411675/151‑1 (MQ=255)                                                                                                       tGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTggg                                  >  1:70794/1‑150 (MQ=255)                                                                                                        gAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgag                                <  2:554069/151‑1 (MQ=255)                                                                                                         aaCAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTggg                                  >  2:12632/1‑148 (MQ=255)                                                                                                          aCAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgag                                >  1:482339/1‑149 (MQ=255)                                                                                                          aCAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGagag                              <  2:69027/151‑1 (MQ=255)                                                                                                            aGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGATGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTATCCCTGTGGGagag                              >  2:607955/1‑149 (MQ=255)                                                                                                                  agcttCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGACATGTCCACGCTGGAATTAAGGGGAAAATGTACCTGTTTCGGGATGAAAAATCAATAAAATGGAGCTGGTGTTCTCCCTCTCCCTGTGGGAGAGGGcc                          <  1:517318/143‑1 (MQ=255)                                                                                                                       cTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                 >  2:583983/1‑151 (MQ=255)                                                                                                                               cACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAc         >  1:408777/1‑151 (MQ=255)                                                                                                                                 cGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAcgc       >  2:241526/1‑151 (MQ=255)                                                                                                                                   ttCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACgcgc     >  2:398705/1‑151 (MQ=255)                                                                                                                                     cAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACgcgc     >  1:125477/1‑149 (MQ=255)                                                                                                                                      aaCAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGc  >  2:490122/1‑151 (MQ=255)                                                                                                                                      aaCAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGc  >  1:564543/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                           AGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGC  >  NC_000913/289009‑289292

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