SEQanswers

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
Human Illumina Paired-end RNA-Seq remove duplication. fabrice Bioinformatics 8 10-15-2012 09:10 PM
RNA-Seq: Single Read and Paired End mRNA-Seq Illumina Libraries from 10 Nanograms Tot Newsbot! Literature Watch 0 11-09-2011 02:10 AM
RNA-seq: Replicates, single-end, paired-end story pasta Bioinformatics 2 07-04-2011 11:51 PM

Reply
 
Thread Tools
Old 06-19-2011, 03:11 PM   #1
fabrice
Member
 
Location: paris

Join Date: Oct 2009
Posts: 86
Default Paired-end Illumina RNA-seq adapter trimming

Does anyone can recommend a tool for Paired-end Illumina RNA-seq adapter trimming? Most of tools are for single-end sequence. I think if there a tools used for paired-end data.
The shortcoming of the single end tools:
1, The single end tools sometimes will make the paired-end data mess.
2, It does not consider paired-end information, if there is a tool for paired-end, such remove will be accurately.
fabrice is offline   Reply With Quote
Old 06-20-2011, 12:14 AM   #2
Simon Anders
Senior Member
 
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 991
Default

HTSeq has facilities useful for this.

You can write a little Python script like this:

Code:
import itertools
import HTSeq

in1 = iter( HTSeq.FastqReader( "mydata_1.fastq" ) )
in2 = iter( HTSeq.FastqReader( "mydata_2.fastq" ) )
out1 = open( "trimmed_1.fastq", "w" )
out2 = open( "trimmed_2.fastq", "w" )

for read1, read2 in itertools.izip( in1, in2 ):
   read1.trim_right_end( "ACGGTC" )
   read2.trim_left_end( "TTCGAC" )
   read1.write_to_fastq_file( out1 )
   read2.write_to_fastq_file( out2 )
      
out1.close()
out2.close()
I haven't tested the script, and you will have to adjust it to your needs, of course.

Note that you can make the trimming tolerant to read errors by either adding a second argument, giving the allowed fraction of mismatches, or using trim_right_end_with_quals, which accepts a maximum sum of mismatch qualities. See here and here for an explanation.
Simon Anders is offline   Reply With Quote
Old 06-20-2011, 12:47 AM   #3
fabrice
Member
 
Location: paris

Join Date: Oct 2009
Posts: 86
Default

Simon Anders,
Thanks for your reply.
Could you please also give an example how to trim 3' reads which have low quality? For example, set a cutoff value. Bellow this cutoff, the nucleotide will be trimmed.

Thanks.
fabrice is offline   Reply With Quote
Old 07-20-2011, 05:31 PM   #4
saharnm
Junior Member
 
Location: United States

Join Date: Jul 2011
Posts: 1
Default

Take a look at "http://solexaqa.sourceforge.net/". It works on both single and paired-end reads.
saharnm is offline   Reply With Quote
Old 11-17-2011, 05:40 AM   #5
vbiaudet
Member
 
Location: Paris

Join Date: Apr 2011
Posts: 13
Default

We have the same problem for paired-end of 100bp and sizing of 200 (fragment size) from Illumina GAIIx RNAseq. In many case we think we got an adapter in 3' that have been sequenced... we think that because the mapping is very bad in 3'. BUT we know that the Quality score is very good (Q> 30). The problem is the 3'end can contains just some bases from the adapter and not the complete sequence... when the list of adapter is known did you find a tool to solve this problem ?

Thanks, vb
vbiaudet is offline   Reply With Quote
Old 03-28-2012, 08:22 PM   #6
ataheri
Junior Member
 
Location: Nashville, TN

Join Date: Nov 2011
Posts: 7
Default Trimmomatic

We are using Trimmomatic as it is able to do trimming and adaptor removal for paired-end reads. At first we tried fastx-toolkit but then we switched to trimmomatic as a better choice for paied-end samples.
ataheri is offline   Reply With Quote
Old 09-23-2013, 05:21 PM   #7
relipmoc
Member
 
Location: Los Angeles, CA

Join Date: Jul 2011
Posts: 58
Default skewer

You may also try skewer which is an adapter trimmer dedicated to Illumina paired-end data

Last edited by relipmoc; 09-23-2013 at 05:23 PM. Reason: title
relipmoc is offline   Reply With Quote
Old 11-17-2013, 11:54 PM   #8
saran.ela123
Junior Member
 
Location: Singapore

Join Date: Aug 2013
Posts: 5
Default

May i know from where to download skewer ?
saran.ela123 is offline   Reply With Quote
Old 01-05-2015, 07:48 AM   #9
FlorDV
Junior Member
 
Location: Uruguay

Join Date: Jan 2015
Posts: 1
Default

Hi saran.ela123

you can download it from here: http://sourceforge.net/projects/skewer/files/
FlorDV is offline   Reply With Quote
Reply

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:41 AM.


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