cmetindex missing in gmapr? Not possible to do bisulfite alignments
1
0
Entering edit mode
@karlnordstroem-14357
Last seen 5.8 years ago

I wanted to generate bisulfite alignments with gmapR, but cannot get it to work. I can generate a normal alignment similar to how it is described in the documentation.

Running gsnap on the command line, you normally have to run cmetindex in order to prepare a normal genomic index for bisulfite alignments. Do you need to do this within gmapR as well?

Looking to the code-base, there are empty wrappers for both cmetindex and atoiindex. They are not visible within R. Is there a plan to add this functionality? Or am I just missing something obvious?

 

gmapr • 900 views
ADD COMMENT
0
Entering edit mode
@michael-lawrence-3846
Last seen 2.3 years ago
United States

I guess we were waiting for a post like this to motivate support for those. They look like they've been implemented, but they've never been tested. You could try calling them, passing the GmapGenome object. Let me know how it goes. Eventually we can export and document them.

ADD COMMENT
0
Entering edit mode

Hey!

Well I'm working on my master thesis under supervision of Dr. Nordström, and I got your reply by email two days ago.

Actually I tried calling gmapR:::cmetindex(genome) and the result was like this:

Error in as.character.default(X[[i]], ...) :
   no method for coercing this S4 class to a vector

Any advice?
Thanks in advance

Best,

Amer

ADD REPLY
0
Entering edit mode

gmapR 1.21.1 (devel) should resolve this issue by exporting cmetindex and adding the missing coercion.

ADD REPLY

Login before adding your answer.

Traffic: 670 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6