Hi I have a question about my RNA seq data. I use hisat2 to align the reads to hg19 genome, and use featureCounts to create gene count matrix. I just realize if no multimapping is allowed, all of my data has no unique gene mapping for TNF which should be one of the most common genes for tumor research. I assume it should have expressions but somehow I don't understand why there is none. If I choose multimapping in featureCounts, TNF has counts, which means all TNF reads have multimapping, which is not for sure belong to TNF. I wonder if anyone here have similar experience or provide me some suggestions. Thanks.
What is the read length, is this single or paired-end data?
Hi, the read length is 50, it is single end data.