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?
•
written
by
mengqingying •
0 answers
No answers yet.
Log in to answer this question.
-
the defination of insert size in RNA-seq, including the length of read, or not?
written by zju.whw •
-
Repetitive kmer profile in RAD seq data
written by markravinet •
-
Picard-Tools Or Trimmomatic Flip Read From Fr Or Rf??
written by mengqingying •
-
Merging Bam Files
written by hlsz.laszlo •
-
Picard Markduplicate Flags Reads
written by SHI Quan •
-
Soft Clipping Position
written by Nicolas Rosewick
-
3' Rna Seq Pipeline
written by k.nirmalraman
-
Strange Behavior In Picard Downsamplesam.Jar Algorithm?
written by jeffrey •
-
Haloplex & Allele Calling
written by Pierre Lindenbaum
-
How To Determine Overlaps From Coordinates
written by k.nirmalraman
Can't see the figure 2. Can you repost it?
It was just a formatting issue, it should be fixed now.