SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
snpeff canonical transcript annotations: multiple annotations problem Nilaksha Bioinformatics 1 03-21-2017 10:16 AM
C elegans known SNPs id0 Bioinformatics 1 06-10-2014 05:10 AM
pdb files for c. elegans bio_informatics Literature Watch 1 03-18-2014 04:03 AM
C.elegans rRNA reduction kits slambert Sample Prep / Library Generation 1 09-30-2013 04:47 PM
RNA-Seq: The Landscape of C. elegans 3'UTRs. Newsbot! Literature Watch 0 06-05-2010 03:00 AM

Reply
 
Thread Tools
Old 08-10-2015, 03:29 PM   #1
virg4l
Junior Member
 
Location: Midwest

Join Date: May 2015
Posts: 7
Default Intron Annotations C. Elegans

I'm working with a mRNA-seq data set where I expect to find non-WT expression that includes bits of introns in mature transcripts. I'd like to be able to use HTSeq-count to add up hits in introns, but my Ensembl GTF (from Illumina's C. elegans iGenomes page) doesn't have introns as a class of annotation.

Does anyone have suggestions about where I can find such a file (or how to generate one)?
virg4l is offline   Reply With Quote
Old 08-11-2015, 01:46 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,479
Default

No one bothers to put introns in annotation files, since the exons are already there and the introns are just between them. One of the more convenient methods to get bounds of introns is to just make/load the TxDb for your organism in R and use the intronsByTranscript() function in GenomicRanges. Note that you'll likely want to post-process that a bit (e.g., make disjoint intronic bins, find which intersect exons and then remove them).
dpryan is offline   Reply With Quote
Reply

Tags
annotation, c. elegans, introns, rna-seq

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 10:27 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO