[Bioc-devel] readGAlignmentsList and read groups

Hi,

You're right, currently the only 'automatic' grouping is by pairs.

For single-end, grouping by tag would straightforward (assuming this is 
a "match on value" not a procedure of checking bit flags etc.).

For paired-end I'm assuming you want the records first grouped by pairs 
then grouped by tag? Maybe it's worth thinking about putting tags common 
to pairs on the outer list elements vs inner. For data are that aren't 
'too large' this would allow easier sorting/splitting.

Valerie

The function readGAlignmentsList[FromBam] sounds general, but it's purpose
currently seems specific to grouping by read pair. What about other types
of grouping, in particular the RG tag? This would be helpful when dealing
with sequencing data for many amplicons.

Just an idea.

Michael

	[[alternative HTML version deleted]]

_______________________________________________
Bioc-devel at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/bioc-devel

Valerie Obenchain

Program in Computational Biology
Division of Public Health Sciences
Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N, M1-B155
P.O. Box 19024
Seattle, WA 98109-1024

E-mail: vobencha at fhcrc.org
Phone:  (206) 667-3158
Fax:    (206) 667-1319