[Bioperl-l] Bio ::seqIO ::tigr

Hilmar Lapp hlapp at gmx.net
Tue Jan 27 12:31:39 EST 2004


A question aside: why do you want to convert to swissprot in order to 
load into biosql? (load_seqdatabase.pl can use any SeqIO reader.)

	-hilmar

On Tuesday, January 27, 2004, at 02:50  AM, matthieu CONTE wrote:

> I currently trying to use the Bio ::seqIO ::tigr module.
> My objective is to download the whole rice genome form Tigr ( adress 
> below)and to integrate it in my BioSQL DB.
> For this I am trying to convert the tigr format in swiss format with 
> the script below
>
>
> use Bio::SeqIO;
>
> my $in = Bio::SeqIO->new(-file 
> =>'</home/conte/pipeline_orthologues/data/orysa_tigr.txt', -format 
> =>'tigr');
>
> my $out = Bio::SeqIO->new(-file => 
> '>/home/conte/pipeline_orthologues/data/orysa_swiss.txt' , 
> -format=>'swiss');
>
> print $out $_ while <$in>;
>
> I obtain:
>
> ------------ EXCEPTION  -------------
> MSG: [19]Required <AUTHOR_LIST> missing
> STACK Bio::SeqIO::tigr::throw 
> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO/tigr.pm:1338
> STACK Bio::SeqIO::tigr::_process_header 
> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO/tigr.pm:700
> STACK Bio::SeqIO::tigr::_process_assembly 
> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO/tigr.pm:535
> STACK Bio::SeqIO::tigr::_process_tigr 
> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO/tigr.pm:453
> STACK Bio::SeqIO::tigr::_process 
> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO/tigr.pm:420
> STACK Bio::SeqIO::tigr::_initialize 
> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO/tigr.pm:90
> STACK Bio::SeqIO::new 
> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO.pm:358
> STACK Bio::SeqIO::new 
> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO.pm:378
> STACK toplevel get_bioseq_tigr.pl:8
>
> Could you please tell me if there is a problem with the parser or with 
> the input data format of Tigr?
>
> Thanks in advance
>
>
>
>
> Matthieu CONTE
> m_conte at hotmail.com
>
> _________________________________________________________________
> MSN Messenger : discutez en direct avec vos amis ! 
> http://www.msn.fr/msger/default.asp
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l
>
>
-- 
-------------------------------------------------------------
Hilmar Lapp                            email: lapp at gnf.org
GNF, San Diego, Ca. 92121              phone: +1-858-812-1757
-------------------------------------------------------------




More information about the Bioperl-l mailing list