perl 引號(hào)

一串兒葡萄皮

部分文件
>Srcc10g063050.1.1 Receptor-like kinase    Try_threonine protein kinase
MKTKHKLTPRPFTPTPSPTHIMGILLICFIFSITNSFVIAQDDAVPAKSFPVFTPEDNFL
IDCGATSSITLPGNKAFQPDQNTAKYLSYTGKDIQACASDKINVPSTLYVNAKIFTTEAI
YTFHASTSGLHWIRLHFFPFKYEEYDLK
">Solyc11g072930.1.1 LRR receptor-like serine_threonine-protein kinase, RLP"
MASSYFLLLVLVLSVFSVSADVFVSLDCGSSEAYTDHETSIDWLGDVDYVANGESHVVPS
NNSISHVSALEVRGLDSTMYSHVDDNSSDPPRITALYLSKFNLSGSLPDFSSMDAL
ETIDLSNNNLDGPIPDFFGTLPNLKELNLANNKFSGPVPASLSNKNGLTLDTSGNSDLCS
SSEESCQNNDSSSPGNDQPTTGSTNNNKKKKKKKNNLPIILGTTISAFLLLWAIVGIFAI
LHYKNKRAATSLINPGQASGGSTPFVDRVQMSEKIEKNPEVTAHDHENSTNV
">Solyc12g00890.1.1 Kinesin  IPR001752  Kinesin, motor region"
ESVINGRDLLGFSLTSPDLVICTGSPDIPARNYGDSPEFLKGCSISLENGI
KGSEEVQAATKLFTDWQGSKDDDLCAPADFELPSPPVEENSSELSVPIVSINVGSTDCIS
SESGIQFSEDKYFCGGNVLSTDTRIEESICASVYQTARVGNFSYHFNNLSAGFYLVDLHF
VEVVLTDGSTGDFSENSPQRNSLEVNGDIKAAGKLQLANVSREK

MMMIX

一串兒葡萄皮 發(fā)表于 2015-08-12 11:09
while(<W>){
   chomp;
    my$i='"\w+"';
        $i=~s/\"//g;

}

你這就是個(gè)空循環(huán), 除了把 W 讀了一遍之外, 啥用都沒有.

另外, 你把 my $i 寫成 my$i 自己看著不糟心么? 多輸入個(gè)空格是能累著你咋滴?

BTW again, 這都不是代碼風(fēng)格問題, 你這完全是審美觀病變壞死的征兆!!!

MMMIX

一串兒葡萄皮 發(fā)表于 2015-08-12 11:09
#!/usr/bin/perl    -w
use warnings;

這兩句完全就是不求甚解的范例: -w 和 use warnings; 的功能是重復(fù)的.

b114213903

回復(fù) 10# 一串兒葡萄皮


    我已經(jīng)在二樓給你提醒過了，ID行有引號(hào)的，不是標(biāo)準(zhǔn)的FASTA格式，BioPerl模塊不能識(shí)別！

1. >Srcc10g063050.1.1 Receptor-like kinase    Try_threonine protein kinase
   
2. MKTKHKLTPRPFTPTPSPTHIMGILLICFIFSITNSFVIAQDDAVPAKSFPVFTPEDNFL
   
3. IDCGATSSITLPGNKAFQPDQNTAKYLSYTGKDIQACASDKINVPSTLYVNAKIFTTEAI
   
4. YTFHASTSGLHWIRLHFFPFKYEEYDLK
   
5. ">Solyc11g072930.1.1 LRR receptor-like serine_threonine-protein kinase, RLP"
   
6. MASSYFLLLVLVLSVFSVSADVFVSLDCGSSEAYTDHETSIDWLGDVDYVANGESHVVPS
   
7. NNSISHVSALEVRGLDSTMYSHVDDNSSDPPRITALYLSKFNLSGSLPDFSSMDAL
   
8. ETIDLSNNNLDGPIPDFFGTLPNLKELNLANNKFSGPVPASLSNKNGLTLDTSGNSDLCS
   
9. SSEESCQNNDSSSPGNDQPTTGSTNNNKKKKKKKNNLPIILGTTISAFLLLWAIVGIFAI
   
10. LHYKNKRAATSLINPGQASGGSTPFVDRVQMSEKIEKNPEVTAHDHENSTNV
    
11. ">Solyc12g00890.1.1 Kinesin  IPR001752  Kinesin, motor region"
    
