Hi!
I am having trouble aligning a DNAStringSet to a BSgenome. The example in
the help works for me....
Could anyone explain to me the error I am getting, please and how I can
make it work? Actually I would like to align these reads to a specific
gene... Is there maybe a better way?
Thanks!
Sonja


> alignShortReads(reads, Hsapiens, seqNames="chr1")
Error in .Call2("ACtree2_build", tb, pp_exclude, base_codes, nodebuf_ptr,
:
  non base DNA letter found in Trusted Band for pattern 248
> class(reads)
[1] "DNAStringSet"
attr(,"package")
[1] "Biostrings"

> class(Hsapiens)
[1] "BSgenome"
attr(,"package")
[1] "BSgenome"

	[[alternative HTML version deleted]]

