IDR on blacklisted region-filtered narrowPeak
1
0
Entering edit mode
4.1 years ago

Hi guys,

I am having a problem performing IDR on my blacklist region-filtered file.

After I called my narrow peaks using MACS2, I used bedtools to filter the blacklisted regions from my narrowPeak file. I then used IDR these filtered files to get the IDR scores and peaks. However, whenever I run my idr code, it gives me an error message:

UnicodeDecodeError: 'utf-8' codec can't decode byte 0x8b in position 1: invalid start byte

I think there is something wrong with the filtered files, since using unfiltered files as my IDR input gives me exactly what I expected.

This is my code for getting rid of the blacklisted region:

for i in Sample*.narrowPeak
do
  bedtools intersect -v -a ${i} -b mm10.blacklist.bed.gz | gzip > filtered_${i}
done

Thanks!

ChIP-Seq ATAC-Seq IDR • 1.6k views
ADD COMMENT
1
Entering edit mode
4.1 years ago
h.mon 34k

It seems you are trying to input a gzipped file to a software that doesn't accept it. Why do't you try:

  bedtools intersect -v -a ${i} -b mm10.blacklist.bed.gz > filtered_${i}
ADD COMMENT
0
Entering edit mode

Thank you so much! This fixed it. Cheers!

ADD REPLY

Login before adding your answer.

Traffic: 570 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6