This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Liftover Gvf From Hg18 To Hg19

Hi,

I'd like to liftover a few gvf files form hg18 to hg19. Do tools exist for this purpose? While the gvf format is supposed to be an extension of gff3, I couldn't get liftover or the gff3ToGenePred tools to accept the files.

Thanks!

liftover hg19
• 6,253 views • link
updated by deanna.church deanna.church • written by hershman • hershman

hey..even i am trying to liftover GVF files from hg19 to hg18...Did u get any solution for that.Is there a GFF to GVF converter otherwise...I tried searching but could not find any as such.Thanks

• 0 views • link
written by gupta.ne • gupta.ne

2 answers

Edit: What kind of errors are you getting?

Edit2: I wrote two perl scripts to make this process go easier. They are not polished code. See below

https://github.com/jewmanchue/zk_tools

grGVF2GFF3 convert your gvf -> gff3

use the UCSC lift over tool

convert back to gvf with LIFTOVERGFF2GVF

I have done this using liftOver from UCSC. Here is the link to the tool:ucsc downloads

GVF should lift over as a GFF3/GFF. If you have any problems let me know.

• 75 views • link
written by Zev.Kronenberg Zev.Kronenberg

I tried using the liftover tool, but it gave me a warning: WARNING: -gff is not recommended. Use 'ldHgGene -out=<file.gp>' and then 'liftOver -genePred <file.gp>'

ldHgGene is no longer supported and the new version of the tool (gff3ToGenePred) gives an error: Input file must conform to the GFF3 specification:

The input file that I'm using came from dgv (ftp://ftp.ebi.ac.uk/pub/databases/dgva/estd59Durbinetal2010/gvf/)

• 0 views • link
written by hershman • hershman

I ignored those warnings. You can spot check that the reference base for any given coordinate matches hg19 on the UCSC genome browser. Everything looked fine in the output.

• 0 views • link
written by Zev.Kronenberg Zev.Kronenberg

I would but the mapped file ended up being just a header and the unmapped file had all of the lines

• 0 views • link
written by hershman • hershman

The NCBI liftover tool natively takes multiple file formats- including gvf and vcf. http://www.ncbi.nlm.nih.gov/genome/tools/remap

• 240 views • link
written by deanna.church deanna.church

Log in to answer this question.

More posts like this
  • How to validate GFF3?
    written by mike.thon • mike.thon
    <p>What tools are available to validate a GFF3 format file? I tried the <a href="http://www.sequenceontology.org/gff3.shtml">one recommended by the Sequence Ontology Project</a> but it has been …
  • BED format - scaffold
    written by 873243 • 873243
    <p>In BED file format, what is the meaning of scaffold instead of chromosome ?</p> <p>Note this format <a href="http://genome.ucsc.edu/FAQ/FAQformat.html#format1">spec</a> .</p>
  • bed files of TopHat2 output
    written by Yizhao • Yizhao
    Hi, I've completed alignment with Tophat2, and got some output files, including deletions.bed, junctions.bed and insertions.bed. My problem is how to tell what info it …
  • Convert .Gff3 File To 12-Column .Bed File
    written by LRStar LRStar
    Hello: I would like to convert .gff3 file to 12-column .bed file, as in this link under "BED Format" (http://genome.ucsc.edu/FAQ/FAQformat.html#format1). I have thus far used …
  • Best Human Genome Reference File For Gatk?
    written by newDNASeqer newDNASeqer
    <p>On GATK website <a href='http://gatkforums.broadinstitute.org/discussion/1204/what-input-files-does-the-gatk-accept'>http://gatkforums.broadinstitute.org/discussion/1204/what-input-files-does-the-gatk-accept</a> and their public FTP server, I found a few difference human genome references in their resource bundle: b37, b36, hg18, …
  • How Are Features With No Strand Encoded In Ucsc Bed Format?
    written by Casey Bergman Casey Bergman
    <p>Hi -</p> <p>I have a set of annotations that have strand information for some of them, but no strand information for others. I'd like to …
  • How To Select Only One Human Genome Build (Hg19) From The Encode Project'S Data
    written by Eric Ho • Eric Ho
    <p>While I looking on the ENCODE project's data in UCSC Genome Browser I discovered both data from hg18 and hg19 are shown even I only …
  • Tool: Bedtools: Analyzing Genomic Features
    written by Istvan Albert Istvan Albert
    All practicing bioinformaticians will face problems that require them to compare, query and select genomic features across an entire genome. As it happens efficient interval …
  • What Does Bed Stand For?
    written by Daniel Standage Daniel Standage
    <p>GFF is for Generic Feature Format, GTF is for Gene transfer format. Is BED (<a href='http://genome.ucsc.edu/FAQ/FAQformat.html#format1'>a la UCSC</a>) also an acronym?</p>
  • Help With Maf'S In Bx-Python
    written by Noah • Noah
    <p>I'm looking to extract some regions from some <a href="http://genome.ucsc.edu/FAQ/FAQformat.html#format5">maf</a> files, and it seems the way to do it in python would be using bx-python. …