I've been using this script for read pair resynchronization.
paired end result for preprocessing
Hi,
I am doing preprocessing analysis for RNA-seq. After I removed the adapter, the total sample reads is still the same as total high quality reads. The worse problem is, when I run the select paired end read, the result obtained was zero bytes. So how can I overcome this problems?
Thank you.
1 answer
By quality-trimming them independently, you have broken read pairing. You should use any program which keeps pairing for trimming (BBDuk, cutadapt, SeqyClean, many others).
You may also use Repair (from BBTools) to fix the pairing from the reads trimmed by DynamicTrim.
I'm quiet confused about why we should trimmed read as pairs. I already check the quality for both reads via fastqc tools and it shows my reads have high quality scores. My problem now is when I want to select the paired end read that remained after I had removed all the contaminant such as adapter, the result shown zero byte means there is nothing in the file.
We should trim reads as pairs to avoid the annoyance of resyncing them again.
• 1 views
•
link
written
by
Devon Ryan
Log in to answer this question.
More posts like this
-
RNA editing sites based on DNA seq and RNA seq
written by Ahiad Chen Zion •
-
barcode of TCR-sequencing
written by yueli7
-
FastQC duplicates questions
written by Beth •
-
How to select the reads with a certain adapter in a fastq file in R
written by Youyy •
-
Only a small portion of my reads have Smart-seq adapters, why?
written by Youyy •
-
Reduce number of reads from 100M to 50M in the sample from the existing RNA-Seq fastq files
written by mohammedtoufiq91
-
adapter trimming PE RNA reads
written by rahmati.razieh83 •
-
Interpreting trim_galore report
written by vaslanzadeh •
-
How can I remove adapter sequence from illumina 2000 paired end data?
written by tcf.hcdg
-
Improving alignment rates for RNA-seq
written by pm2012
"when I run the select paired end read"
Can you clarify what you mean by that? I understood that trimming didn't remove any reads (that's good, it just means you had minimal adapter contamination), but I haven't a clue what you mean by what I quoted above.
What I mean by that quoted is, I want to obtain paired end reads after remove adapter sequence. But in my case, the paired end reads that I obtained give me no result. So, as you know, if we can't identify the paired end reads, we can't proceed to get shuffled paired end sequence.
How are you trimming them? They should be trimmed as pairs, or are you doing that and everything is being removed?
I had trimmed them separately based on base quality (Qphread=20) and sequence length (min 50 bp). I had done these step in previous analysis with different data and it works. However, I am having problem with current analysis. For your additional information, the size of my reads is up to 5gb each while my previous size of reads is only up to 1gb. It is the size of the reads will influence the outcome for paired end reads?
Could you post which tool and the commands you are using?
I only use perl script command line for my preprocessing analysis.
This is example how I trimmed the reads based on base quality (Qphread = 20)
where
reads_R1&reads_R2are referring to forward and reverse file.I have one additional question, now I'm trying to use trimmomatic tool to trim the reads. However, based on the manual that was provided, I need to insert trim log file to run the command. But, I don't know how to create trim log file.
This is the command that I used:
You don't create the log, the tool does and puts it into a file whose name you provide. Note that the log is optional, so you don't need to provide it.
Thank you Devon Ryan for the info. I will try again.