Beefy Boxes and Bandwidth Generously Provided by pair Networks Frank
No such thing as a small change
 
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].


Comment on Re: What is wrong in this code???
Select or Download Code

Log In?
Username:
Password:

What's my password?
Create A New User
Node Status?
node history
Node Type: note [id://1007827]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this? | Other CB clients
Other Users?
Others meditating upon the Monastery: (4)
As of 2014-04-21 01:11 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    April first is:







    Results (489 votes), past polls