Beefy Boxes and Bandwidth Generously Provided by pair Networks
go ahead... be a heretic
 
PerlMonks  

Re^2: find reverse complement of every 2nd line

by bingalee (Acolyte)
on Jul 09, 2013 at 23:12 UTC ( #1043392=note: print w/replies, xml ) Need Help??


in reply to Re: find reverse complement of every 2nd line
in thread find reverse complement of every 2nd line

Hey Dave. For the reverse complement the code will be something like this

my $raw=<IN>; $rev= reverse $raw; $rev=~tr/ATGCatgc/TACGtacg/;

my problem actually is doing this only to every other line and not the entire file

this is what i came up with, but it doesnt work

while(@inf=<IN>) { if(i%2==1) { $inf[i]=reverse $inf[i]; } }

In this i was just trying to find the reverse to see if this code performs the action on the lines I want it to.. but no luck. Im still thinking of any other way. Any suggestions?

Replies are listed 'Best First'.
Re^3: find reverse complement of every 2nd line
by LanX (Chancellor) on Jul 09, 2013 at 23:27 UTC
    variables in Perl need sigils, so it must be $i not just i

    (using strict and warnings should have prevented you)

    and you need to increment it with ++$i within the loop.

    As an alternative you could just test $. % 2 since $. automatically holds the line-number.

    Cheers Rolf

    ( addicted to the Perl Programming Language)

Re^3: find reverse complement of every 2nd line
by syphilis (Chancellor) on Jul 09, 2013 at 23:40 UTC
    my problem actually is doing this only to every other line and not the entire file

    The special variable $. (documented in perldoc perlvar) keeps count of the line numbers, so you can do something like:
    while(<IN>) { next if $. % 2; # skip odd lines # do the processing }
    Cheers,
    Rob

Log In?
Username:
Password:

What's my password?
Create A New User
Node Status?
node history
Node Type: note [id://1043392]
help
Chatterbox?
[Discipulus]: i'v installed ubuntu 16.04lts and i'm quite disappointed: is deadly simmilar to new win osses where you cannot find even simple things...
[Corion]: Hi Discipulus ;)
[Corion]: Mj1234: Maybe Ubuntu didn't package Devel::GDB? In that case, consider using local::lib and install the module from CPAN
[Corion]: But a quick Google search turns up this page, which seems to suggest that it should be packaged
[Discipulus]: iirc active perl has repositories only for last 2 o 3 major release of perl Mj1234. I strongly suggest portable version of strawberry perl and use cpan client to install modules
[choroba]: Also, who needs GDB when there's hdb?

How do I use this? | Other CB clients
Other Users?
Others lurking in the Monastery: (10)
As of 2016-12-07 08:55 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?
    On a regular basis, I'm most likely to spy upon:













    Results (125 votes). Check out past polls.