Question

如何用Perl批量提取指定位置的DNA序列

Answer 1

#!/usr/bin/perl -w
use strict;
open my $fhin, $ARGV[0];
my ($seq,$g,$c);
while(my $line = <$fhin>){
chomp($line);
$seq .= $line;
print $line;
$g += () = $line =~ m/[Gg]/g;
$c += () = $line =~ m/[Cc]/g;
}

my $length = length($seq);

print "Sequence length is: ".$length."\nContains ".$g." Gs and ".$c." Cs\nGC percentage is: ".(($g+$c)/$length).".\n";