#!/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";
}
}
}
####
Fasta file:
>gi|292658763|ref|NM_014143.3| Homo sapiens CD274 molecule (CD274), transcript variant 1, mRNA
GGCGCAACGCTGAGCAGCTGGCGCGTCCCGCGCGGCCCCAGTTCTGCGCAGCTTCCCGAGGCTCCGCACC
AGCCGCGCTTCTGTCCGCCTGCAGGGCATTCCAGAAAGATGAGGATATTTGCTGTCTTTATATTCATGAC
CTACTGGCATTTGCTGAACGCATTTACTGTCACGGTTCCCAAGGACCTATATGTGGTAGAGTATGGTAGC
AATATGACAATTGAATGCAAATTCCCAGTAGAAAAACAATTAGACCTGGCTGCACTAATTGTCTATTGGG
AAATGGAGGATAAGAACATTATTCAATTTGTGCATGGAGAGGAAGACCTGAAGGTTCAGCATAGTAGCTA
CAGACAGAGGGCCCGGCTGTTGAAGGACCAGCTCTCCCTGGGAAATGCTGCACTTCAGATCACAGATGTG
AAATTGCAGGATGCAGGGGTGTACCGCTGCATGATCAGCTATGGTGGTGCCGACTACAAGCGAATTACTG
TGAAAGTCAATGCCCCATACAACAAAATCAACCAAAGAATTTTGGTTGTGGATCCAGTCACCTCTGAACA
TGAACTGACATGTCAGGCTGAGGGCTACCCCAAGGCCGAAGTCATCTGGACAAGCAGTGACCATCAAGTC
CTGAGTGGTAAGACCACCACCACCAATTCCAAGAGAGAGGAGAAGCTTTTCAATGTGACCAGCACACTGA
GAATCAACACAACAACTAATGAGATTTTCTACTGCACTTTTAGGAGATTAGATCCTGAGGAAAACCATAC
AGCTGAATTGGTCATCCCAGAACTACCTCTGGCACATCCTCCAAATGAAAGGACTCACTTGGTAATTCTG
GGAGCCATCTTATTATGCCTTGGTGTAGCACTGACATTCATCTTCCGTTTAAGAAAAGGGAGAATGATGG
ATGTGAAAAAATGTGGCATCCAAGATACAAACTCAAAGAAGCAAAGTGATACACATTTGGAGGAGACGTA
ATCCAGCATTGGAACTTCTGATCTTCAAGCAGGGATTCTCAACCTGTGGTTTAGGGGTTCATCGGGGCTG
AGCGTGACAAGAGGAAGGAATGGGCCCGTGGGATGCAGGCAATGTGGGACTTAAAAGGCCCAAGCACTGA
AAATGGAACCTGGCGAAAGCAGAGGAGGAGAATGAAGAAAGATGGAGTCAAACAGGGAGCCTGGAGGGAG
ACCTTGATACTTTCAAATGCCTGAGGGGCTCATCGACGCCTGTGACAGGGAGAAAGGATACTTCTGAACA
AGGAGCCTCCAAGCAAATCATCCATTGCTCATCCTAGGAAGACGGGTTGAGAATCCCTAATTTGAGGGTC
AGTTCCTGCAGAAGTGCCCTTTGCCTCCACTCAATGCCTCAATTTGTTTTCTGCATGACTGAGAGTCTCA
GTGTTGGAACGGGACAGTATTTATGTATGAGTTTTTCCTATTTATTTTGAGTCTGTGAGGTCTTCTTGTC
ATGTGAGTGTGGTTGTGAATGATTTCTTTTGAAGATATATTGTAGTAGATGTTACAATTTTGTCGCCAAA
CTAAACTTGCTGCTTAATGATTTGCTCACATCTAGTAAAACATGGAGTATTTGTAAGGTGCTTGGTCTCC
TCTATAACTACAAGTATACATTGGAAGCATAAAGATCAAACCGTTGGTTGCATAGGATGTCACCTTTATT
TAACCCATTAATACTCTGGTTGACCTAATCTTATTCTCAGACCTCAAGTGTCTGTGCAGTATCTGTTCCA
TTTAAATATCAGCTTTACAATTATGTGGTAGCCTACACACATAATCTCATTTCATCGCTGTAACCACCCT
GTTGTGATAACCACTATTATTTTACCCATCGTACAGCTGAGGAAGCAAACAGATTAAGTAACTTGCCCAA
ACCAGTAAATAGCAGACCTCAGACTGCCACCCACTGTCCTTTTATAATACAATTTACAGCTATATTTTAC
TTTAAGCAATTCTTTTATTCAAAAACCATTTATTAAGTGCCCTTGCAATATCAATCGCTGTGCCAGGCAT
TGAATCTACAGATGTGAGCAAGACAAAGTACCTGTCCTCAAGGAGCTCATAGTATAATGAGGAGATTAAC
AAGAAAATGTATTATTACAATTTAGTCCAGTGTCATAGCATAAGGATGATGCGAGGGGAAAACCCGAGCA
GTGTTGCCAAGAGGAGGAAATAGGCCAATGTGGTCTGGGACGGTTGGATATACTTAAACATCTTAATAAT
CAGAGTAATTTTCATTTACAAAGAGAGGTCGGTACTTAAAATAACCCTGAAAAATAACACTGGAATTCCT
####
Sorted_position_walk.txt
88
91
92
94
101
113
114
121
122
124
125
126
134
140
146
148
153
159
171
173
183