Hi Monks, I have DNA data for a number of individuals at a few thousand markers. The data is in a matrix format with the rows being individuals and columns being markers. For example:

3851 A A G G T T A A C C ...
3854 A A G G T T A A C C ...
3871 A A G G T T A A G G ...
[download]

The first column is the individual ID, and each column after that is a marker observation. I need to count the number of occurrence of each character per column, and then output a text file such as:

column 2
3 A
column 10
2 C 1 G
[download]

I only know how to do this per column using the unix grep command, but that would be extremely time consuming and necessary. The output does not need to be as specific as I have put, as long as it gives the column number and how many occurrences of each character. Thank you for any advice.

---

---

Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!

Titles consisting of a single word are discouraged, and in most cases are disallowed outright.

Read Where should I post X? if you're not absolutely sure you're posting in the right place.

Please read these before you post! —

• How do I compose an effective node title?
• How do I post a question effectively?
• Markup in the Monastery

Posts may use any of the Perl Monks Approved HTML tags:

a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr

You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)

	For:		Use:
	&		&
	<		&lt;
	>		&gt;
	[		&#91;
	]		&#93;

Link using PerlMonks shortcuts! What shortcuts can I use for linking?

See Writeup Formatting Tips and other pages linked from there for more info.