python - Blasting remotely from biopython -
i'm trying remotely blast 70 200-nt sequences using biopython. i've been trying hours figure out why following python script won't work. i can work read file contains 1 fasta using seqio.read, when try switch seqio.parse, don't in .xml save file create. any ideas? as side note, if knows option syntax excluding organisms results (as possible when using ncbi website, please let me know). thanks help. matt from bio.blast import ncbiwww bio import seqio import tkinter.filedialog tkfd in_file=tkfd.askopenfilename() record = seqio.parse(in_file, format="fasta") out_file = tkfd.asksaveasfilename() save_file = open(out_file, "w") rec in record: print(rec) result_handle = ncbiwww.qblast("blastn", "nt", rec.format("fasta")) save_file.write(result_handle.read()) result_handle.close() else: save_file.close() this content of in_file i'm using test file (the line formatting of file set @ 80 char, might