Entering edit mode
3.2 years ago
ElCascador
▴
30
I'm using this script from ATAqc to compute my TSS enrichment score for ATAC-seq data. But my values are pretty high, around 20+.
I was wondering if this script, using the metaseq package was not counting paired end reads as 2 fragments ? I will run some more tests on single end data to see the difference but I can only do that later.