From Hg18 To Grch38
2
3
Entering edit mode
10.8 years ago
aleka ▴ 110

Does anyone know if it is possible to convert SNP coordinated from Hapmap database to the new reference genome. UCSC doesn't have the liftover yet ready. Any suggestions?

• 8.6k views
ADD COMMENT
1
Entering edit mode
ADD REPLY
1
Entering edit mode
ADD REPLY
6
Entering edit mode
ADD COMMENT
1
Entering edit mode

NCBI's remap is a great alternative to liftOver, and a standalone script can be found here. It doesn't accept HapMap's old format, but will take VCFs. At this FTP link, Ensembl lists HapMap variants in VCFs mapped to GRCh37.p13. And remap-ing GRCh37.p13 to GRCh38 should cause fewer conflicts than hg18 to GRCh38. In the remap tool, be sure to set the Source Assembly to GRCh37.p13 instead of GRCh37.

ADD REPLY
1
Entering edit mode
10.7 years ago

An hg18ToHg38 chain file is now available for liftOver. But since you're interested in HapMap variants, it is best to start with the GRCh37.p13 VCFs from Ensembl. And then liftOver those VCFs to GRCh38 using the hg19ToHg38 chain file. If you have trouble with UCSC's different chromosome names, use NCBI's remap instead, and map GRCh37.p13 to GRCh38.

ADD COMMENT

Login before adding your answer.

Traffic: 1658 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