A (sort-of) workaround would be to count ambiguous reads for both genes:
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
This seems to be an interesting edge-case. So, TUBB3 and MC1R can form fusion transcripts and fusion proteins. So ENSG00000198211 isn't a mistake and htseq-count isn't doing anything incorrect. However, cases like this are really annoying for read counting. The simplest solution would be to just remove the entries for ENSG00000198211 in your GTF file. The actual correct solution would require a different annotation format that can account for non-novel fusion events of genes that are properly annotated and analyzed as separate in many cases. Such a format will probably get created someday, but I imagine it'll take a while.
Leave a comment:
-
Originally posted by dpryan View PostJust use the -o option to track how the reads are being categorized. Then you can simply look and see why they're not being counted.
Thank you for your reply. I have checked the results, It showed htseq classified most of the reads of this region into "ambiguous" condition, and they belong to ENSG00000198211.8 + ENSG00000258947.3. However, i checked these two genes. Both of them are TUBB3 gene. One is Chromosome 16: 89,921,392-89,938,761 and the other is Chromosome 16: 89,919,165-89,936,092. I understand this problem is not caused by the Htseq software. However, if we discard these reads, it will have a huge impact on expression of this TUBB3 gene (Htseq count 7, but there are dozens).
Do you have any ideas about this?
Thanks
Leave a comment:
-
Just use the -o option to track how the reads are being categorized. Then you can simply look and see why they're not being counted.
Leave a comment:
-
Htseq question
Hi, all,
i just got a result from ht-seq. It showed that my interesting gene has 7 counts in the alignment. However, from the view of IGV, i could easily identify much more counts than 7 on this gene. My alignment is from STAR, and i used more stringent parameters to control the multiple alignment, which means there should not be any multiple aligned reads in the output. I am really confused about this.
Any suggestion?Attached FilesTags: None
Latest Articles
Collapse
-
by seqadmin
The first FDA-approved CRISPR-based therapy marked the transition of therapeutic gene editing from a dream to reality1. CRISPR technologies have streamlined gene editing, and CRISPR screens have become an important approach for identifying genes involved in disease processes2. This technique introduces targeted mutations across numerous genes, enabling large-scale identification of gene functions, interactions, and pathways3. Identifying the full range...-
Channel: Articles
08-27-2024, 04:44 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 08:02 AM
|
0 responses
10 views
0 likes
|
Last Post
by seqadmin
Yesterday, 08:02 AM
|
||
Started by seqadmin, 09-03-2024, 08:30 AM
|
0 responses
13 views
0 likes
|
Last Post
by seqadmin
09-03-2024, 08:30 AM
|
||
Started by seqadmin, 08-27-2024, 04:40 AM
|
0 responses
21 views
0 likes
|
Last Post
by seqadmin
08-27-2024, 04:40 AM
|
||
New Single-Molecule Sequencing Platform Introduces Advanced Features for High-Throughput Genomics
by seqadmin
Started by seqadmin, 08-22-2024, 05:00 AM
|
0 responses
358 views
0 likes
|
Last Post
by seqadmin
08-22-2024, 05:00 AM
|
Leave a comment: