SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Any tool can be fed a list of genome locations and output annotations? gooday23 Bioinformatics 4 04-12-2013 04:19 AM
Whole genome alignment (or the opposite of it?) jomaco Bioinformatics 1 02-14-2013 05:57 AM
Please help: imperfect reference genome/get consensus on genome/read alignment? KAP Bioinformatics 1 08-19-2011 08:14 AM
transferring annotations from reference genome to the resequenced genome mike.t Bioinformatics 1 09-17-2010 06:35 AM
Multi-Genome Alignment for QC... james hadfield Bioinformatics 0 08-17-2010 08:51 AM

Reply
 
Thread Tools
Old 04-12-2013, 12:22 AM   #1
oxydeepu
Member
 
Location: bangalore,india

Join Date: Jul 2011
Posts: 41
Default Visualing the alignment on a genome without annotations

Hi All,

I'm working in a lab which uses planaria as a model organism, which have a contig level of genome assembly. I have RNA-Seq data from the same. Now I am trying to visualize the alignment of the reads on the genome. I also want to visualize the transcripts made from the cufflinks assembly of those RNA-Seq reads. The problem is RNA-Seq alignment is using bowtie and Cufflinks transcripts are using blast. Can anyone help me in visualizing the same. I just want to validate that the transcripts alignment on the genome corresponds to the RNA-Seq alignment. Please help. Thank you in advance.

Best Regards,
Deepak
oxydeepu is offline   Reply With Quote
Old 04-12-2013, 01:16 AM   #2
NicoBxl
not just another member
 
Location: Belgium

Join Date: Aug 2010
Posts: 264
Default

I you have a bam file from an aligner (you used bowtie, but you should use a splice-aware aligner like tophat or STAR) and a gtf (or gff) from cufflinks, you can use IGV to visualize the alignment and check the annotation created by cufflinks.
NicoBxl is offline   Reply With Quote
Old 04-12-2013, 01:26 AM   #3
oxydeepu
Member
 
Location: bangalore,india

Join Date: Jul 2011
Posts: 41
Default

Thank You Nico,

I will definitely try the suggestion. One doubt, does IGV, inputs a genome fasta file or does it have a built in database for genomes? What If I want to put in the EST alignment information, which is a blast?.

Thank you for your time once again.
Best Regards,
Deepak
oxydeepu is offline   Reply With Quote
Old 04-12-2013, 01:28 AM   #4
NicoBxl
not just another member
 
Location: Belgium

Join Date: Aug 2010
Posts: 264
Default

IGV has several genomes in its database, but it accept also fasta files.
NicoBxl is offline   Reply With Quote
Old 04-12-2013, 01:41 AM   #5
oxydeepu
Member
 
Location: bangalore,india

Join Date: Jul 2011
Posts: 41
Default

WHat about the second question?
oxydeepu is offline   Reply With Quote
Old 04-12-2013, 01:42 AM   #6
NicoBxl
not just another member
 
Location: Belgium

Join Date: Aug 2010
Posts: 264
Default

you should convert your blast output to gff. Check here http://www.biostars.org/p/277/
NicoBxl is offline   Reply With Quote
Old 04-12-2013, 01:52 AM   #7
oxydeepu
Member
 
Location: bangalore,india

Join Date: Jul 2011
Posts: 41
Default

So now I can include this gff also to the IGV viewer and visualise. So does this take care of the exon informations, any clue about that?

Thanks,
Deepak
oxydeepu is offline   Reply With Quote
Old 04-12-2013, 02:13 AM   #8
NicoBxl
not just another member
 
Location: Belgium

Join Date: Aug 2010
Posts: 264
Default

yes, check the IGV manual to have more infos : http://www.broadinstitute.org/igv/
NicoBxl 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 01:18 PM.


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