I am trying to use IGB for ChIP-Seq analyses and comparisons. I am starting with fastq files. I am grooming the data with the FASTQ Groomer and then aligning the data using Map w/ Bowtie for Illumina. I then filter the data using Filter Sam and convert to a BAM file. When I try to use MACs peak calling I repeatedly get a float division by zero error.
INFO @ Thu, 13 Aug 2015 14:07:40:
# ARGUMENTS LIST:
# name = MACS_in_Galaxy
# format = BAM
# ChIP-seq file = /galaxy-repl/main/files/012/160/dataset_12160005.dat
# control file = None
# effective genome size = 2.70e+09
# tag size = 36
# band width = 200
# model fold = 10
# pvalue cutoff = 1.00e-05
# Ranges for calculating regional lambda are : peak_region,1000,5000,10000
INFO @ Thu, 13 Aug 2015 14:07:40: #1 read tag files...
INFO @ Thu, 13 Aug 2015 14:07:40: #1 read treatment tags...
Traceback (most recent call last):
File "/galaxy/main/deps/macs/1.3.7.1/devteam/package_macs_1_3_7_1/a7ea583a35d2/bin/macs", line 273, in <module>
main()
File "/galaxy/main/deps/macs/1.3.7.1/devteam/package_macs_1_3_7_1/a7ea583a35d2/bin/macs", line 71, in main
options.bg_redundant = bg_r (treat,options.max_dup_tags) # average
File "/galaxy/main/deps/macs/1.3.7.1/devteam/package_macs_1_3_7_1/a7ea583a35d2/bin/macs", line 225, in bg_r
return float(total_duplicates)/(total_tags)
ZeroDivisionError: float division by zero
What do I need to do to resolve this error?
Colin
