Hi Jen,
thanks.
I am a bit confused: on Galaxy the only human genome listed is hg_g1k_v37.
So when I uploaded the new data from "Get data", under "Genome" I selected hg_g1k_v37
Now, all i want is to get cufflinks with gene names: which genome am I supposed to use? the only one I knew was hg19 from iGenomes...but seems i cannot use it. Do I have to select a genome when I upload raw fastq data? I havent stated doign anything so far, its just raw data
Thanks a lot,
Irene
________________________________________
From: Jennifer Jackson [jen@bx.psu.edu]
Sent: Wednesday, November 27, 2013 10:08 PM
To: Irene Bassano
Cc: galaxy-bugs@bx.psu.edu
Subject: Re: [galaxy-bugs] Galaxy tool error report from bsib@leeds.ac.uk
Hello,
iGenomes covers the UCSC build, and this named "Human Feb. 2009
(GRCh37/hg19) (hg19)" in the full name in the UI. The "hg19" key is the
important part - as the name may be abbreviated in some tools, but this
key will be in all. The genome with the "hg_g1k_v37" key is slightly
different and you will have another genome mismatch problem with the
RNA-seq (and most other) tools if you combine this genome with data from
UCSC (the source of hg19) or the wrong iGenomes file.
"hg_g1k_v37" (source: 1000 genomes via GATK) and "hg19" (source: UCSC)
are just about the same, but the identifiers are different. If you want
to examine the differences, both are available on our rsync server and
can be downloaded and compared. On the "Help -> Support" wiki are links
for reference genomes.
The iGenomes GTF file for hg19 is on the public Main server, if that is
more convenient for you, or should you just want to be sure you have the
right one. Look in Shared Data -> Data Libraries -> iGenomes.
Best,
Jen
Galaxy team
ps. please try to send new questions to one of our lists, thanks!
On 11/27/13 12:51 PM, Irene Bassano wrote:
Hi Jen,
I uploaded some fastq files and selected as Genome from the drop down list ""Homo sapiens b37(hg_g1k_v37).
Is this the same as the genome listed in UCSC "February 2009 (GRCh37/hg19)"?
I am using iGenomes to get the gene names rather than annotation such as NM_00xxxx and I fused the UCSC website choosing the newest genome, February 2009
Thanks,
best,
Irene
________________________________________
From: Jennifer Jackson [jen@bx.psu.edu]
Sent: Monday, November 18, 2013 7:26 PM
To: galaxy-bugs@bx.psu.edu; Irene Bassano
Subject: Re: [galaxy-bugs] Galaxy tool error report from bsib@leeds.ac.uk
Hi again,
The database mismatch is a problem here is well for the same reasons. My
guess is that you intended to run dataset #21 against hg19, and that the
run against hg18 was a tool form input mistake?
Good luck with the next runs,
Jen
Galaxy team
On 11/18/13 5:21 AM, galaxy-bugs@bx.psu.edu wrote:
GALAXY TOOL ERROR REPORT
------------------------
This error report was sent from the Galaxy instance hosted on the server
"usegalaxy.org"
-----------------------------------------------------------------------------
This is in reference to dataset id 7081835 from history id 1686699
-----------------------------------------------------------------------------
You should be able to view the history containing the related history item
37: Cufflinks on data 21 and data 34: assembled transcripts
by logging in as a Galaxy admin user to the Galaxy instance referenced above
and pointing your browser to the following link.
usegalaxy.org/history/view?id=d88c1ef77619eb4f
-----------------------------------------------------------------------------
The user 'bsib@leeds.ac.uk' provided the following information:
Same as before but the reference genome is UCSC MAin on Human:refFlat
-----------------------------------------------------------------------------
job id: 6096851
tool id: toolshed.g2.bx.psu.edu/repos/devteam/cufflinks/cufflinks/0.0.6
job pid or drm id: 136948
-----------------------------------------------------------------------------
job command line:
python /galaxy/main/migrated_tools/toolshed.g2.bx.psu.edu/repos/devteam/cufflinks/b01956f26c36/cufflinks/cufflinks_wrapper.py --input=/galaxy-repl/main/files/007/074/dataset_7074022.dat --assembled-isoforms-output=/galaxy-repl/main/files/007/081/dataset_7081835.dat --num-threads="8" -I 300000 -F 0.1 -j 0.15 -G /galaxy-repl/main/psufiles/004/831/dataset_4831015.dat -N -b --ref_file="None" --dbkey=hg18 --index_dir=/galaxy/main/server/tool-data -u
-----------------------------------------------------------------------------
job stderr:
Error running cufflinks.
return code = 1
Command line:
cufflinks -q --no-update-check -I 300000 -F 0.100000 -j 0.150000 -p 8 -G /galaxy-repl/main/psufiles/004/831/dataset_4831015.dat -u -N -b /galaxy/data/hg18/sam_index/hg18.fa /galaxy-repl/main/files/007/074/dataset_7074022.dat
Error: cannot open reference GTF file /galaxy-repl/main/psufiles/004/831/dataset_4831015.dat for reading
-----------------------------------------------------------------------------
job stdout:
cufflinks v2.1.1
cufflinks -q --no-update-check -I 300000 -F 0.100000 -j 0.150000 -p 8 -G /galaxy-repl/main/psufiles/004/831/dataset_4831015.dat -u -N -b /galaxy/data/hg18/sam_index/hg18.fa
-----------------------------------------------------------------------------
job info:
None
-----------------------------------------------------------------------------
job traceback:
None
-----------------------------------------------------------------------------
(This is an automated message).
--
Jennifer Hillman-Jackson
http://galaxyproject.org
--
Jennifer Hillman-Jackson
http://galaxyproject.org
___________________________________________________________
The Galaxy User list should be used for the discussion of
Galaxy analysis and other features on the public server
at usegalaxy.org. Please keep all replies on the list by
using "reply all" in your mail client. For discussion of
local Galaxy instances and the Galaxy source code, please
use the Galaxy Development list:
http://lists.bx.psu.edu/listinfo/galaxy-dev
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/