Question: reasons for htseq-count with zero features
gravatar for shenwei1376
2.3 years ago by
shenwei13760 wrote:

Hi everyone,

I have asked this before, but want specify more this time. I have been using STAR-Htseq analyze the RNASeq data. Everything was find until recently, one of the library gives me 0 counts after the htseq-count. All the other libraries are fine, which means the gtf file or pipeline is find. I am wondering what could be reasons for 0 zero counts?

Many thanks for any discussion!


Wei S

rna-seq htseq-count • 959 views
ADD COMMENTlink modified 2.3 years ago by genomax89k • written 2.3 years ago by shenwei13760

I think auto-correct changes a bunch of fines to finds. Can you edit your post and correct those please?

ADD REPLYlink written 2.3 years ago by RamRS30k

Duplicate of htseq-count with zero features.

ADD REPLYlink written 2.2 years ago by h.mon31k

shenwei1376 , please do not open duplicate questions, especially when queries on the original are yet to be addressed. Doing so is rude to people helping you in the other post and does not make you look good.

ADD REPLYlink modified 2.2 years ago • written 2.2 years ago by RamRS30k
gravatar for genomax
2.3 years ago by
United States
genomax89k wrote:

Have you looked at the alignment report for that sample and if so what does it say. Perhaps you received data that is not from your sample or the sample may have failed during library prep and may need to be repeated. Either possibility will lead to 0 counts. Take a few reads and blast them at NCBI to discount first possibility. For the latter you will need to dig into QC and/or talk with your sequence provider.

Note: This is based on your assertion that rest of the samples have worked well and produce reasonable results with the pipeline you are using.

ADD COMMENTlink modified 2.3 years ago • written 2.3 years ago by genomax89k

Thank you for your suggestions!! The library has been mapped to two indexed genomes. One is indexed based on a gtf file that is re-structured, so all the genes are re-grouped into transcripts (protein informs). Another one is indexed based on regular gtf file. The strange thing is all libraries worked when mapped to first genome, but the last library failed when mapped to the second genome! The pipeline and settings are same!

Thank you for more thoughts!

ADD REPLYlink written 2.3 years ago by shenwei13760
Please log in to add an answer.


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