SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   aligning pre-assembled transcriptomes to genomes (http://seqanswers.com/forums/showthread.php?t=50224)

skeffington 02-11-2015 01:32 AM

aligning pre-assembled transcriptomes to genomes
 
Hi,

I have been given a pre-assembled transcriptome for my organism, which was generated using de novo assembly of RNAseq reads, assembly or reads against a reference genome and by assembly of ESTs. Each method was run using several different assemblers and then the results were clustered to remove redundancy and further processed according to some criteria about reading frames to produce a final transcript set.

I now wish to align this transcript set against the reference genome, so I can view it together with the genome, for example using the IGB. What is the best tool to do this? I presume TopHat is sub-optimal, because it's designed for short reads, not full length transcripts? In addition, my transcript set is in FASTa format - and I think I'm right in saying that you can only run TopHat on FASTq files. I was looking at a tool called LAST - does anyone know if this is suitable?

Please be clear and simple in any answers - I'm an experimentalist learning how to analyse his own sequencing data!

Any help would be greatly appreciated.
Alastair

sarvidsson 02-11-2015 02:32 AM

I've had good results with GMAP (http://research-pub.gene.com/gmap/) for such tasks.
Otherwise BLAT (https://genome.ucsc.edu/FAQ/FAQblat.html) is a quite popular choice.

mastal 02-11-2015 03:28 AM

I would also recommend BLAT.


All times are GMT -8. The time now is 08:33 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.