Question: Trimmomatic : keepBothReads True vs False option
gravatar for candida.vaz
3.6 years ago by
candida.vaz30 wrote:

Hello everyone,

Is anyone familiar with the "keepBothReads" in Trimmomatic :

keepBothReads: After read-through has been detected by palindrome mode, and the adapter sequence removed, the reverse read contains the same sequence information as the forward read, albeit in reverse complement. For this reason, the default behaviour is to entirely drop the reverse read. By specifying "true" for this parameter, the reverse read will also be retained, which may be useful e.g. if the downstream tools cannot handle a combination of paired and unpaired reads

I was wondering if the keepBothReads: "True" option is better than using the default "False" option

  • First, I'm using STAR aligner and Tophat2 for mapping
  • Second, Will there be any compromise in the downstream quantification (by Htseq-count or Cufflinks) if I choose to use "True" and keep both the reads?

Thanks for your time and concern!

rna-seq • 2.1k views
ADD COMMENTlink modified 3.6 years ago by Brian Bushnell17k • written 3.6 years ago by candida.vaz30
gravatar for Brian Bushnell
3.6 years ago by
Walnut Creek, USA
Brian Bushnell17k wrote:

Dropping one read of a pair leaves you with some read pairs and some singletons, which is generally annoying to deal with, so it's better to keep both when they are perfectly good quality-wise. Alternatively, merging the reads so that the singletons represent the consensus of the pair is another good solution. Just dropping one at random is not usually ideal.

ADD COMMENTlink modified 3.6 years ago • written 3.6 years ago by Brian Bushnell17k

Dropping reads for any other reason than quality ... is a great way to ruin your data set.

ADD REPLYlink written 3.6 years ago by YaGalbi1.5k
Please log in to add an answer.


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