Quality trimming in R using ShortRead and Biostrings
[This article was first published on Jermdemo Raised to the Law, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
I wrote an R function to do soft-trimming, right clipping FastQ reads based on quality.
This function has the option of leaving out sequences trimmed to extinction and will do left-side fixed trimming as well.
#softTrim #trim first position lower than minQuality and all subsequent positions #omit sequences that after trimming are shorter than minLength #left trim to firstBase, (1 implies no left trim) #input: ShortReadQ reads # integer minQuality # integer firstBase # integer minLength #output: ShortReadQ trimmed reads library("ShortRead") softTrim
To use:
library("ShortRead") source("softTrimFunction.R") #or whatever you want to name this reads
I strongly recommend reading the excellent UC Riverside HT-Sequencing Wiki cookbook and tutorial if you wish to venture into using R for NGS handling. Among other things, it will explain how to perform casting if you have Solexa scaled (base 64) fastq files. The function should respect that. http://manuals.bioinformatics.ucr.edu/home/ht-seq
To leave a comment for the author, please follow the link and comment on their blog: Jermdemo Raised to the Law.
R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.