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

Output as .cool format instead of .mcool #9

Open
mb6434 opened this issue Jan 26, 2023 · 2 comments
Open

Output as .cool format instead of .mcool #9

mb6434 opened this issue Jan 26, 2023 · 2 comments

Comments

@mb6434
Copy link

mb6434 commented Jan 26, 2023

Hi,

Thank you very much for the development of this tool. I am facing an issue regarding the output format.

Here is the cmd line which I used:

HiCLift --input $file.cool --input-format cooler --out-pre test --output-format cool --out-chromsizes mm10_2.sizes --in-assembly mm9 --out-assembly mm10

I would like to know how could you got a .cool format instead of .mcool.

And also, I am not familiar with the argument --high-res, would you please develop a little bit?!

Thank you for your help and time,

Best regards,

Mathilde

@XiaoTaoWang
Copy link
Owner

Hi Mathilde,

.mcool is just an extension of .cool which stores contact matrices at multiple resolutions, and it's quite straightforward to extract the contact matrices at a specific resolution into a .cool file. Suppose your output .mcool file name is "test.mm10.mcool", and you want to extract the contact matrix at 5kb, just execute the cooler cp command below::

$ cooler cp test.mm10.mcool::resolutions/5000 test.mm10.5kb.cool

You can use the cooler ls command to list all individual cool files (cool URIs) within a .mcool file::

$ cooler ls test.mm10.mcool

By default, the highest resolution of contact matrices reported by HiCLift is 5kb. By specifying "--high-res", the highest reported resolution will be set to 1kb.

Best,
Xiaotao

@mb6434
Copy link
Author

mb6434 commented Jan 27, 2023

Hi Xiaotao,

Thank you very much for your answer and to give to me solutions. This is awesome!

Best reagrds,

Mathilde

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

2 participants