Go Back   SEQanswers > Applications Forums > Epigenetics

Similar Threads
Thread Thread Starter Forum Replies Last Post
SICER (histone ChIP-Seq) rebrendi Bioinformatics 2 02-23-2012 09:46 AM
The difference on wiggle file from SICER and MACS tujchl Bioinformatics 7 01-05-2012 11:58 PM
Remove redundancy from 454 data sklages 454 Pyrosequencing 2 09-23-2010 08:14 AM
threshold alperyilmaz RNA Sequencing 1 02-09-2010 11:09 AM

Thread Tools
Old 06-28-2012, 10:31 AM   #1
Location: Maryland, USA

Join Date: Oct 2009
Posts: 27
Default SICER: redundancy threshold parameter

I am looking for some help in setting the 'redundancy threshold' parameter that is to be specified while running SICER for calling peaks from histone data. The protein of interest here is a high mobility group nucleosome binding protein (HMGN1) which is predicted to have broad histone-like peaks.

As I do not want to remove duplicate reads from the input BED file, what do I specify for the 'redundancy threshold' on the command line. I am not able to leave it blank, as SICER expects something in that field.

I will appreciate any help or pointer in setting this parameter.

ParthavJailwala is offline   Reply With Quote
Old 08-07-2012, 03:55 AM   #2
Location: Czech Republic

Join Date: Aug 2012
Posts: 10

The redundancy threshold parameter is handled by python script src/

if current_count <= cutoff:
     retained += 1
If you really want to retain all the duplicates in your bed file then you have to set this parameter (cutoff) higher than maximum number of duplicates in your bed file. 2147483647 should be in most cases OK

Duplicate is defined like a reed with the same start and stop coordinate, some orientation on the same chromosome. If there will be for example 2147483645 reads of the same specification above it will be then retained.
Kulvait is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 12:41 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO