HiTC problem exporting to my5C
0
0
Entering edit mode
@sergioespeso-gil-6997
Last seen 2.2 years ago
New York

Hi community!

I am not sure if I miss a dependency, but when I am trying to export HiC data mapped to hg38 to a my5C file, even though I add << genome="hg38" >> the headers say: 1|mm9|chr1:1-100000

Is this a bug? Or I am missing a dependency?

x<-importC("hESC100000.matrix", xgi.bed="hESC100000_abs.bed") show(x) export.my5C(x, file="hESC.export5C", genome="hg38", per.chromosome=FALSE)

Thanks!!!!

S.

HiTC HiC-Pro • 369 views
0
Entering edit mode

Hi, Could you please specify which HiTC version you are using ? Thanks. Nicolas

0
Entering edit mode

Hi Nicolas,

I am using R version 3.5.1 and HiTC version '1.26.0'

Should I upgrade to version 1.30?

Much appreciated

0
Entering edit mode

Indeed, I think there is a small bug if the input object is an HTClist ... sorry for that. To move forward, a simple sed on the my5C file should allow you to change all mm9, to hg18. I will try to fix this as soon as possible. Best

0
Entering edit mode

No problem, all good.

Thanks a lot Nicolas!

0
Entering edit mode

Indeed, I think there is a small bug if the input object is an HTClist ... sorry for that. To move forward, a simple sed on the my5C file should allow you to change all mm9, to hg18. I will try to fix this as soon as possible. Best

0
Entering edit mode

Indeed, I think there is a small bug if the input object is an HTClist ... sorry for that. To move forward, a simple sed on the my5C file should allow you to change all mm9, to hg18. I will try to fix this as soon as possible. Best