$$Score_i = Phage\_score_i \cdot GC\_cont\_score_i$$ $$Phage\_score_i = Tr(\mathbf{1}_{pVOG}(gene_i))$$ $$GC\_cont\_score_i = Tr(GC(gene_i)\cdot \mathbf{1}_{pVOG}(gene_i))$$