gravatar for Apprentice

6 hours ago by

Hi.

I'm using Trimmomatic v0.39 to trim low quality base.
I set MINLEN option to 50 as below.

java -Xmx4g -jar trimmomatic-0.39.jar PE 
-threads 1 -phred33 -trimlog ${id}.trimlog 
${id}_R1.fastq.gz 
${id}_R2.fastq.gz 
${id}_R1.paired.fastq.gz ${id}_R1.unpaired.fastq.gz 
${id}_R2.paired.fastq.gz ${id}_R2.unpaired.fastq.gz 
TRAILING:20 MINLEN:50 
ILLUMINACLIP:adapters.fa:2:30:10

However, trimmed fastq file had reads with length<50.
I thought it might be dropped if both reads are less than MINLEN in length in the case of paired end. In other words, if one read length is less than MINLEN, but the other read is longer than MINLEN, it will not be dropped.
Is this idea correct?



Source link