求助：关于核苷酸和氨基酸筛选的问题（已解决，感谢！）

56836430

在这里真诚感谢各位的帮助，解决了我一个大问题！

sunzhiguolu

@hztj2005, @RE_HASH
佩服, 厉害啊... 向大神致敬!

hztj2005

回复 12# sunzhiguolu

我多年用c++写代码，思路容易理顺。但对Perl并不熟悉，这段代码应该是很笨拙的，至少我自己是这样看。


   

sunzhiguolu

回复 13# hztj2005
大神谦虚了, 我感觉解决问题最重要的就是解决问题的思路.
至于使用什么语言, 都是工具, 次要的.

   

rubyish

回复 10# RE_HASH
nide code, very youyisi. ( nage gee wo kan le 1 tian ~ )
3Q ~ youyisi ~
perltidy ~

1. #!/usr/bin/perl
   
2. my ( $T, $R, $C );
   
3. open C, "c.fasta";
   
4. map { $T .= ( $_ > 5 ) ? "Y" : "N"; } <C>;
   
5. close C;
   
6. 
   
7. $T =~ s/(N+)/{ sprintf "(?:.{%d})", length $1 }/ge;
   
8. $T =~ s/(Y+)/{ $C++; $R .= "\$C"; sprintf "(.{%d})", length $1 }/ge;
   
9. 
   
10. sub CutFile {
    
11.     my ( $F, $T, $R ) = @_;
    
12.     open A, $F;
    
13.     local $_ = do { $/ = undef; <A> };
    
14.     close A;
    
15. 
    
16.     s/(?<=\S{40})\n//g;
    
17.     s/$T/qq(qq($R))/gee;
    
18.     s/(\S{50})/$1\n/g;
    
19.     return $_;
    
20. }
    
21. 
    
22. printf "d.fasta\n%s", CutFile( "a.fasta", $T, $R );
    
23. print "\n\n";
    
24. $T =~ s/(\d+)/{$1 * 3}/eg;
    
25. printf "e.fasta\n%s", CutFile( "b.fasta", $T, $R );
    

复制代码