プライマー設計する時などに、
塩基配列とアミノ酸配列を上下に並行してアラインメントした資料があると便利なのですが(参考画像、文章末に説明有)、
塩基配列を入力するだけで自動でそのようなファイルを作製できるツール(PCソフト・ウェブツールなど)などありませんでしょうか?
また、
・一定の塩基数間隔で改行
・一部の配列をラベル(アンダーライン・フォント・注釈など)
・ファイルエクスポート
・プリントアウト
などの機能付きだと、なお助かります。
現状、そのような資料を手打ちですべて作製しています。
ご存じの方おられましたら、
ご教示いただくことできないでしょうか?
よろしくお願いいたします。
(添付した参考画像)
目的遺伝子(仮)のDNA配列の下に、
コドンに対応するアミノ酸の一文字表記を
テキストファイルに打ち込んだもの。
・100塩基毎に改行
・イントロン →グレー、エキソン →ブラック
・プライマー →アンダーライン+イタリック
など手を加えています。
No.1ベストアンサー
- 回答日時:
回答がないので、しゃしゃり出てきました。
多分探してもぴったりのものは難しいと思います。
こんな場合は作ってしまうのが早いと思います。頭の体操のつもりでちょっとやってみました。
Rubyという言語です。もしMacをお使いでしたら、bioと言うライブラリを追加すれば動きます。
sudo gem install bio
で追加できます。Windowsの場合はRubyとrubygemsをインストールしたあとにbioを追加してください。
文字の修飾なしですが出力はほぼ同じ形式になります。文字の修飾は手でやるか、rtfというライブラリを使うとわりと簡単にできます。
参考まで
require 'rubygems'
require 'bio'
primer_seq="CCGTTCAGAG"
primer_name="primer-1"
na_seq="AATCTTAGAATAAAAAATGGGTACCGTTCAGAGACCTTTAGAGATTGCAAGGCATCACAGATGATAAAAAGCTCCATCTCTAGACGTGTTCAGGAGTGGGTTGGGGCTTTGACCTTGACTAGCTGCATCAACTTGGACAAGTCACTTCGCTTCCCTGTGCCTCAGTTTCCTCATCCATAT"
abort unless primer_pos=/#{primer_seq}/=~na_seq
exon_pos=primer_pos+primer_seq.length
exon=Bio::Sequence::NA.new(na_seq[exon_pos..-1])
aa_seq=exon.translate
# 出力
aa_str=" "*exon_pos+aa_seq.scan(/./).collect{|aa| "#{aa} "}.join
primer_str=" "*primer_pos+primer_name
primer_str=primer_str+" "*(na_seq.length-primer_str.length)
na_seq.scan(/.{1,100}/).zip(primer_str.scan(/.{1,100}/), aa_str.scan(/.{1,100}/)) do |ns, ps, as|
puts ps
puts ns
puts as
end
ありがとうございます。
bioruby
br_biofetch.rb
br_bioflat.rb
br_biogetseq.rb
br_pmfetch.rb
はインストールできたのですが、
素人過ぎてそこから先に進めませんでした。
いい機会ですので、
Rubyを少し勉強して
上記のソースコード(?)を実行してみたいと思います。
No.2
- 回答日時:
ANo.1のお礼について
どこかでインストール時にエラーが出ているのでしょうか?
私の環境は、rubyは1.8.7です。(Mac OSX 10.6.8, 10.7.2)
ruby -v
で確認できます。
biorubyのdocumentに
REQUIREMENTS
Ruby 1.8.6 or later (except Ruby 1.9.0) – www.ruby-lang.org/
Ruby 1.8.7-p352 or later is recommended.
Not yet fully ready with Ruby 1.9, although many components can now work in Ruby 1.9.1 and Ruby 1.9.2.
とありますので、rubyが古いようでしたら、新しくrubyをインストールされることをお勧めします。
bio-rubyのインストールでエラーが出ているのでしたら
sudo gem install bio
でもだめでしょうか?
(Windowsの場合はsudoが要らないかも)
古いバージョンではありますが、
環境はki073さんと同じだと思うので
使用には問題はないはずです > ruby 1.8.7 (P249) (Mac OSX 10.7.2)
ruby -vで確認できました。
お恥ずかしい話、
これまでbioinformaticsと呼べるものはmegaなどguiしか使っておらず、
cuiはrの参考書を読んで少しだけいじっただけで、
インストールから先の操作を全く知らないのです。
下記のようなweb siteでもbiorubyが勉強できそうなので、
少しずつですが解読して使ってみようと思います。
BioRuby http://bioruby.open-bio.org/
BioRuby Wiki(日本語) http://dev.bioruby.org/ja/
BioRubyの使い方 http://kanehisa.hgc.jp/~k/lecture/20040824JST/bi …
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 高校 生物の問題で下の塩基配列は、DNAを構成する2本のヌクレオチド鎖のうち、一方の塩基配列の一部を示した 2 2022/05/08 17:01
- 生物学 塩基配列を解析する意味は? 1 2023/02/02 00:51
- 生物学 高校生物の質問です。 タンパク質の平均分子量をアミノ酸の平均分子量で割ってタンパク質を構成するアミノ 1 2023/01/07 16:10
- 生物学 イントロンについて 1 2022/07/24 21:27
- Visual Basic(VBA) 列と行の名前(重複あり)が交差するセルに、データを入力したい 1 2022/06/18 21:20
- 生物学 生物の問題で分からないものがあるので教えてください。 1.次のmRNAを基に作られるアミノ酸の配列を 2 2022/10/19 16:23
- 生物学 アミノ酸配列頑張っは分かってるのですがmRNAがまだ分からないのですが、アミノ酸配列からmRNAに変 1 2023/06/28 15:41
- 環境学・エコロジー アミノ酸について 3 2022/03/29 00:59
- Excel(エクセル) Excelにて、フォルダ内のTextファイルをマクロで統合すると文字化けしてしまう時の解消コード 4 2023/01/01 07:32
- Excel(エクセル) エクセルの条件付き書式 個人シートを参照して集計シートに色付けしたい 1 2023/06/22 00:39
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
DNA
-
PCRの 5’・3’ が分からず困っ...
-
遺伝子操作
-
ハイブリダイゼーションとアニール
-
系統樹の読み方について
-
突然変異(転座・逆位・欠失・...
-
2つのシークエンスによる結果の...
-
コドンとトリプレットの違いに...
-
アミノ酸配列からの逆翻訳プロ...
-
フィボナッチ数列と生物
-
真核生物のコザック配列について
-
HBVの「マイナス鎖」と「プラス...
-
ポリヌクレオチド鎖とヌクレオ...
-
DNA塩基配列から推定されるアミ...
-
PCRプライマーへの制限酵素配列...
-
ヒトのGC含有率はいくらですか?
-
塩基配列を解析する意味は?
-
生物の問題です。 「一本鎖切断...
-
DNA
-
DNAとRNAのどのような性質に違...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
cos末端とは何ですか?DNAの端...
-
系統樹の読み方について
-
2つのシークエンスによる結果の...
-
DNA断片のシークエンス結果の読...
-
Poly(IC)とは
-
DNA
-
遺伝子操作
-
Primer3のLeftとRightの方向
-
DNA塩基配列から推定されるアミ...
-
DNAとRNAのどのような性質に違...
-
poly-Aとは
-
アリル特異的PCR法の原理を教え...
-
生化学に関して
-
一本鎖核酸のモル吸光係数の算出
-
卒論で遺伝子配列とアミノ酸配...
-
HBVの「マイナス鎖」と「プラス...
-
cDNA産物は1本鎖ですか?
-
Blastにおける用語
-
Shine-Dalgarno...
-
PCRの 5’・3’ が分からず困っ...
おすすめ情報