[Bioperl-l] indexing and accessing local database
Qunfeng Dong
qfdong@iastate.edu
Tue, 9 Jul 2002 13:28:36 -0500 (CDT)
No,it didn't print anything, even though AF468672 is one of the accession number
in the indexed fasta file.
Qunfeng
> Qunfeng,
>
> But did the "print $seq->seq" print anything printable?
>
> Brian O.
>
> -----Original Message-----
> From: Qunfeng Dong [mailto:qfdong@iastate.edu]
> Sent: Tuesday, July 09, 2002 2:18 PM
> To: Brian Osborne
> Cc: bioperl-l@bioperl.org
> Subject: RE: [Bioperl-l] indexing and accessing local database
>
> Brian, Thanks for the suggestion, I didn't use Bio::DB::Fasta because I
> couldn't
> find any Doc about how to reuse the index instead of generating it every
> time.
> Anyway, I tried the following code
> ****************************************
> #!/usr/bin/perl -w
> #name: DBindex.pl
> use Bio::DB::Fasta;
> $File_Name = shift;
> $inx = Bio::DB::Fasta->new($File_Name);
> $seq = $inx->get_Seq_by_id("AF468672");
> print $seq->seq;
> *****************************************
> then ran it as
> % DBindex.pl fastaFile
> then I got the following error msg
>
> Use of uninitialized value in numeric lt (<) at
> /usr/lib/perl5/site_perl/5.6.1/Bio/DB/Fasta.pm line 604.
> Odd number of elements in hash assignment at
> /usr/lib/perl5/site_perl/5.6.1/Bio/DB/Fasta.pm line 969.
> Use of uninitialized value in numeric gt (>) at
> /usr/lib/perl5/site_perl/5.6.1/Bio/DB/Fasta.pm line 817.
>
>
>
> Something must be wrong. Maybe I shouldn't use BioPerl and stick to the
> blast
> database with fastacmd :-)
>
>
> Qunfeng
>
>
> > Qunfeng,
> >
> > I'm also seeing a segmentation fault on my Cygwin/Windows machine with
> code
> > like yours, using Bio::Index::Fasta. But I'm surprised that I can index
> > using Bio::DB::Fasta. If you don't care which indexing module you're using
> > then I suggest you try something like this and see what happens:
> >
> > use Bio::DB::Fasta;
> > $File_Name = shift;
> > $inx = Bio::DB::Fasta->new($File_Name);
> > $seq = $inx->get_Seq_by_id("P8413");
> > print $seq->seq;
> >
> >
> > Brian O.
> >
> > -----Original Message-----
> > From: bioperl-l-admin@bioperl.org [mailto:bioperl-l-admin@bioperl.org]On
> > Behalf Of Qunfeng Dong
> > Sent: Tuesday, July 09, 2002 12:44 PM
> > To: Jason Stajich
> > Cc: bioperl-l@bioperl.org
> > Subject: Re: [Bioperl-l] indexing and accessing local database
> >
> > It's Linux and perl 5 (and when I install bioperl, perl is also installed)
> >
> >
> > > What operating system are you using and version of Perl?
> > >
> > > On Tue, 9 Jul 2002, Qunfeng Dong wrote:
> > >
> > > > Hi,
> > > >
> > > > I am trying to follow the sample code from
> > > >
> >
> http://bioperl.org/Core/POD/bptutorial.html#III_1_1_Accessing_remote_databas
> > > >
> > > > to build a local fasta sequence database
> > > > #!/usr/bin/perl -w
> > > > #scriptName - BuildIndex.pl
> > > > use Bio::Index::Fasta; # using fasta file format
> > > > $Index_File_Name = shift;
> > > > $inx = Bio::Index::Fasta->new(
> > > > -filename => $Index_File_Name,
> > > > -write_flag => 1);
> > > > $inx->make_index(@ARGV);
> > > >
> > > >
> > > > then I ran the script BuildIndex.pl as
> > > > % BuildIndex.pl indexName FastaSeqName
> > > > and I got Segmentation fault. Any idea why?
> > > >
> > > > Qunfeng Dong
> > > >
> > > >
> > > >
> > > >
> > > > _______________________________________________
> > > > Bioperl-l mailing list
> > > > Bioperl-l@bioperl.org
> > > > http://bioperl.org/mailman/listinfo/bioperl-l
> > > >
> > >
> > > --
> > > Jason Stajich
> > > Duke University
> > > jason at cgt.mc.duke.edu
> > >
> > >
> >
> >
> > Qunfeng Dong, Ph.D
> > ZmDB Manager
> > 2104 Molecular Biology Building
> > Iowa State University
> > Ames, IA 50010
> >
> >
> >
> >
> >
> >
> > _______________________________________________
> > Bioperl-l mailing list
> > Bioperl-l@bioperl.org
> > http://bioperl.org/mailman/listinfo/bioperl-l
> >
> >
> > _______________________________________________
> > Bioperl-l mailing list
> > Bioperl-l@bioperl.org
> > http://bioperl.org/mailman/listinfo/bioperl-l
> >
>
>
> Qunfeng Dong, Ph.D
> ZmDB Manager
> 2104 Molecular Biology Building
> Iowa State University
> Ames, IA 50010
>
>
>
>
>
>
>
>
Qunfeng Dong, Ph.D
ZmDB Manager
2104 Molecular Biology Building
Iowa State University
Ames, IA 50010