[Bioperl-l] Parsing join(x1..x2,x3..x4)!
Hilmar Lapp
hlapp@gnf.org
Wed, 23 Oct 2002 12:22:20 -0700
Johan, can you provide the sequence in question (e.g. by accession),
and the output your piece of code produces?
Also, I presume you're using version 1.0.2 or higher. If not, please
upgrade and let us know if the problem persists.
-hilmar
On Wednesday, October 23, 2002, at 01:57 AM, Johan Geijer wrote:
> Hello!
>
> I'm having problems to parse the join shown below using the code shown
> below. It's possible to collect both of the starts and stops but it's
> not possible to collect the note.. Does anyone know why.. Solution??
>
>
> Piece of Genbak record I want to extract....
>
> misc_feature join(1511..11036,13779..20629)
> /note="regions of repeated incomplete
> cruzipain ORFs;
> probably pseudogenes"
>
>
>
> A piece of the code I'm using...
>
> use Bio::Seq;
> use Bio::SeqIO;
> use Bio::SeqFeatureI;
>
> foreach $feat ( $seq->all_SeqFeatures() )
> {
> if($feat->primary_tag eq "misc_feature")
> {
> print FASTAOUT " ", $feat->start," ",$feat->end ,"
> ",$feat->strand," ","\n";
> if($feat->has_tag('note'))
> {
> @note = $feat->each_tag_value('note') ;
> print FASTAOUT " ", $note[0], " ";
> }
> }
> }
>
> Regards
>
> Johan
>
> #####################################
>> Johan Geijer
>> Center for Genomics and Bioinformatics
>> Karolinska Institutet
>> Berzelius väg 37
>> 171 77 Stockholm, SWEDEN
>> E-mail: Johan.Geijer@cgb.ki.se
>> #####################################
>> _________________________________________
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l@bioperl.org
> http://bioperl.org/mailman/listinfo/bioperl-l
>
--
-------------------------------------------------------------
Hilmar Lapp email: lapp at gnf.org
GNF, San Diego, Ca. 92121 phone: +1-858-812-1757
-------------------------------------------------------------