[BioC] A problem about trimLRPatterns
Harris A. Jaffee
hj at jhu.edu
Thu Mar 1 22:45:04 CET 2012
I'm not sure what you are doing. Let's change your "A" to "C", so we don't
confuse it with the lower case a's. Then we have this example, which looks
ok to me since the "C" is not trimmed.
pattern = "TCGACTG"
subject = "aaaaaaCTCGACTG"
trimLRPatterns(subject=subject, Rpattern=pattern, with.Rindels=TRUE, max.Rmismatch=1)
What happens to you? Please be explicit.
On Mar 1, 2012, at 4:05 PM, wang peter wrote:
> i have found a problem
> if we allow edit distance as 1, with indel=1
> subject aaaaaaATCGACTG
> pattern TCGACTG
> and i found some reads were trimmed ATCGACTG
> not TCGACTG
> that is ok for RNA-seq
> but not ok for microRNA, if we want to keep the whole microRNA sequences
> 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
> Bioconductor mailing list
> Bioconductor at r-project.org
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor