Hi Rsubread developers and R users,
I am using R version 3.1 and the associated "Rsubread" Bioconductor package.
The "subjunc" function allows the user to output junction reads in the
following format,
#Chr, StartLeftBlock, EndRightBlock, Junction_Name, nSupport, Strand,
StartLeftBlock, EndRightBlock, Color, nBlocks, BlockSizes, BlockStarts
chr10 94009 94602 JUNC00000001 251 + 94009 94602 255,0,0
2 46,47 0,546
chr10 94822 95396 JUNC00000002 993 - 94822 95396
0,255,255 2 30,49 0,525
My question is, for the second last column (BlockSizes column), is there a way
to modify the default maximum block size from 49 to other value (e.g. 20)? I
understand this may involve the modification of source C code and
re-compilation.
Thanks,
Xue
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
