Hi,
I tried your code too, but no success. is it because of the editor
display. I am working on Mac will that make a difference. I cannot get
correct tabs working?
On 4/29/07, Bill Luebkert <[EMAIL PROTECTED]> wrote:
> amit hetawal wrote:
> > Hello All,
> > I really dont know whats causing this problem. Its like i have the
> > input file which looks like:
>
> > 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
>
> Looks OK for me:
>
> use strict;
> use warnings;
> my $filename = $ARGV[0] || 'data/test.dat';
>
> open FILE, "< $filename" or die "Cant open $filename : $!";
> while (<FILE>) {
> chomp;
> next if /^(\s)*$/; # skip blank lines
> s/\s+$//g; # remove trailing spaces
>
> my $sequence = $_;
> my $compliment = $sequence;
> $compliment =~ tr/aAtTcCgGuU/TTAAGGCCTT/;
> my $rev_compliment = reverse $compliment;
> print "sequence : $sequence\t";
> print "rev_sequence : $rev_compliment\n";
> }
>
> __END__
> sequence : TGCAGAAGAGAGTGAGCAC rev_sequence : GTGCTCACTCTCTTCTGCA
> sequence : TGCAGAAGAGAGTGAGCAC rev_sequence : GTGCTCACTCTCTTCTGCA
> sequence : TGCAGAAGAGAGTGAGCAC rev_sequence : GTGCTCACTCTCTTCTGCA
> sequence : TGCAGAAGAGAGTGAGCAC rev_sequence : GTGCTCACTCTCTTCTGCA
> sequence : TGCAGAAGAGAGTGAGCAC rev_sequence : GTGCTCACTCTCTTCTGCA
> sequence : TGCAGAAGAGAGTGAGCAC rev_sequence : GTGCTCACTCTCTTCTGCA
> sequence : CGCAGAAGAGAGTGAGCACA rev_sequence : TGTGCTCACTCTCTTCTGCG
> sequence : TTGCAGAAGAAAGAGAGCAC rev_sequence : GTGCTCTCTTTCTTCTGCAA
> sequence : TTGCAGAAGATAGAGAGCAC rev_sequence : GTGCTCTCTATCTTCTGCAA
>
_______________________________________________
ActivePerl mailing list
[email protected]
To unsubscribe: http://listserv.ActiveState.com/mailman/mysubs
