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.
•
written
by
ropni •
2 answers
What makes you think this is ascii-64 encoded? The charts indicate that it is ascii-33, and all modern data is ascii-33, so I would conclude that it is ascii-33.
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?
Log in to answer this question.
-
Error in "Per Tile Sequence Quality"
written by SomeOne •
-
Fastqc of the fastq file
written by sgadila •
-
QC prior to bacterial genome assembly
written by Pegasus •
-
Per base sequence quality and per base N content
written by misterie
-
Why does read quality drop after adapter trimming with cutadapt?
written by Nabeel Ahmed •
-
Quality check of Medip Seq Data
written by gskbioinfo143 •
-
QC of RNA seq
written by Arash •
-
Filtering multisample VCF based on genotype using SnpSift filter
written by bioinforesearchquestions
-
how to choose optimal K in Consensus clustering
written by LJ
-
Read counts based on genes from samtools, feature count and htseq-count
written by bioinforesearchquestions
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.