mismatch nb in Biostrings
1
0
Entering edit mode
carol white ▴ 680
@carol-white-2174
Last seen 7.7 years ago
European Union
Hi, How does it work with mismatch.min, mismatch.max in the functions of Biostrings for ex, vcountPattern? If mismatch.min=0, mismatch.max=2, then if there is no match, the function searches up to 2 matches, that is, if it find matches with 1mismatch, then it searches in the pattern in the next subject, otherwise it tries to find 2 matches in the same subject? Look forward to your reply, Carol [[alternative HTML version deleted]]
• 663 views
ADD COMMENT
0
Entering edit mode
@harris-a-jaffee-3972
Last seen 8.3 years ago
United States
First, see ?`lowlevel-matching`, and if you are really adventurous, the function nedit_at() in Biostrings/src/lowlevel_matching.c. The degree to which your pattern and subjects disagree, in a sense dictated by your with.indels=FALSE and fixed=TRUE arguments, is determined first, although actually only up to your max.mismatch limit. Then that (modified) ?distance" is tested against your min.mismatch constraint, as it says here: ?isMatchingAt?: If ?subject? is an XString object, then return the logical vector defined by: min.mismatch <= neditAt(...) <= max.mismatch On Jul 7, 2014, at 4:37 AM, carol white <wht_crl at="" yahoo.com=""> wrote: > Hi, > How does it work with mismatch.min, mismatch.max in the functions of Biostrings for ex, vcountPattern? If mismatch.min=0, mismatch.max=2, then if there is no match, the function searches up to 2 matches, that is, if it find matches with 1mismatch, then it searches in the pattern in the next subject, otherwise it tries to find 2 matches in the same subject? > > Look forward to your reply, > > Carol > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
ADD COMMENT

Login before adding your answer.

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