SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BLAST+ creating custom blast database and using blast+ filtering features deniz Bioinformatics 3 07-07-2019 09:04 AM
makeblastdb error while using repeat masker locally novino11 Bioinformatics 1 09-08-2011 11:06 AM
BLAST database error - when changing to new BLAST+ local program biobio Bioinformatics 4 06-15-2011 06:20 AM
makeblastdb for hg18 error openchromatin Bioinformatics 0 03-16-2011 06:24 PM
Blast makeblastdb problem Lazy D-Sledge Bioinformatics 0 12-12-2010 11:01 AM

Reply
 
Thread Tools
Old 03-04-2014, 11:45 AM   #1
milo0615
Member
 
Location: Walnut, California

Join Date: Dec 2012
Posts: 39
Default Blast 'makeblastdb' Error: (1431.1)

Hello,

I am creating a database for a 2.7Gb merged sequences that I have, but I get the following error:

command used: /usr/local/bin/makeblastdb -in Merged_Lsat_1_v4.fa -dbtype nucl -parse_seqids


Building a new DB, current time: 03/04/2014 19:26:23
New DB name: Merged_Lsat_1_v4.fa
New DB title: Merged_Lsat_1_v4.fa
Sequence type: Nucleotide
Keep Linkouts: T
Keep MBits: T
Maximum file size: 1000000000B
Error: (1431.1) FASTA-Reader: Warning: FASTA-Reader: First data line in seq is about 50% ambiguous nucleotides (shouldn't be over 40%)

Does anyone have a workaround or any suggestions on what I can do to fix this problem? blast doesnt quit after the error message, it just hangs in there.

Thank you,

-Milo
milo0615 is offline   Reply With Quote
Old 03-04-2014, 02:06 PM   #2
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

Is the warning correct (do you have lots of IUPAC ambiguity codes instead of just A,C,G,T)?

If not, it could be a sign that there was an error interpreting your FASTA file.
maubp is offline   Reply With Quote
Old 03-04-2014, 03:07 PM   #3
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,992
Default

If Merged_Lsat_1_v4.fa was edited on a PC/Mac and then moved to a UNIX server then this could be an indication of the file format (PC/Mac vs UNIX) problem.
GenoMax is offline   Reply With Quote
Old 03-14-2014, 01:11 PM   #4
Naibin Duan
Junior Member
 
Location: ithaca

Join Date: Jan 2014
Posts: 3
Default

I also have this proble ,but I got the answer ,you should check your file :Merged_Lsat_1_v4.fa.

there must be nnnnnnnnnnnnnnnnnn....... in it ,Good luck !
Naibin Duan is offline   Reply With Quote
Old 05-16-2014, 08:56 AM   #5
kbradnam
Member
 
Location: Davis, CA

Join Date: May 2011
Posts: 53
Default

I blogged about this very topic yesterday:

http://www.acgt.me/blog/2014/5/15/fu...rom-ncbi-blast

It is not a great 'error' message'. Actually, it's a warning message that also has 'Error' pre-pended to it. It only looks at the first line of your sequence to determine this. So if you the same amount of Ns at the end of the sequence, you won't see the error.

So it's mostly harmless, and a little pointless. You will see that message if your first line has >40% ambiguous characters even if

a) the first line is only 5 bp
b) the rest of the entire sequence is free from ambiguity codes
kbradnam is offline   Reply With Quote
Old 01-20-2015, 08:37 AM   #6
ebioman
Member
 
Location: Switzerland

Join Date: Aug 2013
Posts: 41
Default

Hi
But does it essentially add the sequence to the database and (real warning) ,
or does it skip this sequence (error)? It isnt really conclusive to me
ebioman is offline   Reply With Quote
Old 01-20-2015, 06:46 PM   #7
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

My interpretation of the message is this is a warning only, and the sequence was still added to the database.

You can check this using blastdbcmd (either look at the total number of sequences, or use this tool to extract FASTA sequences from the database for an in depth check).
maubp is offline   Reply With Quote
Reply

Tags
blast, error:1431.1, fasta-reader, makeblastdb

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 03:39 PM.


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