12. ESVINGRDLLGFSLTSPDLVICTGSPDIPARNYGDSPEFLKGCSISLENGI
    
13. KGSEEVQAATKLFTDWQGSKDDDLCAPADFELPSPPVEENSSELSVPIVSINVGSTDCIS
    
14. SESGIQFSEDKYFCGGNVLSTDTRIEESICASVYQTARVGNFSYHFNNLSAGFYLVDLHF
    
15. VEVVLTDGSTGDFSENSPQRNSLEVNGDIKAAGKLQLANVSREK

復(fù)制代碼

這些示例數(shù)據(jù)，就用BioPerl的Bio::SeqIO模塊來讀只能識(shí)別第一個(gè)序列！

1. #!/usr/bin/perl
   
2. use strict;
   
3. use Bio::SeqIO;
   
4. use Bio::Seq;
   
5. 
   
6. my $fasta=shift @ARGV;
   
7. (my $Out=$fasta)=~s/(\.[^\.]+)$/_out$1/;
   
8. 
   
9. open (IN,"<$fasta") or die "Open $fasta failed!\n";
   
10. my $Out=Bio::SeqIO->new(-file=>">$Out",-format=>'fasta');
    
11. 
    
12. my ($flag,$seq,$id,$desc)=();
    
13. while(my $line=<IN>){
    
14.         chomp($line);
    
15.         if($line=~/\>/){
    
16.                 if($flag){
    
17.                         print "$id\t$desc\n$seq\n";
    
18.                         my $SEQ_OBJ=Bio::Seq->new(-seq=>$seq,-id=>$id,-desc=>$desc,-alphabet=>'protein');
    
19.                         $Out->write_seq($SEQ_OBJ);
    
20.                 }
    
21.                 $line=~s/[\"\>]//g;
    
22.                 ($id,$desc)=split (/\s+/,$line,2);
    
23.                 $seq=undef;
    
24.                 $flag=1;
    
25.         }else{
    
26.                 $seq.=$line;
    
27.         }
    
28. }
    
29. print "$id\t$desc\n$seq\n";
    
30. my $SEQ_OBJ=Bio::Seq->new(-seq=>$seq,-id=>$id,-desc=>$desc,-alphabet=>'protein');
    
31. $Out->write_seq($SEQ_OBJ);
    
32. 
    
33. $Out->close();

復(fù)制代碼

一串兒葡萄皮

受教了，這里果然是個(gè)好地方，謝謝回復(fù) 12# MMMIX


   

一串兒葡萄皮

(my $Out=$fasta)=~s/(\.[^\.]+)$/_out$1/;這一句怎么理解？O(∩_∩)O謝謝









回復(fù) 14# b114213903


   

MMMIX

一串兒葡萄皮 發(fā)表于 2015-08-12 15:20
(my $Out=$fasta)=~s/(\.[^\.]+)$/_out$1/;

等價(jià)于:

1. 
   
2. my $Out = $fasta;
   
3. $Out =~ s/(\.[^\.]+)$/_out$1/;
   

復(fù)制代碼

一串兒葡萄皮

1. 2
   
2. my ($flag,$seq,$id)=();
   
3. 
   
4. --------------------- WARNING ---------------------
   
5. MSG: sequence '1' doesn't validate, mismatch is  ($,,$,,$,),();
   
6. ---------------------------------------------------
   
7. 
   
8. ------------- EXCEPTION -------------
   
9. MSG: Attempting to set the sequence '1' to [my ($flag,$seq,$id)=();] which does not look healthy
   
10. STACK Bio::PrimarySeq::seq /usr/lib64/perl5/vendor_perl/5.16.0/Bio/PrimarySeq.pm:285
    
11. STACK Bio::PrimarySeq::new /usr/lib64/perl5/vendor_perl/5.16.0/Bio/PrimarySeq.pm:239
    
12. STACK Bio::Seq::new /usr/lib64/perl5/vendor_perl/5.16.0/Bio/Seq.pm:497
    
13. STACK toplevel delete_qw-A.pl:18
    
14. -------------------------------------
    

復(fù)制代碼

回復(fù) 14# b114213903

有這樣的報(bào)錯(cuò)？
   

b114213903

我這運(yùn)行正常，不知道什么原因

104359176

文件替換要考慮特殊情況，這種情況加一個(gè)位置限定符就好了：

1. /^\<"|"$/

復(fù)制代碼

特別聲明引號(hào)出現(xiàn)在開始尖括號(hào)的后面和行末。