[BioPython] Running Blast locally
Muthuraman, Manickam
manickam.muthuraman at wur.nl
Fri Jun 16 13:12:08 UTC 2006
Dear Peter
i am not clear about the subtopic running blast locally
let me explain in detail
i have blast executable files in my home directory i.e
/home/manickam/blast/blastall
i have my database files of nr,swissprot,pdb in /usr/junk/
the files which i can see under /usr/junk/ folder are
nr.00.phr
nr.00.ppi
nr.01.phr
nr.01.ppi nr.pal
pdbaa.00.msk
lot in there and there extenstions are *.phr , ppi ,pal,msk,psq
i am not clear from the manual where do i need to provide the input sequences and how to i store the out put after running the local blast.
below is the following code which i tried and it works but b_record is none.
mport os
my_blast_db=os.path.join(os.getcwd(),'at-est','a-cds-10-7.fasta')
my_blast_file=os.path.join(os.getcwd(),'at-est','test_blast','sorghum_est-test.fasta')
my_blast_exe=os.path.join(os.getcwd(),'blast','/home/manickam/blast/blastall')
from Bio.Blast import NCBIStandalone
blast_out,error_info=NCBIStandalone.blastall(my_blast_exe,'blastp',my_blast_db,my_blast_file)
b_parser=NCBIStandalone.BlastParser()
b_iterator=NCBIStandalone.Iterator(blast_out,b_parser)
b_record=b_iterator.next()
while 1:
b_record=b_iterator.next()
if b_record is None:
break
for alignment in b_record.alignments:
for hsp in alignment.hsps:
print 'seq:',alignment.title
from
manickam
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/ms-tnef
Size: 3446 bytes
Desc: not available
Url : http://lists.open-bio.org/pipermail/biopython/attachments/20060616/2de84992/attachment-0001.bin
More information about the Biopython
mailing list