[R] the functionality of outputing junction reads in Rsubread package

[Wang, Xue, Ph.D.]

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