Trimmomatic execution removes the data from my fastq sequence file. What am I doing wrong?
2
0
Entering edit mode
2.6 years ago

I'm trying to use Trimmomatic 0.38 to trim my Illumina fastq sequences. Something must be wrong with my syntax because after execution I get a new empty file named "ILLUMINACLIP:TruSeq3-PE.fa:2:30:10" and one of my fastq sequences goes from 65G or 0G. Here is my command:

[kmmahan@bnode4 KMM2_analysis]$ java -jar ~/projects/algae/ABY2/Trimmomatic/Trimmomatic-0.38/trimmomatic-0.38.jar PE -threads 140 -phred33 -trimlog ~/projects/algae/ABY2/sequencing/KMM2_analysis/KMM2-A_R1_001.fastq /KMM2-A_R2_001.fastq KMM2_output_forward_paired_trimmed.fastq KMM2_output_reverse_paired_trimmed.fastq KMM2_output_forward_unpaired_trimmed.fastq KMM2_output_reverse_unpaired_trimmed.fastq ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36

TrimmomaticPE: Started with arguments:
 -threads 140 -phred33 -trimlog /home/kmmahan/projects/algae/ABY2/sequencing/KMM2_analysis/KMM2-A_R1_001.fastq /home/kmmahan/projects/algae/ABY2/sequencing/KMM2_analysis/KMM2-A_R2_001.fastq KMM2_output_forward_paired_trimmed.fastq KMM2_output_reverse_paired_trimmed.fastq KMM2_output_forward_unpaired_trimmed.fastq KMM2_output_reverse_unpaired_trimmed.fastq ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36
Trimmomatic • 1.2k views
ADD COMMENT
0
Entering edit mode

Added output file for trimlog and put the 4 output files in correct order and it looks like it's working. Thanks!

ADD REPLY
2
Entering edit mode
2.6 years ago

at first sight you did not provide an output file to write the log info the option -trimlog takes a file name as argument (now it will indeed overwrite one of your input files)

you also need to provide the output as <paired out 1> <unpaired out 1> <paired out 2> <unpaired out 2>

And I understand you want to speed up the analysis running it multi-core but 140 is seriously over the top, it will normally run at ok speed using 4-6-10 cores

ADD COMMENT
0
Entering edit mode
2.6 years ago
h.mon 32k

You have to pass the full path to both R1 and R2 files, not just for the R1 file as you are doing. Also, I don't know if the rest of your command is correct, as I don't use Trimmomatic, but I believe it isn't, as the manual states it needs either a --baseout parameter, or 4 file names for output files:

Paired-end mode requires 2 input files (for forward and reverse reads) and 4 output files (for forward paired, forward unpaired, reverse paired and reverse unpaired reads).

So the order of your output files won't be what you expect. Here is what I think is the correct command:

java -jar ~/projects/algae/ABY2/Trimmomatic/Trimmomatic-0.38/trimmomatic-0.38.jar PE 
  -threads 140 -phred33 -trimlog KMM2_log.txt \
  ~/projects/algae/ABY2/sequencing/KMM2_analysis/KMM2-A_R1_001.fastq \
  ~/projects/algae/ABY2/sequencing/KMM2_analysis/KMM2-A_R2_001.fastq \
  KMM2_output_forward_paired_trimmed.fastq KMM2_output_forward_unpaired_trimmed.fastq \
  KMM2_output_reverse_paired_trimmed.fastq KMM2_output_reverse_unpaired_trimmed.fastq \
  ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36

edit: lieven.sterck is right, you also forgot to provide a name to the log file. I edited my suggestion to fix this as well.

ADD COMMENT

Login before adding your answer.

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