I'm running MACS on some .bam files, and I'm coming across this error on one particular .bam file (but not all .bam files) and I'm not sure what to do about it:
INFO @ Fri, 05 Jun 2015 11:14:58:
# ARGUMENTS LIST:
# name = 1487
# format = AUTO
# ChIP-seq file = 1487_filtered.bam
# control file = 1580_input_filtered2.bam
# effective genome size = 2.70e+09
# band width = 300
# model fold = 10,30
# pvalue cutoff = 1.00e-05
# Large dataset will be scaled towards smaller dataset.
# Range for calculating regional lambda is: 1000 bps and 10000 bps
INFO @ Fri, 05 Jun 2015 11:14:58: #1 read tag files...
INFO @ Fri, 05 Jun 2015 11:14:58: #1 read treatment tags...
Traceback (most recent call last):
File "/usr/local/bin/macs14", line 5, in <module>
pkg_resources.run_script('MACS==1.4.2', 'macs14')
File "/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/pkg_resources.py", line 489, in run_script
self.require(requires)[0].run_script(script_name, ns)
File "/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/pkg_resources.py", line 1214, in run_script
exec script_code in namespace, namespace
File "/Library/Python/2.7/site-packages/MACS-1.4.2-py2.7.egg/EGG-INFO/scripts/macs14", line 366, in <module>
File "/Library/Python/2.7/site-packages/MACS-1.4.2-py2.7.egg/EGG-INFO/scripts/macs14", line 60, in main
File "/Library/Python/2.7/site-packages/MACS-1.4.2-py2.7.egg/EGG-INFO/scripts/macs14", line 333, in load_tag_files_options
File "build/bdist.macosx-10.9-intel/egg/MACS14/IO/Parser.py", line 60, in guess_parser
Exception: Can't detect format!
Other forum postings suggest that I should patch Macs14, as specified in http://ashishagarwal.org/2011/03/26/macs14-patch/
However, the error this patch would fix is not the one I'm dealing with. Has anyone else had this error?
For reference, here is the script I'm using to run MACS:
macs14 -t $name1 -c $name2 -g hs -n $name3 -B -S --call-subpeaks
Uchillama
INFO @ Fri, 05 Jun 2015 11:14:58:
# ARGUMENTS LIST:
# name = 1487
# format = AUTO
# ChIP-seq file = 1487_filtered.bam
# control file = 1580_input_filtered2.bam
# effective genome size = 2.70e+09
# band width = 300
# model fold = 10,30
# pvalue cutoff = 1.00e-05
# Large dataset will be scaled towards smaller dataset.
# Range for calculating regional lambda is: 1000 bps and 10000 bps
INFO @ Fri, 05 Jun 2015 11:14:58: #1 read tag files...
INFO @ Fri, 05 Jun 2015 11:14:58: #1 read treatment tags...
Traceback (most recent call last):
File "/usr/local/bin/macs14", line 5, in <module>
pkg_resources.run_script('MACS==1.4.2', 'macs14')
File "/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/pkg_resources.py", line 489, in run_script
self.require(requires)[0].run_script(script_name, ns)
File "/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/pkg_resources.py", line 1214, in run_script
exec script_code in namespace, namespace
File "/Library/Python/2.7/site-packages/MACS-1.4.2-py2.7.egg/EGG-INFO/scripts/macs14", line 366, in <module>
File "/Library/Python/2.7/site-packages/MACS-1.4.2-py2.7.egg/EGG-INFO/scripts/macs14", line 60, in main
File "/Library/Python/2.7/site-packages/MACS-1.4.2-py2.7.egg/EGG-INFO/scripts/macs14", line 333, in load_tag_files_options
File "build/bdist.macosx-10.9-intel/egg/MACS14/IO/Parser.py", line 60, in guess_parser
Exception: Can't detect format!
Other forum postings suggest that I should patch Macs14, as specified in http://ashishagarwal.org/2011/03/26/macs14-patch/
However, the error this patch would fix is not the one I'm dealing with. Has anyone else had this error?
For reference, here is the script I'm using to run MACS:
macs14 -t $name1 -c $name2 -g hs -n $name3 -B -S --call-subpeaks
Uchillama