[Genome] Liftover dies when using BED files with 12 fields?
Hiram Clawson
hiram at soe.ucsc.edu
Wed Aug 15 10:44:44 PDT 2007
John:
Please try out these updated liftOver binaries:
ftp://hgdownload.cse.ucsc.edu/apache/htdocs/admin/exe/liftOver.v151.MacOSX.ppc
ftp://hgdownload.cse.ucsc.edu/apache/htdocs/admin/exe/liftOverMerge.v151.MacOSX.ppc
They will lift these bed files with the r,g,b string in column 9. I see they have
an annoying side effect of turning the r,g,b string into an integer, but that integer
will still result in the desired color. I now need to fix that output so the r,g,b
string will remain as-is and not turn into an integer.
--Hiram
John Major wrote:
> Hello Hiram-
>
> I am using the binary, for the PPC MAC.
>
> John
>
> Hiram Clawson wrote:
>
>> 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