hichipper
hichipper copied to clipboard
Provided hg19 restriction fragment file includes 0bp entry
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?