Trimmomatic not recognising both adapters
1
0
Entering edit mode
3.1 years ago

I received some WES samples back from a commercial sequencing company, without an adapters file, so I created an adapters.fa file that looks like this:

>Novogene adapter,PE, 5'adapter
AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT
>Novogene adapter, PE, 3' adapter
GATCGGAAGAGCACACGTCTGAACTCCAGTCACATCACGATCTCGTATGCCGTCTTCTGCTTG

I put this into my script as ILLUMINACLIP:adapters.fa: 2:30:10.

This works but in my output file I get:

TrimmomaticPE: Started with arguments:
...
Using Long Clipping Sequence: 'GATCGGAAGAGCACACGTCTGAACTCCAGTCACATCACGATCTCGTATGCCGTCTTCTGCTTG'
ILLUMINACLIP: Using 0 prefix pairs, 1 forward/reverse sequences, 0 forward only sequences, 0 reverse only sequences
Quality encoding detected as phred33
Input Read Pairs: 51487885 Both Surviving: 50756708 (98.58%) Forward Only Surviving: 0 (0.00%) Reverse Only Surviving: 0 (0.00%) Dropped: 731177 (1.42%)
TrimmomaticPE: Completed successfully

Does this mean it only looks for one of the adapters, or does it mean it did not find any of the 5' adapter?

trimmomatic adapters • 1.3k views
ADD COMMENT
2
Entering edit mode
3.1 years ago

I will say that the adapter naming is a bit weird, the difference between the names has odd whitespace. Note how both sequences are called the same:

Novogene

because the sequence id is defined as the word up until the first whitespace.

That might make Trimmomatic consider them a single sequence. In general sequences ought to have unique ids for example:

Novogene_adapter_1

and

Novogene_adapter_2

ADD COMMENT
0
Entering edit mode

Thanks, that's really helpful. I'll edit and try again (or try with TrimGalore where you can just list the sequence).

ADD REPLY
0
Entering edit mode

Yep- I changed the name and it works!

ADD REPLY
0
Entering edit mode

Please accept Istvan's answer below (green check mark) to provide closure to this thread.

ADD REPLY

Login before adding your answer.

Traffic: 2631 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6