Hello All,
I really dont know whats causing this problem. Its like i have the
input file which looks like:
:set list on in vi
TGCAGAAGAGAGTGAGCAC$
TGCAGAAGAGAGTGAGCAC $
TGCAGAAGAGAGTGAGCAC $
TGCAGAAGAGAGTGAGCAC $
TGCAGAAGAGAGTGAGCAC $
TGCAGAAGAGAGTGAGCAC $
CGCAGAAGAGAGTGAGCACA$
TTGCAGAAGAAAGAGAGCAC$
TTGCAGAAGATAGAGAGCAC$
(there are some trailing spaces in above sequences.)
And when i use the below script:
$filename = $ARGV[0];
open (FILE, "< $filename") or die "Cant open $filename : $!";
$sequence = '';
while(<FILE>)
{
next if /^(\s)*$/; # to skip blank lines in file
chomp;
$_=~s/\s+$//g;
$sequence=$_;
$compliment=$sequence;
$compliment =~tr/aAtTcCgGuU/TTAAGGCCTT/;
$rev_compliment=reverse $compliment;
print "sequence : $sequence\t";
print "rev_sequence : $rev_compliment\n";
$sequence = ' ';
$compliment=' ';
$rev_compliment=' ';
}
And when i see the output there is some difference in the tabs in
displays. why is that so. Any ideas
sequence : TTCAAGGACTTCTATTCAGAA rev_sequence : TTCTGAATAGAAGTCCTTGAA
sequence : TCTTGATGAAGAGGAATGGGA rev_sequence : TCCCATTCCTCTTCATCAAGA
sequence : GCTTGGCAATCTCATGTCAAA rev_sequence : TTTGACATGAGATTGCCAAGC
sequence : CTAGACCAATGTTTAGTTAGT rev_sequence : ACTAACTAAACATTGGTCTAG
sequence : TAGTTGTCGTCGTATCTTTGT rev_sequence : ACAAAGATACGACGACAACTA
sequence : GAGTAAGGCGCTTATCTTCT rev_sequence : AGAAGATAAGCGCCTTACTC
sequence : TCTGCAACCGGAAAGGGAGC rev_sequence : GCTCCCTTTCCGGTTGCAGA
last 2 outputs dont display my tabs why is that so.
I checked the escape characters in the file but cant find much in them
it look like :
sequence : TGCAGAAGAGAGTGAGCAC^Irev_sequence : GTGCTCACTCTCTTCTGCA$
sequence : TGCAGAAGAGAGTGAGCAC^Irev_sequence : GTGCTCACTCTCTTCTGCA$
sequence : TGCAGAAGAGAGTGAGCAC^Irev_sequence : GTGCTCACTCTCTTCTGCA$
sequence : CGCAGAAGAGAGTGAGCACA^Irev_sequence : TGTGCTCACTCTCTTCTGCG$
sequence : TTGCAGAAGAAAGAGAGCAC^Irev_sequence : GTGCTCTCTTTCTTCTGCAA$
sequence : TTGCAGAAGATAGAGAGCAC^Irev_sequence : GTGCTCTCTATCTTCTGCAA$
Please advice.
Thanks
_______________________________________________
ActivePerl mailing list
[email protected]
To unsubscribe: http://listserv.ActiveState.com/mailman/mysubs