Question: Picard-Tools Or Trimmomatic Flip Read From Fr Or Rf??
gravatar for mengqingying
5.3 years ago by
United States
mengqingying30 wrote:

I have some PE RNAseq reads contaminated with adapter sequences.

Before removing the adapters, i mapped the reads using bowtie2 -> samtools -> picard-tools "CollectInsertSizeMetrics.jar" and got histogram as shown in the first figure.

After removing the adapters with Trimmomatic, i again run bowtie2 -> samtools -> "CollectInsertSizeMetrics.jar" and got the 2nd histogram with many RF reads, as shown in blue color in the plot. Actually here i already set bowtie2 parameter --no-discordant to remove discordant reads.

I cannot figure out the reason there are so many reads flip their orientations from FR -> RF? I wonder if this change are caused by Trimmomatic, or picard? How this happened?! Any one can help?

Figure1_beforeQC Figure2_afterQC

rnaseq • 1.7k views
ADD COMMENTlink modified 5.3 years ago by Devon Ryan88k • written 5.3 years ago by mengqingying30

Can't see the figure 2. Can you repost it?

ADD REPLYlink written 5.3 years ago by Ashutosh Pandey11k

It was just a formatting issue, it should be fixed now.

ADD REPLYlink written 5.3 years ago by Devon Ryan88k
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: 1882 users visited in the last hour