トラブルシューティングと編集のヘルプを探しています。これは宿題です。私の教授はフォーラムの使用を奨励しています。私はまだ Perl 関数や Subs の経験がないので、理解できるように回答を適切なレベルに制限してください。
スクリプトの目的は、DNA の文字列 (または後で追加するコマンド ラインからのファイル) を読み取り、それを RNA に変換してから、タンパク質の値を大文字の 1 文字のアミノ酸名の形式で返すことです。
スクリプトの機能:
最初の文字から 3 文字の「コドン」を取得し、それらに 1 文字のシンボル (ハッシュ テーブルからの大文字の 1 文字のアミノ酸名) を付けます。
AUG (「M」) で始まり、UAG、UAA、または UGA で終わる文字列である RNA タンパク質を出力します。
ギャップが発生した場合、新しい行が開始され、プロセスが繰り返されます。ギャップは 3 の倍数であると想定できます。
私が知る限りの主な問題:
ハッシュ テーブルを介してデータをループさせる場所がわかりません。Foreach ブロックの前後に配置してみました。また、Foreach ブロックを完全に取り除き、While & If を試しました。
Foreach ブロックは @all_codons 配列のすべてを処理しておらず、AUG でのみ停止しているようです。
明らかで最大の問題は、何も返さないことです。どこかで $next_codon 値に「false」が割り当てられています。各行を少しずつコメントアウトしようとしました-何かを返した最後の行は My $start で、そこからはすべてfalseです。
スクリプト:
$^W = 1;
use strict;
my $dna_string = "CCCCAAATGCTGGGATTACAGGCGTGAGCCACCACGCCCGGCCACTTGGCATGAATTTAATTCCCGCCATAAACCTGTGAGATAGGTAATTCTGTTATATCCACTTTACAAATGAAGAGACTGAGGCAAAGAAAGATGATGTAACTTACGCAAAGC";
my %codon_codes = (
"UUU" => "f", "UUC" => "f", "UUA" => "l", "UUG" => "l",
"CUU" => "l", "CUC" => "l", "CUA" => "l", "CUG" => "l",
"AUU" => "i", "AUC" => "i", "AUA" => "i", "AUG" => "m",
"GUU" => "v", "GUC" => "v", "GUA" => "v", "GUG" => "v",
"UCU" => "s", "UCC" => "s", "UCA" => "s", "UCG" => "s",
"CCU" => "p", "CCC" => "p", "CCA" => "p", "CCG" => "p",
"ACU" => "t", "ACC" => "t", "ACA" => "t", "ACG" => "t",
"GCU" => "a", "GCC" => "a", "GCA" => "a", "GCG" => "a",
"UAU" => "y", "UAC" => "y", "UAA" => " ", "UAG" => " ",
"CAU" => "h", "CAC" => "h", "CAA" => "q", "CAG" => "q",
"AAU" => "n", "AAC" => "n", "AAA" => "k", "AAG" => "k"
);
my $rna_string = $dna_string;
$rna_string =~ tr/[tT]/U/;
my @all_codons = ($rna_string =~ m/.../g);
foreach my $next_codon(@all_codons){
while ($next_codon =~ /AUG/gi){
my $start = pos ($next_codon) -3;
last unless $next_codon =~ /U(AA|GA|AG)/gi;
my $stop = pos($next_codon);
my $genelen = $stop - $start;
my $gene = substr ($next_codon, $start, $genelen);
print "\n" . join($start+1, $stop, $gene,) . "\n";
}
}