Should PARs be hard-masked in human genome while using htseq-count for downstream analysis?
Hi all, I'm doing RNAseq analysis, and I have aligned reads to hg38 use HISAT2, as we know that htseq-count would discard multi-mapped reads, GENCODE gene moduls contain annotation in PARs, e.g. ENST00000399012.6 & ENST00000399012.6_PAR_Y, so should I hard mask PARs in chrY before alignment?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Why Is My Uniquely Mapped Read Count Higher in Samtools Compared to Bowtie2?
written by daffodil •
-
Alternative gap introduced when comparing insertion/deletion results between minimap2 and STAR2
written by peterhuang319 •
-
Discrepancy in Alignment Rates: HISAT2 vs FeatureCounts
written by Max •
-
Htseq-count vs samtools idxstats, which one to choose?
written by 2822462298 •
-
Why my Hisat2 result is 0% aligned concordantly and HTseq result is 100% not aligned.
written by guotingfeng •
-
How to select uniquely and concordantly reads from hisat2 alingment for raw read count
written by Bioinfonext
-
Question: Error HISAT2 output for HTseq-count
written by 11yj3312 •
-
Different alignment rates for Hisat2 and STAR, Hisat2 has lower alignment rate and STAR have many m…
written by SMILE
-
Error HISAT2 output for HTseq-count
written by nicolas.merienne •
-
The comparison between HISAT2 and Tophat2
written by Peng Huang
As far as I know that is often the case, already. Sounds like a good idea if it isn't.
Thanks @WouterDeCoster , I just wonder whether anybody has done that way before, and I found that reference used in Ilumina basespace APP is hard-masked in chrY PARs, maybe I should mask it before alignment, too. Thanks!