忍者ブログ
[PR]
×

[PR]上記の広告は3ヶ月以上新規記事投稿のないブログに表示されています。新しい記事を書く事で広告が消えます。


2017/11/22 00:41 |
FASTAを80行にするプログラムを書いてみた。
先ほど書いた
リボソーマルプロテインは一行に何塩基もあってみずらいので、
ちゃんとしたFASTAファイルにする必要がある。
そこで、Perlでプログラムを書いて整える。

手作業で79個の遺伝子を80行に揃えるのは大変。
一つの遺伝子が5行だとしても400行リターンを押さないといけない。
パールなら簡単。

以下がそのプログラムファイル。
使い方はまた今度説明するかも。

以下そのプログラム

############################################################################
#!/usr/bin/perl -w
use strict;
use warnings;
use POSIX;

# -------------------------------------------------------------- #
# MJsamurai
# @McGill, 2007
# Oyaji
# -------------------------------------------------------------- #


my $counter = 0;


########################### enter the fastafile to hash.

my ($fastafile) = @ARGV;
open FASTA, "<$fastafile";

###########################

while (<FASTA>)
{
 if ($_ =~ /^>/)
 {
        my $header = $_;
        my $line = <FASTA>;

    $counter = $counter + 1;
       
        print "$header";
        print_sequence($line, 80);
 }
}
close FASTA;

print "$counter fasta sequences are here\n";


################# sub print_FASTA sequence ##################

sub print_sequence
{
      my($sequence, $length) = @_;


 # Print sequence in lines of $length
for ( my $pos = 0 ; $pos < length($sequence) ; $pos += $length )
   {
     print substr($sequence, $pos, $length), "\n";
   }
}

################################################################

PR

2007/08/07 00:56 | Comments(0) | TrackBack(0) | パールプログラム(Perl)

トラックバック

トラックバックURL:

コメント

コメントを投稿する






Vodafone絵文字 i-mode絵文字 Ezweb絵文字 (絵文字)



<<79個のリボソーマルプロテイン(80行) | HOME | ユウレイボヤの79個のリボソーマルプロテイン>>
サンカラショッピングモール / ダイエット / ビューティー・美容 / 布団・寝具 / 日用品・雑貨 / 家電製品 / 健康・ヘルスケア / 食べ物 / 飲み物 / アクセサリー / 時計 / ファッション / ブランド / ペット用品 / コスプレ / もしもの人気商品100選 / 新着さんから100 / 送料無料ショップ / もちもちロール / 送料無料ショップ / emsダイエット / ホットシャワー3 /
忍者ブログ[PR]