About ATAC-seq call peak (hmmratac) and differential analysis via MACS3 #675
Unanswered
DayTimeMouse
asked this question in
Q&A
Replies: 0 comments
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Uh oh!
There was an error while loading. Please reload this page.
-
Hi,
Thanks for developing this function to call peaks for ATAC-Seq.
I ran
macs3 hmmratac -i final.bam -n treat --cutoff-analysis-only
. Then, I got the report(see the attachment), according to the document, -l may set to45
(that can capture moderate number ( about 10k ) of peaks with normal width ( average length 500-1000bps long)), but I'm confused about setting-u
and-c
parameter.The report:
treat_cutoff_analysis.txt
Another question, I have control and treatment, which are no repeats. So, can I use MASC3 to do differential analysis? Could you give me a detail approach?
Thanks in advance.
Beta Was this translation helpful? Give feedback.
All reactions