STAR GeneCounts for most genes are 0
Hello all, I am relatively new to this field.
I am doing an RNA-seq alignment and expecting a gene count output for a prokaryote genome. I have an annotation GTF file in which I have converted the third column into "exon" for the alignment to work. Everything seemed well but the readspergene.out.tab file has mostly 0 for the gene counts.
Can anyone help diagnose what the problem is? Thank you
• 1,440 views
•
link
updated
by
Ram
•
written
by
Hansel Ivander •
•
written
by
Hansel Ivander •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Converting STAR Gene-level alignment to TPM expression
written by CTLong •
-
Aligning single transcript to sequence
written by bram_danneels •
-
Mapping reads and quantifying genes - Metagenomic workshop
written by arshad1292 •
-
Data filtering by creating dictionaries in python
written by mdgn •
-
bedmap to annotate bed file with GTF that was converted to BED produces way too many results
written by Tawny
-
htseq_counts process "Killed" during SAM alignment record pairs processed
written by bertb •
-
Add new sequence to clustalw2 alignment
written by DanielC
-
List of Exons per gene bad covered
written by windsur •
-
how to filter cuffdiff results by alignment count?
written by mike.thon •
-
GTF2/GFF3 "feature" types and expression analysis
written by ndejay •
What STAR commands did you run for generating the genome index and for aligning? Also, you may want to share both the STAR alignment log and a few lines from your modified GTF file.
It looks like most of the reads aligned outside your annotated “transcripts”. I would double-check that your annotations are in the same coordinate system as your genome.
I guess the problem is somewhere here and agree that you should post an excerpt from your GTF file. Coming from a prokaryote that doesn't have splicing it is likely not in a proper format and lacks exon annotation (and changing only the type from CDS or whatever to exon doesn't cut it). In principle, you won't even need a splicing aligner here at all for that matter.