Perlその2
大きなmultifasta形式のファイルを100コずつのmultifastaのファイルに分割する呪文。
#!/usr/bin/perl my $i = 0; my $file = "sequence"; while() { chomp; my $string = $_; if($string =~ /^\>/) { $i++; if($i % 100 == 1) { close FILE; open(FILE, ">$file.$i") or die "$file.$i $!\n"; } } print FILE "$string\n"; } close FILE;
がしかし、まだまだ自分で呪文を唱えるには至らず。。
先月のPerl記念日の復習もここのところずっと放置していたので呪文を唱えられないのは当然で、今回の分も含め本腰入れて復習・勉強すべし。