#!/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"; } } }