I'm working with Rat transcriptome (mRNA) using HISAT as aligner and featurecounts (subread) to count reads using BAM files from HISAT. Featurecounts has the possibility to count onli unique reads or multi-mapping reads.
What is the best practice, taking into account multi-mapping reads or not?
When I do the counting process with multi-mapping reads it increases by 1% of reads counted in each sample approximately comparing if I only count unique reads. The bad part is that the count table changed dramatically getting lost differential expression genes between control and experimental samples during DESeq2 Analsisis. So, I'm not sure what is better from the point of view of what option I have to choose, using multi-mapping reads or only unique reads.