求助代码

zhagnqiang829

一 把基因名（16191_NC_003284）下I001序列和I004序列互换，I002和I003互换，E开头的序列位置不变。还句话说就是把原来I001，I002，I003，I004序列顺序颠倒一下，降序排列。
二 I开头的序列比较多，其基因名庞大，需要批量处理。序列如下：
16191_NC_003284
I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
E012ATGACAAGGAAGCGGCGGACATTGAA
I002CACTCAAGTCTGCCGTTT
E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA
I003CACCCATTAAAAAAATCCTATATCAAACACATCGTTTTCAGT
E034CGGATCAGTTGTTTCGTCATATTTC
I004CATTCAAAAGCCTTGAAAAAAATTTTAC
E045AGCAAATGAAGGCCGTTTCAACCGAACGAAAACGCGG

16192_NC_003284
I001CACAAGAACTAGGTACTTTG
E012TCTACAACACGTCTGCTCTCTGTG
I002CATGCTAAAAGCTAGTC
E023GCTGCAGAGAAGAAGGTTCTATGCGGACAACT
处理之后结果为：
16191_NC_003284
I004CATTCAAAAGCCTTGAAAAAAATTTTAC
E012ATGACAAGGAAGCGGCGGACATTGAA
I003CACCCATTAAAAAAATCCTATATCAAACACATCGTTTTCAGT
E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA
I002CACTCAAGTCTGCCGTTT
E034CGGATCAGTTGTTTCGTCATATTTC
I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
E045AGCAAATGAAGGCCGTTTCAACCGAACGAAAACGCGG

16192_NC_003284
I002CATGCTAAAAGCTAGTC
E012TCTACAACACGTCTGCTCTCTGTG
I001CACAAGAACTAGGTACTTTG
E023GCTGCAGAGAAGAAGGTTCTATGCGGACAACT

mcshell

回复 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

复制代码

只是一个红薯

关键代码给点注释就好了回复 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
    

复制代码

mcshell

回复 4# rubyish


    还是原来的头像萌啊

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

复制代码

zhagnqiang829

ok ,Thank you very much.

zhagnqiang829

1 如果我要实现下面情况改怎么处理：数据如下
16191_NC_003284
I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
E012ATGACAAGGAAGCGGCGGACATTGAA
I002CACTCAAGTCTGCCGTTT
E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA

16192_NC_003284
I001CACAAGAACTAGGTACTTTG
E012TCTACAACACGTCTGCTCTCTGTG
2 结果为：
16191_NC_003284
I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
E012ATGACAAGGAAGCGGCGGACATTGAA
I001CACTCAATAAATGAAGGAACTTTTGAGAGCTTAA
E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA
I002CACTCAAGTCTGCCGTTT
E012ATGACAAGGAAGCGGCGGACATTGAA
I002CACTCAAGTCTGCCGTTT
E023ATGGAAGAGATGAGAAAAAGGGCGAGGCTTATGA

16192_NC_003284
I001CACAAGAACTAGGTACTTTG
E012TCTACAACACGTCTGCTCTCTGTG

3 说明：在I001后分别把E0**序列放后面，同理I002，I003.......I00*.
回复 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
    

复制代码