HTSeq-count: gene models/transcripts are in FASTA format not gtf/gff
I've generated .SAM files with RNA-seq reads aligned to a reference genome using STAR. I ultimately want to pull out differential expressed genes from my dataset.
I've chosen to use HTSeq-count to generate count files to ultimately use in either DEseq2 or EdgeR. It seems I need to provide a gff/gtf file with transcripts or gene models. I have access to both the gene models and transcripts for the genome, but both files are in FASTA format.
What can I do to solve this issue? Is there another program, besides HTSeq I could use to generate count files?
Thank you so much!
• 1,486 views
•
link
written
by
sbrashidx •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Error with HTseq RNAseq read count
written by Catalina •
-
Proper HTSeq usage on bacterial genome. Don't quite understand --t
written by SushiRoll •
-
TopHat2 to align RNA with or without annotation file
written by mario_de_manner •
-
Merged vs. Ensembl GTF for htseq-count
written by kousi31 •
-
Take RNASeq data and find number of viral transcripts of specific virus
written by shmo2 •
-
STAR: Is it ok to align with an index created from genecode but count features from gtf from ensemb…
written by simplitia •
-
Error on HT-Seq-count while RNA quantification
written by k.kathirvel93
-
HTSEq-Count returns 0 counts for all genes in the counts file!
written by manish.nu •
-
Multibamcov Or Htseq-Count To Count Read Per Feature ?
written by Nicolas Rosewick
-
Gff Format And Rna Seq Alignment
written by User 2889 •
You could use kallisto and then import into DESeq2 via tximport. The DESeq2 vignette at bioconductor has an example for this.
Which reference genome did you use? I am sure you can get the matching gtf/gff files somewhere.