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

Albert Vilella avilella at ub.edu
Thu Dec 22 05:55:27 EST 2005


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.



More information about the Bioperl-l mailing list