Converting 1kHG genomic coordinates from hg37 to hg38 build
0
0
Entering edit mode
3.1 years ago
AVA ▴ 40

I am trying to change genomic coordinates for 1kHG data from hg37 to hg38 using CrossMap. I have Crossmap.py in the folder CrossMap at path /fs/projects/..../CrossMap/bin/CrossMap.py. I have never run a python script before. I understood that I am supposed to set python3 path to run this script. I used the following command to set python path in a folder one path above CrossMap using:

export PATH="$PATH:/usr/local/bin/python3" and then run the following command to change the coordinates:

path/CrossMap.py vcf <chain_file> <input.vcf> <refgenome.fa> <output_file> The command seems to run as I do not get any warning or error but I don't get any log or output for this. Any suggestions or ideas are highly appreciated.

python3 CrossMap 1kHG hg37 to hg38 • 806 views
ADD COMMENT

Login before adding your answer.

Traffic: 3340 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6