CEGRcode/scriptmanager

ChIP-seq TagPileup Data Cutoff and Full Fragment issue

Closed this issue · 2 comments

owlang commented

In testing out ChIP-seq analysis with ENCODE data, we found a couple of issues possibly relating to the nature of the longer insert sizes of the ChIP-seq data.

Pileups of both 5prime and full-fragment settings result in the appearance of missing data:
Screen Shot 2023-02-21 at 3 27 39 PMScreen Shot 2023-02-21 at 3 27 22 PM

Data accessions: ENCFF865UPM and ENCFF014UUB

I suspect the padding window set when TagPileup retrieves SamRecords from the BAM file is not large enough for the 5' end issue. Cannot comment yet on what is causing the full fragment issue and will need to test for differences with known working datasets

owlang commented

Fragment issue: Dataset was single-end 🙄

owlang commented

BED reference file was irregularly sized from the ENCODE liftOver to hg38 🙄

Not a bug. Be sure to add these warnings for single-end/ENCODE analysis in documentation. Closing ticket.