[BioC] a problem with trimLRPatterns

wang peter wng.peter at gmail.com
Fri Mar 9 20:15:37 CET 2012


i donot know why 0.1 mismatch can work to trim the correct adapter, but
if i set the max.Rmismatch from 1 to 9, it cannot work
thanks

> subject = "GGGAGTAAGAAAGGCACTGAAGGCACTATCAATAGATCGGAAGAGCGGTT"
> Rpattern = "CTGTAGGCACCATCAATAGATCGGAAGAGCGGTTCAGAAGGAATGCCGAG"

> trimLRPatterns(Rpattern = Rpattern, subject = subject, max.Rmismatch=1,with.Lindels=TRUE)
[1] "GGGAGTAAGAAAGGCACTGAAGGCACTATCAATAGATCGGAAGAGCGGTT"
> trimLRPatterns(Rpattern = Rpattern, subject = subject, max.Rmismatch=0.1,with.Lindels=TRUE)
[1] "GGGAGTAAGAAAGGCA"
> trimLRPatterns(Rpattern = Rpattern, subject = subject, max.Rmismatch=2,with.Lindels=TRUE)
[1] "GGGAGTAAGAAAGGCACTGAAGGCACTATCAATAGATCGGAAGAGCGGTT"
> trimLRPatterns(Rpattern = Rpattern, subject = subject, max.Rmismatch=3,with.Lindels=TRUE)
[1] "GGGAGTAAGAAAGGCACTGAAGGCACTATCAATAGATCGGAAGAGCGGTT"

-- 
shan gao
Room 231(Dr.Fei lab)
Boyce Thompson Institute
Cornell University
Tower Road, Ithaca, NY 14853-1801
Office phone: 1-607-254-1267(day)
Official email:sg839 at cornell.edu
Facebook:http://www.facebook.com/profile.php?id=100001986532253



More information about the Bioconductor mailing list