chr	pos	id	ref	alt	mat	pat	unphased	as.matrix(counts[sel, ])
chr17	16110466	rs78348579	A	G	FALSE	TRUE	FALSE	32
chr17	16110466	rs78348579	A	G	FALSE	TRUE	FALSE	0
chr17	16110466	rs78348579	A	G	FALSE	TRUE	FALSE	32