按位点提取vcf文件并获取注释的基因

这算是笔记。。。。

  • 前提是总vcf文件已经用snpeff注释过了
vcftools --gzvcf ../chicken.477.chr1Z.indel.zk.tc.zk.zs.vcf.gz --positions specific_position.txt --recode --recode-INFO-all --out sig.indel.vcf

vcftools --gzvcf ../chicken.477.chr1Z.indel.zk.tc.zk.zs.vcf.gz --positions specific_position.txt --recode --recode-INFO-all --stdout |grep -v "#" |cut -f8 |tr "|" "\t" |cut -f5|tr "-" "\n"|sort|uniq >ensebmlID
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容