nservant / HiC-Pro

HiC-Pro: An optimized and flexible pipeline for Hi-C data processing
Other
372 stars 181 forks source link

Question about manually sorting the '.validPairs'. #580

Closed Miracle-Yao closed 3 months ago

Miracle-Yao commented 1 year ago

Hi all, I have a question about manually sorting the '.validPairs'. I found that after running the 'mapped_2hic_fragments.py ' script, there is no sorting of the '.validPairs', and I ran the code in the log I got from the previous process and found that it doesn't seem to react. Is there any way to fix this? The code I used is listed below.

LANG=en; sort -T tmp -k2,2V -k3,3n -k5,5V -k6,6n -o CTCF_HiChIP.valid_hg38.bwt2pairs.sorted.validPairs CTCF_HiChIP.valid_hg38.bwt2pairs.validPairs

The current directory contains all the files or tmp directories needed for this command, and there are roughly 70,000,000 lines in '.validPairs'. I'm using HiCPro version 2.11.4. Because '-s proc_hic' had some problems running, I switched to running all commands manually.

Thanks

Miracle-Yao commented 1 year ago

All right. More patience.