Trimmomatic not trimming the adapter

1

Entering edit mode

5.5 years ago

S ▴ 100

Hello,

I am using Trimmomatic to trim the adapter sequence from the raw small RNA data. For some reason after trimming I still find the adapter in the trimmed sequence. Does anyone have any suggestion?

This is the command I am using:

java -jar trimmomatic SE -phred33 -threads 10 file_in file_out ILLUMINACLIP:adapter.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 CROP:34 MINLEN:18

This is the raw sequence data (as an example) - the adapter sequence is in bold:

@SN741:1069:CCAUHACXX:2:2201:1721:2117 1:N:0:ATCACG
NAGCGTTGGATTGTTCACCC**TGGAATTCTCGGGTG**CCAAGGAACTCCAGTC
+

Here is the trimmed sequence after running Trimmomatic:

@SN741:1069:CCAUHACXX:2:2201:1721:2117 1:N:0:ATCACG
AGCGTTGGATTGTTCACCC**TGGAATTCTCGGGTG**
+

Any help is appreciated.

Trimmomatic small RNA • 2.6k views

ADD COMMENT • link updated 5.3 years ago by Biostar 20 • written 5.5 years ago by S ▴ 100

0

Entering edit mode

Please use the formatting bar to indicate code. I did it for you this time.

bar

ADD REPLY • link 5.5 years ago by ATpoint 81k

0

Entering edit mode

Sorry, I will do next time, I overlook it. Thanks!

ADD REPLY • link 5.5 years ago by S ▴ 100

0

Entering edit mode

You could try a different tool e.g. fastp.

ADD REPLY • link 5.5 years ago by alex.zaccaron ▴ 410

1

Entering edit mode

Hi Alex,

I actually did try a different tool in the meantime. I tried cutadapt and it worked perfectly. However, I am still wondering whether Trimmomatic will still work for this or not.

Thanks!

ADD REPLY • link 5.5 years ago by S ▴ 100

0

Entering edit mode

Could you please show the content of adapter.fa and also a complete fastq entry including the quality values for the read bases?

fin swimmer

ADD REPLY • link 5.5 years ago by finswimmer 16k

0

Entering edit mode

Hi finswimmer,

The contents of the adapter.fa is indicated in my question under the double asterisk. I meant for it to be displayed in bold but not sure why it isn't. I think it is because it is under the code block? Below is the raw sequence with the quality values.

@SN741:1069:CCAUHACXX:2:2201:1721:2117 1:N:0:ATCACG
NAGCGTTGGATTGTTCACCCTGGAATTCTCGGGTGCCAAGGAACTCCAGTC
+
#0<FFFFFBBFFFFIFFFFFIIIIFFFIIIFBF0BBFIFF0000BFFFIFB

Thanks!

ADD REPLY • link 5.5 years ago by S ▴ 100

0

Entering edit mode

are you sure it's using the sequence from your own adapter.fa file (and not the 'build-in' one)? Perhaps try running it with ILLUMINACLIP:./adapter.fa:2:30:10 (specifically point to the 'local' file)