hichipper icon indicating copy to clipboard operation
hichipper copied to clipboard

Provided hg19 restriction fragment file includes 0bp entry

Open afrendeiro opened this issue 6 years ago • 1 comments

Maybe not an actual problem when using hichipper, but I was just using the provided hg19 MboI restriction fragment files and found something weird in it. It includes a 0-bp entry on chrM.

> wget -O hg19_MboI_resfrag.bed.gz https://github.com/aryeelab/hichipper/blob/master/RestrictionFragmentFiles/hg19_MboI_resfrag.bed.gz?raw=true
> gzip -d hg19_MboI_resfrag.bed.gz
> bedtools sort -i hg19_MboI_resfrag.bed > hg19_MboI_resfrag.sorted.bed
Error: malformed BED entry at line 7127610. Start was greater than end. Exiting.

In that line:

chrM    0       0
chrM    0       741
chrM    741     952
chrM    952     1228
chrM    1228    2350

With HiC-Pro 2.10.0 digest_genome.py:

chrM    0       741     HIC_chrM_1      0       +
chrM    741     952     HIC_chrM_2      0       +
chrM    952     1228    HIC_chrM_3      0       +
chrM    1228    2350    HIC_chrM_4      0       +

Interestingly, this is not true for hg38. Were the files produced with HiC-Pro and if so, with which version?

afrendeiro avatar Jun 14 '18 09:06 afrendeiro