[Bioperl-l] simple PrimarySeq question
Niels Larsen
niels at genomics.dk
Mon Jul 2 20:45:07 UTC 2007
I write hoping someone could show me how to create a PrimarySeq
object without parsing features and all first. The lines below
return
"Can't locate object method "next_seq" via package "Bio::PrimarySeq" at ./tst2 line 16."
whereas calling Bio::SeqIO-> gives no error, but a too big object.
The GenBank record after the __END__ is the "1.gb" file. I could not
find out how from the tutorial or the Bio::PrimarySeq description.
Niels L
#!/usr/bin/env perl
use strict;
use warnings FATAL => qw ( all );
use Data::Dumper;
use Bio::Seq;
use Bio::SeqIO;
my ( $seq_h, $seq );
$seq_h = Bio::PrimarySeq->new( -file => "1.gb", -format => 'genbank' );
# $seq_h = Bio::SeqIO->new( -file => "1.gb", -format => 'genbank' );
$seq = $seq_h->next_seq();
# print Dumper( $seq );
__END__
LOCUS X60065 9 bp mRNA linear MAM 14-NOV-2006
DEFINITION B.bovis beta-2-gpI mRNA for beta-2-glycoprotein I.
ACCESSION X60065 REGION: 1..9
VERSION X60065.1 GI:5
KEYWORDS beta-2 glycoprotein I.
SOURCE Bos taurus (cattle)
ORGANISM Bos taurus
Eukaryota; Metazoa; Chordata; Craniata; Vertebrata; Euteleostomi;
Mammalia; Eutheria; Laurasiatheria; Cetartiodactyla; Ruminantia;
Pecora; Bovidae; Bovinae; Bos.
REFERENCE 1
AUTHORS Bendixen,E., Halkier,T., Magnusson,S., Sottrup-Jensen,L. and
Kristensen,T.
TITLE Complete primary structure of bovine beta 2-glycoprotein I:
localization of the disulfide bridges
JOURNAL Biochemistry 31 (14), 3611-3617 (1992)
PUBMED 1567819
REFERENCE 2 (bases 1 to 9)
AUTHORS Kristensen,T.
TITLE Direct Submission
JOURNAL Submitted (11-JUN-1991) T. Kristensen, Dept of Mol Biology,
University of Aarhus, C F Mollers Alle 130, DK-8000 Aarhus C,
DENMARK
FEATURES Location/Qualifiers
source 1..9
/organism="Bos taurus"
/mol_type="mRNA"
/db_xref="taxon:9913"
/clone="pBB2I"
/tissue_type="liver"
gene <1..>9
/gene="beta-2-gpI"
CDS <1..>9
/gene="beta-2-gpI"
/codon_start=1
/product="beta-2-glycoprotein I"
/protein_id="CAA42669.1"
/db_xref="GI:6"
/db_xref="GOA:P17690"
/db_xref="UniProtKB/Swiss-Prot:P17690"
/translation="PALVLLLGFLCHVAIAGRTCPKPDELPFSTVVPLKRTYEPGEQI
VFSCQPGYVSRGGIRRFTCPLTGLWPINTLKCMPRVCPFAGILENGTVRYTTFEYPNT
ISFSCHTGFYLKGASSAKCTEEGKWSPDLPVCAPITCPPPPIPKFASLSVYKPLAGNN
SFYGSKAVFKCLPHHAMFGNDTVTCTEHGNWTQLPECREVRCPFPSRPDNGFVNHPAN
PVLYYKDTATFGCHETYSLDGPEEVECSKFGNWSAQPSCKASCKLSIKRATVIYEGER
VAIQNKFKNGMLHGQKVSFFCKHKEKKCSYTEDAQCIDGTIEIPKCFKEHSSLAFWKT
DASDVKPC"
sig_peptide <1..>9
/gene="beta-2-gpI"
ORIGIN
1 ccagcgctc
//
More information about the Bioperl-l
mailing list