Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Understanding write_interaction_pairs output #334

Open
Biojenifer opened this issue Nov 11, 2020 · 0 comments
Open

Understanding write_interaction_pairs output #334

Biojenifer opened this issue Nov 11, 2020 · 0 comments

Comments

@Biojenifer
Copy link

Biojenifer commented Nov 11, 2020

Hi!

I have generated a file with pair_wise interactions for a HiC experiment with the following command and I would like to ensure that I am understanding the output:

As you can see, I choose to get genomic coordinates and not bin, and to focus just on chromosome 17.

hic_exp = Experiment('Sample1', resolution = 1000, norm_data = hic_data, identifier='sampA10', exp_type='Hi-C')

hic_exp.write_interaction_pairs(fname = 'interaction_pairs_chr17_1kb.tsv', zscored=False, header=True, true_position=True, uniq=True, focus = '17', format='tsv')

These are the first lines of the output:

elt1	elt2	normalized hi-c
2000	33000	2
2000	153000	1
2000	178000	1
2000	225000	1

This means that, the bin corresponding to the genomic position 2000 (which I understand that includes all positions from chr17:2000 to chr17:2999 as the resolution is 1kb), interacts with the following genomic positions:

chr17:33000-33999
chr17:153000-153999
chr17:178000-178999
chr17:225000-225999

Is this interpretation correct?

@Biojenifer Biojenifer changed the title write_interaction_pairs output Understanding write_interaction_pairs output Nov 11, 2020
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant