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

Problem in creating consensus ATAC peaks #236

Open
gauriprasad08 opened this issue Mar 2, 2023 · 3 comments
Open

Problem in creating consensus ATAC peaks #236

gauriprasad08 opened this issue Mar 2, 2023 · 3 comments

Comments

@gauriprasad08
Copy link

Hi,

I used the project pipeline of PEPATAC to create consensus ATAC peaks for my samples.I get the combined summary stats for all samples (PEPATAC_tutorial_stats_summary.tsv). However, the pipeline fails to run the PEPATAC_collator and do not produce the consensus peaks. I have attached a PEPATAC_collator_log.md file for your reference.

The error does not seem to indicate that any component is missing from the installation, so most likely it is due to a certain feature in my data that the PEPATAC "does not like". But What is that feature is not clear to me?

Kindly help me out to solve this issue.

PEPATAC_collator_log.txt

@nsheff
Copy link
Member

nsheff commented Mar 13, 2023

This script is choking when trying to read in a file. Could it be an empty file? Or is it possible one of your sample names has a space in it?

@gauriprasad08
Copy link
Author

Hi. Thanks for raising this point. I checked my files one by one. None of the files is empty or with space in its name.

@nsheff
Copy link
Member

nsheff commented Mar 14, 2023

Ah, I think I see something. Do you have multiple genomes in this PEP?

You cannot call consensus peaks if the reference genome isn't identical across all samples.

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