Beefy Boxes and Bandwidth Generously Provided by pair Networks
There's more than one way to do things
 
PerlMonks  

Re: What is wrong in this code???

by brap (Pilgrim)
on Dec 07, 2012 at 21:08 UTC ( #1007827=note: print w/replies, xml ) Need Help??


in reply to What is wrong in this code???

Hey anonymous,

I could be way off base here, but it almost sounds like you want an array of counters, where each element of the array represents the position of the line you're counting. Something along the order of (untested):

if ($split_seq_to_examine[$j] eq 'A') {$count_A[$j]++;} elsif ($split_seq_to_examine[$j] eq 'T') {$count_T[$j]++;} elsif ($split_seq_to_examine[$j] eq 'C') {$count_C[$j]++;} elsif ($split_seq_to_examine[$j] eq 'G') {$count_G[$j]++;} elsif ($split_seq_to_examine[$j] eq '-') {$count_non[$j]++;}

To find out what you have in, say, position 56, look at $count_A[56], $count_T[56], $count_G[56], and $count_C[56].

Log In?
Username:
Password:

What's my password?
Create A New User
Node Status?
node history
Node Type: note [id://1007827]
help
Chatterbox?
[Corion]: A good morning to everybody!

How do I use this? | Other CB clients
Other Users?
Others wandering the Monastery: (2)
As of 2017-08-17 07:01 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?
    Who is your favorite scientist and why?



























    Results (282 votes). Check out past polls.

    Notices?