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
网友评论