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

plotFingerprint bigwig support #1290

Open
annajbott opened this issue Feb 12, 2024 · 1 comment
Open

plotFingerprint bigwig support #1290

annajbott opened this issue Feb 12, 2024 · 1 comment

Comments

@annajbott
Copy link

Hello,

I was wondering if plotFingerprint supports bigwig input files?

The documentation mentions using bigwigs or indexed BAM files as input: https://deeptools.readthedocs.io/en/develop/content/tools/plotFingerprint.html (plotFingerprint randomly samples genome regions (bins) of a specified length and sums the per-base coverage in indexed [BAM][] (or bigWig) files that overlap with those regions), but I don't think bigwigs are actually supported. Could this be added as a feature?

Deeptools version = deeptools 3.5.1
Python version = Python 3.8.18

Command:

plotFingerprint -b test.bw -plot test.png --blackListFileName test_blacklist.bed --region chr10:25249964:25252083 --outQualityMetrics test_QC.tsv --outRawCounts test_counts.tab

Screen ouput:
[E::hts_hopen] Failed to open file test.bw

The same command works with indexed bams, so I assume its two typos in the documentation that mention bigwig files (unless I'm missing a flag or something). I can use BAMs, but bigwigs would be much preferred for my usage of plotFingerprint. Thanks!

@WardDeb
Copy link
Member

WardDeb commented Mar 27, 2024

Hi, Thanks for the report. Plotfingerprint indeed doesn't support bigwigs, I've taken the ambiguity out of the docs.
That said, support for this could be added, but this will take a bit more time to implement.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants