Open drkmj opened 5 years ago
Hi,
I am using RSEM v1.3.1 to generate count matrix files for DEG analysis. I ran:
RSEM-1.3.1/rsem-calculate-expression \ ~/CV/1_8.fq.gz \ ~/CV_genome/CV \ ~/CV/RSEM/1_8
the resulting count file looks like:
12815220 32705307 881233 46401760 15576149 17129158 17238217 126285194 1
From which I understand there were 32,705,307 alignable reads.
However, when I load the resulting 1_8.genes.results file into R and sum the total number of reads I only have 30,176,704 reads for that sample.
Any idea why this value would be lower than the number of alignable reads?
Thank you, Kevin
Hi,
I am using RSEM v1.3.1 to generate count matrix files for DEG analysis. I ran:
RSEM-1.3.1/rsem-calculate-expression \ ~/CV/1_8.fq.gz \ ~/CV_genome/CV \ ~/CV/RSEM/1_8
the resulting count file looks like:
12815220 32705307 881233 46401760 15576149 17129158 17238217 126285194 1
From which I understand there were 32,705,307 alignable reads.
However, when I load the resulting 1_8.genes.results file into R and sum the total number of reads I only have 30,176,704 reads for that sample.
Any idea why this value would be lower than the number of alignable reads?
Thank you, Kevin