macs2---使用記錄

-- (base) [yangjy@GSCG01 dum]$ macs2 callpeak -c  Control.merge.rmdup.bam -t H2Aub1.merge.rmdup.bam -f BAM -B -g mm -n H2Aub1
INFO  @ Wed, 03 Feb 2021 12:22:13:
# Command line: callpeak -c Control.merge.rmdup.bam -t H2Aub1.merge.rmdup.bam -f BAM -B -g mm -n H2Aub1
# ARGUMENTS LIST:
# name = H2Aub1
# format = BAM
# ChIP-seq file = ['H2Aub1.merge.rmdup.bam']
# control file = ['Control.merge.rmdup.bam']
# effective genome size = 1.87e+09
# band width = 300
# model fold = [5, 50]
# qvalue cutoff = 5.00e-02
# The maximum gap between significant sites is assigned as the read length/tag size.
# The minimum length of peaks is assigned as the predicted fragment length "d".
# Larger dataset will be scaled towards smaller dataset.
# Range for calculating regional lambda is: 1000 bps and 10000 bps
# Broad region calling is off
# Paired-End mode is off

處理進程

INFO  @ Wed, 03 Feb 2021 12:22:13: #1 read tag files...  
INFO  @ Wed, 03 Feb 2021 12:22:13: #1 read treatment tags...    讀取 實驗組
INFO  @ Wed, 03 Feb 2021 12:22:15:  1000000
INFO  @ Wed, 03 Feb 2021 12:22:16:  2000000
INFO  @ Wed, 03 Feb 2021 12:22:18:  3000000
INFO  @ Wed, 03 Feb 2021 12:22:19:  4000000
INFO  @ Wed, 03 Feb 2021 12:22:21:  5000000
INFO  @ Wed, 03 Feb 2021 12:22:22:  6000000
INFO  @ Wed, 03 Feb 2021 12:22:24:  7000000
INFO  @ Wed, 03 Feb 2021 12:22:25:  8000000
INFO  @ Wed, 03 Feb 2021 12:22:27:  9000000
INFO  @ Wed, 03 Feb 2021 12:22:28:  10000000
INFO  @ Wed, 03 Feb 2021 12:22:30:  11000000
INFO  @ Wed, 03 Feb 2021 12:22:32:  12000000
INFO  @ Wed, 03 Feb 2021 12:22:33:  13000000
INFO  @ Wed, 03 Feb 2021 12:22:35:  14000000
INFO  @ Wed, 03 Feb 2021 12:22:36:  15000000
INFO  @ Wed, 03 Feb 2021 12:22:38:  16000000
INFO  @ Wed, 03 Feb 2021 12:22:39:  17000000
INFO  @ Wed, 03 Feb 2021 12:22:41: 17515747 reads have been read.

INFO  @ Wed, 03 Feb 2021 12:22:41: #1.2 read input tags...   *讀取 對照組
INFO  @ Wed, 03 Feb 2021 12:22:42:  1000000
INFO  @ Wed, 03 Feb 2021 12:22:44:  2000000
INFO  @ Wed, 03 Feb 2021 12:22:45:  3000000
INFO  @ Wed, 03 Feb 2021 12:22:47:  4000000
INFO  @ Wed, 03 Feb 2021 12:22:48:  5000000
INFO  @ Wed, 03 Feb 2021 12:22:50:  6000000
INFO  @ Wed, 03 Feb 2021 12:22:51:  7000000
INFO  @ Wed, 03 Feb 2021 12:22:52:  8000000
INFO  @ Wed, 03 Feb 2021 12:22:54:  9000000
INFO  @ Wed, 03 Feb 2021 12:22:55:  10000000
INFO  @ Wed, 03 Feb 2021 12:22:57:  11000000
INFO  @ Wed, 03 Feb 2021 12:22:58:  12000000
INFO  @ Wed, 03 Feb 2021 12:23:01: 12331048 reads have been read.

