Question: Visualize smooth ChIP peaks with Homer makeUCSCfile
gravatar for varsha619
22 months ago by
varsha61970 wrote:

Hello, I am trying to use Homer makeUCSCfile for my ChIP seq data using as follows -

makeUCSCfile sample.tagdir/ auto -style chipseq > sample.bedGraph.gz

Although I am able to visualize the ChIP peaks in UCSC, I would like to get a smoothened plot instead, please let me know if there is an option that can be used to smoothen the peaks, Homer or otherwise. Thank you.

ucsc chip-seq homer • 950 views
ADD COMMENTlink modified 19 months ago by macmath130 • written 22 months ago by varsha61970
gravatar for macmath
19 months ago by
macmath130 wrote:

You can try check the below parameters Modifying Read Coverage You can manually set the fragment lengths that are visualized and shift their positions, both of which can be useful: -fragLength <# | given> : sets the fragment length, default: uses fragmentLengthEstimate in the tagInfo.txt file of the tag directory. If you want to visualize how the signal changes over large regions, it can be useful to set the fragment length to a very large value (i.e. 10000). If you want to visualize the exact length of the reads, use "-fragLength given".

-adjust <#> : adjust the position of the read by this amount from the 5' end. For example, -adjust -10 would start the coverage 10 bp upstream. This useful when the 5' end of the read represents a localized signal, i.e. DNase nicking site, as opposed to a ChIP-Seq fragment, which implies the factor binds downstream from the 5' end.

-tbp <#> : limit the number of reads considered per position, default: no limit. i.e. "-tbp 1" only counts one read per position.

-inputFragLength <#>, -inputAdjust <#>, -inputtbp <#> work the same for input directories if calculating a ratio.

ADD COMMENTlink written 19 months ago by macmath130
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 715 users visited in the last hour