關於原理部分和更詳細的介紹,見hic-pro: hi-c資料預處理高效工具, 這裡只介紹如何快速使用singularity的hic-pro進行資料分析。
關鍵內容就是,config-hicpro.txt 裡的檔案路徑資訊都必須是絕對路徑,否則預設都位於annotation目錄下。切記,切記,切記。
mkdir -p /opt/biosoft/hic-pro
cd /opt/biosoft/hic-pro
wget
# 使用
singularity exec /opt/biosoft/hic-pro/hicpro_latest_ubuntu.img hic-pro -h
第一步:建立輸入資料資料夾, 資料不能是軟連線形式,只能是複製或者移動
mkdir -p fastq/***
cp ***_r1.fastq.gz ***_r2.fastq.gz fastq/***/
# 如果有多個樣本
mkdir -p fastq/yyy
cp yyy_r1.fastq.gz yyy_r2.fastq.gz fastq/yyy/
enzyme=dpnii
# build reference
mkdir reference
mv 你的參考序列.fasta reference/genome.fa
# enzyme site
singularity exec /opt/biosoft/hic-pro/hicpro_latest_ubuntu.img \
/usr/local/bin/hic-pro_2.11.4/bin/utils/digest_genome.py -r $enzyme \
-o reference/genome_$.bed reference/genome.fa
# genome size
seqkit fx2tab -nl reference/genome.fa | awk '' > reference/genome.chrom.size
# bowtie/2.3.4.3
singularity exec /opt/biosoft/hic-pro/hicpro_latest_ubuntu.img \
bowtie2-build --threads 60 reference/genome.fa reference/genome
第三步: 複製配置檔案並修改
# config
singularity exec /opt/biosoft/hic-pro/hicpro_latest_ubuntu.img \
cp /usr/local/bin/hic-pro_2.11.4/config-hicpro.txt config-hicpro.txt
修改其中的如下項,
第四步:執行
singularity exec /opt/biosoft/hic-pro/hicpro_latest_ubuntu.img hic-pro -i fastq -o results -c config-hicpro.txt ```
Singular和SciTE的結合
singular是乙個針對多項式系統的計算機代數系統,特別強調交換代數,非交換代數,代數幾何和奇異值理論.在windows系統上需要借助於cygwin執行.軟體執行速度不錯,使用起來也很方便.這裡簡要給出如何在scite中編輯並執行singular程式.1 安裝singular和scite 2 開啟...
conda 安裝指定路徑 HiCPro安裝與使用
hic pro可以用來處理hi c資料,從原始的fastq檔案 illumina資料 到標準化的互動圖譜。簡單的來說就是將hi c資料比對到拼裝好的參考基因組上,並形成互動檔案去儲存hi c資料。下面我們就來介紹一下hic pro的安裝與使用。python版本為2.7 conda create n ...
container of 的的的原理
另外一篇,同樣精彩,揭開linux核心中container of的神秘面紗 華清遠見嵌入式學院講師。在linux 核心中有乙個大名鼎鼎的巨集container of 這個巨集是用來幹嘛的呢?我們先來看看它在核心中是怎樣定義的。呵呵,乍一看不知道是什麼東東。我們先來分析一下container of p...