Perl遍历文件中的每一行,并追加到另一个文件中每行的末尾 - 跟进

问题描述:

我对ealier文章有跟进问题。 相关的帖子是: Perl iterating through each line in a file and appending to the end of each line in another filePerl遍历文件中的每一行,并追加到另一个文件中每行的末尾 - 跟进

我用:

use warnings; 
use strict; 

open my $animals, '<', 'File1.txt' or die "Can't open animals: $!"; 
open my $payloads, '<', 'File2.txt' or die "Can't open payloads: $!"; 

my @payloads = <$payloads>; #each line of the file into an array 
close $payloads or die "Can't close payloads: $!"; 

while (my $line = <$animals>) { 
    chomp $line; 
    print $line.$_ foreach (@payloads); 
} 
close $animals or die "Can't close animals: $!"; 

这工作得很好,看起来像这样的文件:

file 1:  file 2: 
line1  lineA 
line2  lineB 
line3  lineC 

而不是看起来像这样的文件:

<01 line1 
<02 line2 

所以我想要做的是以下几点:

file 1:    file 2: 
<01 line1   <AA lineAA 
<02 line2   <AB lineAB 

应该变成:

file 3: 
<01_AA line1lineAA 
<01_AB line1lineAB 
<02_AA line2lineAA 
<02_AB line2lineAB 

我曾尝试通过使用标签上分割字符串来解决它,而在while循环(见下文),但我不能让循环它工作。

我的脚本:

#!C:/perl64/bin/perl.exe 

use warnings; 
use strict; 


open my $file1, '<', 'file1.fasta' or die "Can't open file1: $!"; 
open my $file2, '<', 'file2.fasta' or die "Can't open file2:  $!"; 
open(OUT, '>', 'file3.fasta') or die "Cannot write $!"; 



while (<$file2>) 
{ 
    chomp; 
    my ($F2_Id, @SF2_seq) = split (/\t/, $_); 

    while (<$file1>) 
     { 
      chomp; 
      my ($F1_Id, @F1_seq) = split (/\t/, $_); 
      foreach my $seq (@F1_seq) 
       { 
        print OUT $F1_Id,"_",$F2_Id,"\t",$seq.$_ foreach (@F2_seq),"\n"; 
       } 
      close; 
     } 
} 

我开始用Perl就在最近,所以我可以想像,有很多脚本故障。

对于真的很长的帖子我很抱歉,但我会appriciate任何帮助。

+0

是否所有行都看起来像“ fugu

+0

是的所有线条看起来像 Paeseitje

您可以将第一个文件的id和seq存储在数组数组中。

您还必须用_替换第二个文件中的<

#!/usr/bin/perl 
use warnings; 
use strict; 

open my $LEFT, '<', 'file1.fasta' or die "Can't open file1: $!"; 
open my $RIGHT, '<', 'file2.fasta' or die "Can't open file2: $!"; 
open my $OUT, '>', 'file3.fasta' or die "Cannot write: $!"; 

my @left; 
while (<$LEFT>) { 
    chomp; 
    push @left, [ split /\t/ ]; 
} 

while (<$RIGHT>) { 
    chomp; 
    my ($id, $seq) = split /\t/; 
    $id =~ s/</_/; 
    print {$OUT} "$_->[0]$id\t$_->[1]$seq\n" for @left; 
} 
close $OUT or die "Cannot close: $!"; 
+0

非常感谢你!还有1个小故障。尽管有s />/_ /,第二个'>'并没有被'_'取代。可能是什么原因? – Paeseitje

+0

发现它,有一个'',我查看了它。再次感谢! – Paeseitje