zoukankan      html  css  js  c++  java
  • How to use liftover

    First step, installation:

    download from UCSC genome browser

    chmod + x 

    Second, very important: to download the "Chain File". 

    What is "chain file"?  --- to transfer the gene data from one build to another . for example, from hg19 to hg18

    Third, do like this:

    liftOver oldfile chain_file new_file unmapped_file

    for example:

    liftOver SAEC_DNase.bed hg19ToHg18.over.chain SAEC_DNase_hg18.bed SAEC_DNase_hg18_Unmapped.bed

    or you can do the batch work by shell script:

    #!/bin/bash
    for file in *.bed; do
        newfile=${file/%.bed/_hg18.bed}
        unmapp_file=${file/%.bed/_Unmapped.bed}
        liftOver $file hg19ToHg18.over.chain Hg18_Files/$newfile Hg18_Files/$unmapp_file
    done
  • 相关阅读:
    HDU 5247
    HDU 4965
    CodeForces 445B
    HDU 5835
    CodeForces 731C
    HDU 5783
    CodeForces 660D
    POJ 1631
    HDU 6112
    HDU 5860
  • 原文地址:https://www.cnblogs.com/foreverycc/p/3170807.html
Copyright © 2011-2022 走看看