Question: Want read count with respect to geneID instead of transcript IDs
1
gravatar for archana.bioinfo87
12 months ago by
archana.bioinfo87180 wrote:

Hi,

Maybe my question is very basic but I just want to know how I will get the reads count with respect to geneID. I have RNA-seq data and I used htseq-count and got read count in terms of Ensembl Transcript ID.

Thanks

rna-seq forum • 211 views
ADD COMMENTlink modified 12 months ago by lieven.sterck9.0k • written 12 months ago by archana.bioinfo87180
3
gravatar for colin.kern
12 months ago by
colin.kern930
United States
colin.kern930 wrote:

The default arguments for htseq-count will count by Ensembl gene ID. If you got read counts for Ensembl transcript IDs, I'm guessing your command has "-i transcript_id" in it. Remove that. If your GTF file doesn't use the attribute "gene_id" for gene IDs (GTFs downloaded from Ensembl should), you need to specify it with the -i argument, e.g. "-i geneID".

ADD COMMENTlink modified 12 months ago • written 12 months ago by colin.kern930

Thanks. I got the point.

ADD REPLYlink written 12 months ago by archana.bioinfo87180

I have one more query. Now as I have Ensembl gene IDs and I got the gene symbols with each Ensembl gene IDs. Also, I found that some Ensembl gene IDs having the same gene symbol. In this case, if I want read count w.r.t gene symbol then what I have to consider.

Any help is much appriciated.

Thanks

ADD REPLYlink written 12 months ago by archana.bioinfo87180

This might answer your question: Why am I getting different ensembl gene ids for a given gene symbol?

ADD REPLYlink written 12 months ago by colin.kern930
2
gravatar for lieven.sterck
12 months ago by
lieven.sterck9.0k
VIB, Ghent, Belgium
lieven.sterck9.0k wrote:

Have a look in the manual of htseq-count and pay especially attention to the -t and -f parameters. Then look in your GFF file and choose the feature tags that correspond to your question and use those when running htseq-count

ADD COMMENTlink written 12 months ago by lieven.sterck9.0k

Thanks for the help.

ADD REPLYlink written 12 months ago by archana.bioinfo87180
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1636 users visited in the last hour