[Bioperl-l] Bio::Graphics::Glyph::translation how to show the correct sequences
Scott Cain
scott at scottcain.net
Mon Mar 24 14:46:38 UTC 2014
Hi Christine,
Can you provide the sample data and complete script? Nothing looks
particularly wrong, but I can't tell from you image that anything is wrong
either. Presumably there is something in the data making you think it is
wrong though.
Scott
On Tue, Mar 18, 2014 at 2:15 AM, Christine Hoogland <
christine.hoogland at gmail.com> wrote:
> Hi,
>
> I'm trying to display the frame translation sequences into images built
> with Bio::Graphics showing exons, introns, CDS... like the example shown in
> figure 13a and 13b on this page:
> http://chlamycyc.mpimp-golm.mpg.de/gbrowse/tutorial/dbgff/tutorial.html
>
> Currently sequences shown are always the same, whatever the CDS (see image
> attached), although sequences printed in the console for debugging are fine.
>
> my @CDS = sort {$a->start <=> $b->start} $mRNA->CDS;
> my $cdsft=Bio::Graphics::Feature->new();
> for my $cds (@CDS) {
> $cdsft->add_segment(Bio::Graphics::Feature->new(
> -start=> $cds->start,
> -end => $cds->end,
> -seq => Bio::Seq->new(-seq=> $cds->dna,
> -alphabet => 'dna'),
> -strand => $cds->strand,
> -phase => $cds->phase
> ));
>
> }
>
> $panel->add_track($cdsft,
> -glyph => 'translation',
> -show_sequence => 'true',
> -key => 'Sequence',
> -height => 64,
> -bump => +1,
> -translation => '3frame',
> -label => 1,
> -description => 1,
> );
>
> Is anybody able to see what's wrong with this code?
>
> Thanks a lot!
> Christine
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
--
------------------------------------------------------------------------
Scott Cain, Ph. D. scott at scottcain dot
net
GMOD Coordinator (http://gmod.org/) 216-392-3087
Ontario Institute for Cancer Research
More information about the Bioperl-l
mailing list