This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Formatdb For Chromosome Database?

Hi,

I recently downloaded the new Fugu v5 assembly database from IMCB (http://www.fugu-sg.org/downloads/downloads3.htm). This is my first time using a chromosome level database (I am used to contig and scaffold).

I tried to use formatdb as I normally would to get the output files; however, nothing seems to be happening. Is there a way to use formatdb to format chromosome databases or am I supposed to do it in a different way? Am I missing some steps that precede formatdb?

blast database
• 2,777 views • link
updated by Simon Cockell Simon Cockell • written by User 4143 • User 4143

what is the size of the fasta ? what are your parameters for formatdb ?

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

And what is written to the formatdb.log file?

• 0 views • link
written by Daniel Standage Daniel Standage

1 answer

OK, I downloaded this file and used formatdb from the BLAST 2.2.25 package on OS X.

$ gunzip fugu5_chromosomes.fa.gz
$ formatdb -o T -p F -i fug5_chromosomes.fa
$ blastall -p blastn -d fugu5_chromosomes.fa -i testseq.fa

Works just fine for me. Without any more specific information it is really hard to debug your specific problem. I've packaged the output of formatdb and you can download it from here (94MB).

(for the record, the contents of formatdb.log:

========================[ Apr 1, 2011  2:34 PM ]========================
Version 2.2.25 [Feb-01-2011]
Started database file "fugu5_chromosomes.fa"
Formatted 45 sequences in volume 0
SUCCESS: formatted database fugu5_chromosomes.fa

).

• 1 views • link
written by Simon Cockell Simon Cockell

Log in to answer this question.

More posts like this