Znn3bq.jpeg
²é¿´: 538  |  »Ø¸´: 1

houqingzhen

гæ (³õÈëÎÄ̳)

[ÇóÖú] python-×Ô¶¯·ÃÎÊÍøÂ磬½«htm½á¹û´æÅ̵½±¾µØ

ÐÅϢѧÐÂÊÖÇóÖú£º

ÎÒÏÖÔÚÓиöÎĵµ£¬ÏëÖð¸ö¶ÁÈ¡ÎĵµÀïµÄÄÚÈÝ£¬È»ºóÔÚÍøÒ³Éϵã»÷£¬×îºó½«·µ»ØµÄhtml½á¹û´æÅÌ¡£
ÍøÒ³ÍøÖ·ÊÇ
http://www.ncbi.nlm.nih.gov/COG/old/xognitor.html£¬
ÐèҪѡÔñµÄÑ¡ÏîÊÇBeTs to 17 clades£¬µã»÷µÄÑ¡ÏîÊÇcompare to COGs£¬
×îºó½«·µ»ØµÄhtml ½á¹û´æÅÌ¡£
ÍøÉÏÓиö½Å±¾£¬µ«ÊÇÎÒ²»ÓÃNCBIWWWÕâ¸öÄ£¿é£¬Ï£ÍûÄÜÓÃurl


ÔÚÍøÉϲéÕÒµÄÏà¹Ø½Å±¾ÊÇ£º
from Bio import Fasta
from Bio.Blast import NCBIWWW

filename_for_blast = 'seq/all.fasta' # Òª×öblastµÄfastaÎļþ

parser = Fasta.RecordParser()

file_for_blast = file(filename_for_blast, 'r')
iterator = Fasta.Iterator(file_for_blast, parser)
while 1:
cur_record = iterator.next()
if cur_record is None:
break
print u'´ò¿ªÊä³öÎļþ%s' % cur_record.title+'.htm'
file_for_blast_result = file(cur_record.title+'.htm', 'w')
print 'Now BLAST >', cur_record.title
''' NCBIWWW.blastº¯ÊýĬÈϵÄblast²ÎÊýÊÇ:
expect='10', alignments='50', timeout=20
'''
try:
blast_result = NCBIWWW.blast('blastn', 'nr', cur_record, \
expect='0.001', alignments='1', timeout=1)
file_for_blast_result.write(blast_result.read())
print u'OK£¬¸ã¶¨Ò»¸ö'
except:
print u'³¬Ê±ÁË£¬×öÏÂÒ»¸ö'
file_for_blast_result.close()
file_for_blast.close()

print u'È«²¿¸ã¶¨!'
raw_input()
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

jumper158

½ð³æ (ÕýʽдÊÖ)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

