[BioPython] clustalw problem using standalone or Bio.Clustalw

Emanuel Hey jodyhey at yahoo.com
Tue Jul 24 04:23:13 UTC 2007


I am running windows XP,  Python 2.5.1  and BioPython
trying to do a multiple alignment on a fasta file.
The file is not large, with just 10 sequences to try
things out.  I'm familiar with clustalw as a command
line program. 

With clustalw in the path the following works fine and
opens a command prompt window with the clustalw menu
>>>os.system('clustalw')

However this line returns a '1' when it should do an
alignment on the sequences contained in c.faa
>>> os.system('clustalw c.faa')

I also tried following the examples in the cookbook
http://biopython.org/DIST/docs/tutorial/Tutorial.html
and the examples in the bioinformatics course
http://www.pasteur.fr/recherche/unites/sis/formation/python/ch11s06.html#quest_run_clustalw

For example,  the following call to a constructed
command line 
>>> alignment = Clustalw.do_alignment(cline)

generates this error result


Traceback (most recent call last):
  File "<pyshell#82>", line 1, in <module>
    alignment = Clustalw.do_alignment(cline)
  File "C:\Program
Files\Python25\Lib\site-packages\Bio\Clustalw\__init__.py",
line 117, in do_alignment
    % (out_file, command_line))
IOError: Output .aln file .\testalign.out not
produced, commandline: clustalw
.\chimp_hemoglobin.gb.faa -OUTFILE=.\testalign.out


Thanks much for any pointers

jhey


       
____________________________________________________________________________________
Sick sense of humor? Visit Yahoo! TV's 
Comedy with an Edge to see what's on, when. 
http://tv.yahoo.com/collections/222



More information about the Biopython mailing list