2022-04-21fa文件和gff文件提取cds序列并翻译成蛋白

cat *.fa >hg19.fa

#根据hg19.fa提取CDS序列

gffread gencode.v19.annotation.gff3 -g hg19.fa -y tr_cds.fa

#翻译后蛋白序列

gffread gencode.v19.annotation.gff3 -g hg19.fa -x cds.fa

#获得外显子序列

gffread gencode.v19.annotation.gff3 -g hg19.fa -w exons.fa

#格式转换

gffread gencode.v19.annotation.gff3 -T -o gencode.v19.gtf

gffread merged.gtf -o- > merged.gff3

用于记录

详见其他人简书：

https://www.jianshu.com/p/c3f723c895fe