根据基因组fa文件和gff文件提取cds并翻译成pep

cds文件:注释信息里的蛋白编码序列
pep文件:cds对应序列翻译成的氨基酸序列
从NCBI上下载基因组时,有的并没有上传cds文件和pep文件,此时该怎么办呢?


(1)利用脚本

①根据注释文件提取转录本:
生信笔记系列之序列提取--根据GTF提取转录本
从NCBI基因组数据中获得cds,pep和geneID对应表- 薛猫_柳叶 ...
②将cds转换成pep:
从cds到pep

(2)利用cufflinks中的gffread工具

GTF/GFF格式gffread入门使用

# conda上安装cufflinks,使用之前激活环境
source /data1/spider/liupiao/miniconda3/bin/activate
# 提取cds
gffread in.gff3 -g ref.fa -x cds.fa
# 获得pep
gffread in.gff3 -g ref.fa -y pep.fa

你可能感兴趣的:(根据基因组fa文件和gff文件提取cds并翻译成pep)