Bioops

Bioinformatics=(ACGAAG->AK)+(#!/bin/sh)+(P(A|B)=P(B|A)*P(A)/P(B))

利用Biperl获取序列基本统计信息（Bioperl HOWTO翻译9）

| Comments

获取序列基本统计信息

英文原文

除了前面讲的可通过序列对象的不同方法来获取序列文件中已有的序列信息，也可以利用bioperl获取序列的其他一些信息。例如，SeqStats对象可以获取序列的分子质量，单个氨基酸或核苷酸频率，核苷酸序列所包括的密码子频率，等等。

1
2
3
4
5

use Bio::Tools::SeqStats;
$seq_stats  = Bio::Tools::SeqStats->new($seqobj);
$weight = $seq_stats->get_mol_wt();
$monomer_ref = $seq_stats->count_monomers();
$codon_ref = $seq_stats->count_codons();  # for nucleic acid sequence

注：有些序列包含一些有歧义的code，get_mol_wt()会返回两个数组，一个是最大的分子质量，一个是最小的分子质量。

SeqWords对象可以统计氨基酸或核苷酸片段的频率。详见See Bio::Tools::SeqStats和Bio::Tools::SeqWords。

Comments