Question: Extracting VCF information to use in LiftOver
0
gravatar for cookersjs
3.3 years ago by
cookersjs20
Toronto, ON
cookersjs20 wrote:

Good day,

I am interested in using UCSC's LiftOver tool to convert from GRCh37 to 38. Since VCF's only have the single position in the file, is there any reason I couldn't use the position as both the start and end values?

The reason i'm interested in LiftOver particularly, is that they have the file that can be used from the command-line and it's a simple executable.

I briefly looked at NCBI's Remap (no command-line tool?), and Crossmap (don't want to have python as a dependency).

Any suggestions are welcome!

Thanks

liftover vcf • 1.3k views
ADD COMMENTlink modified 3.3 years ago by Pierre Lindenbaum131k • written 3.3 years ago by cookersjs20
3
gravatar for Pierre Lindenbaum
3.3 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum131k wrote:

Since VCF's only have the single position in the file, is there any reason I couldn't use the position as both the start and end values?

UCSC uses a semi-open (0-based) interval (Cheat Sheet For One-Based Vs Zero-Based Coordinate Systems) using the same start/end (1 based) means you're trying to submit a deletion = 0 bases.

Try to submit:

(variant.chrom):(variant.start-1)-(variant.start)

or use something like LiftoverVcf https://broadinstitute.github.io/picard/command-line-overview.html#LiftoverVcf

ADD COMMENTlink written 3.3 years ago by Pierre Lindenbaum131k

Thanks very much, I'll make sure to do that!

ADD REPLYlink written 3.3 years ago by cookersjs20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2012 users visited in the last hour