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

IGVを使って、bowtieのマッピング結果表示してみます。samtoolsを使ったbamへの変換やsortなど、IGVの表示に必要なステップを1つずつ確認しながら説明します。

スポンサーリンク



bowtieでマッピング

bowtieでマッピングします。マッピング結果をSAM形式で取得します。

bowtie -S -p 4 ref_index -1 read1.fastq -2 read2.fastq mapping.sam

SAMtoolsでファイル変換

IGVのサイトによると、SAMよりもBAM形式の方がオススメと書いてあるので、SAMtoolsを使ってBAMファイルを準備します。また、IGVで読込むBAMファイルはソート・インデックスが必要なので、対応したBAMファイルを作成します。SAMtoolsは、BAMフォーマットのアライメント結果を操作するためのツール群です。

1. samtools viewでファイル変換

まずは、SAMtools viewで、SAM形式からBAM形式へのフォーマット変換を行います。

  • SAM -> BAM
    samtools view -bS mapping.sam > mapping.bam
  • BAM -> SAM(ちなみに)
    samtools view -h mapping.bam > mapping.sam

2. samtools sortでソート処理

BAMファイルのアライメントを左端の座標でソートして並べます。

samtools sort mapping.bam mapping.sort

3. samtools indexでindex付与

ソート済みのBAMファイルにインデックス付与します。コマンド実行後、mapping.sort.bam.baiファイルが作成されます。IGVでの読込みの際、名前がbamファイルと関連付けられている必要があるので、安易に名前を変更しないほうが良いです。

samtools index mapping.sort.bam

IGVでマッピング結果を表示する

ソート・インデックス対応したBAMファイルがあれば、IGVでのマッピング表示は簡単です。(操作の詳細についてまとめました。

  • IGVを起動
  • リファレンス配列の読み込み
    メニューの「Genomes」->「Load Genome from File…」で、Fastaファイルを指定します。
  • マッピング結果の読み込み(BAMファイル)
    メニューの「File」->「Load from File…」で、ソート・インデックス対応したBAMファイルを指定します。尚、読み込みフォルダにbaiファイルを置いておきます。

IGV2

スポンサーリンク





ピックアップ記事

  1. R subset関数 データフレームやmatrixからの条件指定による行・列の抽出

    R の subset関数は、データフレームやマトリックスから条件にマッチした部分集合を取り出します…
  2. R plot 重ねる方法3パターン サンプルでわかるRの使い方

    Rでグラフ (plot) を重ねる方法は、「単純な追加」「図に重ねて描画」「濃淡で重なり表現」の3…
  3. awk split サンプルでわかる列の分割とセパレータ(文字・正規表現)の指定方法

    awkのsplit関数は、セパレータを指定して文字列を分割することができます。セパレータには文字や…

人気記事

  1. IGV, 解析ツール

    IGV 使い方 インストール〜便利な使い方まで | リファレンス・マッピングデータ・アノテーションを読み込んで表示しよう
    IGV(Integrative Genomics View…
  2. Excel, その他, 統計

    z-score 計算方法 エクセル(Excel) 編
    統計処理で、大きく変化しているなどの判断基準にも使われる値…
  3. R データ型 - 文字列・ベクター・データフレーム・マトリックス など-, R言語, スクリプト

    R subset関数 データフレームやmatrixからの条件指定による行・列の抽出
    R の subset関数は、データフレームやマトリックスか…

おすすめ記事

  1. bash 応用, シェルスクリプト

    シェル スクリプト ファイル存在チェック・空のファイルチェック
    bashでスクリプトを作成するときに、よく使うのがファイル…
  2. R言語, グラフ

    R 使い方 軸・ラベルの調整(向き・サイズ・色など) グラフの描き方
    Rによるplot(グラフ)の描画は、手軽で大変便利です。た…
  3. awk, bash 文字列操作, シェルスクリプト

    bash 部分文字列・置換・長さ・連結・分割の文字列処理
    bashのよく使う文字列処理、部分文字列・置換・連結・長さ…