SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
quantification of transcript level bogdan RNA Sequencing 6 06-19-2018 08:04 AM
Velvet and Oases - choice of k value? MikeyG Bioinformatics 5 05-08-2014 09:48 PM
de novo assembly using Trinity versus Velvet-Oases Nol De novo discovery 8 10-26-2013 11:56 AM
Velvet-Oases Memory Problem dacotahm Bioinformatics 0 10-10-2011 09:08 AM
gene and isoform expression level estimation with Cufflinks Jane M RNA Sequencing 1 08-31-2011 11:47 PM

Reply
 
Thread Tools
Old 07-05-2010, 10:10 PM   #1
simon.jarman@aad.gov.au
Junior Member
 
Location: Hobart, Australia

Join Date: May 2009
Posts: 5
Default Velvet/Oases transcript expression level?

Hi,
I am in the process of de novo transcriptome assembly from a non-model eukaryote. I have 17 million 65 bp, paired-end reads that I have assembled with Velvet and then Oases. This gives nice results in that I am getting full-length contigs that appear reasonable when compared to orthologues in related model eukaryotes. However, I am really after relative transcript expression levels. Does anyone know if there as a way that I can get information on that from this software combination?
simon.jarman@aad.gov.au is offline   Reply With Quote
Old 08-10-2010, 03:03 AM   #2
konrad98
Member
 
Location: Exeter, UK

Join Date: Jan 2009
Posts: 17
Default

Hi Simon,

I'm not by any means an expert on this, but you probably have two options:

1. Convert the kmer coverage in the stats.txt file to base-coverage and use that as a proxy for expression level

2. Align the reads to the assembled transcripts using something like Bowtie and then extract the coverage of each transcript using SAMTools

Hope that helps!
konrad98 is offline   Reply With Quote
Old 08-11-2010, 05:33 AM   #3
choy
Member
 
Location: Michigan

Join Date: Jul 2009
Posts: 14
Default Velvet question

How can one extract all of the reads that were used in the assembly, in order to do the re-mapping step you described?
choy is offline   Reply With Quote
Old 08-11-2010, 09:47 AM   #4
konrad98
Member
 
Location: Exeter, UK

Join Date: Jan 2009
Posts: 17
Default

Under the contrib/extractContigReads directory in Velvet you'll find a script to do this
konrad98 is offline   Reply With Quote
Reply

Tags
expresison, oases, quantification, rna seq, velvet

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 02:06 AM.


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