1.8.3 released (Glang-users) - G-language

g.

1.8.3リリースしました。

GenomeProjectorまわりの細かいfixが多いですが、３つ
新しい機能が追加されています。

1. shuffleseq()にk-merを維持するオプションを付けました。

　-kオプションで数字を与えると、k-merの個数（つまり 
1<=x<=k全てにおいて）
　を完全に維持し、なおかつ一様なランダムさでランダム配列を作 
成してくれます。

　文字の種類（塩基なら例えば4, アミノ酸なら20）を 
前提とせず、kも任意で
　完全に個数を維持した配列を作成するのはアルゴリズム的には非常に
　複雑です。詳しくはhelp shuffleseqを見て下さい。

2. $gb->relocate_origin($pos)
      配列の1文字目を$posに移動（よって対応するア 
ノテーションも全て移動）
　したG構造体を返します。これにより、例えば大腸菌のようにori 
が1文字目
　ではないゲノムを一瞬でoriからはじまるようにできます。

3. $gb->reverse_strand()
　DNAの逆鎖に対応したG構造体を返します。つまり、配列は 
complement,
    遺伝子の並びは逆順、direct/complementが逆転したもの 
です。

2と3はオブジェクト指向らしく使うと（Perlだと残念 
ながら美しくないですが）

$gb->reverse_strand()->relocate_origin($ori)->output("out.gbk");

のように、大腸菌ゲノムをoriからスタートさせ逆転した配列の
GenBankファイルが作れます。

＃実際は、$oriを求める時に
＃($ori, $ter) = find_ori_ter($gb->reverse_strand())
＃をする必要があるが。

例えば遺伝子のリストを見るだけならシェルで以下のようにできます。
$gb->reverse_strand()->relocate_origin($ori)->find()

changelogは以下の通り。

===== v.1.8.3    2008.04.21  =====
   *removed deprecated G::Seq::Usage, which was previously integrated  
to G::Seq::Codon
   *added error handling for G::Seq::GCskew::gcsi for when genome size  
is too small
   *$gb->find() is now case insensitive
   *bug fix for arrow direction for G::Seq::GenomeMap::circular_map()
   *added opt_list which shows the default options for Odyssay  
functions. This is in alpha state.
   *major update to G::Seq::Primitive::shuffleseq() to support  
preservation of k-mer count.
   *added $gb->reverse_strand(), $gb->relocate_origin()
   *genomicskew() now returns array of references to result arrays
   *added -at, -purine and -keto options for gcsi()
   *updated bundled genomes

Kazuharu Arakawa, Ph.D.
Institute for Advanced Biosciences, Keio University
252-8520 Japan   Tel/Fax: +81-466-47-5099

-------------- next part --------------
HTMLの添付ファイルを保管しました...
ダウンロード 

G-language

[Glang-users] 1.8.3 released