BLAST-like alignment of a DNAString to a BSgenome?
Entering edit mode
Last seen 3 days ago


I am looking for a way to align a sequence stored in a DNAString object to a BSgenome object with similar speed as with BLAST or long-read aligners. The sequence might be matching a repeated region multiple times, and checking for that is the reason why I would like to do the alignment. Since I find the sequence of interest during an interactive R session, I would prefer to use Bioconductor functions to run the alignment.

I found the previous question which suggested that there was only matchPattern available, which more precise and therefore limited to cases with at most ~15 mismatches. I have not found alternative, but if I missed something, please let me know !

-- Charles

Alignment Genome • 19 views

Login before adding your answer.

Traffic: 281 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6