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!