What is Trimmomatic -validatePairs, and how to use it?
Hi,
I use Trimmomatic v0.35 and I notice the option -validatePairs:
Usage:
PE [-threads <threads>] [-phred33|-phred64] [-trimlog <trimLogFile>] [-quiet] [-validatePairs] [-basein <inputBase> | <inputFile1> <inputFile2>] [-baseout <outputBase> | <outputFile1P> <outputFile1U> <outputFile2P> <outputFile2U>] <trimmer1>...
or:
SE [-threads <threads>] [-phred33|-phred64] [-trimlog <trimLogFile>] [-quiet] <inputFile> <outputFile> <trimmer1>...
But I can't find a description of this option in the manual. I want to filter away singletons from two paired fastq-files and I tried to run trimmomatic with this option and the two paired files as input, but it didn't work.
Thanks, Jon
1 answer
My guess is you are using an incorrect Trimmomatic command, and the issue is not related to -validatePairs. I can use it:
TrimmomaticPE -phred33 -validatePairs 1.fastq 2.fastq \
ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36
results in:
TrimmomaticPE: Started with arguments:
-phred33 -validatePairs 1.fastq 2.fastq ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:12
Multiple cores found: Using 2 threads
WARNING: Pair validation failed at record: 0
Forward read: 1/1
Reverse read: 4/2
Input Read Pairs: 3 Both Surviving: 3 (100.00%) Forward Only Surviving: 0 (0.00%) Reverse Only Surviving: 0 (0.00%) Dropped: 0 (0.00%)
TrimmomaticPE: Completed successfully
TrimmomaticPE -phred33 1.fastq 2.fastq \
ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36
results in:
TrimmomaticPE: Started with arguments:
-phred33 1.fastq 2.fastq ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:12
Multiple cores found: Using 2 threads
Input Read Pairs: 3 Both Surviving: 3 (100.00%) Forward Only Surviving: 0 (0.00%) Reverse Only Surviving: 0 (0.00%) Dropped: 0 (0.00%)
TrimmomaticPE: Completed successfully
Log in to answer this question.
More posts like this
-
Error: Trimmomatic Command not found
written by Tundup •
-
MultiQC doesn't produce data for all samples, why? (trimmomatic)
written by john •
-
Problem with Trimmomatic Multiple Paired-End sequencing.
written by john •
-
How to decide which parameter to change for more quality results while trimming with Trimmomatic?
written by sam •
-
single end Trimmomatic
written by gina02 •
-
Trimmomatic errors
written by marion.ryan •
-
Trimmomatic 0.39 Permission denied accessing fastq.gz
written by e.wing •
-
Trimmomatic couldn't generate output file
written by leukippus0116 •
-
Trimmomatic paired-end mode question
written by slin023 •
-
Trimmomatic 0.36 not acknowledging commands in MacOS terminal in SE mode
written by msimmer92
Hi, Please change the post to a question. It is not a job.
Yes, it's certainly not a job! Thanks! :-)
Describe how "it didn't work. Were there errors? The output is not what you expected?
Do your files contain singletons from the start? This only occurs when they are processed incorrectly. Otherwise, Trimmomatic will output properly paired reads.
My questions is really how to use trimmomatic with this option? I can't find any mention of it in the manual. I tried to launch trimmomatic with
-validatePairs input1.fq input2.fq paired1out.fq single1out.fq paired2out.fq single2out.fqbut nothing happened. Only the usage menu showed.And if you use the same command without
-validatePairs, does it work as expected?Yes, I can do the trimming and everything and it spits out the two correct paired files and the two single files. But in this case I wondered whether I could use trimmomatic just for filtering out singletons, without doing any trimming or filtering.
Probably not. You can use
repair.shfrom the BBTools package for this.Ok, I'll check it out. Thanks for the advice :)