perl - 1 行のテキスト形式にして複数行のテキストに変換する方法は?

Question

私は per を使用している初心者なので、助けが必要です。たとえば、ファイルにアミノ酸配列があります。そのシーケンスは 1 行です。したがって、1 行に 60 個のアミノ酸が必要です。どうすればperlを使用してそれを行うことができますか?

Answer 1

結果の幅を設定する小さなプログラムを次に示しますが、アイデアが得られるはずです。

#!/usr/bin/perl
use strict;
use warnings;
use Bio::SeqIO;

my $in  = Bio::SeqIO->new( -file   => "fasta_junk.fasta" ,
                           -format => 'fasta');

my $out = Bio::SeqIO->new( -file   => '>test.dat',
                           -format => 'fasta');

my $lookup = 'GTGCCAGCAGCCGC';
$out->width(20);

while ( my $seq = $in->next_seq() ) {
    my $pos = index $seq->seq, $lookup;


    # if $pos != -1, ($lookup not found),
    # or $pos != 0, (found $lookup at first position, thus
    #   no preceding characters).
    if ($pos > 0) {
        my $trunc = $seq->trunc(1,$pos);
        $out->write_seq($trunc);
    }
}

この出力（幅20）を生成しました。

>LM1
AAGTCTGACGGAGCAACGCC
GCGTGTATGAAGAAGGTTTT
CGGATCGTAAAGTACTGTCC
GTTAGAGAAGAACAAGGATA
AGAGTAACTGCTTGTCCCTT
GACGGTATCTAACCAGAAAG
CCACGGCTAACTAC

fasta_junk.fasta ファイルは

>LM1
AAGTCTGACGGAGCAACGCCGCGTGTATGAAGAAGGTTTTCGGATCGTAA
AGTACTGTCCGTTAGAGAAGAACAAGGATAAGAGTAACTGCTTGTCCCTT
GACGGTATCTAACCAGAAAGCCACGGCTAACTACGTGCCAGCAGCCGCGG
TAATACGTAGGTGGCAAGCGTTGTCCGGATTTATTGGGCGTAAAGCGCGC
GCAGGCGGTCTTTTAAGTCTGATGTGAAAGCCCCCGGCTTAACCGGGGAG
GGTCATTGGAAACTGGAAGACTGGAGTGCAGAAGAGGAGAGTGGAATTCC
ACGTGTAGCGGTGAAATGCGTAGATATGTGGAGGAACACCAGTGGCGAAG
GCGACTCTCTGGTCTGTAACTGACGCTGAGGCGCGAAAGCGTGGGGAGCA
AACAGGATTAGATACCCTGGTAGTCCACGCCGT

結果を確認するために、さまざまな幅で自分で遊ぶことができます。

Answer 2

open my $infile, '<', "/path/to/sequencefile" or die $!;
open my $outfile, '>', "/path/to/newfile" or die $!;

while(my $line = <$infile>) {
    print $outfile join("\n", split(/\s/, $line)) . "\n";
}

close $infile;
close $outfile;