[BioRuby] Blast with file as a query option?

Kevin English kenglish at gmail.com
Mon Apr 6 00:13:37 UTC 2009


Hello,
  I have to very large local fasta files that wish to blast against one
another and parse the results in bio ruby. I'm wondering if there is a way
to mimic the behavior of this blast command:

blastall -p blastn -i Large_list_sequences_1.fasta -d Large_list_sequences_2


where Large_list_sequences_2 is a formatted fasta db. My current
implementation opens Large_list_sequences_1.fasta and goes through it
sequence by sequence. It seems to run pretty slow. I'm wondering if I can in
some way do the above blast command and loop through the results and get a
performance gain.

For any curious, my code is on github:

http://github.com/kenglishhi/bioflexrails/tree/master

The file that is doing the blasts is under app/model/biodatabase.rb.

I'm trying to write a rails app uses biosql db and allows this biologist to
organize his sequences. I'm very new to bioinformatics but have a lot
experience with Ruby on Rails.

Thanks in advance for you help.

Kevin



More information about the BioRuby mailing list