Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,776,304	A→C	E129A (GAG→GCG)	yfjX →	CP4‑57 prophage; putative antirestriction protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,776,304	0	A	C	100.0%	192.5 / NA	56	E129A (GAG→GCG)	yfjX	CP4‑57 prophage; putative antirestriction protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (35/21);  total (35/21)

GACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAA  >  NC_000913/2776165‑2776450                                                                                                                                            |                                                                                                                                                   gACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTTcc                                                                                                                                         <  2:758218/151‑1 (MQ=255)   cGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGc                                                                                                                                       >  2:129627/1‑151 (MQ=255)           aTAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCggg                                                                                                                                >  2:774688/1‑150 (MQ=255)           aTAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCggg                                                                                                                                <  1:574132/150‑1 (MQ=255)             aaaTGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGAt                                                                                                                              <  1:11025/150‑1 (MQ=255)                   cATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGctc                                                                                                                       >  1:718204/1‑151 (MQ=255)                      ttGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTgc                                                                                                                    >  2:232857/1‑151 (MQ=255)                               ggCATGAATGGCAGTGGTGCGGAATTGACCGCGGAGGCCTCAGGTATCGCTGTCTGCCTGATTGAATATGGCCACCACGCCTGCCTCCCCGAATGTGTTGCCATGGCGGCGCCCTATGACCGCCTGCGGGATTACGCTCTGCAGCACCCTg                                                                                                           <  1:548999/151‑1 (MQ=255)                                  aTGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAAt                                                                                                        <  2:745134/151‑1 (MQ=255)                                    gAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATcc                                                                                                      >  1:775933/1‑151 (MQ=255)                                    gAATGGCAATGGTGCGGAAATGACCGCGGAAGCCGCAGGTATCGCTGGCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTACAGCACCCt                                                                                                            <  1:523080/145‑1 (MQ=255)                                     aaTGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTTTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATcc                                                                                                      <  1:753019/150‑1 (MQ=255)                                     aaTGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATcc                                                                                                      >  2:860617/1‑150 (MQ=255)                                      aTGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAg                                                                                                    <  1:521281/151‑1 (MQ=255)                                        ggCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAgcg                                                                                                  >  2:645302/1‑151 (MQ=255)                                            aTGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCAt                                                                                              >  2:578045/1‑151 (MQ=255)                                              ggTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCAtt                                                                                             >  1:548941/1‑150 (MQ=255)                                               gTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTAt                                                                                           >  1:682764/1‑151 (MQ=255)                                                 gCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATgc                                                                                         >  2:289539/1‑151 (MQ=255)                                                   ggAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATgcgc                                                                                       >  2:671558/1‑151 (MQ=255)                                                     aaaTGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCat                                                                                     >  2:227528/1‑151 (MQ=255)                                                      aaTGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCatt                                                                                    >  2:637601/1‑151 (MQ=255)                                                       aTGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCatt                                                                                    >  1:373060/1‑150 (MQ=255)                                                         gAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCattatt                                                                                 <  1:285776/151‑1 (MQ=255)                                                             gcgGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGAt                                                                              >  1:906976/1‑150 (MQ=255)                                                                  aGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAgg                                                                         <  2:704229/150‑1 (MQ=255)                                                                      gCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACt                                                                    >  1:370267/1‑151 (MQ=255)                                                                        aGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTc                                                                   >  1:755304/1‑150 (MQ=255)                                                                         ggTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCat                                                                 >  1:697710/1‑151 (MQ=255)                                                                          gTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCatg                                                                >  2:811066/1‑151 (MQ=255)                                                                          gTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCatg                                                                >  2:823096/1‑151 (MQ=255)                                                                          gTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCatg                                                                <  2:84589/151‑1 (MQ=255)                                                                            aTCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCatgat                                                              >  2:305153/1‑151 (MQ=255)                                                                              cGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATgg                                                            >  1:38003/1‑151 (MQ=255)                                                                                cTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCatgatg                                                             <  1:470207/148‑1 (MQ=255)                                                                                cTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGa                                                           <  2:183473/150‑1 (MQ=255)                                                                                   tCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACa                                                        <  1:500152/150‑1 (MQ=255)                                                                                     tGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTc                                                     >  2:380132/1‑151 (MQ=255)                                                                                               aaTATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATccc                                            >  1:156878/1‑150 (MQ=255)                                                                                                      ccATCACGCCTGCCTCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGacac                                    <  2:314372/151‑1 (MQ=255)                                                                                                              cctgcctCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACtt                            >  2:16655/1‑151 (MQ=255)                                                                                                                  cctCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCAc                        <  1:623393/151‑1 (MQ=255)                                                                                                                   ctCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCCCTGCAGCACCCTGAATACAGCGCAATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACt                       >  1:498951/1‑151 (MQ=255)                                                                                                                    tCACCGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTg                      >  1:270589/1‑151 (MQ=255)                                                                                                                        cGAATGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCg                  >  1:836/1‑151 (MQ=255)                                                                                                                            tGCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGcac               >  2:490509/1‑150 (MQ=255)                                                                                                                             gCGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGcaca              <  1:97544/150‑1 (MQ=255)                                                                                                                              cGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAAc            >  2:732850/1‑151 (MQ=255)                                                                                                                              cGATGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAAAAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAAc            <  2:6918/151‑1 (MQ=255)                                                                                                                                 tGCCATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGc          >  2:83855/1‑150 (MQ=255)                                                                                                                                   ccATGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACg       <  1:891531/151‑1 (MQ=255)                                                                                                                                      tGACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGt     >  2:700933/1‑150 (MQ=255)                                                                                                                                       gACGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTa    <  1:904220/150‑1 (MQ=255)                                                                                                                                        aCGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTaa   >  1:804002/1‑150 (MQ=255)                                                                                                                                         cGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCCCTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTaaa  <  2:370281/150‑1 (MQ=255)                                                                                                                                         cGGCGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTaaa  >  2:907959/1‑150 (MQ=255)                                                                                                                                            |                                                                                                                                                   GACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAA  >  NC_000913/2776165‑2776450

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