Hi Chris, The testing for blast_library can follow the test for a single sequence blast in blast_test.py. Here's the necessary diff, assuming that blast_library is added as a class member to BlastDB
87,88c87 < for id,s in sp.iteritems(): # FOR EVERY SEQUENCE IN SWISSPROT < sp.blast(s,msa,expmax=1e-10, verbose=False) # GET STRONG HOMOLOGS, SAVE ALIGNMENT IN msa --- > sp.blast_library(sp,msa,expmax=1e-10, verbose=False) 120,122c119,120 < for id,s in sp.iteritems(): # FOR EVERY SEQUENCE IN SWISSPROT < sp.blast(s,msa,expmax=1e-10, verbose=False) # GET STRONG HOMOLOGS, SAVE ALIGNMENT IN msa < msa.build(saveSeqDict=True) # DONE CONSTRUCTING THE ALIGNMENT, SO BUILD THE ALIGNMENT DB INDEXES --- > sp.blast_library(sp,msa,expmax=1e-10, verbose=False) # GET STRONG > HOMOLOGS, SAVE ALIGNMENT IN msa > msa.build(saveSeqDict=True) # DONE CONSTRUCTING THE ALIGNMENT, SO BUILD > THE ALIGNMENT DB INDEXES 138c136 < --- > --~--~---------~--~----~------------~-------~--~----~ You received this message because you are subscribed to the Google Groups "pygr-dev" group. To post to this group, send email to pygr-dev@googlegroups.com To unsubscribe from this group, send email to pygr-dev+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/pygr-dev?hl=en -~----------~----~----~----~------~----~------~--~---
