[BioPython] formatdb over python code
Christof Winter
winter at biotec.tu-dresden.de
Tue Jun 10 10:13:06 UTC 2008
Stefanie Lück wrote, On 06/10/08 10:38:
> Hi!
>
> Does someone know's whether it's possible to make a database with formatdb
> (NCBI) via python code (among Windows) and not over the console?
Here is the Python code I use for that:
cmd = "formatdb -i %s -p T -o F" % database
os.system(cmd)
-p T specifies protein sequences, -o T creates indexes, but fails if the fasta
file does not follow the defline format (see
http://en.wikipedia.org/wiki/Fasta_format#Sequence_identifiers). If it fails,
use -o F.
Christof
More information about the Biopython
mailing list