day22 ChIP-seq 用MACS2找peak

這兩日出現(xiàn)了很多次core dump,表現(xiàn)為qstat正常,但是在e文件中出現(xiàn)報(bào)錯(cuò),且出現(xiàn)core文件。很大可能是因?yàn)?$LD_LIBRARY_PATH的錯(cuò)誤。

一、安裝MACS2

打算用pip安裝。pip install MACS2
注釋?zhuān)哼@是第三個(gè)用pip安裝的軟件了,之前還有cutadapt和multiqc。但是失敗了。。
ERROR: Command errored out with exit status 1:
command: /data/software/anaconda3-202007/bin/python /data/software/anaconda3-202007/lib/python3.8/site-packages/pip install --ignore-installed --no-user --prefix /tmp/pip-build-env-5kfedc4c/overlay --no-warn-script-location --no-binary :none: --only-binary :none: -i https://pypi.org/simple -- pip 'setuptools>=41.2' 'numpy>=1.17' 'cython>=0.29' wheel
cwd: None

下載到本地電腦
https://files.pythonhosted.org/packages/e2/61/85d30ecdd34525113e28cb0c5a9f393f93578165f8d848be5925c0208dfb/MACS2-2.2.7.1.tar.gz
再上傳到服務(wù)器
運(yùn)行:
tar -zxf MACS2-2.2.7.1.tar.gz
python setup.py install

image.png

于是增加了安裝目錄:
python setup.py install --prefix=/data/zds209/.local
因?yàn)榘惭bcupadapt的時(shí)候已經(jīng)把這個(gè)路徑:~/.local/bin設(shè)定為環(huán)境變量里啦。所以這次直接就有啦,不需要再添加path了。

二、MACS2使用

project=/data/zds209/ChIP-seqtest
ls project/bam/*.sorted.bam | while read id do macs2 callpeak -tid -c project/bam/control.sorted.bam -f BAM -B -g mm -nid --outdir $project/bed/

-g是基因組可以是mm,或hs
-n輸出文件前綴
-B 輸出bgd文件,下游bigwig文件生成所需
--nodel 這個(gè)參數(shù),用于雙端測(cè)序數(shù)據(jù)
輸出文件有這些:# NAME 是使用 -n 設(shè)置的輸出文件前綴

  1. NAME_control_lambda.bdg 和 NAME_treat_pileup.bdg
  • bdg即為bedGraph格式,可以導(dǎo)入U(xiǎn)CSC或者轉(zhuǎn)換為bigwig格式。
    -兩個(gè)bdg文件:
    • treat_pileup
    • control_lambda
  1. NAME_model.r # R代碼

  2. NAME_peaks.narrowPeak # BED 6+4 格式

  3. NAME_peaks.xls # 包含 peak 信息的 xls 文件,注意 chrStart 是1-based

  4. NAME_peaks_summits.bed

  • BED 格式,包含peak的summits位置,第5列是-log10pvalue。
    -如果想找motif,推薦使用此文件。

大神的教程 https://zhuanlan.zhihu.com/p/90180058,學(xué)習(xí)ing

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡(jiǎn)書(shū)系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容