macs2 install & 使い方

スポンサーリンク

MACS2

2021.12.222018.04.18

MACS2、良く使われている ChIP-Seq 解析用の peak caller です。

とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談

目次

macs2
論文・データ解析の支援
関連記事

macs2

色々なサイトで、macs2 のインストール方法や使い方がまとめられていますが、
バージョンの違いやインストール方法の好みによって、バラバラです。

もっとも簡単だと思われる方法をまとめました。

前提条件 – Prerequisites

インストールの際、以下の環境が必要になります。

・Python version must be equal to 2.7 to run MACS. （推奨バージョン：version 2.7.9.）
・Numpy (>=1.6)
・Cython (>=0.18)

INSTALL Guide For MACS

macs2 install

Ubuntu Server 18.04 LTS でインストールを確認しました。

update & upgrade

sudo apt update
sudo apt upgrade

Python 2.7.15

sudo apt install python

Python pip

sudo apt install python-pip

numpy

pip install numpy --user
#Successfully installed numpy-1.16.5

cython

pip install cython --user
#Successfully installed cython-0.29.13

MACS2

pip install MACS2 --user
# Successfully installed MACS2-2.1.2 numpy-1.16.5

pathの調整

$ export PATH=~/.local/bin/:$PATH

MACS2 使い方

MACS2 でピークコールを行います。ターゲットとなる処理群の BAM ファイルと、コントロールの BAM ファイルを引数で指定します。
他に、ファイル形式やゲノムサイズ、識別名を指定します。オプションで、ピークコールの際のしきい値を指定することができます。

MACS2 使い方 callpeak

$ macs2 callpeak -t target.bam -c control.bam -f BAM -g hs -n name --outdir outdir_name 2> log.txt

$ macs2 callpeak -t target1.bam target2.bam target3.bam -c control1.bam control2.bam control3.bam -f BAM -g hs -n name --outdir outdir_name 2> log.txt

MACS2 使い方 callpeak よく使うオプション

・-t：ターゲットの bam ファイル
・-c：コントロールの bam ファイル
・-f：入力ファイルの形式
・-g：ゲノムサイズを指定（生物名で指定する場合、ヒトゲノムは hs, マウスゲノム mm など）
・-n：任意の識別子（解析の名前）
・-q, -p：q-value や p-value の閾値。（但し、両方同時には指定できない。）

README for MACS (2.1.2) & MACS2 Usage

論文・データ解析の支援

個別に直接相談できるのが良いと、たくさんの方にご利用頂いています。

とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談

とにかく早く問題解決したい人はこちら＞＞直接、論文相談

関連記事

IGV 使い方遺伝子位置や好きな位置への移動

IGVはリードのマッピング状況を視覚的に確認できるツールです。ゲノムデータ（fastaファイル）とアライメントデータ（bamファイル、及びbaiファイル）を準備すれば簡単に表示できますが、更に遺伝子の位置情報を読み込むともっと使いやすくな...

IGV 使い方インストール〜便利な使い方まで | リファレンス・マッピングデータ・アノテーションを読み込んで表示しよう

IGV（Integrative Genomics Viewer）は、マッピングデータを視覚的に確認することができる大変便利なツールです。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談IGVとはIGVは...

IGV 使い方 bowtieのマッピングデータをsamtoolsで変換してIGVに表示するまで

IGVを使って、bowtieのマッピング結果表示してみます。samtoolsを使ったbamへの変換やsortなど、IGVの表示に必要なステップを１つずつ確認しながら説明します。とにかく早く問題解決したい人はこちら＞＞直接、データ解...

マッピング率が低いときの対応方法 - Trimmomatic 除去アダプターの登録方法

マッピング率が極端に低い場合は、使用しているシーケンサ特有のアダプター配列がリードに残っていないか確認するとうまくいくケースがあります。特に、ライブラリーサイズが小さくなった場合は、シーケンサ内部のソフトで特有の配列が除去しきれず...

fastqc memory error - java実行時により大きなメモリを割り当てる方法

bamファイルに対してfastqcを実行したら、メモリが足りないというエラーが発生しました。実行環境は十分にメモリーを積んでいたので、fastqc実行時のjavaのメモリーオプションを強制的に書き換えて実行したら、正常終了しました。...

FastQC 使い方 Fastq リードの評価

バイオインフォ道場、くまぞうです。バイオインフォマティクスの解析を行うにあたって、次世代シーケンサから得られたリード情報の品質をチェックすることはとても大事なことです。代表的なツールFastQCの使い方をまとめます。とにか...

FastQC 使い方コマンドラインでテキスト形式の結果を活用する

バイオインフォマティクスの解析を行う際、一番最初に行うのはリードの品質チェックです。代表的なツールはFastQCで結果を視覚的に確認できて非常に便利ですが、たくさんデータを使う場合はコマンドライン（command line）を活用して結果...

FastQC - インストール

