[Bioperl-l] subsection of genbank file

Chris Fields cjfields at illinois.edu
Tue Sep 22 17:12:10 UTC 2009


That should be out in the latest alpha on CPAN as well (the final  
1.6.1 should be out this week).

chris

On Sep 22, 2009, at 11:44 AM, Roy Chaudhuri wrote:

> Hi Liam,
>
> My mistake, it looks like the bug had already been reported and  
> fixed, which means I get to go home earlier. I've marked your bug as  
> a duplicate of bug 2810.
>
> You can get the patched version by installing bioperl-live (just  
> downloading the bioperl-live SeqUtils.pm and putting it in the  
> correct place on your system would probably also work).
>
> Cheers.
> Roy.
>
> Roy Chaudhuri wrote:
>> Hi Liam,
>> Yes, that is a bug - I think it is to do with the Feature  
>> Annotation rollback from 1.6, it works fine with 1.5.2. Looks like  
>> the tests I wrote don't check for the presence of tags, just the  
>> coordinates of the feature, so this hasn't been picked up. Submit  
>> it to Bugzilla, and I'll take a look when I get a chance.
>> Cheers.
>> Roy.
>> Liam Elbourne wrote:
>>> Hi Roy,
>>>
>>> Thanks for that, works well, but there are no _gsf_tag_hash  
>>> values? I'm particularly interested in the locus id, obviously the  
>>> translation could be problematic if the whole gene is not included  
>>> after truncation, but things like the note, product, protein_id  
>>> would be good. I had a look at the code for the method and  
>>> couldn't see any obvious why those values didn't make it across.  
>>> Should I submit this as a bug, or is there something I'm missing?
>>>
>>>
>>> Regards,
>>> Liam.
>>>
>>>
>>>
>>> On 18/09/2009, at 8:43 PM, Roy Chaudhuri wrote:
>>>
>>>> Hi Liam,
>>>>
>>>> I just discovered your message, which has not yet been replied  
>>>> to. What you require has been discussed in a recent thread:
>>>> http://bioperl.org/pipermail/bioperl-l/2009-August/031071.html
>>>>
>>>> Try using trunc_with_features from Bio::SeqUtils:
>>>>
>>>> my $sub_seqobj=Bio::SeqUtils->trunc_with_features($seqobj, 300,  
>>>> 2000);
>>>> Cheers.
>>>> Roy.
>>>>
>>>> Liam Elbourne wrote:
>>>>> Hi All,
>>>>> Is there a method or methodology that will produce a fully  
>>>>> fledged Seq  object with all the associated metadata given a  
>>>>> start and end  position? To clarify, I create a sequence object  
>>>>> from a genbank file:
>>>>> ****
>>>>> my $io  = Bio::Seqio->new(as per usual);
>>>>> my $seqobj = $io->next_seq();
>>>>> ****
>>>>> I now want:
>>>>> my $sub_seqobj = $seqobj between 300 and 2000
>>>>> where $sub_seqobj is a Seq object (which I appreciate is an   
>>>>> 'aggregate' of objects) too. The "trunc" method only returns a   
>>>>> PrimarySeq object which lacks all the annotation etc. I've  
>>>>> previously  done this task by iterating through feature by  
>>>>> feature and parsing out  what I needed, but thought there might  
>>>>> be a more elegant approach...
>>>>> Regards,
>>>>> Liam Elbourne.
>>>> -- 
>>>> Dr. Roy Chaudhuri
>>>> Department of Veterinary Medicine
>>>> University of Cambridge, U.K.
>>>> _______________________________________________
>>>> Bioperl-l mailing list
>>>> Bioperl-l at lists.open-bio.org
>>>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>>>> ac.uk ([131.111.51.215]:49455)
>>>> by ppsw-7.c
>>> ______________________________
>>>
>>> Dr Liam Elbourne
>>> Research Fellow (Bioinformatics)
>>> Paulsen Laboratory
>>> Macquarie University
>>> Sydney
>>> Australia.
>>>
>>> http://www2.oxfam.org.au/trailwalker/Sydney/team/228
>>>
>>>
>>>
>> _______________________________________________
>> Bioperl-l mailing list
>> Bioperl-l at lists.open-bio.org
>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l




More information about the Bioperl-l mailing list