perl 重命名ID的代码

一串儿葡萄皮

我有n条fasta格式的序列，现在想对序列的ID名称改变为从1开始 的数字，即seq1 , seq2.....seqn,希望大神帮忙用perl写个代码

我文件格式，ID如下如下：
>AT3G56040.1|PACid:19662817
MANPQASPILHHPQNHLSLFHFRTTTSPRSFSSLHFRKPLLFLSSSSSFSSKLQQSEQQCNNHQVRHVSTVPVEYSTPTPPESDDFLSEIDRLKSLLSKLDVSKDLRRKDAVIDADSRVRRFFSENRGGLSKVFGYLGLNSNEMFLVKCVIAAGQEHALCMNYEEAFGEEEEEYTVRSSVKNALYALVEMIERFDVNSSGYKGRREMGTVLDSEEIAHFRKFLTFLEEIEQFYDCIGGIIGYQVMVLELLHQSSKRRNTNRSQLVEESLGCQYLEMHTPSVLDLTQEEDYASQAALWGIEGLPDLGEIYPLGGAADRLGLIDSETGECLPAAMLAHCGRTLLEGLIRDLQAREFLYFKLYGKQCVTPVAIMTSAAKNNHEHVSSLCERLKWFGRGQSNFRLFEQPLVPAVSAEDGQWIVSKPFVPVSKPGGHGVIWKLAYDKGVFNWFYDHGRKGATVRQVSNVVAATDVTLLALAGIGLRYNKKLGFASCKRNAGATEGINVLMEKKNFDGKWEYGISCIEYTEFDKFDISNRSPSSNGLQADFPANTNILYVDLHSAELIGSSSNAKSLPNMVLNTKKRIEYLDQYGDYHSVMGGRLECTMQNIADNFFNKFPSRCHGSLEDKLDTYIVYNERRKVTSSAKKKKPHASAALHQTPDGALLDILRNGYDLLTECDIKLPMIEANDKYVDSPPPYLILLHPALGPLWEVSRQKFKGGSISSCSELQLEIAEFSWNNVQVDGSLIVTAENAMGSTTPNDNGEPILQYGLRCGKCKLHNVNVVNRGIDWNSKSNVYWRNDVNRLETCKIILHGNAEFEASNVTIEGHHVFEVPDGHKLKITSGNAGLSINLEALKEEVMETGSWYWNYQLNGSHIHLQQVEVSQS*
>AT3G03250.1|PACid:19663528
MAATTENLPQLKSAVDGLTEMSESEKSGFISLVSRYLSGEAQHIEWSKIQTPTDEIVVPYEKMTPVSQDVAETKNLLDKLVVLKLNGGLGTTMGCTGPKSVIEVRDGLTFLDLIVIQIENLNNKYGCKVPLVLMNSFNTHDDTHKIVEKYTNSNVDIHTFNQSKYPRVVADEFVPWPSKGKTDKEGWYPPGHGDVFPALMNSGKLDTFLSQGKEYVFVANSDNLGAIVDLTILKHLIQNKNEYCMEVTPKTLADVKGGTLISYEGKVQLLEIAQVPDEHVNEFKSIEKFKIFNTNNLWVNLKAIKKLVEADALKMEIIPNPKEVDGVKVLQLETAAGAAIRFFDNAIGVNVPRSRFLPVKASSDLLLVQSDLYTLVDGFVTRNKA

华小飞_Perl

来一个~

1. #!/usr/bin/perl
   
2. 
   
3. use warnings;
   
4. use strict;
   
5. 
   
6. my $i;
   
7. while (<>) {
   
8.         chomp;
   
9.         if (/^>/m) {
   
10.                 $i++;
    
11.                 print ">seq", $i, "\n"；
    
12.         }
    
13.         print $_, "\n";
    
14. }

复制代码

sunzhiguolu

1. perl -nle 'print /\A>/ ? ++$n : $_' f

复制代码

b114213903

回复 13# MMMIX


    嗯，您说的很对！这个模块主要就是生物信息员需要用的功能，序列格式转换、网络数据库获取数据、数据处理和分析，报告出图一系列的模块包，很给力！但需要时间学习，甚至安装都有问题，呵呵

MMMIX

回复 12# b114213903


    关键是这个东西它有（专业）针对性，其他一些解决方法怎么看都像是 hackish

b114213903

回复 11# MMMIX


    BioPerl模块，功能还是很强大的！

MMMIX

回复 10# b114213903


    我就知道有对应的模块……

b114213903

1. #!/usr/bin/perl
   
2. use strict;
   
3. use Bio::SeqIO;
   
4. 
   
5. my $fasta=shift @ARGV;
   
6. (my $Out=$fasta)=~s/(\.[^\.]+)$/_out$1/;
   
7. 
   
8. my $IN=Bio::SeqIO->new(-file=>"$fasta",-format=>'fasta');
   
9. my $OUT=Bio::SeqIO->new(-file=>">$Out",-format=>'fasta');
   
10. my $n=0;
    
11. while(my $Seq=$IN->next_seq){
    
12.         $n++;
    
13.         print "Now $n:\t",$Seq->id,"\n";
    
14.         $Seq->id($n);
    
15.         $OUT->write_seq($Seq);
    
16. }
    
17. $IN->close();
    
18. $OUT->close();

复制代码

MMMIX

回复 8# 一串儿葡萄皮


    挑你方便的，能表达清楚意思就行了。

一串儿葡萄皮

回复 5# MMMIX


   thank you. should I say this in English? The problem is resuled.Thank you very much.