Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
EBARDenovo - A new RNA-seq do novo assembler for paired-end Illumina data htchu.taiwan RNA Sequencing 2 06-10-2013 12:13 AM
EBARDenovo - A new RNA-seq do novo assembler for paired-end Illumina data htchu.taiwan Illumina/Solexa 9 04-15-2013 11:08 PM
RNA-Seq: De novo transcriptome assembly of RNA-Seq reads with different strategies. Newsbot! Literature Watch 0 01-10-2012 04:00 AM
RNA-Seq: A pipeline for RNA-seq data processing and quality assessment. Newsbot! Literature Watch 0 01-15-2011 02:50 AM
RNA-Seq: MapSplice: Accurate mapping of RNA-seq reads for splice junction discovery. Newsbot! Literature Watch 2 10-14-2010 08:35 AM

Thread Tools
Old 11-21-2011, 07:10 AM   #1
Junior Member
Location: Boston, MA

Join Date: Jan 2011
Posts: 8
Default Mapping paired reads to de novo RNA-seq assemblies for quality assessment?


I'm working with paired end reads (not strand-specific; 100 bp each end with a total fragment size of 400 bp) from Illumina v1.5 data format; 100 million read pairs, but trying out the workflow on 1m, 2m, 5m, and 10m reads. All reads are adapter trimmed and low quality base trimmed.

I'm hoping to get a few suggestions to two questions I have:

1) I am using Trinity to build de novo RNA-seq assemblies. But the mean/median assembly sizes are smaller than I would expect (503 bp/300 bp). Given the data is from mRNA from a vertebrate, what can I do to increase the mean/median contig size?

2) More importantly, as I have paired-end reads that went into the assembly, I'd like to map the pairs back onto the de novo transcripts to assess if the pairs are being used correctly, and to give me some quality assessment for each transcript. What program would folks recommend that would give me some sort of readout? I would expect something like, per transcript, X# pairs map correctly; Y# singles are mapped; Z# are mapped but violate orientation/insert size, etc.

Thanks much,
BobFreemanMA is offline   Reply With Quote
Old 08-02-2012, 11:36 AM   #2
Junior Member
Location: US

Join Date: Aug 2012
Posts: 1

Hi BobFreemanMA

Did you ever get a reply to this?? I am trying to do exactly the same and I do not want to reinvent the wheel.

aprendiz is offline   Reply With Quote

denovo assembly, mapping, paired end reads, rna-seq, trinity

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 03:51 PM.

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