Dave,

The fact that the Tophat wrapper was not setting the genome based on the alignment genome was actually a bug, but I just fixed it in changeset 5570:0c1251f25c6b.

Let us know if you have further questions.

Regards,
Kelly


On May 17, 2011, at 10:57 AM, Dave Walton wrote:

I'd like to get a better understanding of the point of the database/ build attribute, and pose the question of when is the appropriate time to have it
set?

In our case at the Jackson Laboratory, the most common build is NCBI37/MM9.

However, the feeling of many folks here, is that this should not be set on
our fastq files.  The only place we really run into trouble is with
cufflinks. If you haven't set the db when you get to cufflinks you'll get
an error.

Our suggestion is that there should be one of two options:

1) Tophat has the ability to set the database of the output files based on
the genome that was selected for alignment.

2) There should be a module that can be plugged into a workflow that would set the database of the file prior to passing the file to cufflinks (or any
other two that requires the database attribute to be set).

We are curious if anyone else is running into this issue, and how it is
being solved.

We're thinking about hacking the Tophat wrapper, but I wanted to check with
others before I did this.

Thanks,

Dave


___________________________________________________________
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/

___________________________________________________________
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/

Reply via email to