Has anyone else noticed that the GTF files in the Illumina iGenomes collection don't conform to the GTF specification ? The attribute list must begin with gene_id then transcript_id whereas in the iGenomes GTFs it is gene_id then gene_name then either p_id or transcript_id. It would be nice if those extra attributes were at the end of the list so they worked with more programs.
Unconfigured Ad
Collapse
X
-
Yeah, well Illumina has always been known about going off on their own. In any case I question a program that can not parse an arbitrary field. If the order of the attribute list *must* be in a given order then the specification should have just made room for two more exact columns instead of relying on programs to do the parsing.
-
Comment
-
Even though it's described as being ~5G, the actual download is ~9G.Originally posted by fongchun View Post
Do you also get 0 FPKM when you use the GRCh37 from the Tophat site?
Comment
-
Yes. I tried both the iGenome GRCh37 packages downloaded from Tophat and Cufflinks and both genes.gtf file gave me 0 FPKM for all genes when I just used cufflinks with the -G parameter.
It seems that maybe the current archive version (archive-2011-08-30-22-38-04) doesn't work? Anyone else have problems?
Comment
-
It took me a week to drill back through my pipeline and work out these files were messing it up. In the end I got a fresh genome and GTF from ensembl, then used custom scripts to build all the various files req'd for picard, GATK etc. Its now working like a dream.
Lesson - if it looks too good to be true...@sidderb
Comment
-
I actually ended up somewhat solving the issue. It had to do with my .bam files having have my chromosome number appended with "chr". While the GTF files, didn't have the "chr". Thus when it tried to quantify the expression it was as if there were absolutely no alignments to the genome. Was frustrating to figure this out.Originally posted by sidderb View PostIt took me a week to drill back through my pipeline and work out these files were messing it up. In the end I got a fresh genome and GTF from ensembl, then used custom scripts to build all the various files req'd for picard, GATK etc. Its now working like a dream.
Lesson - if it looks too good to be true...
Lesson - check to see if "chr" is appended to the chromosome number or not when dealing with datasets.
Comment
Latest Articles
Collapse
-
by SEQadmin2
Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?
Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating...-
Channel: Articles
10-14-2026, 09:55 AM -
-
by SEQadmin2
Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.
We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing...-
Channel: Articles
-
ad_right_rmr
Collapse
News
Collapse
| Topics | Statistics | Last Post | ||
|---|---|---|---|---|
|
Started by SEQadmin2, Yesterday, 11:13 AM
|
0 responses
21 views
0 reactions
|
Last Post
by SEQadmin2
Yesterday, 11:13 AM
|
||
|
Started by SEQadmin2, Yesterday, 10:48 AM
|
0 responses
18 views
0 reactions
|
Last Post
by SEQadmin2
Yesterday, 10:48 AM
|
||
|
Started by SEQadmin2, 10-15-2026, 10:52 AM
|
0 responses
36 views
0 reactions
|
Last Post
by SEQadmin2
10-15-2026, 10:52 AM
|
||
|
Started by SEQadmin2, 10-14-2026, 05:39 AM
|
0 responses
42 views
0 reactions
|
Last Post
by SEQadmin2
10-14-2026, 05:39 AM
|
Comment