Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	287,977	A→G	M397V (ATG→GTG)	yagH →	CP4‑6 prophage, putative xylosidase/arabinosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	287,977	0	A	G	100.0%	237.5 / NA	66	M397V (ATG→GTG)	yagH	CP4‑6 prophage, putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/35);  total (31/35)

CGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAG  >  NC_000913/287841‑288083                                                                                                                                         |                                                                                                           cgacgaCACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTg                                                                                                          <  1:1400443/139‑1 (MQ=255)       cacCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTtc                                                                                                    >  2:1023442/1‑139 (MQ=255)       cacCCTCGGCTCGCTAACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTAGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTtc                                                                                                    >  1:1225794/1‑139 (MQ=255)              ggCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCgg                                                                                             >  2:1284097/1‑139 (MQ=255)                 tcgctcACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGc                                                                                          >  1:1331219/1‑139 (MQ=255)                    ctcACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACt                                                                                       >  1:139842/1‑139 (MQ=255)                         cgcgcgCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCcag                                                                                  >  1:949154/1‑139 (MQ=255)                         cgcgcgCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGCCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCcag                                                                                  <  2:1671727/139‑1 (MQ=255)                             cgcCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCCGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCCCGCCGGTGCACTTCCAGCaga                                                                              >  1:2236/1‑139 (MQ=255)                                ccGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                           <  1:1175065/139‑1 (MQ=255)                                ccGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                           <  1:1334820/139‑1 (MQ=255)                                     cttcttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCgggg                                                                      <  2:139842/139‑1 (MQ=255)                                       tcttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCt                                                                    >  1:685065/1‑139 (MQ=255)                                         ttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                  <  2:164544/139‑1 (MQ=255)                                           aCGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGAcc                                                                <  2:1435245/139‑1 (MQ=255)                                            cGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                               <  1:1806481/139‑1 (MQ=255)                                              gCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGTACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGTACTTCCAGCAGAGCGCGGGGCTGACCTGc                                                             <  2:1225794/139‑1 (MQ=255)                                              gCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGc                                                             >  2:1649611/1‑139 (MQ=255)                                                  tATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctact                                                         >  1:620415/1‑139 (MQ=255)                                                  tATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctact                                                         >  1:533090/1‑139 (MQ=255)                                                  tATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctact                                                         >  2:157764/1‑139 (MQ=255)                                                  tATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctact                                                         >  2:593391/1‑139 (MQ=255)                                                    tGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTAGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCtactac                                                       >  1:751236/1‑139 (MQ=255)                                                     ggCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTaca                                                      <  1:593391/139‑1 (MQ=255)                                                          cGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                               <  2:580181/125‑1 (MQ=255)                                                          cGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                               >  1:580181/1‑125 (MQ=255)                                                          cGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGc                                                 <  2:620415/139‑1 (MQ=255)                                                            aCTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCaa                                               >  1:1524054/1‑139 (MQ=255)                                                                gCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCgggg                                                                      >  1:1275336/1‑112 (MQ=255)                                                                gCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCgggg                                                                      <  2:1275336/112‑1 (MQ=255)                                                                   cAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTgg                                        <  1:676038/139‑1 (MQ=255)                                                                   cAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTgg                                        <  1:1023442/139‑1 (MQ=255)                                                                     aTTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAg                                      <  2:1524054/139‑1 (MQ=255)                                                                            cTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctac                                                          <  1:1669929/112‑1 (MQ=255)                                                                            cTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctac                                                          >  2:1669929/1‑112 (MQ=255)                                                                              tCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCtt                             >  2:26971/1‑139 (MQ=255)                                                                               cACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCttt                            <  2:949154/139‑1 (MQ=255)                                                                               cACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCttt                            >  2:1716165/1‑139 (MQ=255)                                                                                aCCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGCGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                           >  1:273028/1‑91 (MQ=255)                                                                                aCCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                           <  2:273028/91‑1 (MQ=255)                                                                                aCCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                  >  2:1217906/1‑100 (MQ=255)                                                                                aCCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                  <  1:1217906/100‑1 (MQ=255)                                                                                       cGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTa                    <  1:157764/139‑1 (MQ=255)                                                                                         aCCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                           <  2:296738/82‑1 (MQ=255)                                                                                         aCCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                           >  1:296738/1‑82 (MQ=255)                                                                                         aCCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACg                  <  1:1284097/139‑1 (MQ=255)                                                                                         aCCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACg                  >  1:165379/1‑139 (MQ=255)                                                                                           cGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgag                <  2:533090/139‑1 (MQ=255)                                                                                            gTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGc                                                                                          >  1:545293/1‑64 (MQ=255)                                                                                            gTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGc                                                                                          <  2:545293/64‑1 (MQ=255)                                                                                            gTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCttt                            <  2:837729/126‑1 (MQ=255)                                                                                            gTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCttt                            >  1:837729/1‑126 (MQ=255)                                                                                                 gcgcCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGa          <  1:26971/139‑1 (MQ=255)                                                                                                 gcgcCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGa          <  2:685065/139‑1 (MQ=255)                                                                                                       cTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCcagca                                                                                <  2:933359/63‑1 (MQ=255)                                                                                                       cTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCcagca                                                                                >  1:933359/1‑63 (MQ=255)                                                                                                       cTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                  <  2:395845/77‑1 (MQ=255)                                                                                                       cTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                  >  1:395845/1‑77 (MQ=255)                                                                                                         ggCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAg  <  2:748722/139‑1 (MQ=255)                                                                                                                   cGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGc                                                                                          >  1:669457/1‑41 (MQ=255)                                                                                                                   cGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGc                                                                                          <  2:669457/41‑1 (MQ=255)                                                                                                                         ccGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGa          <  2:517595/115‑1 (MQ=255)                                                                                                                         ccGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGa          >  1:517595/1‑115 (MQ=255)                                                                                                                                 aaaCGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGAcc                                                                >  2:1534217/1‑53 (MQ=255)                                                                                                                                 aaaCGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGAcc                                                                <  1:1534217/53‑1 (MQ=255)                                                                                                                                      cgGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                           <  1:1166115/37‑1 (MQ=255)                                                                                                                                         |                                                                                                           CGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAG  >  NC_000913/287841‑288083

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