SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
normalization of pair-end RNA-seq data with different sequence length yao_licr RNA Sequencing 1 04-21-2015 04:44 AM
Pair-end sequences assembly with Trinity (data mining doubt) GillermoPonz Bioinformatics 0 04-13-2015 02:35 AM
Illumina Nextera Pair-End Sequence Content Bias-Require trimming for DeNovo Assembly? dave1 De novo discovery 5 07-21-2014 12:20 AM
difference between pair end and single end reads assembly jjjscuedu Bioinformatics 0 08-13-2013 06:22 PM
SPRI bead for PCR cleanup in illumina pair end library preparation Berkeley2010 Illumina/Solexa 2 10-20-2010 12:46 PM

Reply
 
Thread Tools
Old 05-28-2015, 08:19 AM   #1
GillermoPonz
Junior Member
 
Location: Madrid

Join Date: Oct 2014
Posts: 6
Default Pair-end sequence data preparation for Trancriptome assembly

Hi everyone,

I'm about to perform a transcriptome assembly with Trinity, but I have a basic question.

I have two RNA-Seq files containing 40 million Pair-End reads from Illumina Hi-Seq. After performing a FastQC Iíve observed that the first 10 bases have a bad result in the "Per Base Sequence Content" plot and I want to trim them (by trimming I mean Trimming by column, because the quality trimming didnít fix the problem) to avoid problems during the assembly. My problem is that I donít know if by trimming the first 10 bases of the 5í end in both files I will mess up the sequence pair mating and completely ruin the downstream analysis.

Could anybody help me with this? I know is a basic question, but Iím a newbie and I havenít been able to get the answer by myself.
GillermoPonz is offline   Reply With Quote
Old 05-28-2015, 08:23 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,698
Default

Do not trim the data at the beginning of the read. That bias is a "known" feature of RNAseq data and it does not affect any downstream analysis.

See post #261 and #263: http://seqanswers.com/forums/showthr...t=4846&page=14
GenoMax is online now   Reply With Quote
Old 05-28-2015, 09:31 AM   #3
GillermoPonz
Junior Member
 
Location: Madrid

Join Date: Oct 2014
Posts: 6
Default

Thanks for your quick answer, this info has been very useful for me.

But just for curiosity, do you know if it would affect the pairing of the sequences (and in consequence the ability of Trinity to use the pair-end information) if I trim them?
GillermoPonz is offline   Reply With Quote
Old 05-28-2015, 09:37 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,698
Default

As long as you use a trimming program (trimmomatic, bbduk etc) that is paired-end aware, pairing of sequences in your data would be maintained.
GenoMax is online now   Reply With Quote
Reply

Tags
assembly, bioinformactics, rna-seq, trimming, 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 08:16 AM.


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