如何使用Bioperl计算fasta文件中的序列

时间:2016-05-04 04:59:19

标签: bioinformatics bioperl

晚上好,我有一个bioperl代码来计算fasta文件中的序列数,但是我试图在任何给定的fasta文件中修改代码以计算短于20和更长的序列。代码在

之下
use strict;

use Bio::SeqIO;

my $seqfile = 'sequences.fa';

my $in = Bio::SeqIO->new
        (
                -format => 'fasta',
                -file => $seqfile
        ) ;

my $count = 0;

while (my $seq = $in->next_seq)
{
        $count++;
}

print "There are $count sequences\n";

1 个答案:

答案 0 :(得分:4)

您可以使用sequence对象的length函数在while循环中构造if语句。

command.AddWithValue("@TeamId", value);
command.Add("@TeamId", SqlDbType.Int).Value = value;