Finding Reads That Suffered Splicing
3
1
Entering edit mode
11.1 years ago
Gabrielw ▴ 60

I've aligned 454 RNA-Seq reads using bwa. Does somebody know if a software that identify variants that suffered splicing exists ?

rna splicing • 2.7k views
ADD COMMENT
6
Entering edit mode
11.1 years ago

TopHat is probably not the "best" in terms of sensitivity or specificity. Have a look here for a more complete list. I use GSNAP quite regularly for this task.

ADD COMMENT
1
Entering edit mode

Yes, I agree. TopHat does not do gapped alignment, while 454 has this particular homopolymer errors. TopHat may not work well for long reads, either. Perhaps even blat/gmap will be better.

ADD REPLY
0
Entering edit mode

Good point. Blat/gmap should do well; I have not seen much 454 data, so I don't have a good sense of how well, though.

ADD REPLY
2
Entering edit mode
11.1 years ago
Rm 8.2k

For RNA-seq data try Tophat; TopHat is a fast splice junction mapper for RNA-Seq reads.

ADD COMMENT
2
Entering edit mode
11.1 years ago
Cjt ▴ 370

I recommend to have a look at PALMA, described here: http://www.fml.tuebingen.mpg.de/raetsch/suppl/palma

Quality estimations provided on their website

ADD COMMENT
2
Entering edit mode

I used to read this paper. This figure is greatly exaggerating the problem with mapping ESTs. The authors are fabricating/selecting ESTs with extremely short exons and claim the alignment is correct only if all splice sites are correct. For most ESTs and when measuring the exon-level accuracy, the difference between PALMA and BLAT should be negligible (a couple percent at most). PALMA relies on BLAST, which is very inefficient for human alignment.

ADD REPLY
0
Entering edit mode

I used to read this paper. This figure is greatly exaggerating the problem in mapping ESTs. The authors are fabricating/selecting ESTs with extremely short exons and claim the alignment is correct only if all splice sites are correct. For most ESTs and when measuring the exon-level accuracy, the difference between PALMA and BLAT should be negligible. PALMA relies on BLAST, which is very inefficient for human alignment.

ADD REPLY

Login before adding your answer.

Traffic: 1364 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6