Ok. I would like to use Picard tools. I have the gene annotation file in GFF3 format. How do I convert it to REF_FLAT (File) format?
I've completed mapping the RNA-seq to the reference genome using HISAT2. Now I would like to assess the quality of the samples. What is the recommended tool to do it? PS: I am only a beginner.
1 answer
FastQC is good for both pre- and post- alignment QC. Picard tools has a bunch of QC related tools you can use. For starters, I'd first look at the % of uniquely mapped reads to get a general idea of how well your sequences aligned to the reference. For human cell lines, a good alignment will be 50%+. Paraffin embedded tissue or other similar preparation methods will likely have more degradation so I'd assume that the alignments might suffer a bit. QC will typically assess the sequencing quality (how well the sequencer sequenced what you gave it to sequence), and not necessarily the quality of the material you gave the sequencer to sequence.
Any possible reference for the statement "For human cell lines, a good alignment will be 50%+. ".
Log in to answer this question.
RNA-SeQC for one. What do alignment stats look like?
HISAT2 summary stats: