This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
chmod a+x liftOver |
The next step is to get the required chain files needed to convert from one genome build to another. All UCSC genome builds are listed here and you can select the desired "LiftOver files" link under the genome build that you want to convert from. This will take you to a downloads page with links to the chain files. Chain files are appropriately named so you know what builds you are converting to and from. For example, hg19ToHg18.over.chain.gz is the chain file needed to convert from hg19 to hg18. Once downloaded, unzip the file for use.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
gunzip hg19ToHg18.over.chain.gz |
To run liftOver, the useage is:
liftOver oldFile map.chain newFile unMapped
where:
oldFile is the file you want to convert from
map.chain is the chain file used to convert from one build to another
newFile is the converted file you want to create
unMapped is a file that contains all the unmapped positions
For more details on usage, just type liftOver in the command line.
No comments:
Post a Comment