Hic-pro的结果文件转化为.hic文件，在juicebox中实现可视化

2024-06-30 12:54:53

hic数据经过Hic-pro处理后，会生成allvalidpairs文件，这是所有有效配对的文件。一般想要可视化的话，比较复杂。这时我们就可以把它转化为.hic文件，放到juicebox中就很好的可视化。

juicer中的pre命令是用来做这个事情的。只要你的数据符合pre命令处理的格式，这个就很简单。具体用法:https://github.com/aidenlab/juicer/wiki/Pre

但有时我们的Hic-pro的结果文件不符合pre的格式，就需要使用hicpro2juicebox.sh这个脚本进行处理。

1.下载hicpro2juicebox.sh及juicetools(必须)

从github中下载脚本

从github下载juicetools

2.处理allvalidpairs文件（必须）

因为我的文件是从GEO数据库下载的，所以要进行下处理。

先需要解压缩，需要注意的是，解完压缩后，是一个 TXT 文件。所以必须修改文件后缀，把TXT后缀去掉。

比如： allvalidpairs.txt -> allvalidpairs

3.修改hicpro2juicebox.sh（非必须）

在脚本的137行和140行中,有一个参数**–parallel=4**是针对与服务器有GNU的，如果服务器没有GNU的话,就必须要删除这个参数。

awk '{$4=$4!="+"; $7=$7!="+"; n1=split($9, frag1, "_"); n2=split($10, frag2, "_"); } $2<=$5{print $1, $4, $2, $3, frag1[n1], $7, $5, $6, frag2[n2], $11, $12 }$5<$2{ print $1, $7, $5, $6, frag2[n2], $4, $2, $3, frag1[n1], $12, $11}' $VALIDPAIRS | LANG=C sort -T ${TEMP} -k3,3d -k7,7d -S 50% --parallel=4 > ${TEMP}/$$_allValidPairs.pre_juicebox_sortedelseawk '{$4=$4!="+"; $7=$7!="+"} $2<=$5{print $1, $4, $2, $3, 0, $7, $5, $6, 1, $11, $12 }$5<$2{ print $1, $7, $5, $6, 0, $4, $2, $3, 1, $12, $11 }' $VALIDPAIRS | sort -T ${TEMP} -k3,3d  -k7,7d -S 50% --parallel=4 > ${TEMP}/$$_allValidPairs.pre_juicebox_sorted

4.运行命令

bash juicebox.sh -i /slst/home/ningwei/data/CRC/FHC/FHC-1_ALL -j /slst/home/ningwei/package/juicer_tools_1.22.jar -g hg19 -o /slst/home/ningwei/data/CRC/FHC/hicdata

-i 需要进行转化的allvalidpairs文件

-j juicetools

-g 参考基因组

-o 输出目录

Hic-pro的结果文件转化为.hic文件，在juicebox中实现可视化相关推荐

pdf文件转化为wmf文件_如何将PDF文件转换为PowerPoint
pdf文件转化为wmf文件 If you have a PDF file that you would like to convert to a PowerPoint presentation, yo ...
Netlib文件转化为mps文件
Netlib文件转化为mps文件下载地址建立emps.c文件内容如下 /* Expand compressed LP programs (in netlib format) to MPS form ...
jupyter命令把.ipynb文件转化为.py文件
在当前文件夹下面按下Shift + 右键打开cmd终端然后输入: jupyter nbconvert --to script *.ipynb 就能把当前文件夹下面的所有的.ipynb文件转化为.p ...
如何用python修改pdf内容_如何利用python将pdf文件转化为txt文件？
https://www.wukong.com/answer/6579491774144708872/?iid=15906422033&app=news_article&share_an ...
如何实行反编译，将.class文件转化为.java文件
#如何实行反编译,将.class文件转化为.java文件 1.使用win+R快捷键,打开"运行"界面. 2.输入cmd,打开命令窗口. 3.因为我的jad.exe软件在e盘里,所以 ...
如何将Word文件转化为Latex文件：word-to-latex-2.56详解
首先推荐大家读一读这篇博文:教你如何在word中像LaTex那样打公式_ibingow的专栏-CSDN博客 Word转Latex软件word-to-latex 2.56破解版下载地址:百度网盘请输入 ...
java scel_使用java将搜狗词库.scel文件转化为.txt文件
需求:批量将.scel文件转化为可视的txt文件(支持1对1,多对1,多对多),并从中提取中文词(去重),支持追加内容. 成果: 使用: package com.hxl.files; import j ...
python 把txt文件转化为pcd文件并用open3D可视化
python 把txt文件转化为pcd文件并用open3D进行可视化最近在弄3D点云的一些东西,遇到一个问题:很多数据集给的点云数据都是txt格式的,但是在模型上用到的数据集基本是pcd格式!找了很 ...
python快速将excel文件转化为txt文件
在进行数据预处理的时候,我们可能习惯于将用于存储数据的其他格式的文件转换为.txt文件,因为这类文件读取.写入和处理起来都比较高效.这篇博客主要介绍一种快速地将excel文件转化为.txt文件的方法, ...
xps数据怎么导出为txt_WFP: 读取XPS文件或将word、txt文件转化为XPS文件
读取XPS格式文件或将doc,txt文件转化为XPS文件,效果图如下: 1.XAML页面代码: xmlns="http://schemas.microsoft.com/winfx/2006/ ...

最新文章

热门文章