Use matchPattern for large DNAStrings
1
0
Entering edit mode
@vinicius-henrique-da-silva-6713
Last seen 11 months ago
Brazil

I would like to compare two DNAStrings that are larger than 20000 basepairs. I tried the matchPattern function from Biostrings package, however, I got this error:

Error in .valid.algos(pattern, max.mismatch, min.mismatch, with.indels,  :
  patterns with more than 20000 letters are not supported

There is a way to analyze these large strings in Bioconductor?

biostrings annotate • 1.0k views
ADD COMMENT
2
Entering edit mode
@herve-pages-1542
Last seen 1 day ago
Seattle, WA, United States

Hi Vinicius,

Maybe you want to try pairwiseAligment() for this. matchPattern() is probably not the right tool for comparing 2 strings. See ?pairwiseAligment for more information. The Biostrings package also has a nice vignette dedicated to Pairwise Sequence Alignments. See bioconductor.org/packages/Biostrings.

H.

ADD COMMENT

Login before adding your answer.

Traffic: 681 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