[Biopython] Primer3Commandline rookie question

Jukka-Pekka Verta jp.verta at gmail.com
Wed Nov 3 14:58:18 EDT 2010


Hi all,

I'm just starting with Python so mine's a rookie question.

I'm trying to write a little program that would do primer design with Primer3Commandline, yet I can't figure out how to access the output. For example, when I try the same example as in the Application.py file under Primer3Commandline function
 
######

from Bio.Emboss.Applications import Primer3Commandline
from Bio.Emboss import Primer3

sequence_handle = open("sequence", "w").write(">test\nCGGTCGCCGGACTATCAAATGGCCAGCGCTGAATTGTAGTACATTTCTCTGTCCTTTAATTCACAGGATTCTGTCATGTAAGATGCGCGATTACAGGTTTAACTGATCCAAACTTCTGCAAAACCGATAAAGAACCTAAATAAAGGAAAGGGACACTGCATAGATACAAACAAACGAAACGTTATCTCAGCAGGATTCATGGCTTCTGAAATTCAACTCATCACTATTGGACCTTCTCATTTTTGTGAAAAGGCACGATGGGGTTTGGACCGTGCTGGAATTTCTTATAAGGAGTCAAAGCATGCTCCTGTATTTCACATGTTATACACAAGGGGTTTGGGACAGGGAACCAGCTGTCCCAAGCTTGTACTAGGTGAGGGCAAAAATAAAGTGGTGTTGCATGAATCATCAGATATCTTGAAGTTTGCTGATGAGAACATCATGTCTGAAGAGGATCGCTTGTATCCATCGAACCTTGAGCAGTCTGTGCAAGAGTGGGGGCTGCATTCAAGGGAAAGGTTGTTAACTTGGTGTTTAATGCCATTGTTAAAGCCTCTCGTATATAGAACAGTTGGTTGCAACATGCCTGCAGCAAAAGATCATTCTCTTGATAAGATTCGTTGTATATTTTCAAGAGGTAGATAATGTGCTTGCTGATGGCCGGCCTTTCCTTTGTGGCAACGAGTTTACTGCAGCTGACTTGACATTTGCTTCACTTGCTGGTCCTGTACTATGCCCAGTGGGATATGGCACCTATCTAGTACCAACAGAGGAACTTCCATCAGAGATGTCTGAAGTTACCCTTTCTTTAAGAGAGACACTAGCAGGGAAGCATGCATTGAGGATGTACGAAACTGAAAGACACAGGAAAGTGGCTTCGAAGGTGCTAGTGACCTGACTCAGTGCATTTACACCGTGGCCTTTGTGTTCTGTCATATTATTGGCAACTAATATTGATTTACGTAAACTGCATACAGCTTTTCTTTCTTGGGAGTGAATCTTTTGAGCAATTTGCTGTCTTCATATAGCAATAATTTGAGAAAGGACCAAGATTTATTAAACTTCGTATTTTGTTTTTTGAATGTTGATATCTTTGTTGTTATTTGAGCAAGCTTAAGCTGCCTAGAACTATAATTCGGTTATGGCATGGGCATGTATCATTCAACCTTGATCTTTGTGCCTTTTTAACTAGAACCCATTTTTTTCAAAAAA\n")

test = Primer3Commandline(sequence = sequence_handle, auto=True, hybridprobe=True)
test.explainflag = True
test.oligosize=20
test.psizeopt=200
test.outfile = "output.pr3"
print test

#######

from here on I've understood that I have to parse the output with Primer3.read(), but an error message tells me that the output.pr3 file doesn't exist.

#######

test1 = Primer3.read(open("output.pr3", "r"))
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
IOError: [Errno 2] No such file or directory: 'output.pr3'

#######

What would be the proper way to proceed?

Thanks in advance!

JP Verta




More information about the Biopython mailing list