[Bioperl-l] half intergenic 5' and 3' regions of a genbank entry

[Albert Vilella]

Hi all,

I'm trying to obtain each of the genes of a genbank genome prokaryote
file with half of the 5' and 3' intergenic sequences with respect to
the corresponding bounding genes.

So if gene 1 has gene 0 at 1200bp 5' and gene2 at 444bp 3', the script
will output the fasta of gene1 with 600bp 5', the gene, and 222bp 3'.

Any script or stub I can reuse? Any suggestion on how to do that?

Thanks,

    Albert.