I wanted to use diffHic just to perform peak calling on Hi-C data. Do I necessarily have to start from fastq alignment with the script included in the package or I can use an alternative input (e.g. HiC-Pro filtered read pairs)? If so, which function can I use and how have I to format the input file(s)? Thanks in advance.