Beefy Boxes and Bandwidth Generously Provided by pair Networks
Your skill will accomplish
what the force of many cannot
 
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 taking refuge in the Monastery: (8)
As of 2014-10-24 11:31 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    For retirement, I am banking on:










    Results (131 votes), past polls