SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
how to convert general fastq to fastq int format? feng Bioinformatics 21 07-04-2014 12:40 AM
denovo2 conversion pseudo-base for base jbdorr SOLiD 2 08-31-2012 10:28 AM
MACS: output of MACS as an control file? jaaker Bioinformatics 0 08-22-2012 10:09 AM
RNA-Seq: Genome-wide characterization of transcriptional start sites in humans by int Newsbot! Literature Watch 0 03-05-2011 03:01 AM
1st Cycle Int (PF) wouter Illumina/Solexa 8 10-07-2009 12:22 AM

Reply
 
Thread Tools
Old 07-13-2013, 03:01 AM   #1
feralBiologist
Member
 
Location: UK

Join Date: Jun 2011
Posts: 61
Default MACS: ValueError: invalid literal for int() with base 10: '1:N:0:'

I got the above error when trying to run MACS 1.4.2 (I'm copying the parameter/output below). I ran previously the same sample/control with mfiold=32 but got no peaks. Than changed mfold to 10 and got this error. Does anyone know what is going on and what to do about it?

MACS input/output (error message is at the end):

========================================================

effective genome size = 1.00e+09

tag size = 40

band width = 300

model fold = 10

pvalue cutoff = 1.00e-05

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL181494.1, minus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL182298.1, plus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL180583.1, minus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL183373.1, minus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL192011.1, minus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL185054.1, minus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL178673.1, minus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL187281.1, plus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL181933.1, minus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL189987.1, plus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL185212.1, minus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL186088.1, plus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL184880.1, minus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL189261.1, minus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL190030.1, plus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL181273.1, plus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL191240.1, plus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL190257.1, minus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL176051.1, minus strand!

WARNING @ Sat, 13 Jul 2013 04:29:05: NO records for chromosome GL183262.1, minus strand!

....

and plenty of similar warnings about missing chromosomes

....

INFO @ Sat, 13 Jul 2013 04:29:31: #1.2 read input tags...

Traceback (most recent call last):

File "/home/g2main/linux3.2-x86_64/bin/macs", line 273, in

main()

File "/home/g2main/linux3.2-x86_64/bin/macs", line 57, in main

(treat, control) = load_tag_files_options (options)

File "/home/g2main/linux3.2-x86_64/bin/macs", line 256, in load_tag_files_options

control = options.build(open2(options.cfile, gzip_flag=options.gzip_flag))

File "/home/g2main/linux3.2-x86_64/lib/python2.7/MACS/IO/__init__.py", line 1480, in build_fwtrack

(chromosome,fpos,strand) = self.__fw_parse_line(thisline)

File "/home/g2main/linux3.2-x86_64/lib/python2.7/MACS/IO/__init__.py", line 1500, in __fw_parse_line

bwflag = int(thisfields[1])

ValueError: invalid literal for int() with base 10: '1:N:0:'
feralBiologist is offline   Reply With Quote
Reply

Tags
macs chip-seq

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 02:20 AM.


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