#! /usr/bin/perl

while (<>) {
    my $id = $_;
    chomp (my $seq = <>);
    if ($seq =~ /ATCAG/) {
        print $id;
        my @bases = split //, $seq;
        foreach (reverse @bases) {
            s/A/T/ or s/T/A/ or s/G/C/ or s/C/G/;
            print "$_";
        }
        print "\n";
    }
}

#!/usr/bin/perl -w
use strict;
use autodie;

open IN,'<',$ARGV[0];
open OUT,'>',$ARGV[1];

while (<IN>){
        next if /^>/;
        if (/^ATCAG/){
                $_ = reverse $_;
                $_ =~ tr/ATCG/TAGC/;
                print OUT "$_\n";
        }
}

5楼: Originally posted by bovine at 2014-06-13 21:22:26
#!/usr/bin/perl -w
use strict;
use autodie;

open IN,'<',$ARGV;
open OUT,'>',$ARGV;

while (<IN>{
        next if /^>/;
        if (/^ATCAG/){
                $_ = reverse $_;
                $_ =~ tr/ATCG/TAGC/;
...

#!/usr/bin/perl
@ARGV = qw{ 1.txt  2.txt };
open IN,'<',$ARGV[0];
open OUT,'>',$ARGV[1];
while(<IN>){
chomp;
if(/>/){
print OUT "$_\n";
}elsif(/^ATCAG/){
$_=reverse $_;
$_=~ tr/ATCG/TAGC/;
print OUT "$_\n";
}
else{
}
}