#!/usr/bin/perl
use strict;
use warnings;
use Term::ANSIColor;
my $file = $ARGV[0];
if (@ARGV < 1){
print STDERR "Usage: $0 input_fasta_file\n";
exit 1;
}
my ($header, $sequence);
open (A, "<", $file) or die "Check the file: $!";
while (my $line = ){
chomp $line;
if ($line =~ /^(>.*)/){
$header = $1;
}
else{
$sequence .= $line;
}
}
close (A);
$sequence =~s/[\n\s]//;
my @sequence = split ("", $sequence);
#print "@sequence\n";
my @position;
my $pos_file = "sorted_position_walk.txt";
open (A, "<", $pos_file) or die "Check the file: $!";
while (my $line = ){
chomp $line;
my $pos = (split /\t/,$line)[0];
push(@position, $pos);
}
for (my $i=0;$i<=scalar(@sequence);$i++){
foreach my $value(@position){
if($value eq $sequence[$i]){
my $colortext = colored (($value+18), 'bold red');
print "$sequence[$i]\n";
}
}
}