INFO  @ Wed, 03 Feb 2021 12:23:01: #1 tag size is determined as 49 bps   
INFO  @ Wed, 03 Feb 2021 12:23:01: #1 tag size = 49.0
INFO  @ Wed, 03 Feb 2021 12:23:01: #1  total tags in treatment: 17515747
INFO  @ Wed, 03 Feb 2021 12:23:01: #1 user defined the maximum tags...
INFO  @ Wed, 03 Feb 2021 12:23:01: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s)
INFO  @ Wed, 03 Feb 2021 12:23:01: #1  tags after filtering in treatment: 17515746
INFO  @ Wed, 03 Feb 2021 12:23:01: #1  Redundant rate of treatment: 0.00
INFO  @ Wed, 03 Feb 2021 12:23:01: #1  total tags in control: 12331048
INFO  @ Wed, 03 Feb 2021 12:23:01: #1 user defined the maximum tags...
INFO  @ Wed, 03 Feb 2021 12:23:01: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s)
INFO  @ Wed, 03 Feb 2021 12:23:01: #1  tags after filtering in control: 12331048
INFO  @ Wed, 03 Feb 2021 12:23:01: #1  Redundant rate of control: 0.00
INFO  @ Wed, 03 Feb 2021 12:23:01: #1 finished!   # 完成
INFO  @ Wed, 03 Feb 2021 12:23:01: #2 Build Peak Model... # 建立模型
INFO  @ Wed, 03 Feb 2021 12:23:01: #2 looking for paired plus/minus strand peaks...
INFO  @ Wed, 03 Feb 2021 12:23:03: #2 number of paired peaks: 4584  # 
INFO  @ Wed, 03 Feb 2021 12:23:03: start model_add_line...
INFO  @ Wed, 03 Feb 2021 12:23:03: start X-correlation...
INFO  @ Wed, 03 Feb 2021 12:23:03: end of X-cor
INFO  @ Wed, 03 Feb 2021 12:23:03: #2 finished!
INFO  @ Wed, 03 Feb 2021 12:23:03: #2 predicted fragment length is 83 bps
INFO  @ Wed, 03 Feb 2021 12:23:03: #2 alternative fragment length(s) may be 83 bps
INFO  @ Wed, 03 Feb 2021 12:23:03: #2.2 Generate R script for model : H2Aub1_model.r
WARNING @ Wed, 03 Feb 2021 12:23:03: #2 Since the d (83) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem!
WARNING @ Wed, 03 Feb 2021 12:23:03: #2 You may need to consider one of the other alternative d(s): 83
WARNING @ Wed, 03 Feb 2021 12:23:03: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing.
INFO  @ Wed, 03 Feb 2021 12:23:03: #3 Call peaks...
INFO  @ Wed, 03 Feb 2021 12:23:03: #3 Pre-compute pvalue-qvalue table...
INFO  @ Wed, 03 Feb 2021 12:23:48: #3 In the peak calling step, the following will be performed simultaneously:
INFO  @ Wed, 03 Feb 2021 12:23:48: #3   Write bedGraph files for treatment pileup (after scaling if necessary)... H2Aub1_treat_pileup.bdg
INFO  @ Wed, 03 Feb 2021 12:23:48: #3   Write bedGraph files for control lambda (after scaling if necessary)... H2Aub1_control_lambda.bdg
INFO  @ Wed, 03 Feb 2021 12:23:48: #3   Pileup will be based on sequencing depth in control.
INFO  @ Wed, 03 Feb 2021 12:23:48: #3 Call peaks for each chromosome...
INFO  @ Wed, 03 Feb 2021 12:24:34: #4 Write output xls file... H2Aub1_peaks.xls
INFO  @ Wed, 03 Feb 2021 12:24:34: #4 Write peak in narrowPeak format file... H2Aub1_peaks.narrowPeak
INFO  @ Wed, 03 Feb 2021 12:24:34: #4 Write summits bed file... H2Aub1_summits.bed
INFO  @ Wed, 03 Feb 2021 12:24:34: Done!
[1]+  Done                    macs2 callpeak -c Control.merge.rmdup.bam -t H2Aub1.merge.rmdup.bam -f BAM -B -g mm -n H2Aub1 2> H2Aub1.macs2.log
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末典徊,一起剝皮案震驚了整個濱河市纬霞,隨后出現(xiàn)的幾起案子弓摘,更是在濱河造成了極大的恐慌,老刑警劉巖腕够,帶你破解...
    沈念sama閱讀 206,839評論 6 482
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異舌劳,居然都是意外死亡帚湘,警方通過查閱死者的電腦和手機,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,543評論 2 382
  • 文/潘曉璐 我一進店門蒿囤,熙熙樓的掌柜王于貴愁眉苦臉地迎上來客们,“玉大人,你說我怎么就攤上這事材诽〉状欤” “怎么了?”我有些...
    開封第一講書人閱讀 153,116評論 0 344
  • 文/不壞的土叔 我叫張陵脸侥,是天一觀的道長建邓。 經(jīng)常有香客問我,道長睁枕,這世上最難降的妖魔是什么官边? 我笑而不...
    開封第一講書人閱讀 55,371評論 1 279
  • 正文 為了忘掉前任,我火速辦了婚禮外遇,結(jié)果婚禮上注簿,老公的妹妹穿的比我還像新娘。我一直安慰自己跳仿,他們只是感情好诡渴,可當我...
    茶點故事閱讀 64,384評論 5 374
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著菲语,像睡著了一般妄辩。 火紅的嫁衣襯著肌膚如雪惑灵。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 49,111評論 1 285
  • 那天眼耀,我揣著相機與錄音英支,去河邊找鬼。 笑死哮伟,一個胖子當著我的面吹牛干花,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播,決...
    沈念sama閱讀 38,416評論 3 400
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼惠窄!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起修赞,我...
    開封第一講書人閱讀 37,053評論 0 259
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎桑阶,沒想到半個月后柏副,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,558評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡蚣录,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 36,007評論 2 325
  • 正文 我和宋清朗相戀三年割择,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片萎河。...
    茶點故事閱讀 38,117評論 1 334
  • 序言:一個原本活蹦亂跳的男人離奇死亡荔泳,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出虐杯,到底是詐尸還是另有隱情玛歌,我是刑警寧澤,帶...
    沈念sama閱讀 33,756評論 4 324
  • 正文 年R本政府宣布擎椰,位于F島的核電站支子,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏达舒。R本人自食惡果不足惜值朋,卻給世界環(huán)境...
    茶點故事閱讀 39,324評論 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望巩搏。 院中可真熱鬧昨登,春花似錦、人聲如沸贯底。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,315評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至糯俗,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間睦擂,已是汗流浹背得湘。 一陣腳步聲響...
    開封第一講書人閱讀 31,539評論 1 262
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留顿仇,地道東北人淘正。 一個月前我還...
    沈念sama閱讀 45,578評論 2 355
  • 正文 我出身青樓,卻偏偏與公主長得像臼闻,于是被迫代替她去往敵國和親鸿吆。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 42,877評論 2 345

推薦閱讀更多精彩內(nèi)容