[Bioperl-l] How to merge mulitple genbank records into one record
hwang at uga.edu
Mon Apr 24 20:59:40 EDT 2006
Thanks a lot for pointing it out. This is exactly what I need and it
works fine for me. Again, your helps are highly appreciated!
Roy Chaudhuri wrote:
>> The locations of the features refer to the individual 1000000 bp
>> sub-sequences. For example, in the second genbank record
>> 'scaffold:FUGU4:scaffold_1:1000001:2000000:1', the location of a gene
>> is 1760..4580. It is supposed to be 1001760..1004580 on the
> Hi Haiming,
> a couple of months ago I wrote a cat function for Bio::Seq objects that
> should do what you want-
> I think Heikki commited it to the Bio::SeqUtils version in cvs. The
> function will automatically adjust the coordinates so should work for
> your situation. The syntax is:
> Bio::SeqUtils->cat($firstseq, @otherseqs)
> Hope this helps.
> Dr. Roy Chaudhuri
> Bioinformatics Research Fellow
> Division of Immunity and Infection
> University of Birmingham, U.K.
More information about the Bioperl-l