Using multimaping reads or unique reads on featurecounts?
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 only 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-4% of reads counted in each sample, approximately, comparing if I only count unique reads. 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.
• 38 views
Read more here: Source link