Skip to content
Hiroyuki MISHIMA edited this page Sep 17, 2015 · 5 revisions

Liftover HGVD allele frequency DB

Protocol

  1. HGVD was converted into a VCF file. See https://github.com/misshie/hgvd2annovar
  2. Use CrossMap at http://crossmap.sourceforge.net/ to perform automatic liftover
  • Using CrossMap, 1,006 loci in 478,228 loci were unmapped.
  1. Validate result in-house scripts using Ruby UCSC API at https://github.com/misshie/bioruby-ucsc-api .
  2. In SNV sites (total 468,025), REF sequences in 5 sites are discordant to hg38.
  3. In unmapped 1,006 loci, 639 loci has dbSNP142-registered with single positions. Thus, 367 loci may need BLAT remapping.