[Bioperl-l] Parsing join(x1..x2,x3..x4)!

Johan Geijer Johan.Geijer@cgb.ki.se
Wed, 23 Oct 2002 10:57:29 +0200


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鋑 37
>  171 77 Stockholm, SWEDEN
>  E-mail: Johan.Geijer@cgb.ki.se
> #####################################
> _________________________________________