[Bioperl-l] half intergenic 5' and 3' regions of a genbank entry
avilella at ub.edu
Thu Dec 22 05:55:27 EST 2005
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?
More information about the Bioperl-l