Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,423,602	T→C	R99R (CGT→CGC)	stfR →	prophage tail fiber N‑terminal domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,423,602	0	T	C	100.0%	280.8 / NA	75	R99R (CGT→CGC)	stfR	prophage tail fiber N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/42);  total (33/42)

ACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGA  >  NZ_CP009273/1423515‑1423687                                                                                        |                                                                                      acaTGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATgcccgccc                                                                                     <  2:1135069/90‑1 (MQ=255) acaTGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATgcccgccc                                                                                     <  1:2084619/90‑1 (MQ=255)       cGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGc                                                                               <  2:1890860/90‑1 (MQ=255)       cGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGc                                                                               >  1:297058/1‑90 (MQ=255)        gggACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCa                                                                              >  1:2699467/1‑90 (MQ=255)        gggACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCa                                                                              >  2:1327053/1‑90 (MQ=255)              aTTACCGTGTATGAAGATTCTCACCCCGGTACGCTGAATGATTTTCTCGGTGCCATGCCGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                        <  1:1854713/90‑1 (MQ=255)              aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                        <  2:1802124/90‑1 (MQ=255)              aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                        <  1:32276/90‑1 (MQ=255)              aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                        <  2:3003930/90‑1 (MQ=255)              aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                        <  1:1240624/90‑1 (MQ=255)              aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                        <  2:656117/90‑1 (MQ=255)               ttACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGcc                                                                       >  1:2733398/1‑90 (MQ=255)                    gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                        <  1:632496/84‑1 (MQ=255)                    gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                        >  2:632496/1‑84 (MQ=255)                    gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGtttt                                                                  >  1:1624643/1‑90 (MQ=255)                    gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGtttt                                                                  >  2:1348541/1‑90 (MQ=255)                    gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGtttt                                                                  >  1:1989496/1‑90 (MQ=255)                    gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGtttt                                                                  >  1:1600644/1‑90 (MQ=255)                    gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGtttt                                                                  >  2:2519354/1‑90 (MQ=255)                    gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGtttt                                                                  >  1:1604569/1‑90 (MQ=255)                      gtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGa                                                                >  2:1219645/1‑90 (MQ=255)                      gtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGa                                                                >  2:1600769/1‑90 (MQ=255)                      gtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGa                                                                >  2:1875788/1‑90 (MQ=255)                          gAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTg                                                            <  2:2376717/90‑1 (MQ=255)                          gAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTg                                                            <  1:2820245/90‑1 (MQ=255)                            aGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAt                                                          <  2:2017274/90‑1 (MQ=255)                            aGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAt                                                          <  1:2519354/90‑1 (MQ=255)                            aGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAt                                                          <  1:2163650/90‑1 (MQ=255)                            aGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAt                                                          <  2:2472647/90‑1 (MQ=255)                                tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtggtg                                                      >  1:2546122/1‑90 (MQ=255)                                tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtggtg                                                      >  1:1473696/1‑90 (MQ=255)                                tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtggtg                                                      >  2:725896/1‑90 (MQ=255)                                tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtggtg                                                      >  1:54732/1‑90 (MQ=255)                                  tcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGa                                                    <  1:1041504/90‑1 (MQ=255)                                    aaCCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTg                                                            >  1:1185430/1‑80 (MQ=255)                                    aaCCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTg                                                            <  2:1185430/80‑1 (MQ=255)                                         ggTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTg                                             <  2:2768062/90‑1 (MQ=255)                                         ggTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTg                                             <  1:213999/90‑1 (MQ=255)                                          gTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAgg                                               >  2:2311433/1‑87 (MQ=255)                                          gTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAgg                                               <  1:2311433/87‑1 (MQ=255)                                               cTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGt                                       >  2:873623/1‑90 (MQ=255)                                               cTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGt                                       >  2:641034/1‑90 (MQ=255)                                               cTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGt                                       >  1:299814/1‑90 (MQ=255)                                               cTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGt                                       >  2:1488631/1‑90 (MQ=255)                                                  aaTGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAAc                                    >  2:1957434/1‑90 (MQ=255)                                                        tttCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTcc                              <  2:926276/90‑1 (MQ=255)                                                        tttCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTcc                              <  2:2460728/90‑1 (MQ=255)                                                        tttCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTcc                              <  2:2217331/90‑1 (MQ=255)                                                        tttCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTcc                              <  1:2992735/90‑1 (MQ=255)                                                          tctcGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCcgc                            <  1:1348541/90‑1 (MQ=255)                                                          tctcGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCcgc                            <  1:1875788/90‑1 (MQ=255)                                                          tctcGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCcgc                            <  2:1473696/90‑1 (MQ=255)                                                           ctcGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCcgcg                           >  2:1885496/1‑90 (MQ=255)                                                             cGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGt                         <  1:50682/90‑1 (MQ=255)                                                             cGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGt                         <  2:2546122/90‑1 (MQ=255)                                                             cGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGt                         <  2:1989496/90‑1 (MQ=255)                                                                tGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGc                      >  1:1314941/1‑90 (MQ=255)                                                                  ccATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGcac                    <  2:444432/90‑1 (MQ=255)                                                                  ccATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGcac                    <  2:1624643/90‑1 (MQ=255)                                                                  ccATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGcac                    <  1:1600769/90‑1 (MQ=255)                                                                  ccATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGcac                    <  2:2699467/90‑1 (MQ=255)                                                                  ccATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGcac                    <  2:2671384/90‑1 (MQ=255)                                                                  ccATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGcac                    <  2:2519916/90‑1 (MQ=255)                                                                  ccATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGcac                    <  2:825324/90‑1 (MQ=255)                                                                  ccATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGcac                    <  1:323628/90‑1 (MQ=255)                                                                  ccATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGcac                    <  2:135437/90‑1 (MQ=255)                                                                         ggaggaTGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAacac             >  1:1535408/1‑90 (MQ=255)                                                                           aggaTGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACgg           >  2:2630322/1‑90 (MQ=255)                                                                            ggaTGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGc          <  1:488880/90‑1 (MQ=255)                                                                            ggaTGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGc          >  2:369427/1‑90 (MQ=255)                                                                                gatgCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTcc                              <  1:818853/66‑1 (MQ=255)                                                                                gatgCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTcc                              >  2:818853/1‑66 (MQ=255)                                                                                gatgCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGcc      >  2:1833075/1‑90 (MQ=255)                                                                                 atgCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGa                 <  2:2965652/78‑1 (MQ=255)                                                                                 atgCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGa                 >  1:2965652/1‑78 (MQ=255)                                                                                 atgCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCcg     <  1:2767451/90‑1 (MQ=255)                                                                                    cccgcccgGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGt                         <  1:1863024/67‑1 (MQ=255)                                                                                    cccgcccgGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGt                         >  2:1863024/1‑67 (MQ=255)                                                                                    cccgcccgGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCga  <  1:369427/90‑1 (MQ=255)                                                                                    cccgcccgGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCga  >  1:2661852/1‑90 (MQ=255)                                                                                        |                                                                                      ACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGA  >  NZ_CP009273/1423515‑1423687

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