FastQCは、NGSのリードの品質をチェックするツールです。解析の直前で使用し、リードの品質に問題がなければ解析に進みます。もし、品質に問題があるようであれば、低品質の情報をトリミングツールで除去したり、対象となるデータを解析から外すな...

blast 使い方完全一致の連続領域を探索するには？

blastn で完全一致の連続領域を探索する方法blastn オプション word_sizeblastnのパラメータ-word_sizeは、完全一致の初期長さを設定します。長さを指定すると、「完全一致が連続する場所」で「指...

blast 使い方エラーと解決 tblastx

バイオインフォマティクスの解析では、様々なツールを使います。たまにはエラーに遭遇してしまい、解決するのに時間がかかったりします。BLASTは、大変使用頻度の高いツールの１つですが、見慣れないエラーに遭遇したので、エラーと対応方法をまとめま...

blast 使い方 makeblastdb

blastを実行するには、事前にデータベースを構築する必要があります。データベースの対象となるfastaファイルを準備し、makeblastdbコマンドを実行してデータベースを作ります。とにかく早く問題解決したい人はこちら＞＞直接...

blast+ 使い方 best hitの算出 awkとoutfmt7

blastはよく使うツールの１つです。色んな機能があるし、オプションもたくさんあるので使いこなすのが大変です。よく苦労するのが「best hitを１個だけ取り出す」ことですが、実は「outfmt7」を使うと簡単に得ることができます！使用し...

blast 使い方インストール blast+ v2.3.0

今回はblast+のインストールに挑戦です。プリコンパイル版です。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談インストールマニュアルに従い、プリコンパイル版をインストールします。場所# c...

シェルスクリプト FASTAファイルの分割 fastx_toolkit, sed, コマンドの活用例

バイオインフォマティクスでは、大きなファイルを扱う機会が多いです。小さなファイルに分割すれば、サーバーで分散処理して大幅な時間短縮を行うことが出来るかもしれません。いくつかのツールを組合せてFASTAファイルを小さく分割してみます。...

bowtie メモリー関連のエラー Warning: Exhausted best-first chunk memory for read

bowtieを--bestモードで実行したところ、メモリーに関するワーニング（警告）が大量に表示されました。エラーではないので処理は進んでいくのですが、「マッピング結果に影響を与えないか？」という点がとても気になりました。そこで、--ch...

bowtie 使い方マッピングとオプション

bowtie は、代表的なマッピングソフトです。高速にマッピング処理を行うことが出来ます。ただし、ギャップは考慮しません。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談Bowtie - 使い方まずは、...

bowtie2 使い方アライメント情報（alignment summary）をファイルに保存する

bowtie2のアライメント情報（alignment summary）をファイル保存する方法について考えます。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談bowtie2 アライメント情報と標準エラー出力...

bowtie v1.1.2 - インストール

今回はbowtie v1.1.2のインストールに挑戦です。プリコンパイル版です。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談インストール場所cd /usr/localダウンロードwget 展...

bowtie2 v2.2.6 使い方インストール

今回はbowtie2 v2.2.6のインストールに挑戦です。プリコンパイル版です。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談インストールマニュアルに従い、プリコンパイル版をインストールします。...

gffread を使った transcripts fasta 転写物の配列取得

ゲノム配列から転写物の配列を取得転写物（transcript）の配列を取得したい場合、染色体上の位置（開始位置〜終了位置）をたよりに配列情報を抜き出します。位置情報を管理する方法としてよく使われるのがGFF3やGTF形式のファイル...

VCFtools 使い方 VCFファイルからの抽出（BEDファイル対応領域）

SNP Callerには、通常、領域指定オプションがあります。しかし、既に算出されたVCFデータが領域外のデータを含む場合、領域を参照しながら整理し直すのは面倒な作業です。そんなとき、vcftoolsを使えば、計算した後からでも簡単に領域...

tophat2 使い方エラーと対策

tophat2でマッピングしていたらエラー終了しました。一応解決したので報告します。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談Tophat2 マッピングで発生したエラーtophat2のマッピングを...

tophat2 使い方 RNA-Seq マッピング

今回は、tophat2を使って、RNA-seqデータのマッピングに挑戦です。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談tophat2とはtophat2は、代表的なマッピングツールの１つです。RNA...

tophat2 v2.1.0 - インストール

今回は、Tophat2 v2.1.0 のインストールに挑戦です。プリコンパイル版です。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談準備Tophat2が内部で利用する解析ツールを事前にインストールしま...

SRA_Toolkit 使い方データ保存フォルダの変更

SRA Toolkiltのprefetchでデータをダウンロードする際の保存フォルダは変更可能です。vdb-configを使って設定します。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談SRA_toolkit...

SRA Toolkit 使い方公開データのダウンロードとsra fastq変換

公開されているFastqデータを活用したい場合は、SRA_toolkitを使うと便利です。データのダウンロードやFastqへの展開がコマンド１つで簡単に実行できます。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談...

UCSCのゲノム配列ファイル .2bitファイルの展開