Õâ¸ö½Å±¾ÓÐÎÊÌ⣬ÄãµÄprint Óï·¨´òµÄ¾ÍÓÐÎÊÌ⣬Õâ¸öÍøÖ·ÊÇ×öʲôÓõģ¿ÄÜÓÐÏêϸµÄÔ­°æÃ´£¿Ëµ²»¶¨ÎÒ¿ÉÒÔ°ïÉÏÄã
ÓÎ×ßÓÚÅ£AÓëÅ£CÖ®¼ä£¬Ë­µÄÇà´º²»¿à±ï
2Â¥2012-06-10 20:49:32
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ houqingzhen µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼²©] »¯Ñ§×¨ÒµÉ격 +3 ÕÔ×ÓÑò 2026-05-23 4/200 2026-05-24 18:10 by ¹¤´óѧ³¤
[½Ìʦ֮¼Ò] ÂÛÎij·¸åÁË +4 bjvtcliu 2026-05-24 7/350 2026-05-24 17:29 by bjvtcliu
[²©ºóÖ®¼Ò] ÊÛSCIÒ»ÇøT0PÎÄÕ£¬ÎÒ:8.O.5.5.1.O.5.4,¿ÆÄ¿ÆëÈ«,¿É+¼± +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 17:10 by 75ui6h7z2t
[¿¼²©] ÊÛSCIÒ»ÇøT0PÎÄÕ£¬ÎÒ:8.O.5.5.1.O.5.4,¿ÆÄ¿ÆëÈ«,¿É+¼± +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 17:01 by 75ui6h7z2t
[¿¼ÑÐ] ÊÛSCIÒ»ÇøT0PÎÄÕ£¬ÎÒ:8.O.5.5.1.O.5.4,¿ÆÄ¿ÆëÈ«,¿É+¼± +3 a2tycdlnq1 2026-05-23 5/250 2026-05-24 16:21 by hhx1yx9evi
[ÂÛÎÄͶ¸å] ÊÛSCIÒ»ÇøT0PÎÄÕ£¬ÎÒ:8.O.5.5.1.O.5.4,¿ÆÄ¿ÆëÈ«,¿É+¼± +3 a2tycdlnq1 2026-05-23 4/200 2026-05-24 16:16 by hhx1yx9evi
[»ù½ðÉêÇë] ºÓ±±Ê¡×ÔÈ»¿ÆÑ§»ù½ð +6 Peterchao 2026-05-18 9/450 2026-05-24 16:02 by 130067131
[˶²©¼ÒÔ°] ÊÛSCIÒ»ÇøT0PÎÄÕ£¬ÎÒ:8.O.5.5.1.O.5.4,¿ÆÄ¿ÆëÈ«,¿É+¼± +4 pmo95bazuy 2026-05-23 8/400 2026-05-24 15:56 by 1uy1ht2y9r
[»ù½ðÉêÇë] Î÷°²½»´óÐÂýѧԺ¸±Ôº³¤Óó·¸åÂÛÎĽáÌâ +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[½Ìʦ֮¼Ò] ij211´óѧ½Ìʦ°Ñ¸öÈ˽Ìʦ¹Ù·½Ö÷Ò³¸Ä³É:ÎÒÅÜÁËÎÒÅÜÁËÎÒÅÜÁË!¹ÙÐûÅÜ·£¡ +4 zju2000 2026-05-21 5/250 2026-05-24 09:35 by songwz
[¿¼²©] 26/27É격×Ô¼ö 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[»ù½ðÉêÇë] ÇàB·¢ËÍÉÏ»á֪ͨÁËÂð +5 chemBioBro 2026-05-22 7/350 2026-05-23 12:35 by zhuifengzhy
[¿¼²©] ²©Ê¿ÉêÇë +3 ½¹ÏþÃ÷ 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[ÂÛÎÄͶ¸å] Ͷ¸åÇóÖú£¬ÆÚ¿¯ +4 Ï£¼½£¬ÓÐÊé¶Á 2026-05-20 8/400 2026-05-22 10:16 by Ï£¼½£¬ÓÐÊé¶Á
[ÎÄѧ·¼²ÝÔ°] Ï×Ѫ¸Ð´¥ +7 ѽѽºÃɵ 2026-05-19 13/650 2026-05-21 20:15 by ѽѽºÃɵ
[»ù½ðÉêÇë] ¹ú×ÔÈ»ÆÀ·Ö +4 ÎÞÃûÕßµÇɽ 2026-05-20 5/250 2026-05-21 16:35 by swuq
[»ù½ðÉêÇë] ¹ú×ÔÈ»ÉÏ»áÒªÇó +7 ÎÞÃûÕßµÇɽ 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[Óлú½»Á÷] ·´Ó¦ºÜ²î£¬´óÁ¿Ô­ÁÏûÓз´Ó¦ 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[¿¼²©] Èç¹û¹¤×÷ÁËÏë¶Á²©£¬¿ÉÒԱ߹¤×÷±ß¶ÁÈ«ÈÕÖÆÂ 30+3 Ìú´ï»ð³µ 2026-05-18 5/250 2026-05-20 09:33 by tfang
[¿¼²©] ²©Ê¿ÉêÇë +5 ÐÇ¡­¡­ 2026-05-18 6/300 2026-05-18 23:49 by ºýºýͿͿºÃ
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û