perl - 用perl查找DNA序列中的核苷酸

我有DNA序列，我想在人们选择的位置找到该序列的核苷酸。下面是示例:

输入序列DNA:
ACTAAAAATACAAAAATTAGCCAGGCGTGGTGGCAC(序列的长度为33)
输入位置:(12)

我希望结果是12位核苷酸是AAA。

我没有发现该职位的氨基酸的问题。下面是我当前的代码。

print "ENTER THE FILENAME OF THE DNA SEQUENCE:= ";
$DNAfilename = <STDIN>;
chomp $DNAfilename;
unless ( open(DNAFILE, $DNAfilename) ) {
  print "Cannot open file \"$DNAfilename\"\n\n";
}
@DNA = <DNAFILE>;
close DNAFILE;
$DNA = join( '', @DNA);
print " \nThe original DNA file is:\n$DNA \n";
$DNA =~ s/\s//g;

print" enter the number ";
$po=<STDIN>;

@pos=$DNA;
if ($po>length($DNA))
{
  print" no data";
}

else
{
  print " @pos\n\n";
}

请提出建议，如何找到DNA序列的位置。

最佳答案

my $nucleotide = substr $DNA, $po, 3;

这将从位置$po到$po+2的3个核苷酸进行分配，并将其分配给$nucleotide。

关于perl - 用perl查找DNA序列中的核苷酸，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/7090371/