你好,
我尝试编写一个程序,读取包含多个DNA序列的FASTA格式文件,识别序列中所有重复的4聚体(即,多次出现的所有4聚体),并打印出重复的4聚体以及查找它的序列的标题.k聚体仅仅是k个核苷酸的序列(例如,"aaca","gacg"和"tttt"是4聚体).
这是我的代码:
use strict;
use warnings;
my $count = -1;
my $file = "sequences.fa";
my $seq = '';
my @header = ();
my @sequences = ();
my $line = '';
open (READ, $file) || die "Cannot open $file: $!.\n";
while ($line = <READ>){
chomp $line;
if ($line =~ /^>/){
push @header, $line;
$count++;
unless ($seq eq ''){
push @sequences, $seq;
$seq = '';
}
} else {
$seq .= $line;
}
} push @sequences, $line;
for (my $i = …Run Code Online (Sandbox Code Playgroud) 我是非程序员(除了一些Perl和R技能).我今年开始攻读博士学位,并考虑将org-mode作为实验室笔记本.但是,在我花时间学习它之前,我有一堆问题..
最重要的是,在工作中,我在不同的PC上工作,而不是在家工作,但我使用它们来工作,所以我需要的是无缝连接.所以我的想法是将组织文件存储在云上,但实际上我将在我的组织模式文件中链接文件,但实际上这些文件在服务器上(在工作中)我不允许存储在我的云上.这意味着,假设我在电子表格中有一些数据,并希望在我的组织文件中链接此电子表格.但是,在我的私人电脑上,此文件不存在.org-mode如何处理?它甚至可能吗?我的意思是,它有抱怨吗?另外,我的工作PC上有一些文件夹对我来说是锁定的,因此安装和工作通常由IT服务完成.emacs甚至不起作用吗?我尝试使用R并且它可以工作,甚至可以从R中安装软件包,但是,Perl并不像从终端执行perl脚本需要锁定的文件夹.
此外,如果你可以推荐我开始的好资源我会很高兴:)
提前致谢!