Hi, I have a set of RNA-seq data, I used STAR for mapping then used RNA-SeQC to check the quality. After STAR mapping I checked the Log.final.out file, there are 81% reads uniquely mapped. I then used RNA-SeQC to produce quality metrics. I used collapsed gtf as input. However in the metrics.tsv file only 18% reads uniquely mapped, and the gene_tpm.gct.gz file and other files contain nan or 0 reads for each gene value, like below: ENSG00000227232.5 WASH7P 0 or ENSG00000223972.5 DDX11L1 -nan. It seems something wrong with RNA-SeQC step. I downloaded rnaseqc.v2.3.6.linux, and used following command: rnaseqc ${genes_gtf} ${bam_file} . -s ${sample_id} --stranded rf -vv. No error message showed up, so I'm not sure where went wrong. Any thoughts? Thank you.
•
written
by
koaeb •
1 answer
You probably have some name mismatch in your reference genome vs gtf. If you like you can sign up for a free trial at Truwl and run GTEx v10 there: https://truwl.com/workflows/library/GTEx/v10
Log in to answer this question.
-
RNASeq QC and STAR results confused
written by shafi98 •
-
Asking for unmapped genes in RNA sequencing via Galaxy
written by Atefeh •
-
STAR pair-end gene count discrepancy vs. RSEM
written by anon •
-
STAR alignment speed
written by Gilles •
-
ReadsPerGene.out.tab column one contains gene ID AND gene name
written by lalina1606 •
-
STAR mapping results
written by mxlsherry1992 •
-
Empty BAM file upon mapping with STAR
written by makwana.kd •
-
Very low unique mapped reads - STAR alignment
written by jsl •
-
Mismatch number of reads after bedtools bam2fastq on RNA-seq data
written by aditi.qamra
-
Recommended doc for interpreting STAR alignment stats
written by emblake •
Hi, did you ever find a solution for this?
Similar trouble here via the GTEx pipeline (v10, so using RNA-SeQC v2).
From metrics file:
But still all zeros in the count data...