用perl实现匹配多行的问题

nanshanjin

文件部分如下：
## Alignment 0: score=252.0 e_value=5e-08 N=6 chr1&chr7 minus
  0-  0:        vv1g01013787001        vv7g01028210001          3e-56
  0-  1:        vv1g01013782001        vv7g01028199001          3e-64
  0-  2:        vv1g01013772001        vv7g01028189001         5e-101
  0-  3:        vv1g01013765001        vv7g01028179001         4e-117
  0-  4:        vv1g01013763001        vv7g01028174001         5e-136
  0-  5:        vv1g01013735001        vv7g01028171001          1e-36
## Alignment 1: score=272.0 e_value=4.8e-06 N=6 chr10&chr12 plus
  1-  0:        vv10g01004803001        vv12g01030482001          6e-06
  1-  1:        vv10g01004809001        vv12g01030487001         2e-146
  1-  2:        vv10g01004811001        vv12g01030488001         2e-156
  1-  3:        vv10g01004812001        vv12g01030489001         1e-143
  1-  4:        vv10g01004820001        vv12g01030491001              0
  1-  5:        vv10g01012628001        vv12g01030508001          9e-38
## 。。。。。。。。。。。
然后我想匹配所有##行含有chr1&chr7的行和对应的下面几行，我写的代码如下：
use strict;
open IN,"<$ARGV[0]" or die $!;
open OUT,">$ARGV[1]" or die $!;

$/="##";
while(<IN>){         
           if(/(.*)chr10&chr12(.*)(\d+)$/){
               print OUT"$_\n";
               
}
}
没有任何输出，我知道匹配错了，望大神指教，不知道$/="##"这个用的对不对

MMMIX

nanshanjin 发表于 2015-09-23 21:57
然后我想匹配所有##行含有chr1&chr7的行和对应的下面几行，我写的代码如下：
...
           if(/(.*)chr10&chr12(.*)(\d+)$/){

你倒是给自己解释下，这两个怎么能对上？

nanshanjin

嘿嘿，求解释回复 2# MMMIX


   

sunzhiguolu

1. perl -ane 'if($F[0] eq "##"){$P=$F[6] eq "chr1&chr7"}print if($P)' f

复制代码

## Alignment 0: score=252.0 e_value=5e-08 N=6 chr1&chr7 minus
  0-  0:        vv1g01013787001        vv7g01028210001          3e-56
  0-  1:        vv1g01013782001        vv7g01028199001          3e-64
  0-  2:        vv1g01013772001        vv7g01028189001         5e-101
  0-  3:        vv1g01013765001        vv7g01028179001         4e-117
  0-  4:        vv1g01013763001        vv7g01028174001         5e-136
  0-  5:        vv1g01013735001        vv7g01028171001          1e-36

104359176

1. use 5.020;
   
2. use File::Slurp qw(read_file);
   
3. 
   
4. my $text = read_file("sample.txt");
   
5. my @lines = split('##', $text);
   
6. say grep(/chr1&chr7/, @lines);
   
7. 
   

复制代码