[BioC] report a possible bug of
wang peter
wng.peter at gmail.com
Tue Feb 28 17:20:56 CET 2012
hi all: i used ShortRead to trim some sequences
max.mismatchs <- 0.25*1:nchar(DNAString(PCR2rc))
trimmedCoords <- trimLRPatterns(Rpattern = PCR2rc, subject =
sread(highQuaReads), max.Rmismatch= max.mismatchs,
with.Rindels=T,ranges=T)
trimmedReads <- narrow(highQuaReads, start=start(trimmedCoords),
end=end(trimmedCoords))
but it appear some null lines, but it is the first time to happen in
my exprience
@HWI-ST132:506:D0CNUABXX:3:1101:4456:1871 1:N:0:TTCACA
+
@HWI-ST132:506:D0CNUABXX:3:1101:4331:1931 1:N:0:TTCACA
GGTGGCTGTAGTTTAGTGGTAAGAATTCTACG
+
?<?DDD?D:<CDFFFG<FCFACGGIII at AEEH
> sessionInfo()
R version 2.14.1 (2011-12-22)
Platform: x86_64-redhat-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
[7] LC_PAPER=C LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
>
--
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