STAR mapped more reads to transcriptome
Hi,
I am following ENCODE pipeline to analyze my RNA-Seq data. However, STAR gives me more reads in the Aligned.toTranscriptome.out.bam compared to the Aligned.sortedByCoordinate.bam, which means it aligned more reads to transcriptome rather than genome. This is beyond my understanding. Did anyone meet this problem? Thanks a lot!
Best, Lei
• 5,732 views
•
link
written
by
tianleivv •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Does TCGA data only use host reads to make pre-processed data (CNV, Mutations status, RNA-seq)?
written by yahn •
-
Using Salmon to map RNA-Seq and Ribosome Profiling reads in Yeast and downstream analysis with DESe…
written by kyusikkim •
-
Analyzing RNA-Seq with duplicate ensembl ids using DESeq2: should tximport be used?
written by Ridha •
-
What is isogenic replicate in ENCODE project
written by progistar •
-
Discrepancy in number of pseudoaligned reads for read quantification with Kallisto
written by nicolasgoedert1997 •
-
Differential Peak Analysis from ENCODE Pipeline normalized outputs
written by sjquon •
-
STAR - multi-mapping reads cut off
written by is.birds •
-
reasons for htseq-count with zero features
written by shenwei1376 •
-
STAR ENCODE option
written by tianleivv •
-
All the reads aligned using STAR have low mapping quality (0-3). What is happening?
written by kirannbishwa01
It's hard to say w/o knowing specifically how you ran STAR. Maybe a difference in the total # of alignments (which may filter more reads in the genome alignment vs transcriptome)?
Hi Chris, Thanks for your reply! I am actually following the ENCODE STAR+RSEM pipeline (https://github.com/ENCODE-DCC/long-rna-seq-pipeline/blob/master/DAC/STAR_RSEM.sh). I think no matter how I ran STAR, it should filter the reads equally for these mapped to transcriptome and genome. Thus, reads mapped to genome should always no less than reads mapped to transcriptome. Lei
It is still hard to say, because we don't know what are the numbers you are talking about. Is it about total counts to genes / transcripts provided by the two methods? Or is it information about number of mapped reads from STAR logs? Or did you run something like
samtools flagstat?Besides, you showed the script you run, but not how you run the script. What flags did you pass to the script?