bed vcf 文件 GRCH37转 GRCH38 hg19 转hg38

准备文件:

hg19ToHg38.over.chain.gz
hg38ToHg19.over.chain.gz

GRCh37_to_GRCh38.chain.gz
GRCh37_to_GRCh38.chain.gz

下载网页打开很慢,已经放到百度云盘

链接:https://pan.baidu.com/s/1CGQ4TtaibfiMPlkFXSD4yQ
提取码:oqkd

1. bed 文件 hg19转hg38 (GRCH同理)

工具:liftover

liftover只能处理bed格式

liftOver input_hg19.bed  hg19ToHg38.over.chain.gz output_hg38.bed unmap.bed

2. vcf文件 GRCH37转GRCH38

工具:vcf-liftover

本质上是调用liftover,速度很快

更正:这个出来的结果并不是标准的vcf格式,作者的代码有问题,后续处理这个vcf文件的时候会有错误

github链接:https://github.com/liqg/vcf-liftover

zcat input.vcf.gz | ~/vcf-liftover GRCh37_to_GRCh38.chain.gz | bgzip -c > 
  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值