MacOSXのTerminalでEMBOSSを使わせていただいています。 そこで質問なのですが、EMBOSSのパッケージの中に、 「cDNA配列を入力として与えるとその転写産物を出力するコマンド」 というのはありますでしょうか。どうぞよろしくお願い致します。 質問日 Dec 10 '10 at 19:27 redgrapefruit |
EMBOSS では biosed を使うみたいですね biosed hoge.fa -target T -replace U 回答日 Dec 10 '10 at 19:46 dritoshi ♦ mn3 ♦♦ perlのワンライナーで単純変換を試みる場合、適切な処理をいれないとヘッダ行とかの塩基配列ではない部分まで全部変換してしまうので、biosedつかうのがよいです。
(Dec 10 '10 at 19:59)
mya_ ♦
なるほど, biosedだとFastaファイルなどの配列部分だけが置換されるのですね。勉強になりました。
(Dec 10 '10 at 20:08)
redgrapefruit
|
早速のご回答どうもありがとうございます。 perlのワンライナーでの置換、または、EMBOSSでならbiosedというコマンドを使えばよいのですね。 勉強になりました。どうもありがとうございました。 回答日 Dec 10 '10 at 19:53 redgrapefruit |
その通りです。DNA は G A T C から構成されていますが、RNA は G A U C なので、 そのように変換されるEMBOSSのコマンドがあれば教えていただきたいです。 回答日 Dec 10 '10 at 19:37 redgrapefruit |
EMBOSS はわからないですがPerlのワンライナーで簡単に変換することができます。 echo "ATGC" |perl -lpe "s/T/U/g" 回答日 Dec 10 '10 at 19:42 dritoshi ♦ mn3 ♦♦ |