[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Omaha.pm] perl



The little program is so much better :)

On 4/19/12 4:24 PM, "Jay Hannah" <jay@jays.net> wrote:

>On Apr 19, 2012, at 1:59 PM, Klinkebiel, David L wrote:
>> In this instance the strand does not make a difference. If you wanted
>>to make the change based on a F (forward) or R (reverse) what would the
>>perl look like. Thanks
>
>At that point the command-line version is a little confusing in my
>opinion. It looks like this:
>
>   perl -anE 'say join "\t", @F[0,1], ($F[2] eq 'F' ? $F[1] + 50 : $F[1]
>- 50)' inputfile.txt
>
>So I'd probably write a little program instead:
>
>   use strict;
>   use 5.10.0;
>   while (<>) {
>      chomp;
>      my ($chr, $pos, $strand) = split /\t/;
>      if ($strand eq 'F') {
>         say join "\t", $chr, $pos, $pos + 50;
>      } else {
>         say join "\t", $chr, $pos, $pos - 50;
>      }
>   }
>
>So I could comment it, add more logic, whatever.
>
>HTH,
>
>j
>
>
>
>
>_______________________________________________
>Omaha-pm mailing list
>Omaha-pm@pm.org
>http://mail.pm.org/mailman/listinfo/omaha-pm