trimmomatic error with unknown trimmer
1
1
Entering edit mode
8.0 years ago
xuanyuan ▴ 10

I am running trimmomatic using the following command:

java -jar ~/tools/Trimmomatic-0.32/trimmomatic-0.32.jar PE R1_001.fastq R2_001.fastq pair1.fq unpair1.fq pair2.fq unpair2.fq step 1

and then I got the following wrong messages:

TrimmomaticPE: Started with arguments: R1_001.fastq R2_001.fastq pair1.fq unpair1.fq pair2.fq unpair2.fq step 1
Multiple cores found: Using 12 threads
Exception in thread "main" java.lang.RuntimeException: Unknown trimmer: step
        at org.usadellab.trimmomatic.trim.TrimmerFactory.makeTrimmer(TrimmerFactory.java:60)
        at org.usadellab.trimmomatic.TrimmomaticPE.run(TrimmomaticPE.java:495)
        at org.usadellab.trimmomatic.Trimmomatic.main(Trimmomatic.java:35)

Does anybody know what this means? Is there anything wrong?

Thanks

sequence next-gen-sequencing RNA-Seq • 7.0k views
ADD COMMENT
1
Entering edit mode

"step 1" is not a valid argument for trimmomatic. What are you trying to do?

ADD REPLY
0
Entering edit mode

i am try to remove the adaptor from the illumina rna-seq raw data. so should i replace <STEP 1> with <ILLUMINACLIP>? thanks a lot!

ADD REPLY
2
Entering edit mode

Yes. You should write in the following format:

ILLUMINACLIP:<fastaWithAdaptersEtc>:<seed mismatches>:<palindrome clip threshold>:<simple clip threshold>

In <fastaWithAdaptersEtc> field you should specify the path to a fasta file with adapters sequences. For example:

ILLUMINACLIP:/home/dir/adapters.fa:2:30:10
cat adapters.fa:

>adapter1
AACACACAAGAGGA
>adapter2
GGAGGAGGCATTT
ADD REPLY
0
Entering edit mode

get it. thanks a lot!

ADD REPLY
2
Entering edit mode
8.0 years ago

You don't literally use "step 1". That's just a placeholder for a step in how you want it to trim things. For example, you might use "LEADING:10" there, to specify trimming leading bases with Phred scores below 10.

ADD COMMENT
0
Entering edit mode

thanks a lot!

ADD REPLY

Login before adding your answer.

Traffic: 1427 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