Question

请问用Linux或Perl如何从蛋白质序列文件中取出存储在另一个文件中的10个蛋白质ID的Fasta格式序列

Answer 1

#!/usr/bin/perl -w

use strict;
#contact perlcoder weixin
my %hash;
while(<DATA>)
{
    if($_=~/\S+/)
    {
        
    }
    else
    {
       next;    
    }
    chomp;
    my ($key)=$_=~/>(\S+)/;
    my $value=<DATA>;
    chomp($value);
    $hash{$key}=$value;    
}

print $hash{'Gh_A04G0739'},"\n";


__DATA__
>Gh_A12G2559              
ATGGCTACTTTCTTTGGCTATTTTAC
>Gh_A12G2554              
ATGGCGGGAACTATCCAATCCCTAAT
>Gh_A04G0755              
ATGGCCTCCGATCAGACCCTTTTTCA
>Gh_A04G0739              
ATGGAGAAAGCAAAACCTGAAGCACC
>Gh_A04G0738              
ATGGCGCAGGTTTTAGACGACGCTGA
>Gh_A07G1346              
ATGGTTCATTGTCTGCCAAAGGTTTC
>Gh_A07G1331              
ATGGCAGACAAGGATTCTTCAAGGCC
>Gh_A07G1334              
ATGACGACGCCAACTCGAGATGCAAT
>Gh_A07G1335              
ATGGCCGCAACTAGATTCCTCTCTCA
>Gh_A08G1510              
ATGGCTACTGCACCGATAAAGTCTCA
>Gh_A08G1433              
ATGGGTAAAACACCTACTGGCAAGGA

Answer 2

最好输入输出都有样例，贴出来。

否则无法给出答复