Trimming read length & quality in a FASTQ file
Hi all,
I was having a different computational problem, I had to compare a single end RNA-seq to paired end RNA seq. So my approach was to take the forward strand reads from the paired end sample and trim the read length and quality score of this fastq file and then map it using BWA, followed by regular stuff of estimating RPKM.
But how can I trim the read length and quality score?
any one liners in perl or awk to do this or is their something in picard tools?
Please share your experience and knowledge.
Thank you
1 answer
Hi I'm doing something similar, I'm using to trim reads features with cutadapt, is a good option and easy to use.
• 0 views
•
link
updated
by
Ram
•
written
by
gufernandez10 •
•
written
by
gufernandez10 •
Log in to answer this question.
More posts like this
-
Cross-strand correlation for paired-end ChIP-seq data
written by alison.ludzki •
-
How to extract the mean quality score of RNAseq for a sample?
written by Youyy •
-
Are single-end reads processed the same as paired-end reads regarding strandedness?
written by sg197 •
-
paired end data with varying read length
written by manjumoorthy95 •
-
Cutadapt wrong behaviour
written by dllopezr •
-
CoverageBed for RNA-seq strand-specific and paired-end data
written by anaritagrosso •
-
de novo transcriptome assembly work flow, paired end reads
written by nikelle.petrillo
-
Generating Read Length?
written by nikelle.petrillo
-
Copy Number Variation from paired end RNA-Seq data
written by ChIP
-
Paired-End Sequence Preprocessing Effects On Downstream Processes
written by Travis
If the pairs are in separate files, you could just work on each file and
substr()on awk whenNR%2==0, no?https://github.com/lh3/seqtk
http://www.usadellab.org/cms/?page=trimmomatic
In case, both of your paired end reads are in the same file first separate them and then trim the forward reads.
Suggestion: If your goal is to compare these two set of files for the end results (RPKM) in your case I don't think there is any point in just considering forward reads. You can compare a fragment library with a paired-end library and compare their complexities, contamination etc. And please use Splice aware aligner like TopHat or STAR if you want to count the reads spanning exons.