Skip to content

Instantly share code, notes, and snippets.

@pcomans

pcomans/gist:985332

Created May 22, 2011 10:22
Show Gist options
  • Save pcomans/985332 to your computer and use it in GitHub Desktop.
Save pcomans/985332 to your computer and use it in GitHub Desktop.
Download ZIP
Bio::FastaFormat bug
Raw
gistfile1.rb
require 'rubygems'
require 'bio'
blaster = Bio::Blast.local("tblastx", "/Users/philipp/Dropbox/IDP/ESTs/Ephydatia_muelleri/Ephydatia_muelleri", "", "/Users/philipp/Downloads/blast-2.2.25/bin/blastall")
sequence = ">EM1c1_molpal [Ephydatia_muelleri]\nGTTTTGTACAAGAAGACCATTCATAAGAACAATATTGATGACCACACAGTACACTAGAAG"
fasta = Bio::FastaFormat.new(sequence)
result = fasta.query(blaster)
puts result.hits.length
fasta.entry()
result = fasta.query(blaster)
puts result.hits.length
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment