Hi,
In order to have the transcript_id for each sequence extracted from the
cuffmerge .gtf file I had to change the extract_genomic_dna.py by adding
the following lines after line 153:
attributes = gff_util.parse_gff_attributes( feature[8] )
if ( "transcript_id" in attributes ):
name = attributes.get( "transcript_id", None )
This way the variable name gets the transcript_id if it exists.
If it's correct, I would appreciate this modification in future galaxy
distributions.
Thanks!
Adhemar
___________________________________________________________
Please keep all replies on the list by using "reply all"
in your mail client. To manage your subscriptions to this
and other Galaxy lists, please use the interface at:
http://lists.bx.psu.edu/
To search Galaxy mailing lists use the unified search at:
http://galaxyproject.org/search/mailinglists/
