[Genome] Liftover dies when using BED files with 12 fields?

Hiram Clawson hiram at soe.ucsc.edu
Wed Aug 15 09:39:49 PDT 2007


Good Morning John:

Are you using the WEB site for this liftOver operation, or are
you using the liftOver binary in a local situation there ?

We're missing some functionality in some of our utilities that do
not properly recognize this column nine in its new form of r,g,b
numbers.  It used to only be a 0 at all times.  I have recently fixed
one of these situations, I would like to fix this one you mention.

As a work-around, although you lose your r,g,b specification, turn
that column 9 into a 0 and then it will lift.

--Hiram

John Major wrote:
> Hello UCSC folks-
> 
> I am trying to convert a large # of BED files from DM2 to DM3.
> These files are all BED files using all of the 12 allowed fields ( 
> http://genome.ucsc.edu/goldenPath/help/customTrack.html#BED ).
> 
> And example line from my BED file is:
> (space separated)
> chr4 646375 646462 miR-1000 0 + 646375 646462 255,0,0 2 42,41 0,46
> chrX 12896008 12896104 miR-1026 0 + 12896008 12896104 255,0,0 2 47,44 0,52
> 
> When I try to run liftover on these files, I get the following error:
> 
> --Reading liftover chains
> --Mapping coordinates
> --invalid unsigned number: "255,0,0"
> 
> According to the documentation for liftover, I should be able to use 
> valid BED files.  These BED files are displayed properly via the Genome 
> Browser, so am i missing something obvious here?  Do I need to enclose 
> the last 4 fields in some kind of quotes?  Or does liftover not actually 
> work with the last 4 BED file fields?
> 
> thanks!
> john major


More information about the Genome mailing list