SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Reference guided assembly StopCodon De novo discovery 5 06-25-2013 10:49 AM
Best approach- genome de novo assembly with reference WaltL Bioinformatics 2 01-08-2013 12:05 PM
bacterial consensus sequence reference guided assembly Kerensa Bioinformatics 0 01-10-2011 09:38 PM

Reply
 
Thread Tools
Old 04-03-2014, 07:09 PM   #1
IdoBar
Junior Member
 
Location: Sunshine Coast, QLD, Australia

Join Date: Mar 2014
Posts: 1
Default Compare de-novo transcriptome assembly to genome reference guided assembly

Hi everyone,

I hope my question is not a total repeat of previous ones, but I did have a read through some previous threads and couldn't find a clear answer.

I have some RNA-seq data of a non-model species (fish), from a single sample that was sent for ~5GB sequencing and bioinformatics analysis at BGI.
The data was mapped and assembled to a transcriptome using SOAPdenovo (very limited one, I know because of lack of repetitions and coverage, but it was a preliminary trial to have a basic idea what we're looking at).

That was roughly a year ago. Currently, a draft genome of a closely related species was published (contigs and scaffolds, not annotated).
My question is, how can I use this reference genome to improve my transcriptome assembly, or to re-assemble it using that genome as a reference and compare the results to the de-novo assembly I already got.
I would also like to use the genome as a reference to detect splice junctions and focus on alternative splicing of specific genes.

I have tried mapping the reads to the published genome reference using STAR, but I'm not sure how to analyse the resulting BAM files and how to use them downstream for assembly or comparison with the de-novo assembly that I got

Last, I would like to re-annotate the new assembly that I will get, what would be the best tool for it?
I would prefer to use open source programs, and I have access to a strong linux server. (32 cpu, 1TB ram)

I am pretty new to NGS bioinformatics, but I'm eager to learn and would appreciate your feedback.

Thanks, Ido
IdoBar is offline   Reply With Quote
Old 04-04-2014, 01:28 AM   #2
TiborNagy
Senior Member
 
Location: Budapest

Join Date: Mar 2010
Posts: 329
Default

Trinity is an all-in-one tool. It can use reference to make a de-novo transcriptome and it has a module to annotate the results:
http://trinityrnaseq.sourceforge.net...d_trinity.html
http://trinotate.sourceforge.net/
TiborNagy is offline   Reply With Quote
Reply

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 11:23 AM.


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