Extracting FASTA sequences based on position with perl script !!

Shruti Paniwala — Wed, 01 Mar 2017 17:10:11 -0600

#!/usr/bin/env perl
 
#Uses: perl sub-seq.pl input.txt range
 
use strict;
use warnings;
 
my $end   = pop;
my $start = pop;
local $/ = '>';
 
while (<>) {
    chomp;
    next unless /(.+)/;
    my ($header) = "$/$1_$start-$end\n";
    my $seq = ${^POSTMATCH};
    $seq =~ s/\s//g;
    print $header;
    print +( substr $seq, $start - 1, $end ) . "\n";
}

BOL: Extracting FASTA sequences based on position with perl script !!

Extracting FASTA sequences based on position with perl script !!