Cuffquant / Cuffnorm - Why are some feature counts exactly the same?
0
0
Entering edit mode
7.5 years ago

Hi,

I'm getting some weird counts after aligning RNA-seq to some bacterial genomes with bowtie2 and then counting / normalizing with cuffquant / cuffnorm respectively. I'm getting something like this with some of my output:

Feature A - Species 1   4440.79
Feature B - Species 2   4440.79
Feature C - Species 3   4440.79

What I'm assuming is that Feature A, B, and C are identical so cuffquant is evenly distributing the count to each of the species. Am I correct?

As an aside: I do see that when I use something like htseq-count on this data set that I get a LOT of ambiguous reads (which I know it discards when it can't be sure of a unique alignment to one spot in the (meta)genome). This reinforces my assumption above.

Many thanks!

RNA-Seq cuffquant cuffnorm htseq-count • 2.2k views
ADD COMMENT

Login before adding your answer.

Traffic: 2735 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6