LINUX | liftOver:基因組坐標轉換hg38到hg19

1 準備bed文件(R)

rawID <- data.frame(matrix(unlist(rawID),ncol = 3,byrow = T))
colnames(rawID) <- c("chr","start","end")
write.table(rawID,"input.bed",col.names = F,row.names = F,quote = F, sep = "\t")

2 下載坐標轉換對應文件,hg38TOhg19

wget http://hgdownload.soe.ucsc.edu/goldenPath/hg38/liftOver/hg38ToHg19.over.chain.gz
gunzip hg38ToHg19.over.chain.gz

3 下載轉換工具

wget -c http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/liftOver
chmod 755 liftOver
./liftOver

4 進行轉換

./liftOver input.bed hg38ToHg19.over.chain output_hg19.bed unMapped

5 查看結果

#hg38
cat input.bed | head
chr1    1223243 1223968
chr1    1804418 1806538
chr1    9931890 9934860
chr1    10433841    10442441
chr1    15534236    15536814
chr1    15638306    15643650
chr1    15717892    15721388
chr1    20749722    20773610
chr1    20757165    20773610

#hg19
cat output_hg19.bed | head
chr1    1158623 1159348
chr1    1735857 1737977
chr1    9991948 9994918
chr1    10493898    10502498
chr1    15860731    15863309
chr1    15964801    15970145
chr1    16044387    16047883
chr1    21076215    21100103
chr1    21083658    21100103
?著作權歸作者所有,轉載或內容合作請聯(lián)系作者
【社區(qū)內容提示】社區(qū)部分內容疑似由AI輔助生成,瀏覽時請結合常識與多方信息審慎甄別。
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發(fā)布,文章內容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務。

相關閱讀更多精彩內容

友情鏈接更多精彩內容