[Bioperl-l] SeqIO tag/location question
Tom Keller
kellert at ohsu.edu
Tue Jul 6 21:54:36 UTC 2010
Greetings,
A genbank file I am working with has features:
exon 37167..37302
/gene="SLC12A3"
/gene_synonym="FLJ96318; NCCT; TSC"
/inference="alignment:Splign:1.39.8"
/number=24
exon 39189..39256
/gene="SLC12A3"
/gene_synonym="FLJ96318; NCCT; TSC"
/inference="alignment:Splign:1.39.8"
/number=25
exon 48058..50644
/gene="SLC12A3"
/gene_synonym="FLJ96318; NCCT; TSC"
/inference="alignment:Splign:1.39.8"
/number=26
When I parse it using:
my $db = Bio::DB::RefSeq->new();
my $seq = $db->get_Seq_by_acc($accession);
my @features = $seq->get_SeqFeatures;
for my$feat_obj ( @features ) {
if ($feat_obj->primary_tag eq 'exon') {
print "primary tag: ", $feat_obj->primary_tag, " start: ", $feat_obj->start, " end: ",$feat_obj->end, "\n";
for my $tag ($feat_obj->get_all_tags) {
for my $value ($feat_obj->get_tag_values($tag)) {
print " value: ", $value, "\n";
}
}
}
}
I get the following:
primary tag: exon start: 42167 end: 42302
value: SLC12A3
value: FLJ96318; NCCT; TSC
value: alignment:Splign:1.39.8
value: 24
primary tag: exon start: 44189 end: 44256
value: SLC12A3
value: FLJ96318; NCCT; TSC
value: alignment:Splign:1.39.8
value: 25
primary tag: exon start: 53058 end: 55644
value: SLC12A3
value: FLJ96318; NCCT; TSC
value: alignment:Splign:1.39.8
value: 26
Where does the discrepancy in location (+5000) come from?
thanks for your help,
Tom
MMI DNA Services Core Facility<http://www.ohsu.edu/xd/research/research-cores/dna-analysis/>
4-2442
kellert at ohsu.edu<mailto:kellert at ohsu.edu>
RJH (CROET/BasicScience) 6339b
More information about the Bioperl-l
mailing list