Hi,<br><br>I have fixed the protein glyph so that it automatically translates DNA into protein if it is run on top of a dna database. The same restrictions apply as to the cds glyph -- the phase of the feature's subparts need to be correctly defined, or you will get nonsense. At high power, the protein translation appears, just like the cds glyph (I did this by making protein a subclass of cds).<br>
<br>It is now in CPAN under Bio::Graphics version 1.81.<br><br>Lincoln<br><br><div class="gmail_quote">On Fri, Dec 19, 2008 at 11:57 AM, Scott Cain <span dir="ltr"><<a href="mailto:cain.cshl@gmail.com">cain.cshl@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi Lincoln,<br>
<br>
Before you retire it, consider that people might be using it either<br>
with GBrowse as a protein viewer (like SGD does, though I don't know<br>
if they are using the protein glyph) or in some other application that<br>
uses Bio::Graphics in an application outside of GBrowse for displaying<br>
protein information. This protein glyph is actually pretty nice in<br>
that, when zoomed out, it behaves like the dna glyph and displays a kd<br>
plot, and zoomed in displays the residues.<br>
<font color="#888888"><br>
Scott<br>
</font><div><div></div><div class="Wj3C7c"><br>
<br>
On Fri, Dec 19, 2008 at 11:45 AM, Lincoln Stein <<a href="mailto:lincoln.stein@gmail.com">lincoln.stein@gmail.com</a>> wrote:<br>
> Oh gee, it looks like I should retire the protein glyph. It definitely isn't<br>
> behaving the way I think it should, which is to do the translation itself<br>
> based on the phase and dna sequence.<br>
><br>
> The -draw_protein_target option is only supposed to work if you've got a<br>
> protein defined as the target which aligns to a dna source -- for example, a<br>
> tblastn result.<br>
><br>
> Lincoln<br>
><br>
> On Wed, Dec 17, 2008 at 3:12 PM, Newton Vidal <<a href="mailto:nwvidal@gmail.com">nwvidal@gmail.com</a>> wrote:<br>
>><br>
>> Hi, Scott and Lincoln.<br>
>><br>
>> I have another question, related to my problem, but arising another issue:<br>
>><br>
>> I have configured two different glyphs, one as CDS and the other as<br>
>> protein. When I zoomed in the main Gbrowse screen, the CDS glyph is shown as<br>
>> protein sequence and the protein is shown as nucleotide sequence. Is this<br>
>> behavior expected? Why the CDS, that codes for a protein, but is composed of<br>
>> DNA is shown as protein and the protein glyph is not?<br>
>><br>
>> Representative pictures are attached.<br>
>><br>
>><br>
>> Thanks in advance.<br>
>><br>
>> Newton.<br>
>><br>
>><br>
>> On Tue, Dec 16, 2008 at 8:52 PM, Dave Clements, GMOD Help Desk<br>
>> <<a href="mailto:gmodhelp@googlemail.com">gmodhelp@googlemail.com</a>> wrote:<br>
>>><br>
>>> Hi Newton, Scott, Lincoln,<br>
>>><br>
>>> The draw_protein_target option of the segments glyph did not work as<br>
>>> of October. I don't recall any recent problems (or success either)<br>
>>> with the draw_target option of the transcript glyph, although it<br>
>>> wouldn't surprise me if both options use the same code. (See<br>
>>> <a href="http://thread.gmane.org/gmane.science.biology.gmod.gbrowse/808" target="_blank">http://thread.gmane.org/gmane.science.biology.gmod.gbrowse/808</a> for<br>
>>> more on draw_protein_target problems).<br>
>>><br>
>>> I also loaded the GFF3 file for Tcruzi into a mysql database. I did<br>
>>> not get any load errors (or Apache runtime errors), but nothing<br>
>>> displays for me either. All landmarks are unrecognized. I tried<br>
>>> playing around with the conf file but to no avail.<br>
>>><br>
>>> Dave C.<br>
>>><br>
>>><br>
>>> Was this helpful? Let us know at <a href="http://gmod.org/wiki/Help_Desk_Feedback" target="_blank">http://gmod.org/wiki/Help_Desk_Feedback</a><br>
>>><br>
>>><br>
>>><br>
>>> On Tue, Dec 16, 2008 at 1:02 PM, Scott Cain <<a href="mailto:cain.cshl@gmail.com">cain.cshl@gmail.com</a>> wrote:<br>
>>> > Hi Lincoln and Newton,<br>
>>> ><br>
>>> > I was under the impression that the draw_protein option was<br>
>>> > incompletely implemented, but if it works with other database<br>
>>> > adaptors, perhaps there is a bug in the Chado GBrowse adaptor. Have<br>
>>> > you tried the cds glyph, which should display residues when zoomed in<br>
>>> > far enough? If draw_protein doesn't work but the cds glyph does, it<br>
>>> > indicates a Bio::Graphics problem. If neither work, it is more likely<br>
>>> > a Chado adaptor or data problem.<br>
>>> ><br>
>>> > As for the T. cruzi problem--I don't see anything in the data that<br>
>>> > looks problematic. Where there any error messages while loading or in<br>
>>> > the apache error_log (with a preemptive note: error messages about<br>
>>> > Bio::Graphics::BrowserConfig don't mean anything)?<br>
>>> ><br>
>>> > Scott<br>
>>> ><br>
>>> ><br>
>>> > On Tue, Dec 16, 2008 at 3:35 PM, Lincoln Stein<br>
>>> > <<a href="mailto:lincoln.stein@gmail.com">lincoln.stein@gmail.com</a>> wrote:<br>
>>> >> Hi Newton,<br>
>>> >><br>
>>> >> Have you tried passing the draw_protein option?<br>
>>> >><br>
>>> >> [Polypeptide]<br>
>>> >> feature = polypeptide<br>
>>> >> glyph = transcript<br>
>>> >> draw_protein = 1<br>
>>> >><br>
>>> >> Also check out the "protein" and "translation" glyphs.<br>
>>> >><br>
>>> >> The T. cruzi problem sounds like a chado issue -- I will defer to<br>
>>> >> Scott<br>
>>> >> and/or David on this.<br>
>>> >><br>
>>> >> Best,<br>
>>> >><br>
>>> >> Lincoln<br>
>>> >><br>
>>> >> On Tue, Dec 16, 2008 at 2:01 PM, Newton Vidal <<a href="mailto:nwvidal@gmail.com">nwvidal@gmail.com</a>><br>
>>> >> wrote:<br>
>>> >>><br>
>>> >>> Dear all,<br>
>>> >>><br>
>>> >>> I am starting to work with GMOD and Chado, and I have uploaded<br>
>>> >>> genbank<br>
>>> >>> files from five different protozoan genomes, using bp_genbank2gff3.pl<br>
>>> >>> and<br>
>>> >>> gmod_bulk_load_gff3.pl scripts in our local chado v. 1.0 database.<br>
>>> >>><br>
>>> >>> I have run into two major problems when visualizing data with Gbrowse<br>
>>> >>> v.<br>
>>> >>> 1.69:<br>
>>> >>><br>
>>> >>> · The polypeptide feature shows nucleic acid sequence instead<br>
>>> >>> of<br>
>>> >>> amino acid, when clicked;<br>
>>> >>> · For four of the genomes (Leishmania sp. and Trypanosoma<br>
>>> >>> brucei)<br>
>>> >>> the others features work fine, but for T. cruzi, no data is provided<br>
>>> >>> when I<br>
>>> >>> click its features.<br>
>>> >>><br>
>>> >>> I have edited chado.conf file to work with my database and added the<br>
>>> >>> lines<br>
>>> >>> below:<br>
>>> >>><br>
>>> >>> [Polypeptide]<br>
>>> >>> feature = polypeptide<br>
>>> >>> glyph = transcript<br>
>>> >>> height = 10<br>
>>> >>> key = Polypeptide<br>
>>> >>> bgcolor = green<br>
>>> >>><br>
>>> >>> I'm sending attached T. brucei and T. cruzi's GFF files and gbrowse<br>
>>> >>> conf<br>
>>> >>> file. I have a 80MB dump file and I can share it if needed.<br>
>>> >>><br>
>>> >>> Thanks in advance,<br>
>>> >>><br>
>>> >>> Newton Vidal<br>
>>> >><br>
>>> >><br>
>>> >><br>
>>> >> --<br>
>>> >> Lincoln D. Stein<br>
>>> >><br>
>>> >> Ontario Institute for Cancer Research<br>
>>> >> 101 College St., Suite 800<br>
>>> >> Toronto, ON, Canada M5G0A3<br>
>>> >> 416 673-8514<br>
>>> >> Assistant: Renata Musa <<a href="mailto:Renata.Musa@oicr.on.ca">Renata.Musa@oicr.on.ca</a>><br>
>>> >><br>
>>> ><br>
>>> ><br>
>>> ><br>
>>> > --<br>
>>> ><br>
>>> > ------------------------------------------------------------------------<br>
>>> > Scott Cain, Ph. D. scott at scottcain<br>
>>> > dot net<br>
>>> > GMOD Coordinator (<a href="http://gmod.org/" target="_blank">http://gmod.org/</a>) 216-392-3087<br>
>>> > Ontario Institute for Cancer Research<br>
>>> ><br>
>>> ><br>
>><br>
><br>
><br>
><br>
> --<br>
> Lincoln D. Stein<br>
><br>
> Ontario Institute for Cancer Research<br>
> 101 College St., Suite 800<br>
> Toronto, ON, Canada M5G0A3<br>
> 416 673-8514<br>
> Assistant: Renata Musa <<a href="mailto:Renata.Musa@oicr.on.ca">Renata.Musa@oicr.on.ca</a>><br>
><br>
<br>
<br>
<br>
</div></div>--<br>
<div><div></div><div class="Wj3C7c">------------------------------------------------------------------------<br>
Scott Cain, Ph. D. scott at scottcain dot net<br>
GMOD Coordinator (<a href="http://gmod.org/" target="_blank">http://gmod.org/</a>) 216-392-3087<br>
Ontario Institute for Cancer Research<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Lincoln D. Stein<br><br>Ontario Institute for Cancer Research<br>101 College St., Suite 800<br>Toronto, ON, Canada M5G0A3<br>416 673-8514<br>Assistant: Renata Musa <<a href="mailto:Renata.Musa@oicr.on.ca">Renata.Musa@oicr.on.ca</a>><br>