SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Exon-exon junctions kentsis Bioinformatics 2 12-18-2011 06:32 AM
count reads or count base-pairs yuelics Introductions 3 07-29-2011 05:41 AM
Quantification: count reads or count base pairs? yuelics Bioinformatics 0 07-27-2011 04:48 AM
Finding Exon-Intron Junctions without a reference genome brachysclereid Bioinformatics 3 05-22-2011 06:21 AM
Alignment at exon-exon junctions Boel RNA Sequencing 2 12-09-2010 11:12 AM

Reply
 
Thread Tools
Old 02-22-2011, 03:00 PM   #1
suninsky
Junior Member
 
Location: North Carolina

Join Date: Feb 2011
Posts: 3
Default count reads across exon junctions

Suppose I have a bam file of paired-end reads and a gff file for exon locations. Is there an easy way to count the number of split reads (which can be identified by CIGAR string) and the number of paired-end reads that overlap any exon junctions? Assume we do not have any isoform information here...

Thanks a lot.
suninsky is offline   Reply With Quote
Old 02-22-2011, 07:18 PM   #2
ketan_bnf
Member
 
Location: India

Join Date: Oct 2010
Posts: 59
Default

I am also finding a answer to this question, i am working on single end reads.

If somebody know about this please help.
ketan_bnf is offline   Reply With Quote
Old 02-22-2011, 07:48 PM   #3
apfejes
Senior Member
 
Location: Oakland, California

Join Date: Feb 2008
Posts: 236
Default

I think the question is not necessarily the right one. How do you know you've got your exon-junction-spanning reads spanning exons? I believe we usually align all of our transcriptome reads to reference genome + exon junction template, from which we can then reposition our exon-spanning reads.

If you haven't done that type of alignment, you may not even have your reads aligned correctly - but that would be incredibly dependent upon your aligner.
__________________
The more you know, the more you know you don't know. —Aristotle
apfejes is offline   Reply With Quote
Old 02-22-2011, 09:09 PM   #4
ketan_bnf
Member
 
Location: India

Join Date: Oct 2010
Posts: 59
Default

Quote:
Originally Posted by apfejes View Post
I think the question is not necessarily the right one. How do you know you've got your exon-junction-spanning reads spanning exons? I believe we usually align all of our transcriptome reads to reference genome + exon junction template, from which we can then reposition our exon-spanning reads.

If you haven't done that type of alignment, you may not even have your reads aligned correctly - but that would be incredibly dependent upon your aligner.
Sorry suninsky for posting in your thread,

I am having 454 exome reads aligned to ref seq using bwa to sam file.
Is it possible to find reads spanning exon region?
ketan_bnf is offline   Reply With Quote
Old 02-23-2011, 07:34 AM   #5
apfejes
Senior Member
 
Location: Oakland, California

Join Date: Feb 2008
Posts: 236
Default

Exome reads are usually aligned against the genome because they are captured from DNA, so, for the most part, you will not have exon-exon junctions in your data set. I'm not quite clear what it is you're asking.

You're probably best off trying to align as much as you cat to the genome and do your analysis on that set. If you're looking for reads spanning exon fusions, they'll probably be among the unaligned set, or if you want fusion proteins, you could use the paired end data to look for them.
__________________
The more you know, the more you know you don't know. —Aristotle
apfejes is offline   Reply With Quote
Reply

Tags
exon-junction mapping, rna-seq, split reads

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 04:32 PM.


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