IGVでbamファイルが表示されない

IGVでbamファイルを開こうとしたのですが、開くことが出来ても全く表示されません。
ネットで調べても原因としてあげられるのは、「bamファイルをsortする」ということがほとんどです。samtoolsでbamファイルをsortしても症状はおなじでした。
原因を見つけるまで時間がかかったので、ここに書いておきます。

reference genomeの位置は表示されます。

試しにsortしていないbamファイルを利用してみましたが、今度は開くことが出来ません。
bam, bam.baiファイルには問題なし
開くことは出来るので、きちんと***.bamという表示はありました。

下のGeneも表示されているので、Reference Genomeの設定も問題なさそうです。

原因は…
bam file
>chr01
fasta fale
>Chr01
gff3 file
>chr01
何が違うかお気づきですか？

実は染色体名が大文字と小文字で異なっており、同一のものと認識されなかったことが原因です。

他のパターンとして
bam file
>Chr01
fasta fale
>Chr01
gff3 file
>chr01

の場合は一番下の「Gene」の表示が出来ません。
表示されない場合は、同一のIDかどうか？を確認してみてください。

OpenRefineで列を結合~ Excelのconcatenateに相当する作業を...

Open Refineで列を結合する場合

Excelではconcatenateに相当します

chr列とposition列を結合し、間に"."を入れてみます。

Excelでは

=concatenate(A1, ".", B1)

と関数を書きます。

OpenRefineの場合は

chrの列を選択し、Edit cellsからTransformを選択、


cells["chr"].value + "." + cells["position"].value

を入力します。

cells["列名"].value + "文字列" + cells["列名"].value

で表示されます。

sedコマンドを使わずに、スラッシュを含む文字列の置換 ~ Open Refineを利用

送られてきたデータにスラッシュが含まれていたので、sedコマンドが使えない...

そこでOpen Refineを利用

変換したい列を選択し、

Edit cells から Transform...を選択
次のような画面が出てくるので、

value.replace関数を使います。

"/"から"."への置換は

value.replace("/",".")
と入力。

今回は/で列を分割したかったので

Edit columnから

Split into several columns...

How to Split Columnと聞かれるので

by separator

でスラッシュ(/)を入力して分割しました。

ほかにもデータの形を整えるにはOpen Refineは便利です。

DDBJ pipelineからbam, bam.baiファイルをダウンロードしてIGVで利用

DDBJ pipelineを利用して解析を行っているのですが、結果をIGVで表示させたいと考えました。

IGVの使い方を見るとsorted bam fileとbam index fileが必要と書かれています。

samtoolsで変換しなければいけないと考えていましたが、これらのファイルをDDBJ pipelineの結果からダウンロードすれば良いことに気づきました。

mappingの結果のページから、下記の二つのファイルをダウンロードすればOK

ダウンロードが上手くいっているかを確認するため、ダウンロードしたファイルに対してMD5というツールを使いました。

$hkane:openssl md5 ファイル名

出てきた結果が、DDBJ pipelineのページの一番右のカラムのMD5の部分をクリックした結果と同じであれば正常です。

MD5の使い方を教えてくださったK様にこの場を借りてお礼申し上げます。

サイズの大きいzip fileを解凍

bam.zip ファイルが解凍できないというトラブルがあったので、解決法をメモしておきます。

当初は解凍先を外付けHDにしていたためだと考えていました。
HDのフォーマットを大容量用に変更しました。

しかしHDの設定を変更しても、エラーが出てきました。

Mac User の私が利用したのは「p7zip」インストールは
 hkane$ brew install p7zip

==> Downloading https://downloads.sourceforge.net/project/p7zip/p7zip/9.20.1/p7z
######################################################################## 100.0%

と表示されます。

brewのインストール方法はいくつか参考になるページがあるので探してみてください。

使うときは

hkane$ 7za x　zipファイル名

大きなbamファイルも無事解凍出来ました。