[Bioperl-l] problem with Bio::SeqIO->write_seq

Prachi Shah prachi_shroff at yahoo.com
Fri May 2 12:28:06 EDT 2003


Hi,

I was playing around with the Bio::Index::Fasta
module, but got the following error. I index a Fasta
formatted sequence file and then try to search for
sequences with the fetch function.
I copied the code from the documentation of
Bio::Index::Fasta. The only difference is that the
example code in the documentation prints to STDIO and
I wanted it to write to a file. I have a very silly
mistake. I will really appreciate any help.

thanks,
Prachi.

------------- EXCEPTION  -------------
MSG: Did not provide a valid Bio::PrimarySeqI object
STACK Bio::SeqIO::fasta::write_seq
C:/Perl/site/lib/Bio\SeqIO\fasta.pm:166
STACK toplevel intersectionsequence-cmd.pl:367

--------------------------------------


### snippet of code ###
my $inx1 = Bio::Index::Fasta->new(
'-filename' => $Index_File_Name1,
'-write_flag' => 1);
my $out = Bio::SeqIO->new('-format' => 'Fasta','-file'
=> ">$clusterfile");
my $seq = $inx1->fetch($r);
$out->write_seq($seq);
################





__________________________________
Do you Yahoo!?
The New Yahoo! Search - Faster. Easier. Bingo.
http://search.yahoo.com


More information about the Bioperl-l mailing list