in reply to Re: BioPerl StandAloneBlast is returning unexpected undefined SearchIO object
in thread BioPerl StandAloneBlast is returning unexpected undefined SearchIO object

If you tried running this code for a single sequence, translating it into the 6 frames and then tried blasting it what would come??
Unfortunately I get the exact same error :(

Ok I might have been a little unclear on this. I translate one nucleotide sequence into six amino acid sequences. A stop codon in the nucleotide sequence will be translated into an asterisk *. It is therefore easier to sort out the sequences after translating them.

I'm sorry I don't understand your concern regarding my use of blastp. I am blasting protein sequences (amino acids).
  • Comment on Re^2: BioPerl StandAloneBlast is returning unexpected undefined SearchIO object

Replies are listed 'Best First'.
Re^3: BioPerl StandAloneBlast is returning unexpected undefined SearchIO object
by biohisham (Priest) on Jun 25, 2009 at 22:16 UTC
    get the six different reading frames from the nucleotides, decide on which ones to remove because they got any one of the three stop codons, do that first, then translate the rest of the sequences into proteins and do the blasting. that would save you from using the * to mark where a stop codon is...I know I am not able to give you a solution but I am just suggesting something that could land you anywhere near it... I will really advice you to use the * only in Typeglob notations in Perl.. I had a feeling that you used blastp on a protein sequence cuz I did not read your code in a clear way. Be a bit more neat instead of repeating the same if statement over and over, your program could be a way shorter and effective too...
    Excellence is an Endeavor of Persistence. Chance Favors a Prepared Mind