liftover 转换不同版本的GRCh

1.下载liftover : wget http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/liftOver

2.下载坐标注释文件：Index of /goldenPath/hg38/liftOver

3.input文件，bed格式：只定义前三列：chr start end，无表头。start和end不同，若是相同，end可以取start+1

4.坐标转换：chmod +x ./liftover #将liftover变成可执行文件
./liftover inputfile over.chain.gz outputfile unmapfile（对应未对应上的行）

inputfile:输入文件 bed格式

over.chain.gz:坐标注释文件（无需解压）

outputfile：输出文件

unmapfile:输出为对应上的行

#常见报错：Expecting integer field 3 line 225125 of chr19.bed, got 2.9e+07

原因：出现此问题的原因是第三列的bp位置出现了科学计数法的形式，因此无法识别。

解决办法：在利用R对原bp位置执行加1操作时，先设置不使用科学计数法：options(scipen = 100)

代码记录：

options(scipen = 100)
library(dplyr)
df = read.table('AD.txt',sep='\t',header = T)
df = mutate(df,end=BP+1)#增加一列
options(max.print = 30000000)
df = df %>% select(CHR, BP, end, SNP)#改换列的位置
write.table(df,file = 'AD.bed',sep = '\t',row.names = F,col.names = F,quote = F)

df1=df %>% filter(df[,1] %in% "chr1") #提取df中第一列为chr1的列

liftover 转换不同版本的GRCh相关推荐

java pdf版本转换_Java 转换 PDF 版本
在日常工作或生活中,我们使用的某些设备可能对PDF版本有严格要求,此时就需要更改PDF版本来适应设备.本文就将展示如何通过Free Spire.PDF for Java来转换PDF版本.(支持转换PD ...
CAD文件低版本怎么转换高版本
我们在工作当中需要CAD格式的素材,但是在网上下载的CAD文件版本过低不兼容,使得这些素材我们不好使用,这时候我们需要将它转换成高一点的版本,那么CAD文件低版本怎么转换高版本呢? 这里小编就用在线转 ...
cad批量转换低版本如何实现？
在CAD制图工作中,我们可能会遇到一些比较麻烦的问题.其中就有CAD图纸由于版本过高导致无法打开查看的问题,而且有时候CAD图纸过多,若是每一张单独查看也会很麻烦.这时候该如何解决呢?cad批量转换 ...
使用babel转换高版本的js 更好的兼容一些设备
前段时间一直在做移动端的东西我很多都是使用es6的语法最后在实际环境测试的时候发现人家的设备根本就不识别我的那些语法这个时候你就有俩个选择要不转换成低版本要不重新使用低版本js重新书写 ...
如何使用cad版本转换器？转换cad版本的技巧
如何使用cad版本转换器?经常使用cad的小伙伴就知道cad制图软件导出文件分为很多种,其中较为常见的有dwg和dxf两种.然而即使是cad文件同为dwg或dxf格式也被分为多个版本,在遇到cad版本 ...
怎么用cad版本转换器？转换cad版本的方法
怎么用cad版本转换器?cad格式是常见的一种制图工具,然而cad工具也在不断的更新,当你获得cad素材时经常会因为版本不一致而导致文件出现打不开的情况.在这类情况出现时我们通常可以用cad版本转换器 ...
java版 pdf转换_Java 转换 PDF 版本
在日常工作或生活中,我们使用的某些设备可能对PDF版本有严格要求,此时就需要更改PDF版本来适应设备.本文就将展示如何通过Free Spire.PDF for Java来转换PDF版本.(支持转换PD ...
cad在线转换低版本_别再傻瓜式操作了，工作效率上不去？这6个小技巧带你玩转CAD...
CAD这个行业现在发展的越来越好了,许多小伙伴们纷纷加入了这个行业中.对于刚接触的小伙伴们来说,学习CAD的时候是一个比较枯燥而且乏味的过程,很多人在学习CAD的时候都没能坚持到最后,所以今天给大家分 ...
教您MAYA 高版本文件转换低版本方法
近日不停的有会员说为什么下载的模型打不开!有的因为模型错误!有的因为版本不对! 请大家善于看论坛的新手帮助!而不是一味的傻问! 现在将网上搜集的方法公布给大家!(含插件) 方法一: 这个版本转换的问题 ...

liftover 转换不同版本的GRCh

liftover 转换不同版本的GRCh相关推荐

最新文章

热门文章