Could you be a little more clear in what you are trying to accomplish here? You are trying to identify contigs from RNA-Seq data after you used Trinity and/or Oases to align to genomic regions? And you are trying to quantify their numbers based on hits to a reference?
I'm trying to align the contigs back to the reference genome, and align the reads onto the contigs, and then use the alignment information to estimate the expression level of each contig.
<p>I use trinity/Oases for de novo transcriptome assembly. In my pipeline, I remove exact duplicate reads(forward and reverse strands) because I believe that duplicates don't …
<p>Hi</p> <p>Are there existing tools for visualizing RNA-Seq splicing graphs from de novo assembly without a reference genome?</p> <p>What I'd like is something that can …
<p>Hi,</p> <p>Is there publicly available strand-specific RNA-Seq data for the MAQC samples? </p> <p>It seems that there are people who have done this (<a href='http://www.nugeninc.com/nugen/index.cfm/products/cs/ngs/encore-complete-rna-seq-library-systems/'>http://www.nugeninc.com/nugen/index.cfm/products/cs/ngs/encore-complete-rna-seq-library-systems/</a>), …
<p>Hi,</p> <p>In sequencing, the asymptotic numbers of gaps and contigs have been give the Lander-Waterman model. But the results from the Lander-Waterman model is only …
<p>After recently assembling a transcriptome, I took the reads used to generate the assembly and aligned them to the assembled contigs using <a href='http://bio-bwa.sourceforge.net/'>BWA</a>. The …
<p>Hi,</p> <p>Is there any publicly available documentation about Velvet's source code? There are a few papers (<a href="http://www.cs.ucr.edu/~jiang/nsf/olga-wabi.pdf">Separating Metagenomic Short Reads into Genomes via Clustering</a>, …
Could you be a little more clear in what you are trying to accomplish here? You are trying to identify contigs from RNA-Seq data after you used Trinity and/or Oases to align to genomic regions? And you are trying to quantify their numbers based on hits to a reference?
I'm trying to align the contigs back to the reference genome, and align the reads onto the contigs, and then use the alignment information to estimate the expression level of each contig.
Why aren't you just aligning the reads directly to the genome, and counting the number of reads that align to each gene?