Go Back   SEQanswers > Bioinformatics > Bioinformatics

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

Thread Tools
Old 12-08-2012, 12:14 PM   #1
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 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

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 10:41 PM.

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