脚本路径:
/TJPROJ6/RNA_SH/personal_dir/yangdan/new/gongdan/TF/TF.py
动植物转录因子数据库路径:
/TJPROJ6/RNA_SH/personal_dir/yangdan/new/gongdan/TF/AnimalTFDB3 /TJPROJ6/RNA_SH/personal_dir/yangdan/new/gongdan/TF/PlantTFDB
参数设置:
# 需要比对的fa parser.add_argument('-F', '--fa', help='input fasta file', required=True) #可选折植物或动物转录因子注释数据库 parser.add_argument('-T', '--type', help='type of speceis Animal or Plant', required=True) # 输出结果路径 parser.add_argument('-O', '--outdir', help='output dir', required=True)
gene_id | species | tf_family | protein_family | score | evalue |
---|---|---|---|---|---|
AVESA.00001b.r3.5Dg0002536 | Phyllostachys heterocycla | C3H | C3H family protein | 67 | 4.9e-69 |
AVESA.00001b.r3.6Ag0003127 | Aegilops tauschii | WRKY | WRKY family protein | 673 | 5.8e-109 |
AVESA.00001b.r3.4Cg0000978 | Zoysia pacific | TALE | TALE family protein | 326 | 3.0e-29 |
AVESA.00001b.r3.4Dg0003041 | Oryza punctata | NF-YB | NF-YB family protein | 1013 | 5.8E-12 |