はいどーも皆さんこんにちは、adachitです。 以前、bedtoolsを使って、2つのファイルの共通領域を書き出そうとしたところ、染色体名が異なるためエラーが出力されました。 ファイルによって、接頭に「chr」が付いていたり付いていなかったりとバラバラです。 こういう場合、Linuxのsedコマンドを使用することで染色体の名前を変更することができます。 それでは、サンプルファイルを使用して実際にやってみます。 今回は、2つのサンプルファイル、1_test_sorted.bedと2_test_sorted.bedを用意しました。 それぞれのファイルの中身は以下のようになっています。 $ les…