文件提取求助

披星戴月无悔

大家好，我有这样一些数据
# BLASTP 2.2.28+
# Query: AGL000001  locus=Scaffold1:332:538:+
# Database: /home/xuefeil/blastdata/db/lac_pro
# Fields: query acc., subject acc., evalue, alignment length, % identity
# 1 hits found
AGL000001        NC_015428.1_cdsid_YP_004397837.1        5e-16        68        50.00
# BLASTP 2.2.28+
# Query: AGL000002  locus=Scaffold1:660:842:+
# Database: /home/xuefeil/blastdata/db/lac_pro
# 0 hits found
# BLASTP 2.2.28+
# Query: AGL000003  locus=Scaffold1:1026:1379:+
# Database: /home/xuefeil/blastdata/db/lac_pro
# Fields: query acc., subject acc., evalue, alignment length, % identity
# 5 hits found
AGL000003        NC_015428.1_cdsid_YP_004397490.1        3e-05        57        43.86
AGL000003        NC_008530.1_cdsid_YP_815609.1        4e-05        74        29.73
AGL000003        NC_013198.1_cdsid_YP_003172141.1        1e-04        57        40.35
AGL000003        NC_017470.1_cdsid_YP_005854679.1        1e-04        59        35.59
AGL000003        NC_015428.1_cdsid_YP_004397954.1        0.002        81        32.10
# BLASTP 2.2.28+
# Query: AGL000004  locus=Scaffold1:1475:2236:+
# Database: /home/xuefeil/blastdata/db/lac_pro
# Fields: query acc., subject acc., evalue, alignment length, % identity
# 5 hits found
AGL000004        NC_015428.1_cdsid_YP_004399387.1        4e-17        206        31.07
AGL000004        NC_013198.1_cdsid_YP_003172141.1        6e-09        154        32.47
AGL000004        NC_006814.3_cdsid_YP_193465.1        1e-07        114        37.72
AGL000004        NC_015428.1_cdsid_YP_004397490.1        1e-07        234        27.35
AGL000004        NC_015428.1_cdsid_YP_004397954.1        3e-07        147        29.25

想提取含有0 hits found的数据，即得到这样的文件
# BLASTP 2.2.28+
# Query: AGL000002  locus=Scaffold1:660:842:+
# Database: /home/xuefeil/blastdata/db/lac_pro
# 0 hits found
我写了这样的脚本
#!/usr/bin/perl
use strict;
my @ff;
my $item;
open(F1,"< seq1.out") or die "false 1";
open(OUTPUT,'>',"out_file") or die "false 1";
my $f = <F1>;
my $ff = chomp $f;
while ( $ff ) {
push @ff, split /\# BLASTP 2.2.28\+/;
}
foreach $item(@ff) {
chomp;
   if ($item =~ m/# 0 hits/)  {
  print OUTPUT $item;
       }
}
电脑跑不出来，数据有几万条，求大神帮助

q1208c

第一条, 永远不要使用 $_, 这是我在知道perl有这个默认变量并且被坑了之后, 告诉我自己的.

旋转小马

没看明白，LZ你说的输出文件中不也前几行不也没有“0 hits found”的字符串嘛？？？

披星戴月无悔

谢谢你的回复，文件里面是像这样一个一个的记录
# BLASTP 2.2.28+
# Query: AGL000001  locus=Scaffold1:332:538:+
# Database: /home/xuefeil/blastdata/db/lac_pro
# Fields: query acc., subject acc., evalue, alignment length, % identity
# 1 hits found
AGL000001        NC_015428.1_cdsid_YP_004397837.1        5e-16        68        50.00这是一个记录

# BLASTP 2.2.28+
# Query: AGL000002  locus=Scaffold1:660:842:+
# Database: /home/xuefeil/blastdata/db/lac_pro
# 0 hits found
这是第二个记录，每个记录开头为# BLASTP 2.2.28+

想把记录里面含有# 0 hits found的提取出来，如上文第二个记录所示
回复 3# 旋转小马


   

披星戴月无悔

大神能看下我这该怎么解决么，我的代码有什么问题呢回复 2# q1208c


   

substr函数

1. my $f = <F1>;
   
2. my $ff = chomp $f;
   
3. while ( $ff ) {
   
4.     push @ff, split /\# BLASTP 2.2.28\+/;
   
5. }

复制代码

$ff 1
while (1) { forever }

披星戴月无悔

确实forever了，求指导，咋搞啊回复 6# substr函数


   

substr函数

1. #!/usr/bin/perl
   
2. 
   
3. my @print;
   
4. while (<DATA>) {
   
5.     next unless /^#/;
   
6.     push @print, $_;
   
7.     if (/^#\s(\d+)\shits\sfound/) {
   
8.         print @print if $1 == 0;
   
9.         @print = ();
   
10.     }
    
11. 
    
12. }
    
13. 
    
14. __DATA__
    
15. # BLASTP 2.2.28+
    
16. # Query: AGL000001  locus=Scaffold1:332:538:+
    
17. # Database: /home/xuefeil/blastdata/db/lac_pro
    
18. # Fields: query acc., subject acc., evalue, alignment length, % identity
    
19. # 1 hits found
    
20. AGL000001        NC_015428.1_cdsid_YP_004397837.1        5e-16        68        50.00
    
21. # BLASTP 2.2.28+
    
22. # Query: AGL000002  locus=Scaffold1:660:842:+
    
23. # Database: /home/xuefeil/blastdata/db/lac_pro
    
24. # 0 hits found
    
25. # BLASTP 2.2.28+
    
26. # Query: AGL000003  locus=Scaffold1:1026:1379:+
    
27. # Database: /home/xuefeil/blastdata/db/lac_pro
    
28. # Fields: query acc., subject acc., evalue, alignment length, % identity
    
29. # 5 hits found
    
30. AGL000003        NC_015428.1_cdsid_YP_004397490.1        3e-05        57        43.86
    
31. AGL000003        NC_008530.1_cdsid_YP_815609.1        4e-05        74        29.73
    
32. AGL000003        NC_013198.1_cdsid_YP_003172141.1        1e-04        57        40.35
    
33. AGL000003        NC_017470.1_cdsid_YP_005854679.1        1e-04        59        35.59
    
34. AGL000003        NC_015428.1_cdsid_YP_004397954.1        0.002        81        32.10
    
35. # BLASTP 2.2.28+
    
36. # Query: AGL000004  locus=Scaffold1:1475:2236:+
    
37. # Database: /home/xuefeil/blastdata/db/lac_pro
    
38. # Fields: query acc., subject acc., evalue, alignment length, % identity
    
39. # 5 hits found
    
40. AGL000004        NC_015428.1_cdsid_YP_004399387.1        4e-17        206        31.07
    
41. AGL000004        NC_013198.1_cdsid_YP_003172141.1        6e-09        154        32.47
    
42. AGL000004        NC_006814.3_cdsid_YP_193465.1        1e-07        114        37.72
    
43. AGL000004        NC_015428.1_cdsid_YP_004397490.1        1e-07        234        27.35
    
44. AGL000004        NC_015428.1_cdsid_YP_004397954.1        3e-07        147        29.25
    

复制代码

披星戴月无悔

真谢谢了，问题解决了 回复 8# substr函数