Hi, I am doing a quality check for my rna-seq data. Attached below are the data I'm working with. I have run a qc with FastQC and i need some advice on all the data like at which position i have to trim the bases and perhaps how to analyse and understand these data. Thank you.
I would use trim galore
eg. trim_galore --phred64 --paired --trim1 *fastqc
BBDuk in the BBMap package is incredibly fast with many, many options for use, highly recommended and great support from the author.
BBDuk is multithreaded and faster than Trimmomatic in my tests.
@ropni -
It appears that there is something seriously wrong with your data, and the files are corrupt. Illumina runs never produce 752bp or 896bp reads, which FastQC indicates you have. Are you sure this is the raw data? Where did the data come from, and what platform generated it?
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
version 2.3.6
Quality tends to fall off usually towards the end of the read, it is not normally a big deal if it isn't too bad and you can just align anyway - Tophat2 should take quality into consideration when aligning. Alternatively you can trim the end off until a threshold is reached, there is a fast tool called trimmomatic that can do that for you.
It looks like this may be PacBio Iso-Seq data? You would want to keep that in mind as you go forward with any analysis. PacBio Iso-Seq wiki would be useful for you.