Question: Trimmomatic execution removes the data from my fastq sequence file. What am I doing wrong?
0
gravatar for kristina.mahan
7 months ago by
kristina.mahan30 wrote:

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 • 352 views
ADD COMMENTlink modified 7 months ago • written 7 months ago by kristina.mahan30

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

ADD REPLYlink written 7 months ago by kristina.mahan30
2
gravatar for lieven.sterck
7 months ago by
lieven.sterck4.5k
VIB, Ghent, Belgium
lieven.sterck4.5k wrote:

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 COMMENTlink modified 7 months ago • written 7 months ago by lieven.sterck4.5k
0
gravatar for h.mon
7 months ago by
h.mon24k
Brazil
h.mon24k wrote:

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 COMMENTlink modified 7 months ago • written 7 months ago by h.mon24k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1641 users visited in the last hour