Is there any legitimate reason a tophat bam file would have reads for a gene whose RPKM was not reported in the output from tophat or cufflinks?
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Originally posted by mgogol View PostIs there any legitimate reason a tophat bam file would have reads for a gene whose RPKM was not reported in the output from tophat or cufflinks?
I thought tophat uses FPKM units and not RKPM.
-
If I look at the list of genes whose FPKMs are reported by cufflinks, a certain group of genes is not found in that list at all.
If I load a bam file generated from the tophat sam file into IGV, and browse around looking at those genes, they have reads. So they have reads, but FPKMs are not reported.
I'm just wondering if this an error or if there is some legitimate explanation?
Comment
-
if the reads are too few, they may be filtered by the F/--min-isoform-fraction parameter where TopHat filters out junctions supported by too few alignments.
May be you can set it to 0.
Also, I don't intuitively understand -a/--min-anchor-length <int> parameter, but it may also be limiting your reads that are counted for the FPKM.
Comment
-
sure, if you find out whats going on, let us know. If the reads you saw on IGV don't span junctions, then its probably a parameter in bowtie that is throwing it- off-low quality or reads mapping to multiple places in the genome or something else. Its a good thing you caught it though
Comment
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
59 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
57 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
51 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
55 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Comment