$$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_{deviation}(gene_i)\cdot \mathbf{1}_{pVOG}(gene_i))$$ $$GC_{deviation}(gene_i) = \left|mean\_gc-GC(gene_i)\right|%2Bmean\_gc$$