求助代碼

zhagnqiang829

一 把基因名（16191_NC_003284）下I001序列和I004序列互換，I002和I003互換，E開頭的序列位置不變。還句話說(shuō)就是把原來(lái)I001，I002，I003，I004序列順序顛倒一下，降序排列。
二 I開頭的序列比較多，其基因名龐大，需要批量處理。序列如下：
16191_NC_003284
I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
E012ATGACAAGGAAGCGGCGGACATTGAA
I002CACTCAAGTCTGCCGTTT
E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA
I003CACCCATTAAAAAAATCCTATATCAAACACATCGTTTTCAGT
E034CGGATCAGTTGTTTCGTCATATTTC
I004CATTCAAAAGCCTTGAAAAAAATTTTAC
E045AGCAAATGAAGGCCGTTTCAACCGAACGAAAACGCGG

16192_NC_003284
I001CACAAGAACTAGGTACTTTG
E012TCTACAACACGTCTGCTCTCTGTG
I002CATGCTAAAAGCTAGTC
E023GCTGCAGAGAAGAAGGTTCTATGCGGACAACT
處理之后結(jié)果為：
16191_NC_003284
I004CATTCAAAAGCCTTGAAAAAAATTTTAC
E012ATGACAAGGAAGCGGCGGACATTGAA
I003CACCCATTAAAAAAATCCTATATCAAACACATCGTTTTCAGT
E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA
I002CACTCAAGTCTGCCGTTT
E034CGGATCAGTTGTTTCGTCATATTTC
I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
E045AGCAAATGAAGGCCGTTTCAACCGAACGAAAACGCGG

16192_NC_003284
I002CATGCTAAAAGCTAGTC
E012TCTACAACACGTCTGCTCTCTGTG
I001CACAAGAACTAGGTACTTTG
E023GCTGCAGAGAAGAAGGTTCTATGCGGACAACT

mcshell

回復(fù) 1# zhagnqiang829

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. use Data::Dumper;
   
5. my %hash;
   
6. while(<DATA>){
   
7.     /(\d+_NC_\d+)/?$hash{$1}=():push @{$hash{$1}},$_;
   
8. }
   
9. map{print "$_\n";
   
10.     my (@arr_i,@arr_e);
    
11.     for (@{$hash{$_}}){
    
12.     /^I\d+.*/?push  @arr_i,$_:push @arr_e,$_;
    
13.     }
    
14.     for (sort{$b cmp $a}@arr_i){
    
15.         print "$_";
    
16.         print shift @arr_e;
    
17.         }
    
18.     print "\n";
    
19.    
    
20.     }sort keys %hash;
    
21. __DATA__
    
22. 16191_NC_003284
    
23. I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
    
24. E012ATGACAAGGAAGCGGCGGACATTGAA
    
25. I002CACTCAAGTCTGCCGTTT
    
26. E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA
    
27. I003CACCCATTAAAAAAATCCTATATCAAACACATCGTTTTCAGT
    
28. E034CGGATCAGTTGTTTCGTCATATTTC
    
29. I004CATTCAAAAGCCTTGAAAAAAATTTTAC
    
30. E045AGCAAATGAAGGCCGTTTCAACCGAACGAAAACGCGG
    
31. 
    
32. 16192_NC_003284
    
33. I001CACAAGAACTAGGTACTTTG
    
34. E012TCTACAACACGTCTGCTCTCTGTG
    
35. I002CATGCTAAAAGCTAGTC
    
36. E023GCTGCAGAGAAGAAGGTTCTATGCGGACAACT

復(fù)制代碼

只是一個(gè)紅薯

關(guān)鍵代碼給點(diǎn)注釋就好了回復(fù) 2# mcshell


   

rubyish

如下：

1. #!/usr/bin/perl
   
2. local $/ = "\n\n";
   
3. while (<DATA>) {
   
4.     my ( @i, @e );
   
5.     my @a = split;
   
6.     my $t = shift @a;
   
7.     /^I/ ? unshift( @i, $_ ) : push( @e, $_ ) for @a;
   
8.     print "$t\n", map { "$i[$_]\n$e[$_]\n" } 0 .. $#i;
   
9.     print "\n";
   
10. }
    
11. __DATA__
    
12. 16191_NC_003284
    
13. I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
    
14. E012ATGACAAGGAAGCGGCGGACATTGAA
    