UCSCのゲノム配列ファイルは、.2bit形式で圧縮されているものがあります。「twoBitToFa」というツールを使って展開します。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談.2bit形式について...

fastx_toolkit 使い方 FastqからFastaへ変換

Fastq・Fastaファイルともに、バイオインフォマティクスで良く使うファイル形式です。FastqからFastaへ変換についてまとめます。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談Fas...

FASTX-Toolkit 使い方インストール

FASTX-Toolkitは、次世代シーケンサのFASTA・FASTQの前処理に関連するツールを集めたものです。トリミングやクリッピング、クオリティ処理に関するツールが提供されています。インストール方法について手順を整理し...

samtools 使い方 faidx Fastaから部分配列

塩基配列の変異に注目したSNP解析や配列パターンに注目した解析などでは、部分配列を抜き出して詳しく調べる場合があります。今回は、開始位置-終了位置を指定してFASTAファイルから注目する塩基配列の一部（部分配列）を取り出します。...

samtools 使い方 mpileup ( calling SNPs ) & annotation

samtoolsを使ったVariant callingについてまとめます。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談解析作業の流れマッピングCalling SNPsフィルタリングアノテーション付与...

samtools 使い方インストール

SAM tools のインストールについて、バージョン毎にまとめました。手順は同じ。指定バージョンを探すのが面倒なので、それぞれ整理しました。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談インストール...

hisat2 使い方 | index 作成

HISAT2のindexファイルを、hisat2-buildコマンドで自ら生成します。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談HISAT2 pre-built index以下の生物種については、...

HISAT2 使い方インストールとマッピング

HISAT2 インストールHISAT2のサイトにインストールが簡単なバイナリー版が準備してあるので、ダウンロードして使います。使用したのは、「Linux x86_64 binary」版です。HISAT2HI...

SnpEff 使い方アノテーション

samtoolsなどで予測されたVariantに、snpEffを使ってアノテーションをつけてみます。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談snpEffsnpEffとは、予測されたVariant...

Miniconda 使い方 | インストールと環境作成・環境利用

Miniconda install公式サイトから、最新版 Miniconda installer（各OS毎）のシェルファイルをダウンロードします。Miniconda$ wget $ bash Miniconda3...

Trinity 使い方 RNA-seq de novo アセンブル

Trinityは、トランスクリプトーム解析で用いられる「de novo アセンブルツール」です。たくさんのメモリとマルチコアの環境を要求します。100万のペアエンドリードに対して1GB程度のメモリが推奨されています。とにかく早く問...

Trinity Docker 面倒な「install作業なし」で一発起動！

Trinityを、一番「簡単に動かす方法」だと思います。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談Trinity in Docker 使い方アセンブルツール Trinity（de novo tra...

Trinity v2.2.0, v2.4.0 インストール

Trinity v2.2.0、v2.4.0のインストールに挑戦です。ソースダウンロードして、make作業を行なってインストールします。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談インストール...

Trinotate インストール | conda/biocondaを使って簡単インストール

trinotateは、trinityのアセンブル結果にアノテーションをつけるツールです。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談trinotate installTrinotate便利な...

fastp インストール・簡単な使い方・よく使うオプション

Fastqファイル用の高速前処理ツールです。クオリティの評価（トリミング前後）やトリミング処理を行います。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談fastp とは？fastp : fas...

freebayes インストール簡単な使い方

freebayesは、代表的なSNP検出ツールの１つです。ベイズ的アプローチでIndelを検出します。インストールと簡単な使い方についてまとめます。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談fre...

htseq-count 使い方 gene単位・transcript単位の数え方

RNA-Seq解析では、ゲノムに張りついたリードの数を数えます。gene単位・transcript単位の数え方は、オプションで切り替えることができます。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談h...

htseq-count 使い方よく使うオプション

htseq-countは、RNA-Seq解析などでゲノムに張りついたリードの数を数えるツールです。よく使うオプションをまとめます。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談htseq-count ...

HTSeq htseq-count 使い方インストール v0.7.2

htseq-countは、RNA-Seq解析などでリードをゲノムにマッピングしたとき、どこにどのくらい張り付いたかをカウントするプログラムです。HTSeqのインストールには、いくつか事前にインストールしておくべきものがあります。...

STAR RNA-seq aligner 使い方インストール・index・マッピング

STARは、RNA-Seq用のマッピングソフトです。非常に高速なマッピングが可能です。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談STARとは？STAR インストールダウンロ...

塩基配列の表示幅を調整する - Fastaファイルの整形 -

Fastaファイルの整形についてまとめます。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談フォーマットを整えるFASTX_Toolkitは、Fastq・Fasta形式のファイルについて前処理をしたり、...

Fastq Fasta ファイル形式

Fastq・Fastaファイルともに、良く使うファイル形式です。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談ファイル形式Fastq形式1本の配列は4行構成で記述されます。１行目：「@...

タイトルとURLをコピーしました