Genome vs transcriptome alignement in STAR
Hi All,
I suspect that my sample might have dna contamination. Is there a way of telling the fraction of reads mapped to the transcriptome vs genome when using the STAR aligner?
I have used the --quantMode TranscriptomeSAM flag so I guess that can use samtools to count the number of reads mapped to the trascriptome
Thanks!
• 2,965 views
•
link
updated
by
grant.hovhannisyan
•
written
by
roy.granit
•
written
by
roy.granit
1 answer
If you used --quantMode GeneCounts, then in count file there will be N_noFeature, i.e. reads that were not assigned to any feature. It can be a proxy of how much "genomic" reads you have. However if you work a non-model organism it can be that these reads are coming from not yet annotated genes.
• 0 views
•
link
written
by
grant.hovhannisyan
Log in to answer this question.
More posts like this
-
Higher number of reads mapped to transcriptome vs genome in star mapping
written by paguirreazorin •
-
STAR pair-end gene count discrepancy vs. RSEM
written by anon •
-
Best Transcriptome file for Salmon after STAR alignment
written by rayanelkholdi •
-
Mitochondrial RNA as scRNAseq quality check
written by gregoire.destreel •
-
Error after STAR mapping
written by potter •
-
running salmon on bam files
written by inah •
-
STAR mapping statistics don't add up with samtools count
written by Zee_S
-
snakemake errors and threads
written by Assa Yeroslaviz
-
Interpreting ATAC-seq alignment results
written by hara •
-
BWA or STAR for RNAseq?
written by agata88
