I have a problem with liftOver. I have 300,000 snps to convert from hg18 to hg19,so I use liftOver. I download it in my own computer and run it.I also download the chain which is hg18 to hg19. I found that it can handle most of snps, but will fail about 10,000 snps. I search the failed snps in NCBI, and found they are actually in NCBI, and have position of hg18 and hg19. How can I do do convert the remainder snps? I want to handle them in batches.
1 answer
If the failed SNPs all have rs IDS, you should be able to retrieve their current position from Ensembl using either BioMart or the Perl API, provided Ensembl has been able to map them to GRCh37 (a.k.a. hg19).
In that case, I would start with BioMart. Please have a look at the BioMart videos on YouTube to get an idea how to work with it: http://www.youtube.com/user/EnsemblHelpdesk#p/u/4/DXPaBdPM2vs and http://www.youtube.com/user/EnsemblHelpdesk#p/u/6/paC3sOANSJA. I hope this helps. Questions regarding Ensembl BioMart can be directed to the Ensembl Helpdesk (helpdesk@ensembl.org).
•
written
by
Bert Overduin
Log in to answer this question.
-
Converting genome build from hg38 to hg19 using liftOver via bigsnpr - error fetching file
written by dam4l
-
How to change "CompressedGRangesList" to "GRangesList"
written by Aki •
-
Lifting over results in coordinates that do not agree with what the rsid suggests the coordinates s…
written by curious
-
rs number conversion when performing liftover from hg18tohg19
written by giulia.pagani •
-
SNP coordinate systems in UCSC LiftOver (hg19), NCBI Remap (GRCh37(hg19)) and dbSNP (GRCh38(hg38))
written by curious •
-
Reversing direction of chain files (E.g. from hg18-->hg19 to hg19-->hg18)
written by dylkot •
-
LiftOver in R (error)
written by viniciushs.z •
-
How Do I Convert Hg19 Gvf To Hg18 Gvf Using Ucsc Liftover Tool
written by gupta.ne •
-
Convert Vcf From Hg18 To Hg19 Using Liftovervcf.Pl
written by Raony Guimaraes •
-
Does Liftover Have Single-Base Resolution?
written by Ryan Thompson
