[BioC] pairwiseAlignment and issue mapping

James Bullard bullard at berkeley.edu
Sat Aug 15 21:30:37 CEST 2009


Hi All,

this should be an easy question and answer, it is mostly about the  
best way to do this with the Bioconductor S4 classes. I have a  
pairwise alignment: A and I want a function which maps from a position  
in the original sequence (either subject: s, or pattern: p) to the  
position in the aignment. Essentially, I want a function that looks  
like:

subject  alignment
1		1
2		2
3		4
4		5

alignmentPosition(i, which = c("subject", "pattern")), i in  
1:length(s) => j in 1:length(alignment)

If I had just characters, it would be easy, but since I have objects  
of type: PairwiseAlignedFixedSubject, I am slightly overwhelmed by the  
class hierarchy, and the best way to do this.

thanks, jim



More information about the Bioconductor mailing list