#!/usr/bin/perl

use strict;
use warnings;

if ( scalar @ARGV < 2 ) {
    die("Usage:\n\t$0 file1 file2\n\n");
}

my %data;

my $fn = $ARGV[0];
{
    my $sequence;
    open DF, $fn or die $!;
    while ( my $line = <DF> ) {
        chomp $line;
        next if $line =~ m/^\s*$/;
        if ( substr( $line, 0, 1 ) eq '>' ) {
            ( $sequence, undef ) = split /\s+/, $line;
        }
        else {
            @{ $data{$sequence} } = split /\s+/, $line;
        }
    }
    close DF;
}

$fn = $ARGV[1];
{
    my $sequence;
    my $side;
    my $len;
    open DF, $fn or die $!;
    while ( my $line = <DF> ) {
        chomp $line;
        next if $line =~ m/^\s*$/;

        if ( substr( $line, 0, 1 ) eq '>' ) {
            ( $sequence, $side, undef, $len ) =
              split /(?:_|=|\s)+/, $line;
            if ( defined @{ $data{$sequence} } ) {
                my @range = ( 0, $len - 1 );
                if ( $side ne "left" ) {
                    @range = ( -1 * $len, -1 );
                }
                print $line, "\n";
                print join( ' ',
                    @{ $data{$sequence} }
                      [ $range[0] .. $range[1] ] ),
                  "\n";
            }
            else {
                warn "No data read from "
                  . $ARGV[0]
                  . " for sequence $sequence\n";
            }
        }
        else {
            print join( "  ", split //, $line ), "\n";
        }
    }
    close DF;
}
[download]