how is obtained the BED score values (the 5th column of the bed files) in BED files
2
0
Entering edit mode
5.9 years ago
fusion.slope ▴ 250

Hello,

am wondering if anyone of you know how is obtained the BED score values (the 5th column of the bed files) of BED files. In UCSC is written:

"score - A score between 0 and 1000. If the track line useScore attribute is set to 1 for this annotation data set, the score value will determine the level of gray in which this feature is displayed (higher numbers = darker gray). This table shows the Genome Browser's translation of BED score values into shades of gray.."

But how is it obtained?

Thanks in advance for any help.

ChIP-Seq BED Files Score • 7.4k views
ADD COMMENT
0
Entering edit mode

Typically, the score is usually stored in the fifth column, with the strand in the sixth column. The score value can really be any numerical value related to the associated interval.

ADD REPLY
0
Entering edit mode

Oh yes the 5th column, thanks for the correction. But i do not see the point of this value if we do not know how it is computed... am I wrong? I need to use it to evaluate the quality of the peaks in my scoring function but i do not know at this point how much it is informative.

ADD REPLY
1
Entering edit mode

It depends on the table you are pulling your BED data from. Each table should have an informational or schema page associated with it on the Genome Browser, which describes how the score field is calculated.

ADD REPLY
1
Entering edit mode

From ENCODE methods for peak calling of these cell lines:

https://www.encodeproject.org/experiments/ENCSR000AKB/

file "Description excerpt: Track description for UCSC Genome Browser composite trackā€¦"

SCORE: "Regions of statistically significant signal enrichment. The score associated with each enriched interval is the mean signal value across the interval. (Note that a broad region with moderate enrichment may deviate from the background more significantly than a short region with high signal.)"

is this score the one of the bed file? Am not sure..

ADD REPLY
0
Entering edit mode

Am working with the CTCF files of ENCODE, for the dataset GM12878 (lymphoblastoid cell lines). I tried to search how the table was obtained but i did not find anything. If you have any idea about encode ChIP-seq dataset this will help me a lot. Thanks in advance Alex.

ADD REPLY
0
Entering edit mode

I have edited 6th column to 5th column after Alex notification.

ADD REPLY
1
Entering edit mode
5.9 years ago

The score can mean anything and be computed in any way. For ChIPseq it's usually some sort of peak score that's computed in a peak caller-specific manner.

ADD COMMENT
0
Entering edit mode

So there is not any technical foundation for this score? So why is it used?

ADD REPLY
1
Entering edit mode

There's no predefined definition of the score, but the score in a given file from a given source will have a definition.

ADD REPLY
0
Entering edit mode

Hi, this is not clear to me. If the file is a bam file. Is the score the number of reads in the interval?

ADD REPLY
0
Entering edit mode

There's no possible answer to this, it can be anything.

ADD REPLY
0
Entering edit mode
5.2 years ago

I found this on the TopHat manual: in the output section junctions.bed. A UCSC BED track of junctions reported by TopHat. Each junction consists of two connected BED blocks, where each block is as long as the maximal overhang of any read spanning the junction. The score is the number of alignments spanning the junction. I understand that if you use Tophat, the score you see on each junction is the number of reads that spanning the junction.

ADD COMMENT

Login before adding your answer.

Traffic: 2180 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6