SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA Soft Clipping Bio.X2Y Bioinformatics 11 03-09-2015 09:08 AM
Chimeric reads hard clipping fanx Bioinformatics 2 05-29-2013 02:10 AM
unmapped vs soft-clipping CNVboy Bioinformatics 2 04-08-2012 10:41 PM
Clipping adapters nsl RNA Sequencing 6 05-27-2011 02:08 PM
what's clipping point? jordi 454 Pyrosequencing 1 05-26-2009 01:47 AM

Reply
 
Thread Tools
Old 12-08-2012, 11:14 AM   #1
MeganS
Member
 
Location: US

Join Date: Sep 2010
Posts: 14
Default Trimmomatic palindrome clipping (MIN_PREFIX=8)

I have Illumina RNA-Seq data with extremely short fragment sizes (150bp for 100bpPE reads). I am using Trimmomatic because I like how the palindrome clipping handles read-through (and bonus points for handling orphans for me). However, after running Trimmomatic I still see residual adapter contamination in the last 7bp or less of the 3' end of some of the reads.

Looking at the source code, I see there is a parameter MIN_PREFIX=8 (in IlluminaClippingTrimmer.java). When I change it to MIN_PREFIX=1, I get the desired result. I think I understand the code, but I don't code in Java. Is there a reason that I am missing that makes it a bad idea to lower this value? My reasoning is that if the left and right reads overlap in such a way that indicates read-through, even if only 1 base of the adapter was sequenced, then we can safely trim the small bit of adapter and throw out the second read.
MeganS 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 11:40 PM.


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