Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Forum
Planet
All »
Home
Log In
Question: Intersect tool in Galaxy
0
gravatar for katherinedixon4
4.5 years ago by
katherinedixon4 • 10
Canada
katherinedixon4 • 10 wrote:

Hi there,

I am using the Intersect tool in Galaxy to find overlapping intervals of peak calls from ChIP-seq data. I haven't been able to determine what the values in the fifth column in the output interval file represent, and was hoping you could help? The values generally range from about 50-1000. Thank you!
galaxy chip-seq • 1.4k views
ADD COMMENT • link •
modified 4.4 years ago by Jennifer Hillman Jackson ♦ 25k • written 4.5 years ago by katherinedixon4 • 10
1
gravatar for Jennifer Hillman Jackson
4.4 years ago by
Jennifer Hillman Jackson ♦ 25k
United States
Jennifer Hillman Jackson ♦ 25k wrote:

Hi,

Interval files can vary in the arrangement of columns - this is how they differ from other structured formats such as BED. But in general, if the column order is "BED-like", then this is the score column. How to interpret these depend on the tool used to generate the data. I'd review the manual for the tool to see what the ranges represent. If you click on the pencil icon for the dataset, you'll also maybe see that the 5th column is assigned as "score", or you can assign it yourself, to have it interpreted by the tools that use it (for example, in visualization, such as Trackster). 

Take care, Jen, Galaxy team
ADD COMMENT • link written 4.4 years ago by Jennifer Hillman Jackson ♦ 25k

Hi Jen,

Thank you for your detailed answer! I have also been trying the new MACS2 on the Galaxy testing site, particularly to use the bdgbroadcall function that is available there. However I'm not sure what the last few columns in the output file represent, and haven't been able to find information about them. I was wondering if you might know what these values are? The file looks like this:

chr1 4999196 5006411 macs2_bdgbroadcall.bdg_broadRegion654 0 . 5005546 5005711 0 1 165 6350 0 0 0
chr1 6740397 6741770 macs2_bdgbroadcall.bdg_broadRegion1178 0 . 6741607 6741762 0 1 155 1210 0 0 0

 Thank you!

ADD REPLY • link written 4.4 years ago by katherinedixon4 • 10

Glad that helped! And next time, please ask a new question as a new post, to help us stay organized.

First, I should state that this wrapper is still under development, and on the Test instance, so I wouldn't expect all to be necessarily perfect (due to the wrapper/tool itself but especially configuration of the instance which is under constant revision).

The output file is called:
bdgbroadcall_on_callpeak_treatment.tabular

This is the MACS2 documentation where outputs are defined:
https://github.com/taoliu/MACS/

Best, Jen, Galaxy team

 

ADD REPLY • link written 4.4 years ago by Jennifer Hillman Jackson ♦ 25k

Thanks Jen!

ADD REPLY • link written 4.4 years ago by katherinedixon4 • 10
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 171 users visited in the last hour