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

Issue Visualizing Control Data #83

Open
RomeroMatt opened this issue Apr 13, 2022 · 0 comments
Open

Issue Visualizing Control Data #83

RomeroMatt opened this issue Apr 13, 2022 · 0 comments

Comments

@RomeroMatt
Copy link

Hello,
I have been attempting to optimize CUT&Tag for our cells of interest and I am currently analyzing our dataset along with running the controls found from the protocol found online: https://www.protocols.io/view/cut-amp-tag-data-processing-and-analysis-tutorial-e6nvw93x7gmk/v1?step=18

I have been running the script and been getting the same output as the protocol except for visualizing the data in IGV. I used all the same parameters for SEACR and took my output.bed file and created a bedgraph file with the following command: cut -f1-3,5 /Users/matthewa/Data/scripts/cut_tag/peakCalling/SEACR/${histName}_seacr_control.peaks.stringent.bed > /Users/matthewa/Data/scripts/cut_tag/peakCalling/SEACR/${histName}_seacr_control.peaks.stringent.bedgraph
Once putting into IGV, I get the following peaks (pasted below).
H3K27me3_rep1
Protocol pic:
image

This looks different from what is in the protocol and many of my samples have simple blocks at their genomic sites (also pasted below).
H3K27me3_H3K4me3_Block

Not sure exactly what I am doing wrong - I am still new to CUT&Tag/bioinformatics so it is possible that I am simply doing the wrong thing. Any insight would be great! Thanks!
-Matt

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