Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	287,351	C→T	A188V (GCG→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,351	0	C	T	100.0%	159.2 / NA	46	A188V (GCG→GTG)	yagH	CP4‑6 prophage; putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/25);  total (21/25)

CCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGAT  >  NC_000913/287214‑287486                                                                                                                                          |                                                                                                                                        ccGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGgcgtgc                                                                                                                                       >  2:786748/1‑140 (MQ=255)   gCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCAc                                                                                                                                     >  2:821930/1‑140 (MQ=255)   gCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCAc                                                                                                                                     >  1:931685/1‑140 (MQ=255)     ccGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCt                                                                                                                                   >  2:780136/1‑140 (MQ=255)     ccGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCt                                                                                                                                   >  1:89308/1‑140 (MQ=255)               cgcgCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCAc                                                                                                                         >  1:526657/1‑140 (MQ=255)               cgcgCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCAc                                                                                                                         >  2:873412/1‑140 (MQ=255)               cgcgCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCAc                                                                                                                         >  2:960881/1‑140 (MQ=255)                gcgcCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCAcg                                                                                                                        <  1:821930/140‑1 (MQ=255)                    caccacAGCAACCCGCACAACCCCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCggg                                                                                                                    <  2:949432/140‑1 (MQ=255)                    caccacAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCggg                                                                                                                    <  1:477928/140‑1 (MQ=255)                    caccacAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCggg                                                                                                                    <  1:616275/140‑1 (MQ=255)                    caccacAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCggg                                                                                                                    <  1:73523/140‑1 (MQ=255)                     accacAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGa                                                                                                                   <  2:53393/140‑1 (MQ=255)                     accacAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGa                                                                                                                   <  2:163049/140‑1 (MQ=255)                     accacAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGa                                                                                                                   <  1:960881/140‑1 (MQ=255)                       cacaGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCAcg                                                                                                                        <  1:565176/133‑1 (MQ=255)                       cacaGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCAcg                                                                                                                        >  2:565176/1‑133 (MQ=255)                                 cGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTg                                                                                                       >  1:419151/1‑140 (MQ=255)                                 cGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTg                                                                                                       >  1:208136/1‑140 (MQ=255)                                           cATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGccgccg                                                                                             <  2:208136/140‑1 (MQ=255)                                           cATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGccgccg                                                                                             <  1:873412/140‑1 (MQ=255)                                             tCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAg                                                                                           >  2:868674/1‑140 (MQ=255)                                             tCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAg                                                                                           >  2:796376/1‑140 (MQ=255)                                                tgtTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGggc                                                                                        >  1:936995/1‑140 (MQ=255)                                                tgtTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGggc                                                                                        <  1:780136/140‑1 (MQ=255)                                                  ttACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGggcgg                                                                                      >  1:921670/1‑140 (MQ=255)                                                        gCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAg                                                                                >  2:383389/1‑140 (MQ=255)                                                          gTTTGACCCGCAGGCCGGCACGCTGTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCGGAGGGCGGCACCAGCt                                                                              <  2:175701/140‑1 (MQ=255)                                                          gTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCt                                                                              >  1:456706/1‑140 (MQ=255)                                                            ttGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTa                                                                             >  1:892983/1‑139 (MQ=255)                                                            ttGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTa                                                                             <  2:892983/139‑1 (MQ=255)                                                             tGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACg                                                                           >  1:805307/1‑140 (MQ=255)                                                                 ccGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCa                                                                       <  1:954363/140‑1 (MQ=255)                                                                  cGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCAc                                                                      <  2:89308/140‑1 (MQ=255)                                                                    cAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGc                                                                    <  2:805307/140‑1 (MQ=255)                                                                      gACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCcg                                                                  <  2:678861/140‑1 (MQ=255)                                                                         cGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCcgtcg                                                               <  1:796376/140‑1 (MQ=255)                                                                                    cGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGTGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGt                                                    <  2:931685/140‑1 (MQ=255)                                                                                                 aaCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCg                                       >  2:824394/1‑140 (MQ=255)                                                                                                  aCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGa                                      <  2:598348/140‑1 (MQ=255)                                                                                                  aCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGa                                      <  1:868674/140‑1 (MQ=255)                                                                                                               gCACGCCGCTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAg                         <  2:456706/140‑1 (MQ=255)                                                                                                                  cgccgcTCTGCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCCGCGTTCCAAAAATATCGACGGGCCGTACGAGCTg                      >  2:842596/1‑140 (MQ=255)                                                                                                                           gCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGAc             <  1:824394/140‑1 (MQ=255)                                                                                                                           gCTACACCGAAGGCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGAc             <  2:526657/140‑1 (MQ=255)                                                                                                                                      ggCGTGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgatgat  >  2:612008/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        CCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGAT  >  NC_000913/287214‑287486

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