15. I002CACTCAAGTCTGCCGTTT
    
16. E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA
    
17. I003CACCCATTAAAAAAATCCTATATCAAACACATCGTTTTCAGT
    
18. E034CGGATCAGTTGTTTCGTCATATTTC
    
19. I004CATTCAAAAGCCTTGAAAAAAATTTTAC
    
20. E045AGCAAATGAAGGCCGTTTCAACCGAACGAAAACGCGG
    
21. 
    
22. 16192_NC_003284
    
23. I001CACAAGAACTAGGTACTTTG
    
24. E012TCTACAACACGTCTGCTCTCTGTG
    
25. I002CATGCTAAAAGCTAGTC
    
26. E023GCTGCAGAGAAGAAGGTTCTATGCGGACAACT
    

復(fù)制代碼

mcshell

回復(fù) 4# rubyish


    還是原來(lái)的頭像萌啊

yinyuemi

1. #!/usr/bin/env perl
   
2. use strict;
   
3. my @data;
   
4. while(<DATA>){
   
5.     if (/^\d/){print;undef(@data)};
   
6.     push @data,$_ if /^[IE]/;
   
7.     f() if /^$/ ;
   
8. }
   
9. f();
   
10. sub f(){
    
11.      map{print $data[$#data-2*$_-1],$data[2*$_+1]}0..$#data/2;
    
12.      print "\n";
    
13. }
    
14. __DATA__
    
15. 16191_NC_003284
    
16. I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
    
17. E012ATGACAAGGAAGCGGCGGACATTGAA
    
18. I002CACTCAAGTCTGCCGTTT
    
19. E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA
    
20. I003CACCCATTAAAAAAATCCTATATCAAACACATCGTTTTCAGT
    
21. E034CGGATCAGTTGTTTCGTCATATTTC
    
22. I004CATTCAAAAGCCTTGAAAAAAATTTTAC
    
23. E045AGCAAATGAAGGCCGTTTCAACCGAACGAAAACGCGG
    
24. 
    
25. 16192_NC_003284
    
26. I001CACAAGAACTAGGTACTTTG
    
27. E012TCTACAACACGTCTGCTCTCTGTG
    
28. I002CATGCTAAAAGCTAGTC
    
29. E023GCTGCAGAGAAGAAGGTTCTATGCGGACAACT

復(fù)制代碼

zhagnqiang829

ok ,Thank you very much.

zhagnqiang829

1 如果我要實(shí)現(xiàn)下面情況改怎么處理：數(shù)據(jù)如下
16191_NC_003284
I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
E012ATGACAAGGAAGCGGCGGACATTGAA
I002CACTCAAGTCTGCCGTTT
E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA

16192_NC_003284
I001CACAAGAACTAGGTACTTTG
E012TCTACAACACGTCTGCTCTCTGTG
2 結(jié)果為：
16191_NC_003284
I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
E012ATGACAAGGAAGCGGCGGACATTGAA
I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA
I002CACTCAAGTCTGCCGTTT
E012ATGACAAGGAAGCGGCGGACATTGAA
I002CACTCAAGTCTGCCGTTT
E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA

16192_NC_003284
I001CACAAGAACTAGGTACTTTG
E012TCTACAACACGTCTGCTCTCTGTG

3 說(shuō)明：在I001后分別把E0**序列放后面，同理I002，I003.......I00*.
回復(fù) 4# rubyish


   

rubyish

在I001后分別把E0**序列放后面，同理I002，I003.......I00*.
如下:

1. #!/usr/bin/perl
   
2. $/ = "\n\n";
   
3. while (<DATA>) {
   
4.     my @i = ( [], [] );
   
5.     my @a  = /^(\w+\n)/mg;
   
6.     my $t  = shift @a;
   
7.     push $i[/^E/], $_ for @a;
   
8.     print $t, map { $a = $_; map "$a$_", @{$i[1]} } @{$i[0]};
   
9.     print "\n";
   
10. }
    
11. __DATA__
    
12. 16191_NC_003284
    
13. I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
    
14. E012ATGACAAGGAAGCGGCGGACATTGAA
    
15. I002CACTCAAGTCTGCCGTTT
    
16. E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA
    
17. 
    
18. 16192_NC_003284
    
19. I001CACAAGAACTAGGTACTTTG
    
20. E012TCTACAACACGTCTGCTCTCTGTG
    

復(fù)制代碼