[BioC] Blast a BSgenome

Ugo Borello ugo.borello at inserm.fr
Fri Jul 5 13:18:22 CEST 2013


Good morning,
I have a custom made Bsgenome data package and I would like to align a cDNA
sequence (~1kb) with the genome sequence of my Bsgenome, and retrieve the
aligned sequences.
Any suggestions? Any BLAST-like function in Bioconductor?
If I understand correctly Biostrings::matchPattern works with small
sequences, right? And I need in addition to account for gaps in the
alignment!
Does annotate::blastSequences could accept my Bsgenome genomic sequence as
argument? 
Thank you in advance for your help.

Ugo



More information about the Bioconductor mailing list