I haven't yet worked out what you are doing with the @two array but here is what I have working so far.

#!perl use strict; use warnings; use autodie; use Bio::SeqIO; STDOUT->autoflush(1); my $gbffpath = 'your path'; opendir DIR, $gbffpath; my @files = grep {/\.gbff$/ } readdir DIR; closedir DIR; foreach my $gbfffile (@files){ my $in = Bio::SeqIO->new(-file=>$gbffpath.'/'.$gbfffile, -format=>"genbank"); my $obj = $in->next_seq(); my @feat = $obj->get_SeqFeatures(); my $start = 0; my $length = 9; foreach my $n (0..$#feat){ my $pt = $feat[$n]->primary_tag(); $start += $feat[$n]->start(); my $end = $start + $length; my $subseq = $obj->subseq($start,$end); printf "%2d %-10s %2d %2d %s\n",$n,$pt,$start,$end,$subseq; $start = $end; } }
poj

In reply to Re: Trying to get intergenic sequences using BioPerl by poj
in thread Trying to get intergenic sequences using BioPerl by Sosi

Title:
Use:  <p> text here (a paragraph) </p>
and:  <code> code here </code>
to format your post, it's "PerlMonks-approved HTML":



  • Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
  • Titles consisting of a single word are discouraged, and in most cases are disallowed outright.
  • Read Where should I post X? if you're not absolutely sure you're posting in the right place.
  • Please read these before you post! —
  • Posts may use any of the Perl Monks Approved HTML tags:
    a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
  • You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)
            For:     Use:
    & &amp;
    < &lt;
    > &gt;
    [ &#91;
    ] &#93;
  • Link using PerlMonks shortcuts! What shortcuts can I use for linking?
  • See Writeup Formatting Tips and other pages linked from there for more info.