STAR: Is it ok to align with an index created from genecode but count features from gtf from ensembl?
I'm wondering if there maybe some coordinate issues as with UCSC and ensembl. I happen to have bam files which were aligned with STAR using indices that were created from gencode however, I would like to use htseq count but with gtf donwloaded from ensembl instead?
thanks! A
• 2,069 views
•
link
updated
by
Kristoffer Vitting-Seerup
•
written
by
simplitia •
•
written
by
simplitia •
1 answer
You will most likely have a problem. Gencode uses chromosome names: "chr 1", "chr 2" etc. Ensembl uses chromosome names "1", "2" etc. Why not use gencode gtf file?
• 0 views
•
link
written
by
Kristoffer Vitting-Seerup
Log in to answer this question.
More posts like this
-
Consistency in the GTF files used for the index building of STAR alignment and RSEM
written by Xunzhi.Zhang •
-
Merged vs. Ensembl GTF for htseq-count
written by kousi31 •
-
HTSeq errors in the command line
written by Morris_Chair
-
HTSEq-Count returns 0 counts for all genes in the counts file!
written by manish.nu •
-
mm9 gtf file gencode
written by skylinesky •
-
GAPDH missing from htseq-counts file
written by melia_g •
-
Ensembl GTF for Mouse that includes intron, exons, utrs, etc.
written by achamess
-
High number of no feature while using htseq-count on Tophat2 aligned as well as STAR aligned paired…
written by fibi_prog_omics •
-
Results differ when using different ensembl versions
written by gudraephouto •
-
RNAseq: Mapping to ORFs
written by chimeric •
Why not keeping things simple and download GTF from Gencode? There are a lot of pitfalls in bioinformatics due to different formatting of files and gene names etc. Don't mix annotation files!