gevro

HI @ahl27, Thanks! What is your email address so I can send you a link and info for a docker that exactly reproduces the bug?

Also, how do I pass ASCII (raw) value 0 to BStringSet? This works: ``` > BStringSet(rawToChar(as.raw(70))) BStringSet object of length 1: width seq [1] 1 F ``` This doesn't work...

These ascii values are encoding kinetic information for PACbio sequencing data. they are a standard part of the Pacbio data format. So it would be useful to have ability to...

Sorry, I by mistake wrote ASCII. What I mean is simply integer values from 0 to 255. Here is an example tag from a read in a BAM file: ip:B:C,30,29,31,34,34,34,27,28,27,30,28,31,29,46,36,29,29,38,27,33,33,30,34,27,39,26,37,46,34,27,29,29,34,32,33,36,34,29,32,30,25,28,30,30,34,31,28,31,33,25,37,34,22,27,32,45,26,23,25,26,25,30,35,30,28,31,33,37,23,28,31,41,30,27,25,30,29,30,30,40,30,29,31,28,31,45,67,29,34,34,36,23,27,25,33,36,35,33,27,33,30,30,26,27,26,38,82,26,36,33,47,28,25,30,34,32,31,38,27,24,42,29,30,27,28,25,31,26,30,36,25,26,27,33,31,34,49,32,32,28,24,24,33,27,26,26,34,31,30,29,29,22,25,29,30,25,39,32,29,28,26,26,27,42,28,32,28,28,32,34,23,28,28,31,34,32,34,31,39,28,39,26,30,28,26,33,41,32,24,24,34,25,33,29,27,42,31,40,33,28,27,34,27,22,22,27,25,29,29,36,35,47,29,26,29,26,24,25,28,33,25,26,40,26,25,23,29,26,68,49,29,27,25,33,27,31,31,23,38,31,35,36,31,34,42,24,25,26,44,26,30,33,30,29,30,31,24,28,26,25,43,30,27,32,26,25,30,25,30,28,30,31,43,28,27,24,32,29,23,26,25,23,22,31,24,30,27,42,36,23,26,26,28,28,36,34,26,27,36,27,34,29,35,45,22,32,38,32,32,28,36,38,31,27,31,20,25,28,31,31,27,33,37,69,25,28,29,32,28,28,42,35,29,29,37,73,68,49,25,25,31,38,26,26,35,23,40,53,30,37,30,25,30,26,31,59,31,66,32,30,40,28,26,37,23,25,34,44,23,33,32...

That would be great - thanks!

Ok. I'm assuming per @hpages that you all will be able to create a better way to do this soon.

Hi all, As the original poster of this issue - thank you for the efforts! I also wanted to point out another situation / consideration: doing sequenceLayer on a BStringSet,...

Just to add - I think that sequenceLayer also would need to be changed to allow assignment of characters for special CIGAR operations outside the range of 0 to 255,...

Thanks for your thoughts. I defer to @hpages of course on potential future ways of addressing this additional issue. Or if there is a work-around (I can't think of one).

Thanks very much for trying to help here. Hopefully I can clarify why I'm stuck a bit better. I'm analyzing kinetic PacBio data, from the tags 'ip' and 'pw' in...