’求助�。�！求助�。�！文件匹配輸出特定的內(nèi)容‘

stanley_tam

聽你這么說，貌似我就是這么處理的。。。要不驗證一下結(jié)果
anyway我就是給個代碼思路，畢竟需求會變的 回復(fù) 30# little_joe


   

sunzhiguolu

回復(fù) 29# little_joe
之前那個, 在判斷邏輯上有點偏差, 再試下,

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. sub collect_data{
   
6.     my (%ha, %hFilt, $id);
   
7.     open (my $fhA, '<', shift);
   
8.     while (defined (local $_ = <$fhA>)){
   
9.         if (/\A(\D\S+)/){
   
10.             %hFilt = () if (!exists $ha{$1});
    
11.             $id = $1;
    
12.             next;
    
13.         }
    
14.         s/\s+\z//;
    
15.         push (@{$ha{$id}}, $_) if (!$hFilt{$_}++);
    
16.     }
    
17.     close ($fhA);
    
18.     \%ha;
    
19. }
    
20. 
    
21. sub compare{
    
22.     my ($id, $rh, @aData) = (@_);
    
23.     my %hStat;
    
24.     foreach my $v (@{$rh->{$id}}){
    
25.         foreach my $V (grep {$v >= $_->[0] && $v <= $_->[2]} @aData){
    
26.             my @aChars = split (//, $V->[4]);
    
27.             next if ($aChars[$v - $V->[0]] eq '-');
    
28.             $hStat{$aChars[$v - $V->[0]]}++;
    
29.         }
    
30.     }
    
31.     if (%hStat){
    
32.         print "$id\t";
    
33.         print "$_:$hStat{$_}\t" for keys %hStat;
    
34.         print "\n";
    
35.     }
    
36. }
    
37. 
    
38. my $rha = collect_data (shift);
    
39. 
    
40. open (my $fh, '<', shift);
    
41. my ($Flag, $Offset, $ID, @aData) = 0;
    
42. while (defined (local $_ = <$fh>)){
    
43.     next if (/\A\s*\z/ or !/\A(?:Query|Sbjct)/);
    
44.     if (/\AQuery=\h+(\S+)/){
    
45.         compare ($ID, $rha, splice (@aData)) if (@aData);
    
46.         exists ($rha->{$1}) ? do {($ID, $Flag, $Offset) = ($1, 1, -1)} : ($Flag = 0);
    
47.         next;
    
48.     }
    
49.     next if (!$Flag);
    
50.     if (/\A(Query|Sbjct)\h+(\d+)\h+(\H+)\h+(\d+)/){
    
51.         push (@{$aData[$1 eq "Query" ? ++$Offset : $Offset]}, $2, $3, $4);
    
52.     }
    
53. }
    
54. close ($fh);
    
55. compare ($ID, $rha, splice (@aData)) if ($Flag);
    

復(fù)制代碼

perl abc.pl a b
-------------------------------
NP_415088.1-1   W:4
NP_415560.1-1   S:1
YP_025310.1-1   T:1     V:1
YP_588459.1-1   N:1     K:1
YP_002791252.1-1        S:1
NP_418697.1-1   E:7     W:7


   

little_joe

再次感謝，試了您剛才的代碼，結(jié)果如圖 但是YP_025310.1-1   T:1     V:1結(jié)果不太對，原因可能是將Query中的“-”當(dāng)成字母讀取了，文件中當(dāng)Query行中出現(xiàn)“-”時不會影響序號區(qū)間例如：Query：1 ASDFGHJKL 9和Query：1  ASD-FGHJKL 9它們兩個都是包含9個元素，“-”不會影響字母的序號例如：Query：1 ASD-FGHJKL 9
                 Sbjct：1 ASDFGHHLKJ 9假設(shè)文件A數(shù)值為4，則找到了Query行的F此時Sbjct行對應(yīng)G的回復(fù) 32# sunzhiguolu


   

little_joe

試了您剛才的代碼，結(jié)果是這樣的 不知道什么原因。請教請教！回復(fù) 31# stanley_tam


   

stanley_tam

我的代碼直接跑perl fii.pl（如果文件名字是fii.pl的話）
pl文件跟 A.word   B.word放在同一個目錄下。。。估計你得把文件名字重命名，或者修改main函數(shù)中的文件名字。
回復(fù) 34# little_joe


   

sunzhiguolu

回復(fù) 33# little_joe
對于 Query 行的 '-' 字符 在什么情況下 會對 Query 行產(chǎn)生影響? (針對 B 文件)
舉例如下:
A 文件:
YP_025310.1-1
73
78

B 文件:
Query= YP_025310.1-1
...
Query  59   FSGDMLPDGGINRRWLTPA-SKLRQPVVCLNRSGNMVGQSVLPMAG-LSGKKKATKNSLF  116
              G       IN R  T A S++ +P+          G+++    G    KK+  KNSL
Sbjct  61   APGWWDKPPLINARVETAATSRMFKPLW-------QHGRAICFADGWFEWKKEGDKNSLI  113
...

符合 A 文件中數(shù)值 73, 78 條件的只有這兩條記錄. (Query 59, Sbjct 61)

1.> 假設(shè) Query 59 這一行, 在 73 之前的位置出現(xiàn)了兩個 '-' 字符, 那么 Query 59 這行的 73 數(shù)值應(yīng)該落在那個字符上面?  (你在說明時, 請使用同一示例文件進(jìn)行說明)
2.> 還是只考慮 當(dāng) A 中 數(shù)值 (如: 78 ) 在 B Query 行正好處于 '-' 位置時才順延到下一個字符, 其他情形忽略.

   

little_joe

試了，是剛才運行的方式不對，謝謝��！回復(fù) 35# stanley_tam


   

little_joe

對于您的疑問
Query= YP_025310.1-1
...
Query  59   FSGDML--PDGGINRRWLTPA-SKLRQPVVCLNRSGNMVGQSVLPMAG-LSGKKKATKNSLF  116
              G       IN R  T A S++ +P+          G+++    G    KK+  KNSL
Sbjct  61   APGWWDKPPLINARVETAATSRMFKPLW-------QHGRAICFADGWFEWKKEGDKNSLI  113
...
如果Query行在73之前出現(xiàn)兩個“--”則73位還是W，但是對應(yīng)的Sbjct是不同的。因為-在Query行中出現(xiàn)其實是運行軟件加注的Gap，不是序列中的內(nèi)容，而Query和Sbjct前后的數(shù)字是指序列序號，也就是字母的序號，與“-”無關(guān)。
例如：
Query  58   VSKWVWKRGWPSSSSNRRKPAVTSRKLTCWDARLNGWRGSTATVRPATRQT-----LIRT  112
            V  W  + GW S +   R       +LT    +     G    +    RQ      + R
Sbjct  15   VQSWKQRDGWDSVAPISRVEMSLEARLTQLIIKPQKTGGDFKEIDLLGRQIERLARVNRY  74

Query  113  SLTATKAG----VANR--------KRIFSVTRPSKSWSRFSLSSLSTISCTGIAPGLSTA  160
            S T  +A     VANR        K+ FSVTRPSKSWSRFSLSSLS ISCTGIAPGLSTA
Sbjct  75   SQTGNEADLNPNVANRNKGGRRKPKKNFSVTRPSKSWSRFSLSSLSIISCTGIAPGLSTA  134

Query  161  SAIS  164
            SAIS
Sbjct  135  SAIS  138
雖然第一行Query和第3行Query均有“-”出現(xiàn)，但其序列編號并未將其包含在內(nèi)。回復(fù) 36# sunzhiguolu


   

sunzhiguolu

回復(fù) 38# little_joe
對于這種情況, B 文件 如何取值? (A 文件的數(shù)值:81)
Query= NP_418697.1-1
...
Query  34   NNICERAIRPVVMGRKAWLFAGSLVAGNRAAQIMSLLGT---------AVWSHMLGWRTS  84
            NNI E A+R V +GRK +LF GS   G R A + SL+GT         +   H+LG   
Sbjct  57   NNIAENALRAVSLGRKNFLFFGSDHGGERGALLYSLIGTCKLNDVDPESYLRHVLGVIAD  116

Query 34 行 應(yīng)該落在那個字符上面, Sbjct 57 行 應(yīng)該落在那個字符上面

   

little_joe

Query= NP_418697.1-1
...
Query  34   NNICERAIRPVVMGRKAWLFAGSLVAGNRAAQIMSLLGT---------AVWSHMLGWRTS  84
            NNI E A+R V +GRK +LF GS   G R A + SL+GT         +   H+LG   
Sbjct  57   NNIAENALRAVSLGRKNFLFFGSDHGGERGALLYSLIGTCKLNDVDPESYLRHVLGVIAD  116
文件A是81的話，Query行取值應(yīng)該為W，Sbjct行則取對應(yīng)的V，謝謝您。
也就是說當(dāng)Query行出現(xiàn)“-”不會影響數(shù)值匹配，但是如果Sbjct行出現(xiàn)“-”則Query對應(yīng)的Sbjct則會發(fā)生變化回復(fù) 39# sunzhiguolu