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
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 -
-
by seqadmin
Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...-
Channel: Articles
03-22-2024, 06:39 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
27 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
31 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
27 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
52 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Comment