Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Low coverage after alignment Rangika Bioinformatics 0 09-28-2016 07:57 PM
Subread low percentage of assigned fragments in ENCODE data kajot Bioinformatics 2 03-07-2016 05:40 AM
Low overall alignment rate maivantan Bioinformatics 7 03-01-2016 12:38 PM
Wang, 'Low Cost Library Construction...': Low Concentration MAdkisson RNA Sequencing 2 02-19-2014 02:33 PM
Superscript "Stitching" low complexity fragments? rflrob RNA Sequencing 0 01-15-2014 03:59 PM

Thread Tools
Old 03-15-2017, 03:18 PM   #1
Junior Member
Location: Canada

Join Date: Sep 2016
Posts: 4
Default Low GC, fragments, low read alignment

I am new to Seqanswers (Hi!) but so far have not yet found a helpful thread for my particular problem.

I am trying to do two de novo assemblies of two species of sponges using Trinity. It is an RNA-seq study comparing different tissues and life stages within the respective species, with 2x150 bp reads sequenced to a depth of 124x. There are two whole transcriptomes already published. The GC content of my transcriptome assemblies are ~3-7% lower than the other ones, and GC content falls only after assembly (ie. not after read trimming). Furthermore, my sequences are short (eg. N50 800 vs 1800bp, median length 300 vs 800bp, mean length 600vs 1200bp). The nail on the coffin is that very few reads align as proper pairs (vs individually) - ~50%.

As far as I can tell, there isn't anything about the reads themselves that stand out as red flags, and raw read quality was high.

There are many places where this could have gone wrong, but I would guess it is during assembly. Has anyone received this particular profile before (low GC, many fragments)? If so, what did you do to achieve a better assembly?

jmah is offline   Reply With Quote

assembly improvements, gc content, short contigs, 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 04:09 AM.

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