GO注释步骤-整合eggnog和interpro的GO结果

# 7. GO 注释
mkdir -p /home/train/13.functional_annotation/07.GO
cd /home/train/13.functional_annotation/07.GO

# 整合eggNOG和InterPro中的GO注释结果
go_from_eggNOG_and_interpro.pl ../04.eggNOG/eggNOG.emapper.annotations ../05.InterPro/interpro.tsv > go.annot
go_reducing_go_number_para.pl /opt/biosoft/go_class/bin/go-basic.obo go.annot 8 > go_reduced.annot
sort go_reduced.annot > go.annot; rm go_reduced.annot
gene_annotation_from_table.pl go.annot > GO.txt

# GO wego分类图（按照goslim_agr.obo方法分了53类）
/opt/biosoft/go_class/bin/annot2wego.pl go.annot > go.wego
get_Genes_From_GO.pl /opt/biosoft/go_class/bin/go-basic.obo go.wego > go_class.tab
go_svg.pl --outdir ./ --name out --color "green" --mark "Whole Genome Genes" --note "GO Class of whole genome genes of Malassezia sympodialis" go.wego
perl -p -i -e 's/MovePer:0.125/MovePer:0.5/; s/FontSize:.*/FontSize:24/;' out.lst
/opt/biosoft/go_class/svg/distributing_svg.pl out.lst out.svg
/opt/biosoft/go_class/bin/changsvgsize.pl out.svg 150 -100
convert out.svg out.png

# GO分类（按照goslim_agr.obo方法分了53类）
go_class.pl /opt/biosoft/go_class/bin/go-basic.obo go.annot --go_class_method /opt/biosoft/go_class/bin/goslim_agr.obo > go_class.tab

cp go.annot ../functional_annotation.GO.tab
cp GO.txt ../functional_annotation.GO.txt
cp out.svg ../functional_annotation.GO_class.svg
cp out.png ../functional_annotation.GO_class.png
cp go_class.tab ../functional_annotation.GO_class.tab
cd ..