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.
0 answers
No answers yet.
Log in to answer this question.
-
Converting GRCh38 30× 1000G data to GRCh37
written by devenvyas
-
Using a Docker image with ensembl-vep
written by langziv •
-
CrossMap VCF liftover (GRCh37 > GRCh38): END info field for indel calls not lifted over
written by mrodrigues.fernanda •
-
output path issue on Busco (python)
written by Darrill •
-
Failed to liftover from hg38 to hg19
written by BAGeno
-
Liftover errors with CrossMap and Picard
written by rshoobs •
-
How to change the path so that screen can take the default python path?
written by saadleeshehreen
-
CrossMap Error With BigWig Files
written by skasowitz •
-
command line python script run on a file in different directory
written by ahmedakhokhar
-
script to run CrossMap for a set of bam files
written by ifudontmind_plzz
