小编ic2*_*luk的帖子

找到DNA序列中所有重复的4聚体 - Perl

你好,

我尝试编写一个程序,读取包含多个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)

regex perl dna-sequence fasta

4
推荐指数
1
解决办法
367
查看次数

Org-mode作为非程序员的笔记本

我是非程序员(除了一些Perl和R技能).我今年开始攻读博士学位,并考虑将org-mode作为实验室笔记本.但是,在我花时间学习它之前,我有一堆问题..

最重要的是,在工作中,我在不同的PC上工作,而不是在家工作,但我使用它们来工作,所以我需要的是无缝连接.所以我的想法是将组织文件存储在云上,但实际上我将在我的组织模式文件中链接文件,但实际上这些文件在服务器上(在工作中)我不允许存储在我的云上.这意味着,假设我在电子表格中有一些数据,并希望在我的组织文件中链接此电子表格.但是,在我的私人电脑上,此文件不存在.org-mode如何处理?它甚至可能吗?我的意思是,它有抱怨吗?另外,我的工作PC上有一些文件夹对我来说是锁定的,因此安装和工作通常由IT服务完成.emacs甚至不起作用吗?我尝试使用R并且它可以工作,甚至可以从R中安装软件包,但是,Perl并不像从终端执行perl脚本需要锁定的文件夹.

此外,如果你可以推荐我开始的好资源我会很高兴:)

提前致谢!

cloud emacs org-mode

2
推荐指数
1
解决办法
674
查看次数

标签 统计

cloud ×1

dna-sequence ×1

emacs ×1

fasta ×1

org-mode ×1

perl ×1

regex ×1