本文へスキップ

Bedtools

bedtoolsはbam形式のファイルをbed形式のファイルに変更したり,得られたbedを用いて遺伝子範囲にどのくらいのリードがあるのかカウントできたりと,様々な役割がある。

インストール方法

ここでダウンロードする。
詳細なマニュアル

使い方

1. bamからbedに変換

 bamToBed -i filename.bam > filename.bed

2.bedからbamに変換

 bedToBam -i <BED/GFF/VCF> -g <GENOME> > filename.bam

3. bedからfastaに変換

 fastaFromBed [OPTIONS] -fi <input FASTA> -bed <BED/GFF/VCF> -fo <output FASTA>
[options]
-name: bed fileの三列目が,fastaの名前に変更できる。

4. 複数のbedをまとめる

 mergeBed -i <BED/GFF/VCF>

5.レファレンスのbed範囲内に,どのくらいのリードが張り付いたかカウントする

 intersectBed -a reference.bed -b sample.bed -c > sample.count

6. もっとも近い領域を探す

 closestBed [OPTIONS] -a <BED/GFF/VCF> -b <BED/GFF/VCF>
[options]
-d :距離が明記される




inserted by FC2 system