Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	285,648	C→G	A85G (GCC→GGC)	yagG →	CP4‑6 prophage; putative sugar transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	285,648	0	C	G	100.0%	173.8 / NA	48	A85G (GCC→GGC)	yagG	CP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (23/25);  total (23/25)

CCTGTCGGCGGGGATTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGC  >  NC_000913/285511‑285793                                                                                                                                          |                                                                                                                                                  ccTGTCGGCGGGGATTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGc                                                                                                                                      <  1:67372/151‑1 (MQ=255) ccTGTCGGCGGGGATTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGc                                                                                                                                      >  1:592605/1‑151 (MQ=255)     tCGGCGGGGATTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCCCCCGGTCCTGCTGTGGGGGGGCATCCCGTTCGGCATc                                                                                                                                   >  2:110755/1‑150 (MQ=255)          ggggATTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGc                                                                                                                             >  2:565153/1‑151 (MQ=255)              aTTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGc                                                                                                                         >  1:573214/1‑151 (MQ=255)                 aTGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGCCCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCt                                                                                                                        <  2:343971/149‑1 (MQ=255)                 aTGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTg                                                                                                                       <  1:133581/150‑1 (MQ=255)                  tGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGAc                                                                                                                     >  1:252226/1‑151 (MQ=255)                            gTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTaca                                                                                                             <  1:5126/149‑1 (MQ=255)                             tttttGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACg                                                                                                           <  1:304311/150‑1 (MQ=255)                                  ggTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGAc                                                                                                     >  1:546107/1‑151 (MQ=255)                                   gTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACt                                                                                                    <  2:437310/151‑1 (MQ=255)                                      tcCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTtct                                                                                                 >  2:419180/1‑151 (MQ=255)                                      tcCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTtc                                                                                                  >  2:483145/1‑150 (MQ=255)                                        ccGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTcc                                                                                               <  1:96766/151‑1 (MQ=255)                                          gcgTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGc                                                                                             >  1:49547/1‑151 (MQ=255)                                                  gACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGc                                                                                      >  2:628119/1‑150 (MQ=255)                                                    cGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCa                                                                                     >  1:580592/1‑149 (MQ=255)                                                      ccGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGat                                                                                 <  1:525581/151‑1 (MQ=255)                                                      ccGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGat                                                                                 <  1:535285/151‑1 (MQ=255)                                                          cACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGatcatc                                                                             >  1:435726/1‑151 (MQ=255)                                                             cGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTAc                                                                          >  1:476659/1‑151 (MQ=255)                                                               aCCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGc                                                                        >  2:240253/1‑151 (MQ=255)                                                                 ccGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCt                                                                      >  1:381942/1‑151 (MQ=255)                                                                      gATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCg                                                                   <  1:5708/149‑1 (MQ=255)                                                                              tgctgGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACAt                                                         >  1:615088/1‑151 (MQ=255)                                                                               gctgGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACAtt                                                        >  1:7427/1‑151 (MQ=255)                                                                               gctgGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACAtt                                                        <  1:97846/151‑1 (MQ=255)                                                                               gcgggTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACa                                                          <  2:449990/146‑1 (MQ=255)                                                                                  ggTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATtct                                                      <  2:7427/150‑1 (MQ=255)                                                                                         cgcacccgcacGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGa                                                                                  <  1:234018/115‑1 (MQ=255)                                                                                         cgcacccgcacGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGa                                                                                  >  2:234013/1‑115 (MQ=255)                                                                                                  acGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGCGACCTACATTCTCCTGACCCTGGTCTac                                      >  1:303664/1‑150 (MQ=255)                                                                                                     cgGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCt                                  <  1:573032/151‑1 (MQ=255)                                                                                                       gCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTc                                <  1:664337/151‑1 (MQ=255)                                                                                                          cGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGt                              <  2:647401/150‑1 (MQ=255)                                                                                                           ggCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGtt                             >  2:332013/1‑150 (MQ=255)                                                                                                             ccAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAAc                          <  2:476287/151‑1 (MQ=255)                                                                                                              cAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACg                         <  1:483157/151‑1 (MQ=255)                                                                                                                  tCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGc                      <  1:299234/150‑1 (MQ=255)                                                                                                                  tCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGc                      <  2:646541/150‑1 (MQ=255)                                                                                                                   ccGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCg                    <  2:263856/151‑1 (MQ=255)                                                                                                                      cccGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTAc                 <  2:127968/151‑1 (MQ=255)                                                                                                                      cccGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTAc                 >  2:563269/1‑151 (MQ=255)                                                                                                                          ttCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTAAGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTgcg             >  2:291412/1‑151 (MQ=255)                                                                                                                           tCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTgcg             <  2:595878/150‑1 (MQ=255)                                                                                                                                    ggggggGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCggg   <  1:62552/151‑1 (MQ=255)                                                                                                                                    ggggggGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCggg   >  2:328017/1‑151 (MQ=255)                                                                                                                                     ggggggCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGc  <  2:605609/151‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                  CCTGTCGGCGGGGATTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGC  >  NC_000913/285511‑285793

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