Question: What is a good Trimmomatic parameter?
0
gravatar for Arindam Ghosh
16 months ago by
Arindam Ghosh160
India
Arindam Ghosh160 wrote:

I am trying to clean RNAseq reads with TRIMMOMATIC using the parameters as:

java -jar trimmomatic-0.36.jar PE -trimlog file_trim_log input_1.fastq.gz input_2.fastq.gz output_1P_clean.fq.gz output_1U_clean.fq.gz output_2P_clean.fq.gz output_2U_clean.fq.gz ILLUMINACLIP:TruSeq3-PE-2.fa:2:30:10 SLIDINGWINDOW:10:30 LEADING:28 TRAILING:28 MINLEN:75

ILLUMINACLIP:TruSeq3-PE-2.fa:2:30:10 : My reads contains TruSeq Universal adapters as indicated by FASTQC.

SLIDINGWINDOW:10:30 : As per manual

LEADING:28 : Willing to remove bases whose quality falls below 28 in FASTQC per base quality module

TRAILING:28 : Willing to remove bases whose quality falls below 28 in FASTQC per base quality module

MINLEN:75 : Minimum length should not fall below 75bp

Are the parameters reasonable or too strict?

For Leading and Trailing the quality score in several manuals were given 3. Isn't that too low?

I am a bit concerned because only 70-80% of the paired reads survived. The rest either fall of Forward/Reverse Only Surviving or Dropped.

ADD COMMENTlink modified 16 months ago • written 16 months ago by Arindam Ghosh160
5

This is totally up to you ! It's more a biological question.

AFAIK trimming reads in RNAseq is not a mandatory step anymore. People say that if your reads are not good enought to be aligned they will not. (If you want to do an alignment after that of course)

But if you really want to trim these reads, you will have to check the 5' and 3' quality of your data using tools like FastQC

If you have a really low quality on the first 5 bases (let's say), you can trim them, but if you have good quality on these bases, why trim them ?

Here, the good post you need : Trimming single end reads for STAR?

And the publication : https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4766705/

ADD REPLYlink modified 16 months ago • written 16 months ago by Bastien HervĂ©4.4k
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: 890 users visited in the last hour