Question: Macs Problem
0
gravatar for Christopher Scharer
8.1 years ago by
Christopher Scharer10 wrote:
Hi, I recently mapped some ChIP-seq data to the mm8 version of the mouse genome with Bowtie to create a SAM file. However, I can not run MACS on the mapped data because I keep getting the error below. Any suggestions?? Thanks, Chris Messages from MACS: INFO @ Thu, 11 Nov 2010 10:07:30: # ARGUMENTS LIST: # name = MACS_in_Galaxy # format = SAM # ChIP-seq file = /galaxy/home/g2main/galaxy_main/database/files/001/727/dataset_1727796 .dat # control file = None # effective genome size = 2.70e+09 # tag size = 25 # band width = 300 # model fold = 32 # pvalue cutoff = 1.00e-05 # Ranges for calculating regional lambda are : peak_region,1000,5000,10000 INFO @ Thu, 11 Nov 2010 10:07:30: #1 read tag files... INFO @ Thu, 11 Nov 2010 10:07:30: #1 read treatment tags... INFO @ Thu, 11 Nov 2010 10:07:43: 1000000 INFO @ Thu, 11 Nov 2010 10:07:57: 2000000 INFO @ Thu, 11 Nov 2010 10:08:08: 3000000 INFO @ Thu, 11 Nov 2010 10:08:21: 4000000 INFO @ Thu, 11 Nov 2010 10:08:33: 5000000 Traceback (most recent call last): File "/home/g2main/linux2.6-x86_64/bin/macs", line 273, in main() File "/home/g2main/linux2.6-x86_64/bin/macs", line 57, in main (treat, control) = load_tag_files_options (options) File "/home/g2main/linux2.6-x86_64/bin/macs", line 252, in load_tag_files_options treat = options.build(open2(options.tfile, gzip_flag=options.gzip_flag)) File "/home/g2main/linux2.6-x86_64/lib/python2.6/MACS/IO/__init__.py", line 1480, in build_fwtrack (chromosome,fpos,strand) = self.__fw_parse_line(thisline) File "/home/g2main/linux2.6-x86_64/lib/python2.6/MACS/IO/__init__.py", line 1500, in __fw_parse_line bwflag = int(thisfields[1]) ValueError: invalid literal for int() with base 10: 'CTCF:7:185:443:687' -- Chris Scharer, PhD Post-doctoral Fellow Laboratory of Dr. Jeremy Boss Dept of Immunology and Microbiology Emory University Atlanta GA 30322 Ph: 404-727-5959
alignment bowtie • 1.1k views
ADD COMMENTlink modified 8.0 years ago by Jennifer Hillman Jackson25k • written 8.1 years ago by Christopher Scharer10
0
gravatar for Jennifer Hillman Jackson
8.0 years ago by
United States
Jennifer Hillman Jackson25k wrote:
Hi Chris, This is the second report of a recent problem with MACS. If you have time, please send us a history with the problem in it. Use: Options -> Share or Publish -> Make History Accessible via Link and email the link back to me. Also please send along any special settings made when using the tool. Thanks for reporting the issue and we look forward to your reply. Hopefully we will be able to quickly determine & correct any issues, Best! Jen Galaxy team -- Jennifer Jackson http://usegalaxy.org
ADD COMMENTlink written 8.0 years ago by Jennifer Hillman Jackson25k
0
gravatar for Jennifer Hillman Jackson
8.0 years ago by
United States
Jennifer Hillman Jackson25k wrote:
Hello Chris, Are you still having problems? Would you have time send a link to your history? Use Options -> Share or Publish -> click to share the link and email back (to me) and we will provide feedback. Thanks! Jen Galaxy team -- Jennifer Jackson http://usegalaxy.org
ADD COMMENTlink written 8.0 years ago by Jennifer Hillman Jackson25k
Hi Chris, It appears that MACS does not properly handle whitespace characters in sequence identifiers. We'll look into a more streamlined way to handle this issue, but in the mean time, the easiest way to fix this issue is to use the FASTQ manipulation tool before mapping your data, similar to this screenshot (the 'From' box contains a space):
ADD REPLYlink written 8.0 years ago by Daniel Blankenberg ♦♦ 1.7k
Please log in to add an answer.

Help
Access

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